CROSS REFERENCE TO RELATED PROVISIONAL APPLICATION 
   This application is based on Provisional Application Ser. No. 60/231,788 filed on Sep. 12, 2000. 


   BACKGROUND OF THE INVENTION 
   1. Field of the Invention 
   The invention is generally related to methods and processes for analyzing well production data and maximizing efficiency of reservoir production therefrom and is specifically directed to the evaluation of multilayer commingled reservoirs using commingled production data and production logging information. 
   2. Discussion of the Prior Art 
   Field production performance data and multiple pressure transient tests over a period of time for oil and gas wells in geopressured reservoirs have been found to often exhibit marked changes in reservoir effective permeability over the producing life of the wells. Similarly, the use of quantitative fractured well diagnostics to evaluate the production performance of hydraulically fracture wells have clearly shown that effective fracture half-length and conductivity can be dramatically reduced over the producing life of the wells. A thorough investigation of this topic may be found in the paper presented by Bobby D. Poe, the inventor of the subject application, entitled: “Evaluation of Reservoir and Hydraulic Fracture Properties in Geopressure Reservoir,” Society of Petroleum Engineers, SPE 64732. 
   Some of the earliest references to the fact that subterranean reservoirs do not always behave as rigid and non-deformable bodies of porous media may be found in the groundwater literature, see for example, “Compressibility and Elasticity of Artesian Aquifers,” by O. E. Meinzer, Econ. Geol. (1928) 23, 263–271. and “Engineering Hydraulics,” by C. E. Jacob, John Wiley and Sons, Inc. New York (1950) 321–386. 
   The observations of early experimental and numerical studies of the effects of stress-dependent reservoir properties demonstrated that low permeability formations exhibit a proportionally greater reduction in permeability than high permeability formations. The stress-dependence of reservoir permeability and fracture conductivity over the practical producing life of low permeability geopressured reservoirs has resulted in the following observations:
         1. Field evidence of reservoir effective permeability degradation with even short production time can often be observed in geopressured reservoirs.   2. Quantitative evaluation of the field production performance of hydraulic fractures in both normal and geopressured reservoirs have resulted in the observation that the fracture conductivity of hydraulically fractured wells commonly decreases with production time.   3. Multiphase fracture flow has been demonstrated to dramatically reduce the effective conductivity of fractures.   4. Pre-fracture estimates of formation effective permeability derived from pressure transient tests or production analyses are often not representative of the reservoir effective permeability exhibited in the post-fracture production performance.       

   The analysis of production data of wells to determine productivity has been used for almost fifty years in an effort to determine in advance what the response of a well will be to production-stimulation treatment. A discourse on early techniques may be found in the paper presented by R. E. Gladfelter, entitled “Selecting Wells Which Will Respond to Production-Simulation Treatment,” Drilling and Production Procedures, API (American Petroleum Institute), Dallas, Tex., 117–129 (1955). The pressure-transient solution of the diffusivity equation describing oil and gas flow in the reservoir is commonly used, in which the flow rate normalized pressure drops are given by:
 
( P   i   −P   wf )/ q   o , and
 
{ P   p ( P   i )− P   p ( P   wf )}/ q   g ,
 
