US7389185B2ExpiredUtilityA1

Methods and systems for determining reservoir properties of subterranean formations with pre-existing fractures

Assignee: HALLIBURTON ENERGY SERV INCPriority: Oct 7, 2005Filed: Oct 7, 2005Granted: Jun 17, 2008
Est. expiryOct 7, 2025(expired)· nominal 20-yr term from priority
Inventors:David Craig
E21B 49/008
93
PatentIndex Score
52
Cited by
96
References
21
Claims

Abstract

Methods and systems are provided for evaluating subsurface earth oil and gas formations. More particularly, methods and systems are provided for determining reservoir properties such as reservoir transmissibilities and average reservoir pressures of formation layer(s) using quantitative refracture-candidate diagnostic methods. The methods herein may use pressure falloff data from the introduction of an injection fluid at a pressure above the formation fracture pressure to analyze reservoir properties. The model recognizes that a new induced fracture creates additional storage volume in the formation and that a quantitative refracture-candidate diagnostic test in a layer may exhibit variable storage during the pressure falloff, and a change in storage may be observed at hydraulic fracture closure. From the estimated formation properties, the methods may be useful for, among other things, determining whether a pre-existing fracture is damaged and evaluating the effectiveness of a previous fracturing treatment to determine whether a formation requires restimulation.

Claims

