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

7.2  FRACTURE LENGTH DISTRIBUTION                                   135

            Another change of variables, t =sin f/J/d.,  in  the equation (7.12) reduces it to

                                          w/2
                                  /(d.)= IF (si~f/J) df/J                   (7.13)
                                          0
            The equation  (7.13)  is  one  of the  SchlOmilch  integral  equations.  The solu-
         tion of this equation can be obtained as follows.  After introducing the function
         F 1(sin¢J/d.)  =  F(d./sinf/J)  into  the  equation  (7.13)  and  then  differentiating  it
         with respect to d., we obtain
                                         w/2

                             d~f'(ds) =-IF{ (sinf/J/d.)sinf/Jdf/J           (7.14)
                                         0
            By substituting xf sin 1/J fords in (7.14), integrating it with respect to 1/J from 0
         to 1rj2, and changing the order of integration in the right side, the equation (7.14)
         can be reduced to
                w/2                      w/2   w/2

                I   _j:_ !' (~) d,P  =  2 I  df/J  IF,' (sin,Psinf/J) sinf/Jd'I/J   (7.15)
                  SID 2 1/J
                           SID 1/J
                                                  l
                                                        X
                0                        0    0
            Having introduced a  new  variable sin A =  sin f/Jsin 1/J,  after changing the order
         of integration in the right side, we obtain
              w/2                      w/2                  w/2
                 .x:  !' (~) d,P =-IF{ (sin A) cosAdA I
                                                                  sinf/Jdf/J
              I   SID  1/J   SID 1/J           x              J cos2 A - cos2 ¢J
              0                        0                    .\
            Taking account of the fact that
                                  w/2
                                         sin ¢J df/J   1r
                                  I   Jcos2A- cos2¢J  =  2'
                                  .\
         after integrating with respect to A,  we obtain the resultant expression

                                           w/2
                          F(rt) = F(oo)- ~ J -4:- f' (  .rt.,,) d,P         (7.16)
                                         7r   SID  1/J   SID 'I'
                                           0
            Formula (7.16)  permits to find  the length distribution function  F(rt)  of disk
         fractures if the length distribution function  /(d.) of fracture traces is known.
            In particular, if the distribution  f(d 8 )  is defined as an exponential relation

                                    f(d.) =  aexp(-cd.)
   136   137   138   139   140   141   142   143   144   145   146