for oil and gas reservoir analyses, respectively, wherein:
         P i  is the initial reservoir pressure (psia),   P wf  is the sandface flowing pressure (psia)   q o  is the oil flow rate (STB/D)   P p  is the pseudopressure function, psia 2 /cp and   q g  is the gas flow rate (Mcsf/D).       

   While analysis of production data using flow rate normalized pressures and the pressure transient solutions work reasonably well during the infinite-acting radial flow regime of unfractured wells, boundary flow results have indicated that the production normalization follows an exponential trend rather than the logarithmic unit slope exhibited during the pseudosteady state flow regime of the pressure-transient solution. 
   Throughout most of the production history of a well, a terminal pressure is imposed on the operating system, whether it is the separator operating pressure, sales line pressure, or even atmospheric pressure at the stock tank. In any of these cases, the inner boundary condition is a Dirichlet condition (specified terminal pressure). Whether the terminal pressure inner boundary condition is specified at some point in the surface facilities or at the sandface, the inner boundary condition is Dirichlet and the rate-transient solutions are typically used. It is also well known that at late production times the inner boundary condition at the bottom of the well bore is generally more closely approximated with a constant bottomhole flowing pressure rather than a constant rate inner boundary condition. 
   An additional problem that arises in the use of pressure-transient solutions as the basis for the analysis of production data is the quantity of noise inherent in the data. The use of pressure derivative functions to reduce the uniqueness problems associated with production data analysis of fractured wells during the early fracture transient behavior even further magnifies the effects of noise in the data, commonly requiring smoothing of the derivatives necessary at the least or making the data uninterpretable at the worst. 
   There have been numerous attempts to develop more meaningful production data analyses in an effort to maximize the production level of fractured wells. One such example is shown and described in U.S. Pat. No. 5,960,369 issued to B. H. Samaroo, describing a production profile predictor method for a well having more than one completion wherein the process is applied to each completion provided that the well can produce from any of a plurality of zones or in the event of multiple zone production, the production is commingled. 
   From the foregoing, it can be determined that production of fractured wells could be enhanced if production performance could be properly utilized to determine fracture efficiency. However, to date no reliable method for generating meaningful data has been devised. The examples of the prior art are at best speculative and have produced unpredictable and inaccurate results. 
   SUMMARY OF THE INVENTION 
   The subject invention is a method of and process for evaluating reservoir intrinsic properties, such as reservoir effective permeability, radial flow steady-state skin effect, reservoir drainage area, and dual porosity reservoir parameters omega (dimensionless fissure to total system storativity) and lambda (matrix to fissure crossflow parameter) of the individual unfractured reservoir layers in a multilayer commingled reservoir system using commingled reservoir production data, such as wellhead flowing pressures, temperatures and flow rates and/or cumulatives of the oil, gas, and water phases, and production log information (or pressure gauge and spinner survey measurements). The method and process of the invention also permit the evaluation of the hydraulic fracture properties of the fractured reservoir layers in the commingled multilayer system, i.e., the effective fracture half-length, effective fracture conductivity, permeability anisotropy, reservoir drainage area, and the dual porosity reservoir parameters omega and lambda. The effects of multiphase and non-Darcy fracture flow are also considered in the analysis of fractured reservoir layers. 
   The subject invention is directed to a method of and process for fractured well diagnostics for production data analysis for providing production optimization of reservoir completions via available production analysis and production logging data. The method of the invention is a quantitative analysis procedure for reservoir and fracture properties using commingled reservoir production data, production logs and radial flow and fractured interval analyses. This permits the in situ determination of reservoir and fracture properties for permitting proper and optimum treatment placement and design of the reservoir. The invention provides a rigorous analysis procedure for multilayer commingled reservoir production performance. Production logging data is used to correctly allocate production to each completed interval and defined reservoir zone. This improves the stimulation and completion design and identifies zones to improve stimulation. 
   The subject invention is a computational method and procedure for computing the individual zone production histories of a commingled multi-layered reservoir. The data used in the analysis are the commingled well production data, the wellhead flowing temperatures and pressures, the complete wellbore and tubular goods description, and production log information. This data is used to construct the equivalent individual layer production histories. The computed individual completed interval production histories that are generated are the individual layer hydrocarbon liquid, gas, and water flow rates and cumulative production values, and the mid-completed interval wellbore flowing pressures as a function of time. These individual completed interval production histories can then be evaluated as simply drawdown transients to obtain reliable estimates of the in situ reservoir effective permeability, drainage area, apparent radial flow steady-state skin effect and the effective hydraulic fracture properties, namely, half-length and conductivity. 
   Typically, an initial production log is run soon after a well is put on production and the completion fluids have been produced back from the formation. Depending on the formation, the stimulation/completion operations performed on the well and the size and productive capacity of the reservoir, a second production log is run after a measurable amount of stabilized production has been obtained from the well. Usually, additional production logs are run at periodic intervals to monitor how the layer flow contributions and wellbore pressures continue to vary with respect to production time. The use of production logs in this manner provides the only viable means of interpreting commingled reservoir production performance without the use of permanent downhole instrumentation. 
   The subject invention is directed to the development of a computational model that performs the production allocation of the individual completed intervals in a commingled reservoir system using the fractional flow rates of the individual completed intervals, determined from production logs and the commingled system total well fluid phase flow rates. The individual completed interval flow rate histories generated include the individual completed interval fluid phase flow rates and cumulative production values as a function of production time, as well as the mid-zone wellbore flowing pressures. The computed mid-zone flowing wellbore pressures at the production time levels of the production log runs are then compared with the actual measured wellbore pressures at those depths and time level to ascertain which wellbore pressure traverse model most closely matches the measured pressures. 
   The identified wellbore pressure traverse model is then used to model the bottom hole wellbore flowing pressures for all of the rest of the production time levels for which there are not production log measurements available. This use of the identified pressure traverse model to generate the unmeasured wellbore flowing pressure is the only assumption required in the entire analysis. It is fundamentally sound unless there are dramatic changes in the character of the produced well fluids or in the stimulation/damage of the completed intervals which is not reflected in the composite production log history, primarily due to inadequate sampling of the changes in the completed intervals producing fractional flow rates. With an adequate sampling of the changing fractional flow rate contributions of the individual completed intervals in a commingled reservoir, this analysis technique is superior to other multi-layer testing and analysis procedures. 
   The method and process of the subject invention provide a fully-coupled commingled reservoir system analysis model for allocating the commingled system production data to the individual completed intervals in the well and constructing wellbore flowing pressure histories for the individual completed intervals in the well. No assumptions are required to be made as to the stimulation/damage steady-state skin effect, effective permeability (or formation conductivity), initial pore pressure level, drainage area extent, or intrinsic formation properties of the completed intervals in a commingled reservoir system. The method of the invention considers only the actual measured response of the commingled system using production logs and industry accepted wellbore pressure traverse computational models. 
   The fundamental basis for the invention is a computationally rigorous technique of computing the wellbore pressure traverses to the midpoints (or other desired points) of each completed interval using one or more of a number of petroleum industry accepted wellbore pressure traverse computational methods in combination with the wellbore tubular configuration and geometry, wellbore deviation survey information, completed interval depths and perforation information, wellhead measured production rates (or cumulatives) and the wellhead pressures and temperatures of the commingled multilayer reservoir system performance. The computed pressure traverse wellbore pressures are compared with the measured wellbore pressures of either a production log or a wellbore pressure survey. This permits the identification of the pressure traverse computational method that results in the best agreement with the physical measurements made. 
   The invention permits the use of information from multiple production logs run at various periods of time over the producing life of the well. The invention also permits the specification of crossflow between the commingled system reservoir layers in the wellbore. The invention evaluates the pressure traverse in each wellbore segment using the fluid flow rates in that wellbore section, the wellbore pressure at the top of that wellbore section, and the temperature and fluid density distributions in that section of the wellbore traverse. The method and process of the invention actually uses downhole physical measurements of the wellbore flowing pressures, temperatures, fluid densities, and the individual reservoir layer flow contributions to accurately determine the production histories of each of the individual layers in a commingled multilayer reservoir system. The results of the analysis of the individual reservoir layers can be used with the commingled reservoir algorithm to reconstruct a synthetic production log to match with the actual recorded production logs that are measured in the well. The invention has an automatic Levenberg-Marquardt non-linear minimization procedure that can be used to invert these production history records to determine the individual completed interval fracture and reservoir properties. The invention also has the option to automatically re-evaluate the initially specified unfractured completed intervals that indicate negative radial flow steady-state skin effects as finite-conductivity vertically fractured completed intervals. 
   The method and process of the subject invention permits for the first time a reliable, accurate, verifiable computationally rigorous analysis of the production performance of a well completed in a multilayer commingled reservoir system using physically measured wellbore flow rates, pressures, temperatures, and fluid densities from the production logs or spinner surveys and pressure gauges to accomplish the allocation of the flow rates in each of the completed reservoir intervals. The combination of the production log information and the wellbore traverse calculation procedures results in a reliable, accurate continuous representation of the wellbore pressure histories of each of the completed intervals in a multilayer commingled reservoir system. The results may then be used in quantitative analyses to identify unstimulated, under-stimulated, or simply poorly performing completed intervals in the wellbore that can be stimulated or otherwise re-worked to improve productivity. The invention may include a full reservoir and wellbore fluids PVT (Pressure-Volume-Temperature) analysis module. 


   
     BRIEF DESCRIPTION OF THE DRAWINGS 
       FIG. 1  is a flow chart of the process of the subject invention. 
       FIG. 2  is an illustration of the systematic and sequential computational procedure in accordance with the subject invention. 
   


   DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
   The subject invention is directed to a computational model for computing the wellbore pressure traverses and individual layer production contributions of the individual completed intervals in a commingled reservoir. Direct physical measurements of the individual layer flow contributions to the total well production and the actual wellbore flowing pressures are recorded and included in the analysis. There are numerous wellbore pressure traverse models available for computing the bottom hole flowing and static wellbore pressures from surface pressures, temperatures and flow rates, as will be well known to those skilled in the art. The selection of the appropriate pressure traverse model is determined by comparison with the actual wellbore pressure measurements. In a commingled reservoir the layer fractional flow contribution to the total well production rate also commonly varies with respect to time. There are many factors that govern the individual layer contributions to the total well production rate with respect to time. Among these are differences in the layer initial pressures, effective permeability, stimulation or damage steady-state skin effect, drainage area, net pay thickness, and the diffusivity and storativity of the different layers. Other factors that are not directly reservoir-controlled that affect the contribution of each of the layers to the commingled reservoir well production are the changing wellbore pressures, completion losses and changing gas and liquid produced fluid ratios with respect to time. 
   Production logs (PLs) provide a direct means of measuring the wellbore flowing pressures, temperatures, and actual reservoir layer flow contributions at specific points in time, with which to calibrate the computed pressure traverse models. It is preferable to run multiple production logs on wells producing commingled reservoirs to track the variation in the individual completed interval contributions with respect to production time. 
   It is known that the commingled system total production rate commonly does not equal or even come close to equaling the sum of the individual completed interval isolated flow rates when each interval is tested in isolation from the other completed intervals in the well. There are several factors causing this, including but not limited to (1) invariably higher flowing wellbore pressures present in the commingled system across each of the completed intervals than when they were measured individually, and (2) possible crossflow between the completed intervals. 
   As more particularly shown in the flowchart of  FIG. 1 , the subject invention is directed to a computational model that performs the production allocation of the individual completed intervals in a commingled reservoir system using the fractional flow rates of the individual completed intervals, determined from the production logs and the commingled system total well fluid phase flow rates. This depicts the analysis process for a reservoir with three completed reservoir layers in which the upper and lower reservoir layers have been hydraulically fractured. The middle reservoir completed interval has not been fracture stimulated. The wellbore pressure traverse is computed using the total well commingled production flow rates to the midpoint of the top completed interval. Then the fluid flow rates in the wellbore between the midpoint of the top and middle completed intervals are evaluated using the total fluid phase flow rates of the commingled system minus the flow rates from the top completed interval. The pressure traverse in the wellbore between the midpoints of the middle and lower completed intervals is evaluated using the fluid phase flow rates that are the difference between the commingled system total fluid phase flow rates and the sum of the phase flow rates from the top and middle completed intervals. The individual completed interval flow rate histories generated in this analysis include the individual completed interval fluid flow rates and cumulative production values as a function of production time, as well as the mid-zone wellbore flowing pressures. The computed mid-zone flowing wellbore pressures at the production time levels of the production log runs are then compared with the actual measured wellbore pressures at those depths and time level to ascertain which wellbore pressure traverse model most closely matches the measured pressures. 
   The identified wellbore pressure traverse model is then used to model the bottomhole wellbore flowing pressure for all of the rest of the production time levels for which there are not production log measurements available. This use of the identified pressure traverse model to generate the unmeasured wellbore flowing pressures is the only major assumption made in the process. It is fundamentally sound unless there are dramatic changes in the character of the produced well fluids or in the stimulation/damage of the completed intervals which is not reflected in composite production log history, primarily due to inadequate sampling of the changes in the completed intervals producing fractional flow rates. With an adequate sampling of the changing fractional flow rate contributions of the individual completed intervals in a commingled reservoir, this analysis technique produces accurate results. 
     FIG. 2  is an illustration of the systematic and sequential computational procedure in accordance with the subject invention. Beginning at the wellhead  10 , the pressure traverses to the midpoint of each completed interval are computed in a sequential manner. The fluid flow rates in each successively deeper segment of the wellbore are decreased from the previous wellbore segment by the production from the completed intervals above that segment of the wellbore. The mathematical relationships that describe the fluid phase flow rates (into or out) of each of the completed intervals in the wellbore are given as follows for oil, gas, and water production of the j th  completed interval, respectively:
   q   oj ( t )= q   ot ( t ) f   oj ( t ),   q   gi ( t )= q   gt ( t ) f   gj ( t ),   q   wj ( t )= q   wt ( t ) f   wj ( t ), 
