AMP 数学物理的发展 1687 - 9139 1687 - 9120 Hindawi出版公司 429835年 10.1155 / 2013/429835 429835年 研究文章 在媒体断裂数值Fractional-Calculus两相流模型 了雯雯 1 Changpin 1 扛鼎之作 2 斯利瓦斯塔瓦 H。 1 数学系 上海大学 上海200444 中国 shu.edu.cn 2 学校的数学和统计数据 湖北工程大学 孝感 湖北432000年 中国 hbeu.cn 2013年 6 8 2013年 2013年 17 05年 2013年 13 07年 2013年 13 07年 2013年 2013年 版权©2013了雯雯钟等。 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。

两相流的数值模拟裂隙多孔介质在地下流是一个重要的话题,环境问题,和石油油藏工程。传统的模型并不适用于很多情况下,因为它缺乏记忆裂隙介质的性质。在本文中,我们开发一个新的数值与部分时间导数公式裂隙多孔介质中两相流。在拟议的配方,不同的部分时间衍生品应用于骨折和矩阵区域,因为他们有不同的内存特性。我们进一步发展二级时间离散方法,它使用一个压力大时间步和一个小的时间步长饱和度。压力方程的隐式地在每一个大的时间步长,而饱和度是由一个显式的更新部分时间在每次分步计划。最后,进行数值试验验证了该数值模型的有效性。 1。介绍</t我tle> <p>多相流的数值模拟裂隙多孔介质在地下流动非常重要,环境问题,和石油油藏工程。与通常的异构多孔介质相比,媒体有两种截然不同的两个空间尺度:骨折骨折和矩阵。比矩阵骨折具有较高的渗透率,但他们的体积相比非常小的矩阵。几个概念模型(<gydF4y2Baxref ref-type="bibr" rid="B1"> 1</gydF4y2Baxref>- - - - - -<gydF4y2Baxref ref-type="bibr" rid="B23"> 11</gydF4y2Baxref>)开发了模拟多孔介质中的多相流动,例如,single-porosity模型,双孔隙度/双重渗透率模型,discrete-fracture模型。可以找到这些模型的评估(<gydF4y2Baxref ref-type="bibr" rid="B13"> 4</gydF4y2Baxref>]。</pgydF4y2Ba> <p>这些模型可以用常数处理好两相流的物理参数。然而,流体可能扰乱多孔形成导致粒子迁移导致毛孔堵塞或化学反应中扩大毛孔或减少毛孔的大小(<gydF4y2Baxref ref-type="bibr" rid="B7"> 12</gydF4y2Baxref>]。事实上,孔隙度和渗透率通常取决于流体压力和饱和度。不可压缩的两相流在石油油藏工程、孔隙度变化小的压力。然而,注入流体孔隙度在很大程度上可以改变。位移过程中石油通过水、注入水会使媒体湿,然后将增大孔隙度,渗透率的变化。因此,应该考虑物理参数的可变性建模现实的两相流。在多孔介质流,卡普托(<gydF4y2Baxref ref-type="bibr" rid="B3"> 13</gydF4y2Baxref>]研究了通量的行为在多孔介质使用记忆形式,取而代之的是一个普通的时间导数的分数阶导数。”(<gydF4y2Baxref ref-type="bibr" rid="B12"> 14</gydF4y2Baxref>]研究了部分时间导数的泛化先前Natale-Salusti模型非线性温度和压力波,流体多孔岩石中传播。</pgydF4y2Ba> <p>在本文中,我们将研究一个部分时间导数经典的两相流模型的泛化。普通的时间导数取代卡普托分数导数积分变量的下限。因为卡普托分数导数由卷积幂律的内核和普通函数的导数,它是一个有用的工具来描述一个幂律的频率变化的物理系数。由于孔隙度和渗透率的对比矩阵和骨折,建议制定使用不同的部分时间衍生品流的断裂和矩阵区域代表不同的内存性能。</pgydF4y2Ba> <p>研制了许多数值方法在文献中,例如,(<gydF4y2Baxref ref-type="bibr" rid="B8"> 15</gydF4y2Baxref>- - - - - -<gydF4y2Baxref ref-type="bibr" rid="B18"> 17</gydF4y2Baxref>]。为了模拟分数两相流,我们提出一种两级时间离散方法基于物理属性的压力快速变化小于饱和(<gydF4y2Baxref ref-type="bibr" rid="B4"> 18</gydF4y2Baxref>];也就是一个大时间步用于压力,以及一个小的时间步长饱和度。当地的分数微分方程广泛研究[<gydF4y2Baxref ref-type="bibr" rid="B24"> 19</gydF4y2Baxref>]。在我们的方法中,内存的饱和度限制在每个大时间步自压力取决于饱和度没有历史记忆。该方法也被视为经典的部分概括隐式压力显式饱和度(IMPES)方法(<gydF4y2Baxref ref-type="bibr" rid="B6"> 20.</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B10"> 21</gydF4y2Baxref>),这是一个受欢迎的时域方法用于多相流模拟。像IMPES,我们将耦合系统分为压力方程和饱和度方程基于两相流的物理过程的性质和治疗的饱和度和毛细管压力压力方程显式地消除其非线性。每个大时间步进一步分为几个子步骤,饱和度是由一个显式的更新部分时间计划。cell-centered有限差分方法(<gydF4y2Baxref ref-type="bibr" rid="B11"> 22</gydF4y2Baxref>)是用于空间离散化。最后,给出数值结果验证了该数值模型的有效性。</pgydF4y2Ba> </sec> <sec id="sec2"> <title>2。分两阶段模型不可压缩流</t我tle> <sec id="sec2.1"> <title>2.1。分级模型</t我tle> <p>卡普托分数阶导数给出的卷积幂律内核和普通函数的导数。所以这是一个有用的工具来考虑一个幂律的频率变异系数的一个简单的卷积。我们现在介绍卡普托分数阶导数的定义(<gydF4y2Baxref ref-type="bibr" rid="B18"> 17</gydF4y2Baxref>]<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M1"> <mml:mtable> <mml:mlabeledtr id="EEq1"> <mml:mtd> <mml:mtext> (1)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msubsup> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mi> y</米米l:mi> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> =</米米l:mo> <mml:mfrac> <mml:mrow> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> Γ</米米l:mi> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:mn mathvariant="normal"> 1</米米l:mn> <mml:mo> - - - - - -</米米l:mo> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> </mml:mrow> </mml:mfrac> <mml:mrow> <mml:msubsup> <mml:mo stretchy="false"> ∫</米米l:mo> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mrow> <mml:mi></mml:mi> </mml:mrow> </mml:mrow> <mml:msup> <mml:mrow> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:mi> t</米米l:mi> <mml:mo> - - - - - -</米米l:mo> <mml:mi> τ</米米l:mi> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mi> γ</米米l:mi> </mml:mrow> </mml:msup> <mml:mi> y</米米l:mi> <mml:mi> ′</米米l:mi> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:mi> τ</米米l:mi> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mi> d</米米l:mi> <mml:mi> τ</米米l:mi> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M2"> <mml:mi> y</米米l:mi> <mml:mo stretchy="false"> (</米米l:mo> <mml:mi> t</米米l:mi> <mml:mo stretchy="false"> )</米米l:mo> </mml:math> </inline-formula>是一个函数,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M3"> <mml:mn mathvariant="normal"> 0</米米l:mn> <mml:mo> <</米米l:mo> <mml:mi> γ</米米l:mi> <mml:mo> <</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:math> </inline-formula>。在这里,我们使用修改卡普托分数阶导数;也就是说,积分下限的(<gydF4y2Baxref ref-type="disp-formula" rid="EEq1"> 1</gydF4y2Baxref>)的函数<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M4"> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:math> </inline-formula>而不是像往常一样不断;也就是说,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M5"> <mml:mi> 一个</米米l:mi> <mml:mo> =</米米l:mo> <mml:mi> 一个</米米l:mi> <mml:mo stretchy="false"> (</米米l:mo> <mml:mi> t</米米l:mi> <mml:mo stretchy="false"> )</米米l:mo> </mml:math> </inline-formula>。分数阶导数可以用来描述复杂的问题,涉及到记忆的时间,因为外地的财产。从这个,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M6"> <mml:mi> 一个</米米l:mi> <mml:mo stretchy="false"> (</米米l:mo> <mml:mi> t</米米l:mi> <mml:mo stretchy="false"> )</米米l:mo> </mml:math> </inline-formula>表示时间的内存范围<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M7"> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:math> </inline-formula>,这记忆属性随时间。</pgydF4y2Ba> <p>使用卡普托分数导数,我们引入一个记忆形式为两阶段不可压缩和非混相流体在多孔介质中。由一个下标表示润湿阶段<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M8"> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:math> </inline-formula>和非润湿相<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M9"> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:math> </inline-formula>。让<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M10"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>的饱和阶段<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M11"> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:math> </inline-formula>。两阶段饱和度受到以下限制:<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M12"> <mml:mtable> <mml:mlabeledtr id="EEq2"> <mml:mtd> <mml:mtext> (2)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula></p> <p>在多孔介质流动速度<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M13"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>的每个阶段<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M14"> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:math> </inline-formula>被达西定律描述为<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M15"> <mml:mtable> <mml:mlabeledtr id="EEq3"> <mml:mtd> <mml:mtext> (3)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mo> - - - - - -</米米l:mo> <mml:mfrac> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> k</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> r</米米l:mi> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> μ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfrac> <mml:mi mathvariant="bold-italic"> K</米米l:mi> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:mo> ∇</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mi> g</米米l:mi> <mml:mo> ∇</米米l:mo> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> ,</米米l:mo> <mml:mo> </mml:mo> <mml:mi> α</米米l:mi> <mml:mo> =</米米l:mo> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> n</米米l:mi> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M16"> <mml:mrow> <mml:mi mathvariant="bold-italic"> K</米米l:mi> </mml:mrow> </mml:math> </inline-formula>多孔介质的绝对渗透率张量,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M17"> <mml:mrow> <mml:mi> g</米米l:mi> </mml:mrow> </mml:math> </inline-formula>重力加速度,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M18"> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> </mml:math> </inline-formula>是深度,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M19"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> k</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> r</米米l:mi> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M20"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> μ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M21"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M22"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>相对磁导率、粘度、压力,和密度分别为每个阶段。</pgydF4y2Ba> <p>每个阶段的质量守恒方程给出<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M23"> <mml:mtable> <mml:mlabeledtr id="EEq4"> <mml:mtd> <mml:mtext> (4)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mi> ϕ</米米l:mi> <mml:msubsup> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> </mml:msubsup> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> <mml:mo> </mml:mo> <mml:mi> α</米米l:mi> <mml:mo> =</米米l:mo> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> n</米米l:mi> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M24"> <mml:mrow> <mml:mi> ϕ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>介质的孔隙度和吗<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M25"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>是外部的质量流率。</pgydF4y2Ba> <p>非润湿相之间的差异和润湿相压力毛细管压力所描述的:<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M26"> <mml:mtable> <mml:mlabeledtr id="EEq5"> <mml:mtd> <mml:mtext> (5)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> </mml:mrow> </mml:msub> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> - - - - - -</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula></p> <p>表示<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M27"> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> k</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> r</米米l:mi> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> /</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> μ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> </mml:math> </inline-formula>和<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M28"> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:mi> ρ</米米l:mi> <mml:mi> g</米米l:mi> <mml:mi> z</米米l:mi> </mml:math> </inline-formula>。用(<gydF4y2Baxref ref-type="disp-formula" rid="EEq3"> 3</gydF4y2Baxref>)(<gydF4y2Baxref ref-type="disp-formula" rid="EEq4"> 4</gydF4y2Baxref>),我们得到<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M29"> <mml:mtable> <mml:mlabeledtr id="EEq6"> <mml:mtd> <mml:mtext> (6)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mi> ϕ</米米l:mi> <mml:msubsup> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> </mml:msubsup> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mi mathvariant="bold-italic"> K</米米l:mi> <mml:mo> ∇</米米l:mo> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> <mml:mo> </mml:mo> <mml:mi> α</米米l:mi> <mml:mo> =</米米l:mo> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> n</米米l:mi> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>求和的形式(<gydF4y2Baxref ref-type="disp-formula" rid="EEq6"> 6</gydF4y2Baxref>)两个阶段,考虑到(<gydF4y2Baxref ref-type="disp-formula" rid="EEq2"> 2</gydF4y2Baxref>),我们可以达到<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M30"> <mml:mtable> <mml:mlabeledtr id="EEq7"> <mml:mtd> <mml:mtext> (7)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mo> - - - - - -</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> <mml:mi mathvariant="bold-italic"> K</米米l:mi> <mml:mo> ∇</米米l:mo> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> - - - - - -</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:msub> <mml:mi mathvariant="bold-italic"> K</米米l:mi> <mml:mo> ∇</米米l:mo> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M31"> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:msub> </mml:math> </inline-formula>,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M32"> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:mo stretchy="false"> (</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> - - - - - -</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo stretchy="false"> )</米米l:mo> <mml:mi> g</米米l:mi> <mml:mi> z</米米l:mi> </mml:math> </inline-formula>和<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M33"> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:msub> </mml:math> </inline-formula>。此外,定义<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M34"> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mo> - - - - - -</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> <mml:mi mathvariant="bold-italic"> K</米米l:mi> <mml:mo> ∇</米米l:mo> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> </mml:math> </inline-formula>和<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M35"> <mml:msub> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> /</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> </mml:math> </inline-formula>。然后,我们有<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M36"> <mml:mtable> <mml:mlabeledtr id="EEq8"> <mml:mtd> <mml:mtext> (8)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>因此,润湿相饱和度的方程<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M37"> <mml:mtable> <mml:mlabeledtr id="EEq9"> <mml:mtd> <mml:mtext> (9)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mi> ϕ</米米l:mi> <mml:msubsup> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> </mml:msubsup> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula></p> <p>我们注意到上述两相流公式是一个分级扩展的配方<gydF4y2Baxref ref-type="bibr" rid="B13"> 4</gydF4y2Baxref>]。我们认为不可压缩两相流的质量古典保守配方:<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M38"> <mml:mtable> <mml:mlabeledtr id="EEq10"> <mml:mtd> <mml:mtext> (10)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ϕ</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> +</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> <mml:mo> </mml:mo> <mml:mi> α</米米l:mi> <mml:mo> =</米米l:mo> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> n</米米l:mi> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>孔隙度往往取决于流体压力和饱和度。在一些不可压缩的两相流问题,石油油藏工程,例如,孔隙度变化小的压力。然而,注入流体孔隙度在很大程度上可以改变。位移过程中石油通过水、注入水会使媒体湿,然后是孔隙度将会扩大。总结两个阶段的大规模保守的配方,我们获得<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M39"> <mml:mtable> <mml:mlabeledtr id="EEq11"> <mml:mtd> <mml:mtext> (11)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ϕ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> - - - - - -</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> <mml:mi mathvariant="bold-italic"> K</米米l:mi> <mml:mo> ∇</米米l:mo> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> - - - - - -</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:msub> <mml:mi mathvariant="bold-italic"> K</米米l:mi> <mml:mo> ∇</米米l:mo> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>由于孔隙度与压力变化小,我们假设<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M40"> <mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ϕ</米米l:mi> </mml:mrow> <mml:mo> /</米米l:mo> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mrow> <mml:mo> ≃</米米l:mo> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:math> </inline-formula>在(<gydF4y2Baxref ref-type="disp-formula" rid="EEq11"> 11</gydF4y2Baxref>),因此,压力方程得到相同的形式(<gydF4y2Baxref ref-type="disp-formula" rid="EEq9"> 9</gydF4y2Baxref>)。润湿相饱和度方程可以表示为<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M41"> <mml:mtable> <mml:mlabeledtr id="EEq12"> <mml:mtd> <mml:mtext> (12)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ϕ</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> +</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>由于孔隙度取决于<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M42"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>由链式法则,我们得到<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M43"> <mml:mtable> <mml:mlabeledtr id="EEq13"> <mml:mtd> <mml:mtext> (13)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ϕ</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> =</米米l:mo> <mml:mi> Υ</米米l:mi> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M44"> <mml:mi> Υ</米米l:mi> <mml:mo stretchy="false"> (</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo stretchy="false"> )</米米l:mo> <mml:mo> =</米米l:mo> <mml:mi> ϕ</米米l:mi> <mml:mo> +</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:msup> <mml:mrow> <mml:mi> ϕ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> ′</米米l:mi> </mml:mrow> </mml:msup> <mml:mo stretchy="false"> (</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo stretchy="false"> )</米米l:mo> </mml:math> </inline-formula>。一般来说,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M45"> <mml:mi> Υ</米米l:mi> <mml:mo stretchy="false"> (</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo stretchy="false"> )</米米l:mo> </mml:math> </inline-formula>可以用幂指数定律描述。有了这个论点,我们可以推断出部分经典的两阶段模型的泛化。这清楚地解释了使用分级模型的物理原因。</pgydF4y2Ba> <p>最后,我们完成了我们的模型的边界和初始条件。我们把边界<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M46"> <mml:mo> ∂</米米l:mo> <mml:mi> Ω</米米l:mi> </mml:math> </inline-formula>计算域的<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M47"> <mml:mrow> <mml:mi> Ω</米米l:mi> </mml:mrow> </mml:math> </inline-formula>分成两个不重叠的部分:Dirichelt部分<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M48"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi mathvariant="normal"> Γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>和牛曼部分<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M49"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi mathvariant="normal"> Γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>,在那里<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M50"> <mml:mo> ∂</米米l:mo> <mml:mi> Ω</米米l:mi> <mml:mo> =</米米l:mo> <mml:msup> <mml:mrow> <mml:mi mathvariant="normal"> Γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> ∪</米米l:mo> <mml:msup> <mml:mrow> <mml:mi mathvariant="normal"> Γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> </mml:msup> </mml:math> </inline-formula>。压力方程(<gydF4y2Baxref ref-type="disp-formula" rid="EEq6"> 6</gydF4y2Baxref>)是受如下边界条件:<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M51"> <mml:mtable> <mml:mlabeledtr id="EEq14"> <mml:mtd rowspan="2"> <mml:mtext> (14)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:mtext> 或</米米l:mtext> <mml:msub> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> =</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> </mml:msup> <mml:mtext> 在</米米l:mtext> <mml:msup> <mml:mrow> <mml:mi mathvariant="normal"> Γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> <mml:mtr> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ·</米米l:mo> <mml:mi mathvariant="bold-italic"> n</米米l:mi> <mml:mo> =</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mtext> 在</米米l:mtext> <mml:msup> <mml:mrow> <mml:mi mathvariant="normal"> Γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M52"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>是压力<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M53"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi mathvariant="normal"> Γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M54"> <mml:mrow> <mml:mi mathvariant="bold-italic"> n</米米l:mi> </mml:mrow> </mml:math> </inline-formula>向外单位法向量<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M55"> <mml:mo> ∂</米米l:mo> <mml:mi> Ω</米米l:mi> </mml:math> </inline-formula>,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M56"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>对流入率吗<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M57"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi mathvariant="normal"> Γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>。边界条件的饱和度<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M58"> <mml:mtable> <mml:mlabeledtr id="EEq16"> <mml:mtd> <mml:mtext> (15)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> </mml:msup> <mml:mtext> 在</米米l:mtext> <mml:msup> <mml:mrow> <mml:mi mathvariant="normal"> Γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>给出的初始饱和度润湿相<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M59"> <mml:mtable> <mml:mlabeledtr id="EEq17"> <mml:mtd> <mml:mtext> (16)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mtext> 在</米米l:mtext> <mml:mi> Ω</米米l:mi> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula></p> </sec> <sec id="sec2.2"> <title>2.2。Discrete-Fracture模型部分衍生品</t我tle> <p>在这里,discrete-fracture模型(<gydF4y2Baxref ref-type="bibr" rid="B13"> 4</gydF4y2Baxref>)扩展到的部分时间导数。将矩阵和骨折gridcell的discrete-fracture模型不同的几何尺寸;也就是说,如果域<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M60"> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:math> </inline-formula>维空间,矩阵区域<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M61"> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:math> </inline-formula>维空间,但骨折是简化为矩阵gridcell接口(<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M62"> <mml:mi> n</米米l:mi> <mml:mo> - - - - - -</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:math> </inline-formula>维。这种治疗将产生的长度尺度的对比的显式表示骨折孔径single-porosity模型,所以它能够大大提高计算效率,方便在实际实现。</pgydF4y2Ba> <p>我们现在整个域分解为两部分:矩阵<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M63"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> Ω</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>和断裂<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M64"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> Ω</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。骨折是矩阵块包围。我们使用下标<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M65"> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:math> </inline-formula>表示矩阵的下标<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M66"> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:math> </inline-formula>代表断裂系统。在矩阵域是由的压力<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M67"> <mml:mtable> <mml:mlabeledtr id="EEq18"> <mml:mtd> <mml:mtext> (17)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mo> - - - - - -</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ∇</米米l:mo> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> - - - - - -</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> n</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ∇</米米l:mo> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>这是受matrix-fracture接口条件:<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M68"> <mml:mtable> <mml:mlabeledtr id="EEq19"> <mml:mtd rowspan="2"> <mml:mtext> (18)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> <mml:mtr> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> <mml:mo> </mml:mo> <mml:mtext> 在</米米l:mtext> <mml:mo> ∂</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> Ω</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ∩</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> Ω</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mtr> </mml:mtable> </mml:math> </disp-formula>在断裂系统,我们表示的裂缝宽度<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M69"> <mml:mrow> <mml:mi> ε</米米l:mi> </mml:mrow> </mml:math> </inline-formula>和假设势常数沿裂缝宽度,然后获得骨折的压力方程<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M70"> <mml:mtable> <mml:mlabeledtr id="EEq20"> <mml:mtd> <mml:mtext> (19)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mo> - - - - - -</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ∇</米米l:mo> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> - - - - - -</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> n</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ∇</米米l:mo> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M71"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>是整个matrix-fracture界面传质。上面的公式类似于经典的模型。</pgydF4y2Ba> <p>如前所述,部分属性代表了孔隙度和渗透率的变化,有非凡的对比矩阵和骨折。结果,不同的部分时间衍生品应该用于流矩阵地区和断裂系统。利用部分时间导数,我们可以表达区域的饱和度方程矩阵<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M72"> <mml:mtable> <mml:mlabeledtr id="EEq21"> <mml:mtd> <mml:mtext> (20)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ϕ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:msubsup> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:msubsup> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>随着界面条件<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M73"> <mml:mtable> <mml:mlabeledtr id="EEq22"> <mml:mtd> <mml:mtext> (21)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> <mml:mo> </mml:mo> <mml:mtext> 在</米米l:mtext> <mml:mo> ∂</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> Ω</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ∩</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> Ω</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>同样,在断裂系统是由饱和度方程<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M74"> <mml:mtable> <mml:mlabeledtr id="EEq23"> <mml:mtd> <mml:mtext> (22)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ϕ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:msubsup> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:msubsup> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M75"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>代表了整个matrix-fracture界面传质。</pgydF4y2Ba> </sec> </sec> <sec id="sec3"> <title>3所示。数值方法</t我tle> <p>在本节中,我们将分两阶段模型的数值方法不可压缩流。在下面,我们专注于时间离散化方案。</pgydF4y2Ba> <p>我们首先把总时间间隔<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M76"> <mml:mo stretchy="false"> (</米米l:mo> <mml:mn mathvariant="normal"> 0</米米l:mn> <mml:mo> ,</米米l:mo> <mml:mi> T</米米l:mi> <mml:mo stretchy="false"> ]</米米l:mo> </mml:math> </inline-formula>成<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M77"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>相等的时间步骤<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M78"> <mml:mn mathvariant="normal"> 0</米米l:mn> <mml:mo> =</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> </mml:msup> <mml:mo> <</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> <mml:mo> <</米米l:mo> <mml:mo> ⋯</米米l:mo> <mml:mo> <</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:msup> <mml:mo> =</米米l:mo> <mml:mi> T</米米l:mi> </mml:math> </inline-formula>和表示时间步长<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M79"> <mml:msub> <mml:mrow> <mml:mi> h</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mrow> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> <mml:mo> /</米米l:mo> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mrow> </mml:math> </inline-formula>。这一次部门用于压力。自饱和比压力变化更快,我们使用一个较小的时间步长饱和度。每个子区间<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M80"> <mml:mo stretchy="false"> (</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> ,</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> <mml:mo stretchy="false"> ]</米米l:mo> </mml:math> </inline-formula>被划分到<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M81"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>sub-subintervals作为<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M82"> <mml:mo stretchy="false"> (</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> ,</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> <mml:mo stretchy="false"> ]</米米l:mo> <mml:mo> =</米米l:mo> <mml:mrow> <mml:msubsup> <mml:mo stretchy="false"> ⋃</米米l:mo> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> - - - - - -</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mrow> <mml:mi></mml:mi> </mml:mrow> </mml:mrow> <mml:mo stretchy="false"> (</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> ,</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> <mml:mo stretchy="false"> ]</米米l:mo> </mml:math> </inline-formula>,在那里<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M83"> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> </mml:msup> <mml:mo> =</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msup> </mml:math> </inline-formula>和<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M84"> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:msup> <mml:mo> =</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>,表示sub-subinterval长度<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M85"> <mml:msub> <mml:mrow> <mml:mi> h</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mrow> <mml:mrow> <mml:mrow> <mml:mo stretchy="false"> (</米米l:mo> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> <mml:mo> - - - - - -</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> <mml:mo stretchy="false"> )</米米l:mo> </mml:mrow> </mml:mrow> <mml:mo> /</米米l:mo> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mrow> </mml:math> </inline-formula>。表示一个变量的值<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M86"> <mml:mrow> <mml:mi> v</米米l:mi> </mml:mrow> </mml:math> </inline-formula>在<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M87"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>时间点,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M88"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> v</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>和一个<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M89"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>通过<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M90"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> v</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>。</pgydF4y2Ba> <p>对压力方程,饱和度值之前的时间步骤,和毛细管的潜力<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M91"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>在每一个细胞通过细胞显式计算饱和度从以前的时间步长和毛细管压力的功能。的变量<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M92"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M93"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M94"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>压力方程也通过细胞显式计算饱和度的前一个时间步。由此,我们获得的压力方程矩阵域:<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M95"> <mml:mtable> <mml:mlabeledtr id="EEq24"> <mml:mtd> <mml:mtext> (23)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mo> - - - - - -</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msubsup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msubsup> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ∇</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo> - - - - - -</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:msub> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msubsup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msubsup> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ∇</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mo> =</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>的上标<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M96"> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:math> </inline-formula>代表了时间步。方程(<gydF4y2Baxref ref-type="disp-formula" rid="EEq18"> 17</gydF4y2Baxref>)受matrix-fracture接口条件:<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M97"> <mml:mtable> <mml:mlabeledtr id="EEq25"> <mml:mtd> <mml:mtext> (24)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo> =</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo> ,</米米l:mo> <mml:mo> </mml:mo> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mo> =</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mo> ,</米米l:mo> <mml:mo> </mml:mo> <mml:mtext> 在</米米l:mtext> <mml:mo> ∂</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> Ω</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ∩</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> Ω</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>相似的表达形式在骨折(引用下标<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M98"> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:math> </inline-formula>),<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M99"> <mml:mtable> <mml:mlabeledtr id="EEq26"> <mml:mtd> <mml:mtext> (25)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:maligngroup></mml:maligngroup> <mml:malignmark></mml:malignmark> <mml:mo> - - - - - -</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msubsup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msubsup> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ∇</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo linebreak="newline" indentalign="left"></mml:mo> <mml:mo> </mml:mo> <mml:mo> </mml:mo> <mml:mo> - - - - - -</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:msub> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msubsup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msubsup> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:msub> <mml:mrow> <mml:mi mathvariant="bold-italic"> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ∇</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mo> =</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo> +</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula></p> <p>一旦压力<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M100"> <mml:mrow> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> </mml:mrow> </mml:math> </inline-formula>和<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M101"> <mml:mrow> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> </mml:mrow> </mml:math> </inline-formula>计算,速度可以评估<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M102"> <mml:mtable> <mml:mlabeledtr id="EEq27"> <mml:mtd> <mml:mtext> (26)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo> =</米米l:mo> <mml:mo> - - - - - -</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mi mathvariant="bold-italic"> K</米米l:mi> <mml:mo> ∇</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="normal"> Φ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula></p> <p>如前所述,下限<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M103"> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:math> </inline-formula>的积分(<gydF4y2Baxref ref-type="disp-formula" rid="EEq1"> 1</gydF4y2Baxref>随着时间的推移)是一个函数。这个函数可以选择根据实际问题。在多孔介质两相流压力快速变化小于饱和的时间(<gydF4y2Baxref ref-type="bibr" rid="B4"> 18</gydF4y2Baxref>),因此一个大时间步的压力。我们也可以看到,压力由饱和度但没有任何历史记忆。因此,内存的饱和度可以限制在每个时间步的压力。描述该内存的财产,我们定义<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M104"> <mml:mtable> <mml:mlabeledtr id="EEq28"> <mml:mtd> <mml:mtext> (27)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mi> 一个</米米l:mi> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> =</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> ,</米米l:mo> <mml:mo> </mml:mo> <mml:mi> t</米米l:mi> <mml:mo> ∈</米米l:mo> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> ,</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> <mml:mo class="right"> ]</米米l:mo> </mml:mrow> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>因此,卡普托分数导数的饱和度<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M105"> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> </mml:math> </inline-formula>被定义为<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M106"> <mml:mtable> <mml:mlabeledtr id="EEq29"> <mml:mtd> <mml:mtext> (28)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:maligngroup></mml:maligngroup> <mml:malignmark></mml:malignmark> <mml:msubsup> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> ,</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mi> 年代</米米l:mi> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> =</米米l:mo> <mml:mfrac> <mml:mrow> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> Γ</米米l:mi> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:mn mathvariant="normal"> 1</米米l:mn> <mml:mo> - - - - - -</米米l:mo> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> </mml:mrow> </mml:mfrac> <mml:mrow> <mml:msubsup> <mml:mo stretchy="false"> ∫</米米l:mo> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mrow> <mml:mi></mml:mi> </mml:mrow> </mml:mrow> <mml:msup> <mml:mrow> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:mi> t</米米l:mi> <mml:mo> - - - - - -</米米l:mo> <mml:mi> τ</米米l:mi> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mi> γ</米米l:mi> </mml:mrow> </mml:msup> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:mi> τ</米米l:mi> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mi> d</米米l:mi> <mml:mi> τ</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mo linebreak="newline" indentalign="left"></mml:mo> <mml:mo> </mml:mo> <mml:mo> </mml:mo> <mml:mo> </mml:mo> <mml:mo> </mml:mo> <mml:mo> </mml:mo> <mml:mo> </mml:mo> <mml:mo> </mml:mo> <mml:mo> </mml:mo> <mml:mo> </mml:mo> <mml:mo> </mml:mo> <mml:mo> </mml:mo> <mml:mo> </mml:mo> <mml:mi> t</米米l:mi> <mml:mo> ∈</米米l:mo> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> ,</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> <mml:mo class="right"> ]</米米l:mo> </mml:mrow> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>我们现在介绍显式时间近似的离散化方案<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M107"> <mml:msubsup> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mi> 年代</米米l:mi> </mml:math> </inline-formula>在<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M108"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>作为<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M109"> <mml:mtable> <mml:mlabeledtr id="EEq30"> <mml:mtd> <mml:mtext> (29)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:maligngroup></mml:maligngroup> <mml:malignmark></mml:malignmark> <mml:msubsup> <mml:mrow> <mml:mover accent="true"> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> <mml:mo> ~</米米l:mo> </mml:mover> </mml:mrow> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> ,</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> </mml:msubsup> <mml:msup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> <mml:mo linebreak="newline" indentalign="left"></mml:mo> <mml:mo> </mml:mo> <mml:mo> =</米米l:mo> <mml:mrow> <mml:mo symmetric="false" class="cases"> {</米米l:mo> <mml:mrow> <mml:mtable class="cases"> <mml:mtr> <mml:mtd columnalign="left"> <mml:msub> <mml:mrow> <mml:mi> b</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> <mml:mo> - - - - - -</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> ,</米米l:mo> </mml:mtd> <mml:mtd columnalign="left"> <mml:mi> j</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0</米米l:mn> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mtr> <mml:mtr> <mml:mtd columnalign="left"> <mml:msub> <mml:mrow> <mml:mi> b</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> <mml:mo> - - - - - -</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> </mml:mtd> <mml:mtd columnalign="left"></mml:mtd> </mml:mtr> <mml:mtr> <mml:mtd columnalign="left"> <mml:mo> </mml:mo> <mml:mo> +</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> b</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mrow> <mml:munderover> <mml:mo stretchy="false"> ∑</米米l:mo> <mml:mrow> <mml:mi> k</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:munderover> <mml:mrow> <mml:mi></mml:mi> </mml:mrow> </mml:mrow> <mml:msub> <mml:mrow> <mml:mi> b</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mo> - - - - - -</米米l:mo> <mml:mi> k</米米l:mi> </mml:mrow> </mml:msub> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> <mml:mo> - - - - - -</米米l:mo> <mml:mi> k</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> <mml:mo> - - - - - -</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> <mml:mo> - - - - - -</米米l:mo> <mml:mi> k</米米l:mi> </mml:mrow> </mml:msup> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> ,</米米l:mo> </mml:mtd> <mml:mtd columnalign="left"> <mml:mi> j</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> <mml:mo> ,</米米l:mo> <mml:mo> …</米米l:mo> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mtr> </mml:mtable> </mml:mrow> </mml:mrow> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M110"> <mml:mtable> <mml:mlabeledtr id="EEq31"> <mml:mtd rowspan="2"> <mml:mtext> (30)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> b</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mfrac> <mml:mrow> <mml:msubsup> <mml:mrow> <mml:mi> h</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mi> γ</米米l:mi> </mml:mrow> </mml:msubsup> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> Γ</米米l:mi> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:mn mathvariant="normal"> 2</米米l:mn> <mml:mo> - - - - - -</米米l:mo> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> </mml:mrow> </mml:mfrac> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> <mml:mtr> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> b</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> k</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msup> <mml:mrow> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:mi> k</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 1</米米l:mn> <mml:mo> - - - - - -</米米l:mo> <mml:mi> γ</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> - - - - - -</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> k</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 1</米米l:mn> <mml:mo> - - - - - -</米米l:mo> <mml:mi> γ</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> ,</米米l:mo> <mml:mo> </mml:mo> <mml:mi> k</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> <mml:mo> ,</米米l:mo> <mml:mo> …</米米l:mo> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> N</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> - - - - - -</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mtr> </mml:mtable> </mml:math> </disp-formula></p> <p>显式方案用于饱和度方程的矩阵域<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M111"> <mml:mtable> <mml:mlabeledtr id="EEq32"> <mml:mtd> <mml:mtext> (31)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ϕ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:msubsup> <mml:mrow> <mml:mover accent="true"> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> <mml:mo> ~</米米l:mo> </mml:mover> </mml:mrow> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> ,</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:msubsup> <mml:msubsup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo> +</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msubsup> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msubsup> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> =</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> 米</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>和裂缝网络<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M112"> <mml:mtable> <mml:mlabeledtr id="EEq33"> <mml:mtd> <mml:mtext> (32)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ϕ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:msubsup> <mml:mrow> <mml:mover accent="true"> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> <mml:mo> ~</米米l:mo> </mml:mover> </mml:mrow> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msup> <mml:mo> ,</米米l:mo> <mml:msup> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:msubsup> <mml:msubsup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo> +</米米l:mo> <mml:mo> ∇</米米l:mo> <mml:mo> ·</米米l:mo> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msubsup> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msubsup> <mml:msubsup> <mml:mrow> <mml:mi mathvariant="bold-italic"> u</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> =</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo> +</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi> 问</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> ,</米米l:mo> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula></p> <p>空间离散化方案,cell-centered有限差分方法用于压力方程,而逆风有限体积方法对饱和度方程。对空间离散化方案的详细描述,我们指的是(<gydF4y2Baxref ref-type="bibr" rid="B15"> 23</gydF4y2Baxref>]。</pgydF4y2Ba> </sec> <sec id="sec4"> <title>4所示。数值测试</t我tle> <p>在本节中,提供了两个数值例子证明拟议的两相流数值模型随着时间部分衍生品。</pgydF4y2Ba> <p>在所有的测试中,绝对渗透率是斜张量和多孔介质是各向同性的。我们使用下面的毛细管压力函数(<gydF4y2Baxref ref-type="bibr" rid="B14"> 24</gydF4y2Baxref>]:<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M113"> <mml:mtable> <mml:mlabeledtr id="EEq34"> <mml:mtd> <mml:mtext> (33)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> </mml:mrow> </mml:msub> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> =</米米l:mo> <mml:mo> - - - - - -</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> B</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> </mml:mrow> </mml:msub> <mml:mi> 日志</米米l:mi> <mml:mo></mml:mo> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M114"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> B</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>是一个积极的绝对渗透率参数相关。两个阶段的相对渗透率计算<d我年代p- - - - - -formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M115"> <mml:mtable> <mml:mlabeledtr id="EEq35"> <mml:mtd> <mml:mtext> (34)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> k</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> r</米米l:mi> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 3</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo> ,</米米l:mo> <mml:mo> </mml:mo> <mml:mo> </mml:mo> <mml:msub> <mml:mrow> <mml:mi> k</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> r</米米l:mi> <mml:mi> n</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msup> <mml:mrow> <mml:mrow> <mml:mo class="left"> (</米米l:mo> <mml:mrow> <mml:mn mathvariant="normal"> 1</米米l:mn> <mml:mo> - - - - - -</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> 年代</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mo class="right"> )</米米l:mo> </mml:mrow> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 3</米米l:mn> </mml:mrow> </mml:msup> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula></p> <p>我们考虑一个水平多孔介质<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M116"> <mml:mn> 20.</米米l:mn> <mml:mo> </mml:mo> <mml:mtext> 米</米米l:mtext> <mml:mo> ×</米米l:mo> <mml:mn> 15</米米l:mn> <mml:mo> </mml:mo> <mml:mtext> 米</米米l:mtext> <mml:mo> ×</米米l:mo> <mml:mn> 1</米米l:mn> <mml:mo> </mml:mo> <mml:mtext> 米</米米l:mtext> </mml:math> </inline-formula>与多个相互关联的骨折(<gydF4y2Baxref ref-type="bibr" rid="B15"> 23</gydF4y2Baxref>),如图所示<gydF4y2Baxref ref-type="fig" rid="fig1"> 1</gydF4y2Baxref>。骨折的宽度是0.01米。骨折的疏密度矩阵和媒体0.15和1,分别。渗透率的矩阵块和50个md和骨折<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M117"> <mml:mn mathvariant="normal"> 1</米米l:mn> <mml:msup> <mml:mrow> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 5</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>医学博士。水和油的粘度都等于1 cP。注射速率为0.2 PV /年。</pgydF4y2Ba> <fig id="fig1"> <label>图1</gydF4y2Balabel> <p>骨折的分布。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2013/429835.fig.001"></graphic> </fig> <p>因为中等水平,它是合理的忽略重力的影响。我们注入水的左端介质,最初是谁的空隙完全饱和油,生产石油在右边。没有其他的注入和没有提取领域的内部。通量向外面的其他边界消失。</pgydF4y2Ba> <p>图<gydF4y2Baxref ref-type="fig" rid="fig2"> 2</gydF4y2Baxref>显示部分的影响时间泛化平均含水饱和度元太的不同点。的传奇人物<gydF4y2Baxref ref-type="fig" rid="fig2"> 2</gydF4y2Baxref>,这个案子<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M118"> <mml:mo stretchy="false"> (</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo stretchy="false"> ]</米米l:mo> </mml:math> </inline-formula>代表了部分的订单时间导数的含水饱和度断裂和矩阵区域。从图<gydF4y2Baxref ref-type="fig" rid="fig2"> 2</gydF4y2Baxref>,我们观察时间衰减较慢,相比普通的情况。</pgydF4y2Ba> <fig id="fig2"> <label>图2</gydF4y2Balabel> <p>含水饱和度与不同部分衍生品。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2013/429835.fig.002"></graphic> </fig> <p>数据<gydF4y2Baxref ref-type="fig" rid="fig3"> 3</gydF4y2Baxref>,<gydF4y2Baxref ref-type="fig" rid="fig4"> 4</gydF4y2Baxref>,<gydF4y2Baxref ref-type="fig" rid="fig5"> 5</gydF4y2Baxref>,<gydF4y2Baxref ref-type="fig" rid="fig6"> 6</gydF4y2Baxref>,<gydF4y2Baxref ref-type="fig" rid="fig7"> 7</gydF4y2Baxref>,<gydF4y2Baxref ref-type="fig" rid="fig8"> 8</gydF4y2Baxref>,<gydF4y2Baxref ref-type="fig" rid="fig9"> 9</gydF4y2Baxref>,<gydF4y2Baxref ref-type="fig" rid="fig10"> 10</gydF4y2Baxref>,<gydF4y2Baxref ref-type="fig" rid="fig11"> 11</gydF4y2Baxref>显示了含水饱和度轮廓在不同时间与间歇时间导数的三双。从这些数据,我们可以看到的存在时间延迟效应当分数阶小于1。这表明,孔隙度和渗透率变化由于注入水润湿媒体,使流体流动缓慢。</pgydF4y2Ba> <fig id="fig3"> <label>图3</gydF4y2Balabel> <p>水饱和度等值线为1.25年<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M119"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:math> </inline-formula>和<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M120"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:math> </inline-formula>。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2013/429835.fig.003"></graphic> </fig> <fig id="fig4"> <label>图4</gydF4y2Balabel> <p>水饱和度等值在5年<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M121"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:math> </inline-formula>和<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M122"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:math> </inline-formula>。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2013/429835.fig.004"></graphic> </fig> <fig id="fig5"> <label>图5</gydF4y2Balabel> <p>水饱和度等值在10年<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M123"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:math> </inline-formula>和<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M124"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:math> </inline-formula>。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2013/429835.fig.005"></graphic> </fig> <fig id="fig6"> <label>图6</gydF4y2Balabel> <p>水饱和度等值线为1.25年<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M125"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0.9</米米l:mn> </mml:math> </inline-formula>和<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M126"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0.95</米米l:mn> </mml:math> </inline-formula>。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2013/429835.fig.006"></graphic> </fig> <fig id="fig7"> <label>图7</gydF4y2Balabel> <p>水饱和度等值在5年<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M127"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0.9</米米l:mn> </mml:math> </inline-formula>和<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M128"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0.95</米米l:mn> </mml:math> </inline-formula>。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2013/429835.fig.007"></graphic> </fig> <fig id="fig8"> <label>图8</gydF4y2Balabel> <p>水饱和度等值在10年<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M129"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0.9</米米l:mn> </mml:math> </inline-formula>和<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M130"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0.95</米米l:mn> </mml:math> </inline-formula>。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2013/429835.fig.008"></graphic> </fig> <fig id="fig9"> <label>图9</gydF4y2Balabel> <p>水饱和度等值线为1.25年<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M131"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0.85</米米l:mn> </mml:math> </inline-formula>和<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M132"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0.9</米米l:mn> </mml:math> </inline-formula>。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2013/429835.fig.009"></graphic> </fig> <fig id="fig10"> <label>图10</gydF4y2Balabel> <p>水饱和度等值在5年<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M133"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0.85</米米l:mn> </mml:math> </inline-formula>和<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M134"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0.9</米米l:mn> </mml:math> </inline-formula>。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2013/429835.fig.0010"></graphic> </fig> <fig id="fig11"> <label>图11</gydF4y2Balabel> <p>水饱和度等值在10年<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M135"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0.85</米米l:mn> </mml:math> </inline-formula>和<我nline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M136"> <mml:msub> <mml:mrow> <mml:mi> γ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0.9</米米l:mn> </mml:math> </inline-formula>。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2013/429835.fig.0011"></graphic> </fig> </sec> <back> <ack> <title>确认</t我tle> <p>这项工作是由上海市教育委员会的关键程序下格兰特(没有。12 zz084)和中国教育部(没有的关键项目。212109)。</pgydF4y2Ba> </ack> <ref-list> <ref id="B1" content-type="article"> <label>1</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 秋雨</年代urn一个米e> <given-names> r·G。</g我ven-names> </name> <name> <surname> 阿奈特</年代urn一个米e> <given-names> r . C。</g我ven-names> </name> <name> <surname> 朗格弗德</年代urn一个米e> <given-names> d . W。</g我ven-names> </name> </person-group> <article-title> 通过有限元技术模拟流体在fractured-porous岩体</一个rt我cle-title> <source> <italic> 国际期刊液体中的数值方法</我t一个lic> <year> 1984年</ygydF4y2Baear> <volume> 4</vgydF4y2Baolume> <issue> 4</我年代年代ue> <fpage> 337年</fp一个ge> <lpage> 348年</gydF4y2Balpage> <pub-id pub-id-type="other"> 2 - s2.0 - 0021639798</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="doi"> 10.1002 / fld.1650040404</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="other"> ZBL0579.76095</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B2" content-type="article"> <label>2</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> Barenblatt</年代urn一个米e> <given-names> g . I。</g我ven-names> </name> <name> <surname> Zheltov</年代urn一个米e> <given-names> i . P。</g我ven-names> </name> <name> <surname> Kochina</年代urn一个米e> <given-names> i . N。</g我ven-names> </name> </person-group> <article-title> 基本概念均匀液体在裂缝性渗流理论的岩石(地层)</一个rt我cle-title> <source> <italic> 应用数学和力学》期刊上</我t一个lic> <year> 1960年</ygydF4y2Baear> <volume> 24</vgydF4y2Baolume> <issue> 5</我年代年代ue> <fpage> 1286年</fp一个ge> <lpage> 1303年</gydF4y2Balpage> <pub-id pub-id-type="other"> 2 - s2.0 - 36248938881</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="doi"> 10.1016 / 0021 - 8928 (60)90107 - 6</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="other"> ZBL0104.21702</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B5" content-type="book"> <label>3</gydF4y2Balabel> <nlm-citation publication-type="book"> <person-group person-group-type="author"> <name> <surname> 陈</年代urn一个米e> <given-names> Z。</g我ven-names> </name> <name> <surname> 欢</年代urn一个米e> <given-names> G。</g我ven-names> </name> <name> <surname> 马</年代urn一个米e> <given-names> Y。</g我ven-names> </name> </person-group> <source> <italic> 计算方法在多孔介质中多相流</我t一个lic> <year> 2006年</ygydF4y2Baear> <publisher-loc> 美国费城,宾夕法尼亚州</pubgydF4y2Balisher-loc> <publisher-name> 暹罗</pubgydF4y2Balisher-name> <pub-id pub-id-type="doi"> 10.1137/1.9780898718942</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="other"> MR2217767</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B13" content-type="article"> <label>4</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> Hoteit</年代urn一个米e> <given-names> H。</g我ven-names> </name> <name> <surname> Firoozabadi</年代urn一个米e> <given-names> 一个。</g我ven-names> </name> </person-group> <article-title> 一个有效的媒体不可压缩两相流数值模型的断裂</一个rt我cle-title> <source> <italic> 水资源的进步</我t一个lic> <year> 2008年</ygydF4y2Baear> <volume> 31日</vgydF4y2Baolume> <issue> 6</我年代年代ue> <fpage> 891年</fp一个ge> <lpage> 905年</gydF4y2Balpage> <pub-id pub-id-type="other"> 2 - s2.0 - 43649093294</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="doi"> 10.1016 / j.advwatres.2008.02.004</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B16" content-type="article"> <label>5</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 齐米。</年代urn一个米e> <given-names> H。</g我ven-names> </name> </person-group> <article-title> 压力瞬态分析天然裂缝性储层的均匀断裂分布</一个rt我cle-title> <source> <italic> SPE杂志</我t一个lic> <year> 1969年</ygydF4y2Baear> <volume> 9</vgydF4y2Baolume> <issue> 4</我年代年代ue> <fpage> 451年</fp一个ge> <lpage> 462年</gydF4y2Balpage> <pub-id pub-id-type="other"> 2 - s2.0 - 0014618506</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="doi"> 10.2118 / 2156 - a</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B17" content-type="article"> <label>6</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 李</年代urn一个米e> <given-names> s . H。</g我ven-names> </name> <name> <surname> 詹森</年代urn一个米e> <given-names> c . L。</g我ven-names> </name> <name> <surname> 湖</年代urn一个米e> <given-names> m F。</g我ven-names> </name> </person-group> <article-title> 高效的有限差分模型与多个长度尺度骨折流在一个水库</一个rt我cle-title> <source> <italic> SPE杂志</我t一个lic> <year> 2000年</ygydF4y2Baear> <volume> 5</vgydF4y2Baolume> <issue> 3</我年代年代ue> <fpage> 268年</fp一个ge> <lpage> 275年</gydF4y2Balpage> <pub-id pub-id-type="other"> 2 - s2.0 - 0034261209</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B19" content-type="article"> <label>7</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> Noorishad</年代urn一个米e> <given-names> J。</g我ven-names> </name> <name> <surname> 麦哈</年代urn一个米e> <given-names> M。</g我ven-names> </name> </person-group> <article-title> 一个上游有限元方法解决裂隙多孔介质中瞬态输运方程</一个rt我cle-title> <source> <italic> 水资源研究</我t一个lic> <year> 1982年</ygydF4y2Baear> <volume> 18</vgydF4y2Baolume> <issue> 3</我年代年代ue> <fpage> 588年</fp一个ge> <lpage> 596年</gydF4y2Balpage> <pub-id pub-id-type="other"> 2 - s2.0 - 0020332120</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="doi"> 10.1029 / WR018i003p00588</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B20" content-type="article"> <label>8</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> Pruess</年代urn一个米e> <given-names> K。</g我ven-names> </name> <name> <surname> 纳史木汗</年代urn一个米e> <given-names> t . N。</g我ven-names> </name> </person-group> <article-title> 一个实际的方法建模在裂隙多孔介质流体和热流</一个rt我cle-title> <source> <italic> SPE杂志</我t一个lic> <year> 1985年</ygydF4y2Baear> <volume> 25</vgydF4y2Baolume> <issue> 1</我年代年代ue> <fpage> 14</fp一个ge> <lpage> 26</gydF4y2Balpage> <pub-id pub-id-type="other"> 2 - s2.0 - 0022011972</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B21" content-type="article"> <label>9</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 太阳</年代urn一个米e> <given-names> 年代。</g我ven-names> </name> <name> <surname> Firoozabadi</年代urn一个米e> <given-names> 一个。</g我ven-names> </name> <name> <surname> 口</年代urn一个米e> <given-names> J。</g我ven-names> </name> </person-group> <article-title> 数值模拟两相二元流体混合使用混合有限元素</一个rt我cle-title> <source> <italic> 计算地球科学</我t一个lic> <year> 2012年</ygydF4y2Baear> <volume> 16</vgydF4y2Baolume> <issue> 4</我年代年代ue> <fpage> 1101年</fp一个ge> <lpage> 1124年</gydF4y2Balpage> <pub-id pub-id-type="doi"> 10.1007 / s10596 - 012 - 9306 - 2</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B22" content-type="article"> <label>10</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 托马斯。</年代urn一个米e> <given-names> l·K。</g我ven-names> </name> <name> <surname> 迪克森</年代urn一个米e> <given-names> t . N。</g我ven-names> </name> <name> <surname> 皮尔森</年代urn一个米e> <given-names> r·G。</g我ven-names> </name> </person-group> <article-title> 裂缝性储层模拟</一个rt我cle-title> <source> <italic> SPE杂志</我t一个lic> <year> 1983年</ygydF4y2Baear> <volume> 23</vgydF4y2Baolume> <issue> 1</我年代年代ue> <fpage> 42</fp一个ge> <lpage> 54</gydF4y2Balpage> <pub-id pub-id-type="other"> 2 - s2.0 - 0020635993</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="doi"> 10.2118 / 9305 - pa</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B23" content-type="article"> <label>11</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 沃伦</年代urn一个米e> <given-names> j·E。</g我ven-names> </name> <name> <surname> 根</年代urn一个米e> <given-names> p . J。</g我ven-names> </name> </person-group> <article-title> 天然裂缝性储层的行为</一个rt我cle-title> <source> <italic> SPE杂志</我t一个lic> <year> 1963年</ygydF4y2Baear> <volume> 3</vgydF4y2Baolume> <issue> 3</我年代年代ue> <fpage> 245年</fp一个ge> <lpage> 255年</gydF4y2Balpage> <pub-id pub-id-type="doi"> 10.2118 / 426 - pa</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B7" content-type="article"> <label>12</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> di朱塞佩</年代urn一个米e> <given-names> E。</g我ven-names> </name> <name> <surname> 莫洛尼</年代urn一个米e> <given-names> M。</g我ven-names> </name> <name> <surname> 卡普托</年代urn一个米e> <given-names> M。</g我ven-names> </name> </person-group> <article-title> 通量与记忆:多孔介质模型和实验</一个rt我cle-title> <source> <italic> 多孔介质中传输</我t一个lic> <year> 2010年</ygydF4y2Baear> <volume> 83年</vgydF4y2Baolume> <issue> 3</我年代年代ue> <fpage> 479年</fp一个ge> <lpage> 500年</gydF4y2Balpage> <pub-id pub-id-type="other"> 2 - s2.0 - 77953871762</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="doi"> 10.1007 / s11242 - 009 - 9456 - 4</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B3" content-type="article"> <label>13</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 卡普托</年代urn一个米e> <given-names> M。</g我ven-names> </name> </person-group> <article-title> 通量在多孔介质模型的记忆</一个rt我cle-title> <source> <italic> 水资源研究</我t一个lic> <year> 2000年</ygydF4y2Baear> <volume> 36</vgydF4y2Baolume> <issue> 3</我年代年代ue> <fpage> 693年</fp一个ge> <lpage> 705年</gydF4y2Balpage> <pub-id pub-id-type="other"> 2 - s2.0 - 0034019492</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="doi"> 10.1029/1999 wr900299</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B12" content-type="article"> <label>14</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 食</年代urn一个米e> <given-names> R。</g我ven-names> </name> </person-group> <article-title> Fractional-calculus模型流体的温度和压力波多孔岩石</一个rt我cle-title> <source> <italic> 物理评论E</我t一个lic> <year> 2011年</ygydF4y2Baear> <volume> 84年</vgydF4y2Baolume> <issue> 3</我年代年代ue> <lpage> 6</gydF4y2Balpage> <pub-id pub-id-type="publisher-id"> 036605年</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="doi"> 10.1103 / PhysRevE.84.036605</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B8" content-type="article"> <label>15</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> Diethelm</年代urn一个米e> <given-names> K。</g我ven-names> </name> <name> <surname> 福特</年代urn一个米e> <given-names> n . J。</g我ven-names> </name> <name> <surname> 释放</年代urn一个米e> <given-names> 答:D。</g我ven-names> </name> <name> <surname> Luchko</年代urn一个米e> <given-names> 余。</g我ven-names> </name> </person-group> <article-title> 算法的分数微积分:数值方法的选择</一个rt我cle-title> <source> <italic> 计算机在应用力学和工程方法</我t一个lic> <year> 2005年</ygydF4y2Baear> <volume> 194年</vgydF4y2Baolume> <issue> 6 - 8</我年代年代ue> <fpage> 743年</fp一个ge> <lpage> 773年</gydF4y2Balpage> <pub-id pub-id-type="other"> 2 - s2.0 - 10644238068</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="doi"> 10.1016 / j.cma.2004.06.006</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B9" content-type="article"> <label>16</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 丁</年代urn一个米e> <given-names> H。</g我ven-names> </name> <name> <surname> 李</年代urn一个米e> <given-names> C。</g我ven-names> </name> </person-group> <article-title> 混合样条函数方法reaction-subdiffusion方程</一个rt我cle-title> <source> <italic> 计算物理学杂志</我t一个lic> <year> 2013年</ygydF4y2Baear> <volume> 242年</vgydF4y2Baolume> <fpage> 103年</fp一个ge> <lpage> 123年</gydF4y2Balpage> <pub-id pub-id-type="doi"> 10.1016 / j.jcp.2013.02.014</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B18" content-type="article"> <label>17</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 李</年代urn一个米e> <given-names> C。</g我ven-names> </name> <name> <surname> 陈</年代urn一个米e> <given-names> 一个。</g我ven-names> </name> <name> <surname> 叶</年代urn一个米e> <given-names> J。</g我ven-names> </name> </person-group> <article-title> 数值方法分数微积分和部分常微分方程</一个rt我cle-title> <source> <italic> 计算物理学杂志</我t一个lic> <year> 2011年</ygydF4y2Baear> <volume> 230年</vgydF4y2Baolume> <issue> 9</我年代年代ue> <fpage> 3352年</fp一个ge> <lpage> 3368年</gydF4y2Balpage> <pub-id pub-id-type="doi"> 10.1016 / j.jcp.2011.01.030</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="other"> MR2780467</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="other"> ZBL1218.65070</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B4" content-type="article"> <label>18</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 陈</年代urn一个米e> <given-names> Z。</g我ven-names> </name> <name> <surname> 欢</年代urn一个米e> <given-names> G。</g我ven-names> </name> <name> <surname> 李</年代urn一个米e> <given-names> B。</g我ven-names> </name> </person-group> <article-title> 一种改进的多孔介质中两相流IMPES方法</一个rt我cle-title> <source> <italic> 多孔介质中传输</我t一个lic> <year> 2004年</ygydF4y2Baear> <volume> 54</vgydF4y2Baolume> <issue> 3</我年代年代ue> <fpage> 361年</fp一个ge> <lpage> 376年</gydF4y2Balpage> <pub-id pub-id-type="doi"> 10.1023 / B: TIPM.0000003667.86625.15</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="other"> MR2031831</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B24" content-type="article"> <label>19</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 杨</年代urn一个米e> <given-names> X.-J。</g我ven-names> </name> <name> <surname> 斯利瓦斯塔瓦</年代urn一个米e> <given-names> h . M。</g我ven-names> </name> <name> <surname> 他</年代urn一个米e> <given-names> 黄永发。</g我ven-names> </name> <name> <surname> Baleanu</年代urn一个米e> <given-names> D。</g我ven-names> </name> </person-group> <article-title> Cantor-type圆柱坐标方法,微分方程与当地部分衍生品</一个rt我cle-title> <source> <italic> 物理信</我t一个lic> <year> 2013年</ygydF4y2Baear> <volume> 377年</vgydF4y2Baolume> <issue> 28 - 30</我年代年代ue> <fpage> 1696年</fp一个ge> <lpage> 1700年</gydF4y2Balpage> <pub-id pub-id-type="doi"> 10.1016 / j.physleta.2013.04.012</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B6" content-type="article"> <label>20.</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 外套</年代urn一个米e> <given-names> k . H。</g我ven-names> </name> </person-group> <article-title> 储层模拟:最先进的</一个rt我cle-title> <source> <italic> 《石油技术</我t一个lic> <year> 1982年</ygydF4y2Baear> <volume> 34</vgydF4y2Baolume> <issue> 8</我年代年代ue> <fpage> 1633年</fp一个ge> <lpage> 1642年</gydF4y2Balpage> <pub-id pub-id-type="doi"> 10.2118 / 10020 - pa</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B10" content-type="article"> <label>21</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 教唆犯</年代urn一个米e> <given-names> r·G。</g我ven-names> </name> <name> <surname> 斯图尔特</年代urn一个米e> <given-names> c . H。</g我ven-names> <suffix> Jr。</年代uff我x> </name> </person-group> <article-title> 二维多相水库模拟器的新方法</一个rt我cle-title> <source> <italic> SPE杂志</我t一个lic> <year> 1966年</ygydF4y2Baear> <volume> 6</vgydF4y2Baolume> <issue> 2</我年代年代ue> <fpage> 175年</fp一个ge> <lpage> 182年</gydF4y2Balpage> <pub-id pub-id-type="doi"> 10.2118 / 1188 - pa</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B11" content-type="article"> <label>22</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 福赛斯</年代urn一个米e> <given-names> p。</g我ven-names> <suffix> Jr。</年代uff我x> </name> <name> <surname> 马斯</年代urn一个米e> <given-names> p . H。</g我ven-names> </name> </person-group> <article-title> 二次收敛cell-centered网格</一个rt我cle-title> <source> <italic> 应用数值数学</我t一个lic> <year> 1988年</ygydF4y2Baear> <volume> 4</vgydF4y2Baolume> <issue> 5</我年代年代ue> <fpage> 377年</fp一个ge> <lpage> 394年</gydF4y2Balpage> <pub-id pub-id-type="doi"> 10.1016 / 0168 - 9274 (88)90016 - 5</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="other"> MR948505</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="other"> ZBL0651.65086</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B15" content-type="article"> <label>23</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 口</年代urn一个米e> <given-names> J。</g我ven-names> </name> <name> <surname> 太阳</年代urn一个米e> <given-names> 年代。</g我ven-names> </name> <name> <surname> 余</年代urn一个米e> <given-names> B。</g我ven-names> </name> </person-group> <article-title> 多尺度多尺度多重物理量time-splitting战略过程的两相流媒体骨折</一个rt我cle-title> <source> <italic> 应用数学学报</我t一个lic> <year> 2011年</ygydF4y2Baear> <volume> 2011年</vgydF4y2Baolume> <lpage> 24</gydF4y2Balpage> <pub-id pub-id-type="publisher-id"> 861905年</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="doi"> 10.1155 / 2011/861905</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="other"> MR2777043</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> <ref id="B14" content-type="article"> <label>24</gydF4y2Balabel> <nlm-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> Hoteit</年代urn一个米e> <given-names> H。</g我ven-names> </name> <name> <surname> Firoozabadi</年代urn一个米e> <given-names> 一个。</g我ven-names> </name> </person-group> <article-title> 两相流的数值模拟非均匀渗透媒体不同毛细现象的压力</一个rt我cle-title> <source> <italic> 水资源的进步</我t一个lic> <year> 2008年</ygydF4y2Baear> <volume> 31日</vgydF4y2Baolume> <issue> 1</我年代年代ue> <fpage> 56</fp一个ge> <lpage> 73年</gydF4y2Balpage> <pub-id pub-id-type="other"> 2 - s2.0 - 36049041288</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> <pub-id pub-id-type="doi"> 10.1016 / j.advwatres.2007.06.006</pub- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -- - - - - -我d> </nlm-citation> </ref> </ref-list> </back> </article> </body> </html>