Page 102 - Percolation Models for Transport in Porous Media With
P. 102

5.1  IMMISCIBLE DISPLACEMENT                                          95

         zone is equal to zero.  Finally, when x 1 > x > x(Vt) we have the dynamic stage of
         the flow.
            Now estimate the values of saturation for  the porous medium in  these zones.
         When x < x(Vt) the fraction of the displacing phase in the V -chain at the moment
         of blocking of the V 1-chain by the Yo-chain is t:  = V /Vt, and therefore the fraction
         of the trapped displaced phase in this chain is  1 - t:.
            After averaging over the region of characteristic size  R(Vo),  using {5.1)-(5.3),
         we have
                                     Vt
                                     I {1  - V/Vt) n(V) B(V) dV
                                     0
                             So(x) = .::..-..,..,.-------                    {5.9)
                                         Vm
                                         I  n(V) B(V) dV
                                         0
         where B(V) = 1 in the pore model, and B(V) = V  in the capillary model [26].
            In  the region  x  > x(V 1 )  trees grow  practically without  interaction.  Around
         each  trunk  a  "crown"  is  formed,  i.e.,  a  zone  where  the  first  phase  has  already
         been displaced from  {Fig.30).  Assuming that crowns have a triangular form  and
         grow around the tree trunks with constant rates,  we find  the distribution of the
         saturation of the medium with the displacing phase in  this region.  The displace-
         ment zone  begins to form  around a V -chain after the length of the tree trunk in
         the course of its growth exceeds the considered coordinate x(Vt ).  The transverse
         dimension of the region is proportional to the quantity

                                   l(V) = Xj(V/Vm)- X,
         which characterizes the distance between the tree top and the plane with coordi-
         nate x.  Therefore the saturation by  the displaced  phase in  the region  x1  > x  >
         x(Vt) equals
                                        Vm
                              S = 1 - {3  J n(V) l(V) B(V) dV              {5.10)
                                      V(x2)
            Here the coefficient  {3  is determined from  the condition of matching solutions
         (5.9)  and {5.10)  for  x  =  x(Vt );  x 2  is  the current coordinate, measured from  the
         level  x  = x(Vt)  (Fig.30),  and V(x 2 )  is  the velocity of the phase interface in  the
         chain,  where it manages to move  up  to the levetx2 •  The velocity is  determined
         from  relationship (5.8)

                                                                            (5.11)

            Thus,  the  system  (5.6)  - (5.11)  allows  to  determine  the  distribution  of the
         saturation of the displaced  phase for  an arbitrary position of the front  x  1, if the
         probability density function Q>(V)  is known.
   97   98   99   100   101   102   103   104   105   106   107