where:
         q oj  is the j th  completed interval hydrocarbon liquid flow rate, STB/D,   q of  is the composite system hydicarbon liquid flow rate, STB/D,   f oj  is the j th  completed interval hydrocarbon liquid flow rate liquid contribution of the total well hydrocarbon liquid flow rate, fraction,   q gi  is the j th  interval gas flow rate, Mcsf/D   j is the index of completed intervals,   q gt  is the composite system total well gas flow rate, Mscf/D,   f gj  is the j th  completed interval gas flow rate fraction of total well gas flow rate, fraction,   q wj  is the j th  interval water flow rate, STB/D   q wt  is the composite system total well water flow rate, STB/D   f wj  is the j th  completed interval water flow rate fraction of total well water flow rate, fraction.       
   The corresponding fluid phase flow rates in each segment of the weilbore are also defined mathematically with the relationships as follows for oil, gas and water for the n th  wellbore pressure traverse segment, respectively. 
   
     
       
         
           
             
               
                 
                   
                     q 
                     
                       o 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       n 
                     
                   
                   ⁡ 
                   
                     ( 
                     t 
                     ) 
                   
                 
                 = 
                   
                 ⁢ 
                 
                   
                     
                       q 
                       
                         o 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         t 
                       
                     
                     ⁡ 
                     
                       ( 
                       t 
                       ) 
                     
                   
                   - 
                   
                     
                       ∑ 
                       
                         
                           
                             
                               j 
                               = 
                               1 
                             
                           
                         
                         
                           
                             
                               n 
                               &gt; 
                               1 
                             
                           
                         
                       
                       
                         n 
                         - 
                         1 
                       
                     
                     ⁢ 
                     
                       
                         q 
                         oj 
                       
                       ⁡ 
                       
                         ( 
                         t 
                         ) 
                       
                     
                   
                 
               
             
           
           
             
               
                 
                   
                     q 
                     
                       g 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       n 
                     
                   
                   ⁡ 
                   
                     ( 
                     t 
                     ) 
                   
                 
                 = 
                   
                 ⁢ 
                 
                   
                     
                       q 
                       
                         g 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         t 
                       
                     
                     ⁡ 
                     
                       ( 
                       t 
                       ) 
                     
                   
                   - 
                   
                     
                       ∑ 
                       
                         
                           
                             
                               j 
                               = 
                               1 
                             
                           
                         
                         
                           
                             
                               n 
                               &gt; 
                               1 
                             
                           
                         
                       
                       
                         n 
                         - 
                         1 
                       
                     
                     ⁢ 
                     
                       
                         q 
                         gj 
                       
                       ⁡ 
                       
                         ( 
                         t 
                         ) 
                       
                     
                   
                 
               
             
           
           
             
               
                 
                   
                     q 
                     
                       w 
                       ⁢ 
                       
                           
                       
                       ⁢ 
                       n 
                     
                   
                   ⁡ 
                   
                     ( 
                     t 
                     ) 
                   
                 
                 = 
                   
                 ⁢ 
                 
                   
                     
                       q 
                       
                         w 
                         ⁢ 
                         
                             
                         
                         ⁢ 
                         t 
                       
                     
                     ⁡ 
                     
                       ( 
                       t 
                       ) 
                     
                   
                   - 
                   
                     
                       ∑ 
                       
                         
                           
                             
                               j 
                               = 
                               1 
                             
                           
                         
                         
                           
                             
                               n 
                               &gt; 
                               1 
                             
                           
                         
                       
                       
                         n 
                         - 
                         1 
                       
                     
                     ⁢ 
                     
                       
                         q 
                         wj 
                       
                       ⁡ 
                       
                         ( 
                         t 
                         ) 
                       
                     
                   
                 
               
             
           
         
       
     
   
   The flow rate and pressure traverse computations are performed in a sequential manner for each wellbore segment, starting at the surface or wellhead  10  and ending with the deepest completed interval in the wellbore, for both production and injection scenarios. The wellbore flow rate and pressure traverse calculation procedures employed permit the evaluation of production, injection or shut in wells. 
   The fundamental inflow relationships that govern the transient performance of a commingled multi-layered reservoir are fully honored in the analysis provided by the method of the subject invention. Assuming that accurate production logs are run in a well, when a spinner passes a completed interval without a decrease in wellbore flow rate (comparing wellbore flow rates at the top and bottom of the completed interval, higher or equal flow rate at the top than at the bottom), no fluid is entering the interval from the wellbore (no loss to the completed interval, i.e., no crossflow). Secondly, once the minimum threshold wellbore fluid flow rate is achieved to obtain stable and accurate spinner operation, all higher flow rate measurements are also accurate. Lastly, the sum of all of the completed interval contributions equals the commingled system production flow rates for both production and injection wells. 
   In the preferred embodiment of the invention, two ASCII input data files are used for the analysis. One file is the analysis control file that contains the variable values for defining how the analysis is to be performed (which fluid property and pressure traverse correlations are uses, as well as the wellbore geometry and production log information). The other file contains commingled system wellhead flowing pressures and temperatures, and either the individual fluid phase flow rates or cumulative production values as a function of production time. 
   Upon execution of the analysis two output files are generated. The general output file contains all of the input data specified for the analysis, the intermediate computational results, and the individual completed interval and defined reservoir unit production histories. The dump file contains only the tabular output results for the defined reservoir units that are ready to be imported and used in quantitative analysis models. 
   The analysis control file contains a large number of analysis control parameters that use can be used to tailor the production allocation analysis to match most commonly encountered wellbore and reservoir conditions.
A method for providing production optimization of reservoir completions having a plurality of completed intervals via available production analysis and production logging data provides a quantitative analysis procedure for reservoir and fracture properties of a commingled reservoir system, that includes the steps of measuring pressure for specific zones in a reservoir; selecting a pressure traverse model; computing midzone pressures using the traverse model; comparing the computed midzone pressures with the measured pressures; and modeling the bottomhole pressure of the reservoir based on the traverse model.