exact text as granted — not AI-modified
1. A method for determining a reservoir transmissibility of at least one layer of a subterranean formation having preexisting fractures having a reservoir fluid comprising the steps of:
 (a) isolating the at least one layer of the subterranean formation to be tested; 
 (b) introducing an injection fluid into the at least one layer of the subterranean formation at an injection pressure exceeding the subterranean formation fracture pressure for an injection period; 
 (c) shutting in the wellbore for a shut-in period; 
 (d) measuring pressure falloff data from the subterranean formation during the injection period and during a subsequent shut-in period; and 
 (e) determining quantitatively the reservoir transmissibility of the at least one layer of the subterranean formation by analyzing the pressure falloff data with a quantitative refracture-candidate diagnostic model. 
 
     
     
       2. The method of  claim 1  wherein step (e) is accomplished by transforming the pressure falloff data to equivalent constant-rate pressures and using type curve analysis to match the equivalent constant-rate pressures to a type curve to determine quantitatively the reservoir transmissibility. 
     
     
       3. The method of  claim 1  wherein step (e) is accomplished by:
 transforming the pressure falloff data to obtain equivalent constant-rate pressures; 
 preparing a log-log graph of the equivalent constant-rate pressures versus time; and 
 determine quantitatively the reservoir transmissibility of the at least one layer of the subterranean formation by analyzing the variable-rate pressure falloff data using type-curve analysis according to the quantitative refracture-candidate diagnostic model. 
 
     
     
       4. The method of  claim 2  wherein the reservoir fluid is compressible; and wherein the transforming of the pressure falloff data is based on the properties of the compressible reservoir fluid in the reservoir wherein the transforming step comprises:
 determining a shut-in time relative to the end of the injection period; 
 determining an adjusted time; and 
 determining an adjusted pseudopressure difference. 
 
     
     
       5. The method of  claim 4  wherein the transforming step comprises:
 determining the shut-in time relative to the end of the injection: Δt=t−t ne ; 
 determining the adjusted time: 
 
       
         
           
             
               
                 
                   t 
                   a 
                 
                 = 
                 
                   
                     ( 
                     
                       
                         μ 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         
                           c 
                           t 
                         
                       
                       _ 
                     
                     ) 
                   
                   ⁢ 
                   
                     
                       ∫ 
                       0 
                       
                         Δ 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         t 
                       
                     
                     ⁢ 
                     
                       
                         
                           ⅆ 
                           Δ 
                         
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         t 
                       
                       
                         
                           ( 
                           
                             μ 
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             
                               c 
                               t 
                             
                           
                           ) 
                         
                         w 
                       
                     
                   
                 
               
               ; 
             
           
         
       
       and
 determining the adjusted pseudopressure difference: ΔP a (t)=P aw (t)−P ai  where 
 
       
         
           
             
               
                 
                   p 
                   a 
                 
                 = 
                 
                   
                     
                       
                         
                           μ 
                           _ 
                         
                         g 
                       
                       ⁢ 
                       
                         z 
                         _ 
                       
                     
                     p 
                   
                   ⁢ 
                   
                     
                       ∫ 
                       0 
                       p 
                     
                     ⁢ 
                     
                       
                         p 
                         ⁢ 
                         
                           ⅆ 
                           p 
                         
                       
                       
                         
                           μ 
                           g 
                         
                         ⁢ 
                         z 
                       
                     
                   
                 
               
               ; 
             
           
         
         wherein: 
         t ne  is the time at the end of the injection period; 
           μ  is the viscosity of the reservoir fluid at average reservoir pressure; 
         (μc t ) w  is the viscosity compressibility product of wellbore fluid at time t; 
         (μc t ) 0  is the viscosity compressibility product of wellbore fluid at time t=t ne ; 
         p is the pressure; 
           p  is the average reservoir pressure; 
         P aw (t) is the adjusted pressure at time t; 
         p ai  is the adjusted pressure at time t=t ne ; 
         c t  is the total compressibility; 
           c   t  is the total compressibility at average reservoir pressure; and 
         z is the real gas deviator factor. 
       
     
     
       6. The method of  claim 5  further comprising the step of preparing a log-log graph of a pressure function versus time: I(ΔP a )=f(t a ),
 where 
 
       
         
           
             
               
                 I 
                 ⁡ 
                 
                   ( 
                   
                     Δ 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     
                       p 
                       a 
                     
                   
                   ) 
                 
               
               = 
               
                 
                   ∫ 
                   0 
                   a 
                 
                 ⁢ 
                 
                   Δ 
                   ⁢ 
                   
                       
                   
                   ⁢ 
                   
                     p 
                     a 
                   
                   ⁢ 
                   
                     
                       ⅆ 
                       
                         t 
                         a 
                       
                     
                     . 
                   
                 
               
             
           
         
       
     
     
       7. The method of  claim 5  further comprising the step of preparing a log-log graph of a pressure derivative function versus time: ΔP a ′=f(t a ),
 where 
 
       
         
           
             
               
                 Δ 
                 ⁢ 
                 
                     
                 
                 ⁢ 
                 
                   p 
                   a 
                   ′ 
                 
               
               = 
               
                 
                   
                     ⅆ 
                     
                       ( 
                       
                         Δ 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         
                           p 
                           a 
                         
                       
                       ) 
                     
                   
                   
                     ⅆ 
                     
                       ( 
                       
                         ln 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         
                           t 
                           a 
                         
                       
                       ) 
                     
                   
                 
                 = 
                 
                   Δ 
                   ⁢ 
                   
                       
                   
                   ⁢ 
                   
                     p 
                     a 
                   
                   ⁢ 
                   
                     
                       t 
                       a 
                     
                     . 
                   
                 
               
             
           
         
       
     
     
       8. The method of  claim 2  wherein the reservoir fluid is slightly compressible; and wherein the transforming of the pressure falloff data is based on the properties of the slightly compressible reservoir fluid in the reservoir wherein the transforming step comprise:
 determining a shut-in time relative to the end of the injection period; and 
 determining a pressure difference; 
 wherein: 
 t ne  is the time at the end of the injection period; 
 P w (t) is the pressure at time t; and 
 P i  is the initial pressure at time t=t ne . 
 
     
     
       9. The method of  claim 8  wherein the transforming step comprises:
 determining the shut-in time relative to the end of the injection: Δt=t−t ne ; and 
 determining the pressure difference: ΔP(t)=P w (t)−P i ; and 
 wherein: 
 t ne  is the time at the end of the injection period; 
 P w (t) is the pressure at time t; and 
 P i  is the initial pressure at time t=t ne . 
 
     
     
       10. The method of  claim 8  further comprising the step of plotting a log-log graph of a pressure function versus time: I(Δp)=f(Δt). 
     
     
       11. The method of  claim 9  where 
       
         
           
             
               
                 I 
                 ⁡ 
                 
                   ( 
                   
                     Δ 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     p 
                   
                   ) 
                 
               
               = 
               
                 
                   ∫ 
                   0 
                   
                     Δ 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     t 
                   
                 
                 ⁢ 
                 
                   Δ 
                   ⁢ 
                   
                       
                   
                   ⁢ 
                   p 
                   ⁢ 
                   
                     ⅆ 
                     Δ 
                   
                   ⁢ 
                   
                       
                   
                   ⁢ 
                   t 
                   ⁢ 
                   
                       
                   
                   ⁢ 
                   or 
                   ⁢ 
                   
                       
                   
                   ⁢ 
                   
                     
                       ∫ 
                       0 
                       t 
                     
                     ⁢ 
                     
                       Δ 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       p 
                       ⁢ 
                       
                         
                           ⅆ 
                           t 
                         
                         . 
                       
                     
                   
                 
               
             
           
         
       
     
     
       12. The method of  claim 8  further comprising the step of plotting a log-log graph of a pressure derivatives function versus time: Δp′=f(Δt). 
     
     
       13. The method of  claim 12  where 
       
         
           
             
               
                 Δ 
                 ⁢ 
                 
                     
                 
                 ⁢ 
                 
                   p 
                   ′ 
                 
               
               = 
               
                 
                   
                     ⅆ 
                     
                       ( 
                       
                         Δ 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         p 
                       
                       ) 
                     
                   
                   
                     ⅆ 
                     
                       ( 
                       
                         ln 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         Δ 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         t 
                       
                       ) 
                     
                   
                 
                 = 
                 
                   
                     Δ 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     p 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     Δ 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     t 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     or 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     
                       
                         ⅆ 
                         
                           ( 
                           
                             Δ 
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             p 
                           
                           ) 
                         
                       
                       
                         ⅆ 
                         
                           ( 
                           
                             ln 
                             ⁢ 
                             
                                 
                             
                             ⁢ 
                             t 
                           
                           ) 
                         
                       
                     
                   
                   = 
                   
                     Δ 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     p 
                     ⁢ 
                     
                         
                     
                     ⁢ 
                     
                       t 
                       . 
                     
                   
                 
               
             
           
         
       
     
     
       14. The method of  claim 9  wherein the reservoir transmissibility is determined quantitatively in field units from a before-closure match point as: 
       
         
           
             
               
                 kh 
                 μ 
               
               = 
               
                 141.2 
                 ⁢ 
                 
                   ( 
                   24 
                   ) 
                 
                 ⁢ 
                 
                   
                     p 
                     wsD 
                   
                   ⁡ 
                   
                     ( 
                     0 
                     ) 
                   
                 
                 ⁢ 
                 
                   
                     
                       
                         
                           C 
                           Lfbc 
                         
                         ⁡ 
                         
                           ( 
                           
                             
                               p 
                               0 
                             
                             - 
                             
                               p 
                               i 
                             
                           
                           ) 
                         
                       
                       ⁡ 
                       
                         [ 
                         
                           
                             
                               p 
                               LfbcD 
                             
                             ⁡ 
                             
                               ( 
                               
                                 t 
                                 D 
                               
                               ) 
                             
                           
                           
                             I 
                             ⁡ 
                             
                               ( 
                               
                                 Δ 
                                 ⁢ 
                                 
                                     
                                 
                                 ⁢ 
                                 p 
                               
                               ) 
                             
                           
                         
                         ] 
                       
                     
                     M 
                   
                   . 
                 
               
             
           
         
       
     
     
       15. The method of  claim 9  wherein the reservoir transmissibility is determined quantitatively in field units from an after-closure match point as: 
       
         
           
             
               
                 kh 
                 μ 
               
               = 
               
                 141.2 
                 ⁢ 
                 
                   
                     ( 
                     24 
                     ) 
                   
                   ⁡ 
                   
                     [ 
                     
                       
                         
                           
                             
                               
                                 p 
                                 wsD 
                               
                               ⁡ 
                               
                                 ( 
                                 0 
                                 ) 
                               
                             
                             ⁢ 
                             
                               C 
                               Lfbc 
                             
                           
                         
                       
                       
                         
                           
                             - 
                             
                               
                                 p 
                                 wsD 
                               
                               ( 
                               
                                 
                                   
                                     ( 
                                     
                                       t 
                                       c 
                                     
                                     ) 
                                   
                                   Lfd 
                                 
                                 ⁡ 
                                 
                                   [ 
                                   
                                     
                                       C 
                                       Lfbc 
                                     
                                     - 
                                     
                                       C 
                                       Lfac 
                                     
                                   
                                   ] 
                                 
                               
                             
                           
                         
                       
                     
                     ] 
                   
                 
                 ⁢ 
                 
                   
                     
                       
                         ( 
                         
                           
                             p 
                             0 
                           
                           - 
                           
                             p 
                             i 
                           
                         
                         ) 
                       
                       ⁡ 
                       
                         [ 
                         
                           
                             
                               p 
                               LfacD 
                             
                             ⁡ 
                             
                               ( 
                               
                                 t 
                                 D 
                               
                               ) 
                             
                           
                           
                             I 
                             ⁡ 
                             
                               ( 
                               
                                 Δ 
                                 ⁢ 
                                 
                                     
                                 
                                 ⁢ 
                                 p 
                               
                               ) 
                             
                           
                         
                         ] 
                       
                     
                     M 
                   
                   . 
                 
               
             
           
         
       
     
     
       16. The method of  claim 5  wherein the injection fluid is compressible and contains desirable additives for compatibility with the subterranean formation wherein the reservoir transmissibility is determined quantitatively in field units from a before-closure match point as: 
       
         
           
             
               
                 kh 
                 μ 
               
               = 
               
                 141.2 
                 ⁢ 
                 
                   ( 
                   24 
                   ) 
                 
                 ⁢ 
                 
                   
                     p 
                     awsD 
                   
                   ⁡ 
                   
                     ( 
                     0 
                     ) 
                   
                 
                 ⁢ 
                 
                   
                     
                       
                         
                           C 
                           Lfbc 
                         
                         ⁡ 
                         
                           ( 
                           
                             
                               p 
                               
                                 a 
                                 ⁢ 
                                 
                                     
                                 
                                 ⁢ 
                                 0 
                               
                             
                             - 
                             
                               p 
                               ai 
                             
                           
                           ) 
                         
                       
                       ⁡ 
                       
                         [ 
                         
                           
                             
                               p 
                               LfbcD 
                             
                             ⁡ 
                             
                               ( 
                               
                                 t 
                                 D 
                               
                               ) 
                             
                           
                           
                             I 
                             ⁡ 
                             
                               ( 
                               
                                 Δ 
                                 ⁢ 
                                 
                                     
                                 
                                 ⁢ 
                                 
                                   p 
                                   a 
                                 
                               
                               ) 
                             
                           
                         
                         ] 
                       
                     
                     M 
                   
                   . 
                 
               
             
           
         
       
     
     
       17. The method of  claim 5  wherein the injection fluid is compressible and contains desirable additives for compatibility with the subterranean formation wherein the reservoir transmissibility is determined quantitatively in field units from an after-closure match point as: 
       
         
           
             
               
                 kh 
                 μ 
               
               = 
               
                 141.2 
                 ⁢ 
                 
                   
                     ( 
                     24 
                     ) 
                   
                   ⁡ 
                   
                     [ 
                     
                       
                         
                           
                             
                               
                                 p 
                                 awsD 
                               
                               ⁡ 
                               
                                 ( 
                                 0 
                                 ) 
                               
                             
                             ⁢ 
                             
                               C 
                               Lfbc 
                             
                           
                         
                       
                       
                         
                           
                             - 
                             
                               
                                 
                                   p 
                                   awsD 
                                 
                                 ⁡ 
                                 
                                   ( 
                                   
                                     
                                       ( 
                                       
                                         t 
                                         c 
                                       
                                       ) 
                                     
                                     Lfd 
                                   
                                   ) 
                                 
                               
                               ⁡ 
                               
                                 [ 
                                 
                                   
                                     C 
                                     
                                       
                                           
                                       
                                       ⁢ 
                                       Lfbc 
                                     
                                   
                                   - 
                                   
                                     C 
                                     
                                       
                                           
                                       
                                       ⁢ 
                                       Lfac 
                                     
                                   
                                 
                                 ] 
                               
                             
                           
                         
                       
                     
                     ] 
                   
                 
                 ⁢ 
                 
                   
                     
                       
                         ( 
                         
                           
                             p 
                             
                               a 
                               ⁢ 
                               
                                   
                               
                               ⁢ 
                               0 
                             
                           
                           - 
                           
                             p 
                             ai 
                           
                         
                         ) 
                       
                       ⁡ 
                       
                         [ 
                         
                           
                             
                               p 
                               LfacD 
                             
                             ⁡ 
                             
                               ( 
                               
                                 t 
                                 D 
                               
                               ) 
                             
                           
                           
                             I 
                             ⁡ 
                             
                               ( 
                               
                                 Δ 
                                 ⁢ 
                                 
                                     
                                 
                                 ⁢ 
                                 
                                   p 
                                   a 
                                 
                               
                               ) 
                             
                           
                         
                         ] 
                       
                     
                     M 
                   
                   . 
                 
               
             
           
         
       
     
     
       18. A system for determining a reservoir transmissibility of at least one layer of a subterranean formation by using variable-rate pressure falloff data from the at least one layer of the subterranean formation measured during an injection period and during a subsequent shut-in period, the system comprising:
 a plurality of pressure sensors for measuring pressure falloff data; and 
 a processor operable to transform the pressure falloff data to obtain equivalent constant-rate pressures and to determine quantitatively the reservoir transmissibility of the at least one layer of the subterranean formation by analyzing the variable-rate pressure falloff data using type-curve analysis according to a quantitative refracture-candidate diagnostic model. 
 
     
     
       19. A computer program, stored on a tangible storage medium, for analyzing at least one downhole property, the program comprising executable instructions that cause a computer to:
 determine quantitatively a reservoir transmissibility of the at least one layer of the subterranean formation by analyzing the variable-rate pressure falloff data with a quantitative refracture-candidate diagnostic model. 
 
     
     
       20. The computer program of  claim 19  wherein the determining step is accomplished by transforming the variable-rate pressure falloff data to equivalent constant-rate pressures and using type curve analysis to match the equivalent constant-rate rate pressures to a type curve to determine quantitatively the reservoir transmissibility. 
     
     
       21. The computer program of  claim 19  wherein the determining step is accomplished by transforming the variable-rate pressure falloff data to equivalent constant-rate pressures and using after closure analysis to determine quantitatively the reservoir transmissibility.

Join the waitlist — get patent alerts

Track US7389185B2 — get alerts on status changes and closely related new filings.

We store only your email — no account needed. See our privacy policy.