# Patent

## Title
Method for reconstructing the total organic carbon content from compositional modeling analysis

## Abstract

Methods are provided for utilizing the results of compositional modeling analysis to obtain accurate total organic carbon values without the need for an oxidation step or lengthy sample preparation, and also to calculate the organic carbon value attributable to contaminants, such as drilling additives.

## Background
FIELD OF THE INVENTION
The present invention relates to a method for reconstructing the total organic carbon (TOC) values in geological samples by a combination of pyrolysis, chemical analysis, and compositional modeling in order to determine the quality of unconventional oil resources or source rocks for use during development drilling and/or exploration operations.
BACKGROUND OF THE INVENTION
The use of open system pyrolysis, as used in the commercial Rock EvalTM analysis, and total organic carbon measurement employed by the Rock EvalTM and LECOTM analysis systems are methods that have long been used in the assessment of petroleum source rocks. See, e.g., Peters, K. E., 1986, Guidelines for Evaluating Petroleum Source Rock Using Programmed Pyrolysis, Bulletin of the American Association of Petroleum Geologists, v. 70, p. 318-329; U.S. Pat. No. 5,811,308, Expitalie et al., Method for Determining Petroleum Characteristics of Geologic Sediments Sep. 22, 1998; Langford, F. F. and M.-M. Blank-Valleron, 1990, Interpreting Rock-Eval Pyrolysis Data Using Graphs of Pyrolizable Hydrocarbons vs. Total Organic Carbon, Bulletin of the American Association of Petroleum Geologists, v. 74, p. 799-804; Lafargue, E., J. Expitalié, F. Marquis, and D. Pillot, 2000, Rock-Eval 6 Applications in Hydrocarbon Exploration, Production, and in Soil Contamination Studies: Revue de L'institut Francais du Parole, Vol. 53, No, 4, p. 421-437. These methods, however, rely on bulk measurements of all organic matter present in a sample and can only provide information from empirically-derived cutoff values and parameters. Methods and systems described in U.S. Pat. No. 7,363,206 and WO 2008/100614 (PCT/US2008/002102) demonstrate that pyrolysis data can be used to characterize the relative amounts of organic matter and contaminants in geological samples.
The Pyrolytic Oil-Productivity Index (POPI) technology originally disclosed in U.S. Pat. No. 5,866,814 has been used successfully in exploration and development wells to quickly and accurately assess the reservoir rock for tar mats and other types of organic matter as the drilling progresses through these layers in the reservoir rock. The POPI method has been incorporated in the commercially available software program identified as GC-ROXTM which is an acronym for “Geochemical Residual Oil eXpert-modeling”. The GC-ROXTM program can be used to optimize oil field development by tar mat identification and quantification, and can be used to validate, organize and store field-data. Earlier observations from various field data had suggested a cut-off of productivity exists at about 3% tar. However, relying on the more precise information produced by the GC-ROX software, a target of 1% tar has been used commercially for injectivity purposes. Interest has been shown in linking the POPI method and information to viscosity. It is thought that proxies could be developed, but the main application of the technique is identification and quantification of tar.
The process of petroleum source rock maturation and hydrocarbon generation is well understood, with the effects of burial and temperature causing kerogen to break down into bitumen, free oil that is expelled and migrates, and eventually into gas. Methods using compositional modeling are disclosed in U.S. Pat. No. 7,363,206 which permit the user to identify the relative proportions by percent of a plurality of end member components in a sample of reservoir rock. This methodology represented a significant advance in the art, since prior methods only provided bulk values or parameters.
The compositional modeling method for assessing residual hydrocarbon staining that is taught in U.S. Pat. No. 7,363,206 at column 6, line 56 to column 9, line 48, is preferred for use in the present invention. As used herein and specifically in the claims, the term “compositional modeling” shall mean the method disclosed in U.S. Pat. No. 7,363,206 and as described below. The disclosure of U.S. Pat. No. 7,363,206 is incorporated herein by reference.
These factors are: (1) the amount of total hydrocarbon yield, and (2) the similarity of the hydrocarbon staining to the produced oils. Pyrolysis instruments are useful for quantifying the amount of hydrocarbon staining and the POPI method assesses the similarity to produced oils by subdividing the hydrocarbons into the Light Volatile (LV), Thermally Distillable (TD), and Thermally Crackable (TC) components (FIG. 1). However, it has been discovered that visual inspection of pyrograms also can be useful in assessing the type of hydrocarbons present because oil, tar, pyrobitumen, and other typical organic matter each also have very characteristic appearances.
FIGS. 2a through 2d are examples of pyrograms for samples with a nearly uniform composition of specific hydrocarbon or pyrolytically identifiable type of organic matter end-member components. These plots show the hydrocarbon yield on the y-axis for each data step that is recorded on the x-axis.
The number of data steps for a particular analysis can vary based on the type of pyrolyis instrument used. Two such commercially available instruments are Vinci's ROCK-EVALTM and Humble's SOURCE ROCK ANALYZERTM. The data that is obtained from the instrument and converted into a digital file can also vary. In the case description and examples that follow, the SOURCE ROCK ANALYZERTM was used and the data were output into digital form using a MICROSOFT EXCELTM CSV file that recorded the yield and temperature over 611 data steps. The first 111 steps record the isothermal hold at 195 degrees C. for 3 minutes and the next 500 steps record the programmed temperature run from 195 degrees C. to 630 degrees C. In general, the temperature associated with any specific step is the same from run to run, so that the step number can be associated with the temperature of the oven during the run.
Compositional modeling for a sample includes entering in an appropriately programmed computer the yield at each individual data step as a value that is made up of the aggregate yield of the various end-member components. In this method, a specific and consistent temperature is associated with each step. The difference between the modeled yield calculated in accordance with the algorithm from the theoretical end-member component or components and the actual yield provides the basis for assessing whether a particular solution accurately reflects the actual composition in the reservoir rock sample.
Each such solution that is assessed must sum the difference between the calculated yield and the actual yield over all the data steps for the sample. Any of a number of statistical methods can be used in quantifying the overall error for any proposed solution. The modeling relies on iteratively varying the concentration of the various components until the aggregate error is minimized and the curves appear very similar. In one preferred embodiment, the software utilizes the iterative process of proposing different compositions, calculating a hypothetical curve based on the yield at each data step, assessing the error for each particular solution, and then minimizing this aggregate error. These method steps can advantageously be completed by the use of macros and the Solver add-on application that is a standard component of the Microsoft ExcelTM program, and its use greatly automates and expedites the process. There are other software packages that can also be utilized to facilitate the methods used to model hydrocarbon composition that are commercially available and include Corel Quattro Pro, Lotus 1-2-3, Corel Paradox, Lotus Approach, Microsoft Access and Microsoft Visual FoxPro.
Table 1, shown in FIG. 6, is the output file from a Humble SOURCE ROCK ANALYZERTM that has been converted to CSV format. The output data file records the same information in the same location for each sample tested and facilitates its extraction by spreadsheet data analysis programs. The header information at the very top of the report records the calculated parameters from the instrument and the run parameters. Starting at row 22, the instrument records the curve signal in the first three columns of the file. The first column contains the data step number, the second column records the signal from the flame ionization detector (FID) in milliVolts (mV), and the third column records the temperature of the oven associated with the data step.
In order to convert the output of the instrument into hydrocarbon yield using the instrument software, a known standard compound or composition from the reservoir region is analyzed. With the data from the standard, the instrument can calculate the conversion factor (CF) to relate millivolts from the FID to hydrocarbon yield in units of milligrams per gram of rock (g Rock). From the data file, these conversion factors are calculated for each sample by summing the total signal in column two and then dividing this signal by the total hydrocarbon yield of the sample in accordance with the following mathematical expression:

CFFID=[ΣSignalstep 1-611 (mV)]/[(LV+TD+TC) (mgHC/g Rock)]  (1)

The signal that is taken for any particular data step is then be converted into mgHC/gRock by simply dividing the signal by CF:

Yieldstep X (mgHC/g Rock)=[SignalSTEP X (mV)]/[CFFID (mV/mgHC/g ROCK)]  (2)

In a preferred embodiment, all instrument output is converted into yields for the purpose of making the relevant calculations to combine end-members (EM) and to compute the results based on modeled solutions. This is done because the actual yield that is given for each end-member sample and for each sample that is being modeled will be unique. In order to model the relative composition of end-members that make up a particular sample, the data for each end-member is normalized so that the total hydrocarbon yields of each recalculated end-member is the same as the actual sample. Therefore, the quantity of an end-member component that would be present for a pure end-member having the same yield as the sample can be expressed as follows:

Yield=[YieldEMSTEPX]*([Total Yield(THC)Sample]/[Total Yield End-Member])  (3)

In the above equation (3) and those that follow the “*” notation is used to indicate multiplication.

Equation (3) is used to calculate the aggregate yield that would be found for a hypothetical sample that contained various percentages of end-members that are needed to describe the sample behavior. Thus, the calculated yield for a proposed hydrocarbon composition at any given data step is the sum of the percentages of each end-member (% EM1 to 5) divided by 100 and multiplied by the yield of the end-member at step x (YieldEM1 to 5,X) times the ratio of the total yield of the sample divided by the total yield of the end-member (THCsample/THCEM1 to 5). This step can be expressed as follows:




  
    
      
        
          Calculated
          ⁢
          
            
          
          ⁢
          
            Yield
            
              
                
                  EM
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
                -
                5
              
              ,
              X
            
          
        
        =
        
          
            
              (
              
                %
                ⁢
                
                  
                
                ⁢
                
                  EM
                  1
                
                ⁢
                
                  /
                
                ⁢
                100
              
              )
            
            *
            
              Yield
              
                
                  EM
                  ⁢
                  
                    
                  
                  ⁢
                  1
                
                ,
                X
              
            
            *
            
              (
              
                
                  THC
                  sample
                
                ⁢
                
                  /
                
                ⁢
                
                  THC
                  
                    EM
                    ⁢
                    
                      
                    
                    ⁢
                    1
                  
                
              
              )
            
          
          +
          
            
              (
              
                %
                ⁢
                
                  
                
                ⁢
                
                  EM
                  2
                
                ⁢
                
                  /
                
                ⁢
                100
              
              )
            
            *
            
              Yield
              
                
                  EM
                  ⁢
                  
                    
                  
                  ⁢
                  2
                
                ,
                X
              
            
            *
            
              (
              
                
                  THC
                  sample
                
                ⁢
                
                  /
                
                ⁢
                
                  THC
                  
                    EM
                    ⁢
                    
                      
                    
                    ⁢
                    2
                  
                
              
              )
            
          
          +
          
            
              (
              
                %
                ⁢
                
                  
                
                ⁢
                
                  M
                  3
                
                ⁢
                
                  /
                
                ⁢
                100
              
              )
            
            *
            
              Yield
              
                
                  EM
                  ⁢
                  
                    
                  
                  ⁢
                  3
                
                ,
                X
              
            
            *
            
              (
              
                
                  THC
                  sample
                
                ⁢
                
                  /
                
                ⁢
                
                  THC
                  
                    EM
                    ⁢
                    
                      
                    
                    ⁢
                    3
                  
                
              
              )
            
          
          +
          
            
              (
              
                %
                ⁢
                
                  
                
                ⁢
                
                  EM
                  4
                
                ⁢
                
                  /
                
                ⁢
                100
              
              )
            
            *
            
              Yield
              
                
                  EM
                  ⁢
                  
                    
                  
                  ⁢
                  4
                
                ,
                X
              
            
            *
            
              (
              
                
                  THC
                  sample
                
                ⁢
                
                  /
                
                ⁢
                
                  THC
                  
                    EM
                    ⁢
                    
                      
                    
                    ⁢
                    4
                  
                
              
              )
            
          
          +
          
            
              (
              
                %
                ⁢
                
                  
                
                ⁢
                
                  
                    EM
                    5
                  
                  /
                  100
                
              
              )
            
            *
            
              Yield
              
                
                  EM
                  ⁢
                  
                    
                  
                  ⁢
                  5
                
                ,
                X
              
            
            *
            
              (
              
                
                  THC
                  sample
                
                ⁢
                
                  /
                
                ⁢
                
                  THC
                  
                    EM
                    ⁢
                    
                      
                    
                    ⁢
                    5
                  
                
              
              )
            
          
        
      
    
    
      
        (
        4
        )
      
    
  




The error between a particular modeled solution for Step X and the actual analytical result for Step X can be obtained by simple difference. However, since some of these values will be positive and some negative, the treatment of errors for all values calculated, e.g., Steps 1-611, is easier to accomplish through the application of a root mean squares (RMS) calculation. Other statistical treatments can be used to also achieve the same results if they employ the difference between each modeled yield and actual yield as an absolute value.
In one preferred method, the RMS average difference is calculated in terms of a percentage that relates to the total response of the sample and can be expressed as follows:

%RMSCALC vs ACTUAL=100*((AVERAGESTEP, 1-611(YieldCALC−YieldACTUAL)2)1/2/(AVERAGESTEP, 1-611(YieldACTUAL))  (5)

The modeling process comprises the steps of varying the percentage of the end-members that are present in the system in which (EM1-5 are preferably used) until the calculated curve matches the actual curve and the % RMS error is minimized. Due to the fact that so many calculations must be made to assess any solution, the use of a spreadsheet program to perform these calculations and automatically plot the result that is achieved greatly simplifies the process. Moreover, a software application such as Solver that is present as an add-on in Microsoft ExcelTM, can greatly expedite the data processing capability of iteratively solving problems with multiple variables that seek to converge on a desired solution which in this case is minimizing error.
FIG. 3 shows the graphic interface that was utilized for a five end-member component system at well-site Z. The graphic illustration includes curves for the current sample, the calculated solution based on the percentage of the components, the oil end-member, the tar end-member, the shale end-member (typical of dispersed kerogen found in shaley lithologies), the coal end-member, and the drilling mud end-member (contamination). The parameter listed on the top line as DEVRM is the RMS deviation as a percent of total yield and is the value that is minimized in obtaining a reasonable solution for a given sample. When all samples are analyzed for a particular well, the results can be plotted as in FIG. 4 to reveal how the composition varies throughout the sampled section. The numerical references 1-4 on the graph are to the legends at the top. Plots such as these are very useful in identifying important trends, such as increasing tar, or in identifying individual coal or tar units that may have important implications in reservoir performance. An alternative method of presenting the data as is shown in FIG. 5, is to plot the relative contribution of each EM component by depth with each curve being adjusted for changes in yield in the samples. This type of plot is particularly useful for identifying true tar mats that typically have an associated dramatic increase in hydrocarbon yield as opposed to a change in composition that appears to be tar, but is present in relatively low concentrations and not likely to affect reservoir performance.
As used herein, it will be understood from the description that “compositional modeling” includes the following steps:

    
    
        a. identifying the end-member components known to be present in reservoir rock in the oil field;
        b. preparing individual pyrograms consisting of a pcd for each of the components identified in step (a);
        c. storing the pcd for each component in a digital data file;
        d. conducting a pyrolytic analysis of a sample from the oil field of reservoir rock that contains one or more hydrocarbon and organic matter components of the type identified in step (a) to obtain pcd for the sample;
        e. comparing the pcd for the sample with the pcd obtained in step (b) for each of the components, measuring and recording the difference between the sample pcd and the pcd for each of the components;
        f. applying a statistical analysis to minimize the aggregate differences between the pcd for the sample and a combination of pcd selected from the components;
        g. recording for retention and display for analysis the resulting pcd that constitutes the minimum aggregate error over the temperature range of the pyrolytic analysis; and
        h. analyzing the displayed data to identify the respective components.
    
    


As used herein, “unconventional oil” and “unconventional oil sources” include oil shales; organic-rich fine-grained carbonates, low-porosity low-permeability sandstones/siltstones/carbonates that are adjacent to hydrocarbon source units; oil sands, and heavy crude oils. Unconventional oil reservoirs are often source rocks that are being exploited to produce hydrocarbons that were unable to migrate as in conventional hydrocarbon reservoirs. One problem faced by geoscientists is how to classify the organic matter and hydrocarbon components at specific depths and/or intervals in the well, since it is recognized that such information, if available, would be extremely useful in assessing the reservoir during exploration and development drilling.
Knowledge of the quantity and types of organic matter and hydrocarbons present in unconventional reservoirs is critical to well performance evaluation and geochemical methods are used to provide this information. Traditional petrophysical analysis focuses principally on the rock matrix; however, the storage and flow capacity of unconventional reservoirs is highly dependent on the organic porosity resulting from the transformation of kerogen. Although total organic carbon (TOC) values are important to reservoir evaluation and assessment, traditional TOC measurement is not available at the well site in real-time in order to meaningfully assist in optimizing drilling operations.
The prior art methods require either time-consuming sample processing, such as the demineralization needed for the LECOTM analyzer, e.g., the LECO (EC 12) carbon analyzer; and preparation of samples using the TOC measurement or they require long analytical times by more complex instrumentation that is not suitable for use at well drilling sites where the rock samples are recovered from the drilling fluid. In addition, the prior art methods are bulk analytical methods which do not differentiate between the materials present in mixtures.
A new method that can be deployed at the well site to provide information that can be used to guide the drilling operation for exploratory and development wells is needed.
Definitions
As used herein, it is to be understood that the following terms and designations have the meanings indicated:
Bitumen—generally refers to the organic material that is extractable through the use of organic solvents and encompasses those obtained through the use of strong organic solvents such as methylene chloride, chloroform, and toluene.
End Member (EM)—A consistent type of organic matter or hydrocarbon that can be distinguished by pyrolytic analysis. End members include oil, soluble tar, pyrobitumen (insoluble tar), kerogen, coal, drilling mud and other contaminants. Specific end members are associated with specific fields and reservoirs. A reference to “local” end members means end members that have been determined to be present in the field and/or nearby wells based on analysis and examination of core samples, well logs, and drilling rock samples, also referred to below as “comparative samples”.
EMx Weight OM—the weight attributable to the organic matter (in milligrams per gram of rock) of one component end member (x) in a sample.
EMx Weight NSOs OM—the weight attributable to the elemental nitrogen, sulfur, and oxygen content (in milligrams per gram of rock) of one component end member (x) in a sample.
EMx Weight H—the weight attributable to the elemental hydrogen content (in milligrams per gram of rock) of one component end member (x) in a sample.
FID—Flame Ionization Detector
Free Oil—hydrocarbons that have been generated through the maturation of kerogen that are similar to producible oil and may be expelled from a source rock given a sufficient saturation level.
HC—Abbreviation for hydrocarbons, THC is used for Total Hydrocarbons.
H/COM—the ratio of hydrogen to carbon in the organic matter of an end member component.
Inert Carbon—Organic Carbon, that can be oxidized to CO2 and CO by application of heat as in traditional source rock analysis in the oxidation cycle or by LECO carbon analyzer, but that cannot be liberated by pyrolysis in an inert atmosphere either due to character of the organic matter or the adsorptive properties of the mineral matrix.
Kerogen—the product of biochemical degradation, polycondensation, insolubilization of biologically-derived material due to burial, time and temperature.
LV—Abbreviation for light volatile components. As used herein, LV refers specifically to the weight in milligrams of HC released per gram of rock at the initial static temperature condition of 180° C. (when the crucible containing the rock sample is inserted into the pyrolytic chamber) prior to the temperature-programmed pyrolysis of the sample.
Non-Pyrolizable Hydrocarbons (HCNon-py)—compounds composed primarily of carbon and hydrogen that cannot be liberated through evaporation or decomposition of organic material during pyrolysis and require combustion to analyze the different elemental components.
NOSx%—the weight percent of Nitrogen, Oxygen, and Sulfur in the organic matter of an end member component.
Organic Matter (OM)—generally biologically produced materials that are composed primarily of carbon and hydrogen, such as shale- and coal-like materials.
OM/HCpy—the ratio of organic matter to pyrolizable hydrocarbons in an end member component.
POPI—Abbreviation for the Pyrolytic Oil-Productivity Index. The POPI is calculated from the pyrolytic data by the following equation: POPI=ln(LV+TD+TC)×(TD/TC), where ln is the logarithmic value and TD and TC are as defined below.
Pyrolizable Hydrocarbons (HCpy)—compounds composed primarily of carbon and hydrogen that can be liberated through evaporation or decomposition of organic material in response to the application of heat in an inert atmosphere (i.e., not containing oxygen).
Pyrolytic Characterization Data (pcd)—Data values measured at a predetermined number of data points, each data point corresponding to a prescribed temperature.
Rock-EvalTM S1 Yield—weight in milligrams of HC released per gram of rock at the initial static temperature condition of 300° C. of a standard Rock-EvalTM pyrolysis analysis as described in Peters, K. E., 1986, Guidelines for Evaluating Petroleum Source Rock Using Programmed Pyrolysis, Bulletin of the American Association of Petroleum Geologists, v. 70, p. 318-329.
Rock-EvalTM S2 Yield—weight in milligrams of HC released per gram of rock during the programmed pyrolysis portion of a standard Rock-Eval analysis, where the temperature is raised from 300° C. to 550-600° C. at a rate of 25° C./minute, as described in Peters, K. E., 1986, Guidelines for Evaluating Petroleum Source Rock Using Programmed Pyrolysis, Bulletin of the American Association of Petroleum Geologists, v. 70, p. 318-329.
Source Rock—A sedimentary rock, deposited in a low-energy environment, generally of fine-grained nature (i.e., composed of silt and clay-sized particles), typically fissile, commonly composed of silica, clay, and carbonate minerals, and with sufficient organic matter content and quality to generate and expel hydrocarbons (i.e., petroleum).
TD—Abbreviation for “thermally distillable” components that, as used herein, refers specifically to the weight in milligrams of HC released per gram of rock at a temperature between 180° C. (195° C. on a Humble SR Analyzer) and Tmin(° C.).
TC—Abbreviation for “thermally crackable” components that as used herein, refers specifically to the weight in milligrams of HC released per gram of rock at a temperature between Tmin(° C.) and 600° C. (630° C. on a Humble SR Analyzer).
THC—Abbreviation for total hydrocarbons.
Tmin—The temperature at which the hydrocarbon yield, as measured by a flame ionization detector (FID) during pyrolysis employing the POPI method, reaches a minima between the peaks representing the thermally distilled and thermally cracked hydrocarbon peaks, generally occurring between 380° C. and 420° C.
Total Hydrocarbon Index (THI)—Represents the total HC released, including during the initial heating and programmed pyrolysis from 195° C. and 630° C., relative to Total Organic Carbon in a sample. The equation for calculating THI is: THI=[(LV+TD+TC)/TOC]×100.
THIOM—the Total Hydrocarbon Index of organic matter in an end member component.
Total Organic Carbon (TOC)—The TOC is the weight percent of organic carbon found in a rock sample.
TOCEMx [%]—The total organic carbon in weight percent that is attributable to a specific end member. This is found by taking total weight of the organic matter represented by the end member (EMx Weight OM), subtracting the weight of the end member attributable to elemental Nitrogen, Sulfur, and Oxygen (EMx Weight NSOs), subtracting the weight of the end member attributable to elemental hydrogen (EMx Weight H), then converting it to a decimal by dividing by 1000 mg/g Rock, and finally multiplying by 100 to put the number into percent form.
TOCinert—The portion of total organic carbon in weight percent that cannot be decomposed and detected by pyrolysis, but requires measurement of carbon dioxide and monoxide produced during an oxidation step at high temperature.
TOCRCN—The reconstructed total organic carbon, which is not directly measured, but rather by summing the TOC for each end member and the inert carbon present in the sample (also known as, TOCCoMod) as follows: TOCRCN=TOCEM1+TOCEM2 . . . . TOCEMx+TOCinert.
Wt % HOM—the weight percent of Hydrogen in the organic matter as determined by CHNOS elemental analysis.
YieldEMx—the hydrocarbon yield (in milligrams of hydrocarbon per gram of rock) that is attributable to each end member from applying compositional modeling as described in U.S. Pat. No. 7,363,206.
SUMMARY OF THE INVENTION
In accordance with the present invention, a method is provided for utilizing the results from compositional modeling analysis to provide accurate total organic carbon (TOC) data for assessment without the need for either an oxidation step or lengthy sample preparation as required by LECO TOC apparatus and method. The present method utilizes nonobvious procedures to assess the level of inert carbon present in geological reservoir rock samples and identifies and quantifies differences in the amount of pyrolizable hydrocarbons related to the geological maturity of the sample. The method of the invention also accounts for the presence of organic carbon attributable to known contaminants, such as drilling additives, which can include a large variety.
The associated methods of the present invention are unique because they allow the determination of total organic carbon in geological samples without time-consuming steps for sample preparation, do not subject the sample to combustion, and only require a flame ionization detector (FID). The methods provide superior assessment of the types of organic matter present in samples and also indicate the presence of mobile hydrocarbons.
Although other workers have noted the contribution of different materials present in geological samples to be analyzed for total organic carbon, all of these prior art methods rely on simple bulk measurements of the materials. Because the method of the present invention provides a detailed assessment of the organic materials found in a rock sample, a unique understanding of the organic matter transformation process is also achieved by this new method. None of the prior art methods can provide a detailed picture in the form of a graphic display of the organic matter contributing to the total organic carbon value from a rock sample.
The method includes the following steps:

    
    
        a. collecting a plurality of representative rock samples typical of unconventional oil reservoirs and source rocks;
        b. performing pyrolytic analysis on the rock samples in accordance with the POPI method(s), as described in U.S. Pat. No. 5,866,814;
        c. obtaining and reviewing the pyrolysis data for samples from a reservoir to assess the likely end-member components that appear to be present;
        d. selecting samples that appear to have a nearly uniform composition of a particular end-member, that is, those appearing to be composed mostly of one member, e.g., “free oil”, bitumen, kerogen, or the like, based upon an inspection of the pyrolysis curve or by a visual inspection of the sample by experienced technical or engineering personnel;
        e. performing separation from the rock sample of any moveable hydrocarbon components by extraction with a non-polar solvent, e.g., cyclohexane, and saving both the solvent extract and the extracted rock;
        f. performing separation of soluble tar from the extracted rock sample by extraction with a strong organic solvent (e.g., methylene chloride), saving both the solvent extract and the extracted rock;
        g. performing separation of residual organic matter on a portion of the strong organic solvent-extracted rock sample by demineralization methods;
        h. analyzing all extracts and organic matter separations by elemental analysis to determine the percent of C, H, N, O and S in each sample;
        i. analyzing all solvent extracts, extracted rock samples, and separated organic matter by the pyrolysis-TOC method to determine the total hydrocarbon index (THI) in accordance with prior art method(s), as described in WO 2008/100614 or US 2010/0057409 A1;
        j. analyzing the pyrolysis data using the prior art compositional modeling method, as described in U.S. Pat. No. 7,363,206;
        k. recording the results for the percentage of total yield for end-members for the group of samples;
        l. employing prior art methods, as described in WO 2008/100614 or US 2010/0057409 A1 to obtain parameters for the total hydrogen index (THIOM), the ratio of organic matter to pyrolizable hydrocarbons (OM/HCpy), the ratio of hydrogen to carbon in the organic matter (H/COM), the percent of Nitrogen, Oxygen, and Sulfur in the organic matter (NOSx%), and the weight of Hydrogen in the organic matter (Wt % HOM);
        m. analyzing a plurality of samples from a well that has a source rock or unconventional reservoir for which the reconstructed TOC values are needed;
        n. applying the compositional modeling method described in U.S. Pat. No. 7,363,206 to obtain the hydrocarbon yield attributable to each end member (YieldEMx) in a group of samples from the source rock or unconventional reservoir;
        o. determining the weight of organic matter represented by the hydrocarbon yield by applying the equation: EMx Weight OM=YieldEMx*OM/HCpy;
        p. determining the weight of elemental nitrogen, sulfur and oxygen by the equation: EMx Weight NSOs=NOSx%*EMx Weight OM;
        q. determining weight of Hydrogen for each end member in a sample by the equation: EMx Weight H=Wt % HOM*EMx Weight OM;
        r. calculating the total organic carbon for the end member (TOCEMx) as follows: TOCEMx=(EMx Weight OM−EMx Weight NSOs−EMx Weight H)×100/1000 mg/g Rock; and
        s. summing the TOC for each end member and the inert carbon present in the formation to obtain the value of reconstructed total organic carbon (TOCRCN or TOCCoMod) as follows: TOCRCN=TOCEM1+TOCEM2 . . . TOCEMx+TOCinert.
    
    


The reconstruction of the TOC measurement from compositional modeling in accordance with the above method requires significantly less analytical time to produce data for decision making than the methods of the prior art. Consequently, the results are produced with greater frequency to assist with the identification of potentially productive reservoir intervals. Also, because the individual organic materials present within the rock sample can be differentiated, including organic contaminants from the drilling mud, the results can be corrected to include only TOC that is indigenous to the rocks that have been produced during the drilling of the well.



BRIEF DESCRIPTION OF THE DRAWINGS
The invention will be described in further detail below and with reference to the attached drawings in which:
FIGS. 1-5 are representative of prior art graphic plots, and specifically:
FIG. 1 is a typical output pyrogram from an instrument performing open-system temperature programmed pyrolysis;
FIG. 2a is a typical pyrogram of API 30° oil from the K reservoir region;
FIG. 2b is a typical pyrogram of tar from the K reservoir region;
FIG. 2c is a typical pyrogram of coaly organic matter from the K reservoir region;
FIG. 2d is a typical pyrogram of organic-rich shale from the K reservoir region;
FIG. 3 is a plot showing the compositional modeling interface that was developed to perform the calculations, where the plot shows pyrograms associated with the current sample, the calculated solution and the end-members that are used in the modeling process;
FIG. 4 is the resultant compositional modeling for well Z-95 showing the percentage of oil, tar, shaley OM and coal present in reservoir rock samples;
FIG. 5 is the plot resulting from the compositional modeling of well Z-95 showing the relative pyrolytic yield of oil, tar, shaley OM and coal present in reservoir rock samples and demonstrates how hydrocarbon yield increases significantly in true tar mats or coal beds;
FIG. 6, which is referred to in the specification as Table 1, is an illustrative representative reproduction for a hypothetical sample of the output file of a commercial pyrolysis instrument;
FIG. 7 is a plot of S2 hydrocarbon yield versus TOC for the nearby “AMAD-1” well.
FIG. 8 is a plot of the end member pyrograms used for modeling a well identified as “DBYT-1”;
FIG. 9 shows the results from modeling a typical pyrogram from the same DBYT-1 well;
FIG. 10 is a plot comparing the determination of TOC from modeling POPI data to measured TOC from instrumental data;
FIG. 11 is a composite log generated by the GC-ROX software for the DBYT-1 well;
FIGS. 12A-B, which is referred to in the specification as Table 2, is a calculation of TOCRCN for well DBYT-1;
FIGS. 13A-B, which is referred to in the specification as Table 3, is a calculation of TOCRCN for well Amad_2;
FIG. 14 is a cross-plot of instrument TOC versus reconstructed TOC (TOCRCN) from the DBYT-1 well;
FIG. 15 is a cross-plot of instrument TOC versus reconstructed TOC (TOCRCN) from a well identified as “AMAD-2”;
FIG. 16 is a plot of the pyrogram for a representative sample from well DBYT-1 produced during standard Rock-Eval analysis;
FIG. 17 is the corresponding plot for a representative sample from well MZLJ-29;
FIGS. 18 and 19 are the corresponding plots of the samples of FIGS. 14 and 15 when assessed by the method of the invention; and
FIG. 20 is a composite log of the data obtained from the method of the invention that corresponds to FIG. 11 from the DBYT-1 well.



DETAILED DESCRIPTION OF THE INVENTION
In applying the method for reconstructing the TOC from compositional modeling analysis, rock samples from a well are collected and analyzed in accordance with the methods described in U.S. Pat. Nos. 5,866,814, 6,823,298 and 7,363,206, the disclosures of which are incorporated herein by reference. In the use of the method of U.S. Pat. No. 7,363,206, it is necessary to identify the local end members that are suitable for use in applying compositional modeling. This is accomplished either by analyzing samples previously obtained from wells in the same field where the zone of interest is located as described below, or by relying exclusively on the samples collected from the well that is being assessed during exploratory or developmental drilling. Preferably, samples from other nearby similar wells are obtained. If there is no other analogue that is appropriate, then samples from the well being drilled could be used. However, this would provide less confidence in any interpretations made. If the samples are from the well being assessed, there can still be a need for real-time results.
Identifying appropriate end members can be accomplished by the examination of preexisting core samples, drill cuttings and electric logs that have been collected for the field or nearby wells. These sources will provide geologists and petroleum engineers with knowledge of the local oil fields and a basis for identifying the presence of organic-rich shales that are source rocks or unconventional reservoirs. Thus, typical steps in finding appropriate local organic end members would include collecting rock samples from zones where such organic matter occurrences have been noted in the equivalent geologic intervals from reservoir regions in proximity to the well that is being assessed.
Initial screening of data samples prepared and subjected to the pyrolytic oil-productivity index (POPI) method analysis as disclosed in U.S. Pat. No. 5,866,814 includes comparing resultant pyrograms to identify the similarity of curve shapes and serves to guide the user as to which samples are to be utilized for the organic matter separations and analysis steps that follow.
The basic pyrolysis analytical method as taught in U.S. Pat. No. 5,866,814 at column 3, line 51 to column 4, line 33, is preferred for use in the present invention. As used herein and specifically in the claims, the term “pyrolysis” shall mean the method disclosed in U.S. Pat. No. 5,866,814 and as described below. The disclosure of U.S. Pat. No. 5,866,814 is incorporated herein by reference.
In the pyrolytic analysis process, small samples (e.g., ≦100 mg) of powdered rock are placed in a steel crucible. The crucible is placed in a furnace and the sample is heated in a stream of helium gas to an initial temperature of 180 degrees C. After heating at 180 degrees C. for about three minutes, the temperature is increased. The rate of increase in the temperature is about 25 degrees C./min. or less, and preferably about 10 degrees C./min, and progresses from 180 degrees C. to about 600 degrees C.
The helium gas carries hydrocarbon products released from the rock sample in the furnace to a detector which is sensitive to organic compounds. During the process, the following three types of events occur:

    
    
        1. Hydrocarbons that can be volatilized at or below 180.degree. C. are desorbed and detected while the temperature is held constant during the first 3 minutes of the procedure. These are called light volatile hydrocarbons (LVHC or LV).
        2. At temperatures between 180 degrees C. and about 400 degrees C., thermal desorption of solvent-extractable bitumen, or the light oil fraction, occurs. These are called thermally distilled hydrocarbons or “distillables” (TDHC or TD).
        3. At temperatures above about 400 degrees C., pyrolysis (cracking) of heavier hydrocarbons, or asphaltenes, occurs. The materials that thermally crack are called thermally cracked hydrocarbons or “pyrolyzables” (TCHC or TC).
    
    


These events give rise to three ‘peaks’ on the initial instrument output (referred to as a pyrogram). The peak for the static 180 degrees C. temperature is a standard output parameter of either the Vinci or Humble instruments. It is referred to as either S1 or volatile total petroleum hydrocarbons (VTPH), respectively. In the present invention, the value will be referred to as light volatiles, LV, as described above. Data generated from the temperature-programmed pyrolysis portion of the procedure is reprocessed manually by the operator to determine the quantity of hydrocarbons in milligrams per gram of sample above and below Tmin. This reprocessing is a trivial exercise for an experienced operator and can be accomplished routinely with either the Vinci or Humble instruments. The first peak above 180 degrees C. represents the amount of thermally distillable hydrocarbons in the sample and is referred to as TD, the second peak above 180 degrees represents the amount of pyrolyzables or thermally “cracked” hydrocarbons in the sample and is referred to as TC. In the case of lighter hydrocarbons or the analysis of oil samples directly for calibration, Tmin may not be discernable. In this case, if the sample analysis is repeatable at 400 degrees C., the values of LV, TD, and TC employed in the method of the present invention are with respect to the specific temperature ranges defined above.
The methods described in published application US 2010/0057409 are applied to identify the end members that are potentially present in the samples through screening and the analysis of consistent curve shapes. The samples are then subjected to physical separation of these components by solvent extractions to assess the uniformity of the organic matter in the samples. Cyclohexane is preferably used to separate any mobile, or “free” hydrocarbons present in the sample that could be produced from the well. Other saturated alkane solvents can also be used to provide similar results, but it has been found that cyclohexane produces an extract that often matches fluids actually produced from hydrocarbon-bearing zones.
After extraction of the component characteristic of mobile hydrocarbons, any remaining soluble organic matter is extracted. The soluble material found in the source rock or an unconventional reservoir will typically have a high asphaltene component and is representative of bitumen, which is formed early in the process of kerogen degradation during source rock maturation. A preferred solvent for bituminous compounds is dichloromethane, which is a strong organic solvent that can solubilize most asphaltenes and heavy hydrocarbon components. However, any other strong or highly polar organic solvent such as chloroform, carbon disulfide, and the like, can be used so long as it extracts all or most of the remaining soluble organic matter. The remaining insoluble organic matter in the sample after solvent extraction of the source rock can be representative of several materials, such as kerogen, inert carbon and recycled and oxidized organic matter.
After appropriate end members have been identified, the methods described in published application no. WO 2008/100614 (PCT/US2008/002102) are applied to determine the values for the Total Hydrogen Index (THIOM), the ratio of organic matter-to-pyrolizable hydrocarbons (OM/HCpy), the ratio of hydrogen-to-carbon in the organic matter (H/COM), the percent of nitrogen, oxygen, and sulfur in the organic matter (NOSx%), and the weight of hydrogen in the organic matter (Wt % HOM).
These steps as taught in US 2010/0057409 at paragraphs [0062] to [0075], is preferred for use in the present invention and as described below. The disclosure of US 2010/0057409 is incorporated herein by reference.
In the assessment of reservoir rock samples, the goal is to characterize all of the hydrocarbons, whether they are “free” hydrocarbons or those bound in a complex structure. In addition, the methods differ in that the analytical procedure utilizes a starting temperature that is much lower (180° to 195° C.). Nonetheless, for each end-member, e.g., oil, tar and pyrobitumen, the amount of hydrocarbon by pyrolysis per gram of carbon will be fairly consistent. This parameter is referred to as the Total Hydrocarbon Index (THI) and is calculated as follows:

THI=[(LV+TD+TC)/TOC]×100  (6)

The units for THI are the same as HI, i.e., mg of hydrocarbon per gram of organic carbon. The Rock-Eval 6 or Humble POPI/TOC analyzers can be used to assess the differences in hydrogen for various end-members. THI is a ratio, with both TOC and THC (LV+TD+TC) determined during the same analysis, thus the errors associated with isolation of the OM, weighing, and small sample sizes that can occur in elemental analysis do not affect the data. It is important that sufficient separations of the organic end-members be obtained, so that results are consistent. Additional information that is needed for each organic matter end-member are the results from elemental analysis, i.e., as would be obtained from a CHNOS analyzer. With the percentage amounts of these elements and THI, the average amount of hydrogen present in the pyrolyzable and non-pyrolyzable portions of end-members can be determined in order to calculate the weight of the overall material. As stated above, the presence of heteroatoms in the various end-members and incomplete pyrolysis does result in some under-reporting of the quantity of hydrocarbon structural units based on the FID response. However, these effects are treated as contributing to the same result when assessing the pyrolyzable versus non-pyrolyzable portion of organic matter.
A typical separation of organic material in reservoir rock can be accomplished through a series of extraction/analytical steps. It is desirable to analyze the produced oil, which can be measured by placing a few microliters of oil on silica gel in a crucible and drying it in an oven at 30° C. for about 6 hours to remove the volatile components. The results of this analysis, plus CHNOS on the oil, provide a means for calculating the mass of oil per gram of rock contained in a sample. Another way of assessing the oil fraction is to extract the rock sample with cyclohexane. The extract that is obtained from this procedure will typically resemble the moveable hydrocarbons in the reservoir and the subsequent measurement of THI and CHNOS will provide data that are suitable for VOM analysis. In addition to the measurement of parameters for the extract, the THI and CHNOS data are also measured on the rock extracted with cyclohexane. This data represents the “immovable” bitumen remaining on the rock after cyclohexane extraction.
The next step in the process is to perform a second extraction on the cyclohexane extracted rock using a strong polar solvent, such as methylene chloride. The resulting extract is considered typical of the remaining tar or asphaltene component of the hydrocarbon staining, while the residual organic matter left in the rock consists of either pyrobitumen (tar that has been altered, lost a significant portion of its hydrogen, and become insoluble even in strong organic solvents) or another insoluble material like coal or kerogen. A variety of mechanisms have been proposed for the formation of pyrobitumen, including thermochemical sulfate reduction (TSR) and thermal alteration. However, there is no single accepted pathway. Nonetheless, hydrogen content in pyrobitumen is lower and the yield of hydrocarbon per gram of carbon is also lower. As with the cyclohexane extract and remaining rock, the methylene chloride extract and methylene chloride extracted rock are both analyzed by pyrolysis to determine the THI, and by elemental analysis for the CHNOS composition of the samples.
Determining Pyrolytic Yield and Mass/Volume Relationship for End-Members
The use of pyrolysis data to determine the volume of various organic constituents in the residual hydrocarbon staining is based on the determination of the Total Hydrocarbon Index for various reservoir organic matter (OM) types and the results from the previously described compositional modeling (CoMod) method. Detailed analysis of the end-members present in a reservoir is only needed for a limited number of samples in order to develop a relationship between the weight of the hydrocarbon component and the total weight of an end-member. Table 1 illustrates in tabular form the calculation steps required to develop the conversion factors for oil, tar and pyrobitumen and their application to determine the volume of end-members in conjunction with CoMod results. Since THI is based on the mg of hydrocarbon per gram of TOC, the composition of the end-members is also given relative to one gram of TOC. In addition, since pyrolysis instruments assess the weight of hydrocarbon in terms of mg per gram of rock, the amount of rock matrix assumed in this calculation is also 1 gram.
In this example, the results from Rock Eval 6 for this reservoir provide a THI of 1050, 525, and 250 respectively for the oil, tar, and pyrobitumen end-members. The results from elemental CHNOS analysis indicate that the end-members have H/COM ratios of 1.9, 1.05, and 0.65, respectively, for oil, tar and pyrobitumen. Because the calculations are based on 1 gram of carbon, the weight of hydrogen in the organic matter based on elemental analysis in the sample can be calculated as follows:

Wt. HOM (mg/1 g TOC)=H/COM×(1000 mg C/Mol.Wt.Carbon)  (7)

As stated previously, utilization of the combined pyrolysis and FID method results in the under-reporting of the weight-percent of hydrogen, because some of the hydrogen is associated with the pyrolyzable OM and some is associated with the non-pyrolyzable OM. Because the aliphatic bonds in complex molecules like asphaltenes and kerogen are most readily broken, it is assumed that the stoichiometry of aliphatic structural units, i.e., —CnH2n, are dominant in the pyrolyzable portion of the organic matter. This results in an average weight-percent for hydrogen in the pyrolyzable hydrocarbons (% H CnH2n) of 14.3%. Applying this assumption, the weight of hydrogen in the pyrolyzable hydrocarbons is determined as follows:

Wt. HHCPy (mg/1 g TOC)=%H CnH2n/100×THI(mg HC/g TOC)×1 g TOC  (8)

The weight of hydrogen in non-pyrolyzable OM is determined as follows:

Wt. HNon-Py(mg/1 g TOC)=Wt. HOM(mg/1 g TOC)−Wt. HHCPy(mg/1 g TOC)  (9)

In accordance with the method of the present invention, the weights of hydrogen in the non-pyrolyzable portion of oil, tar and pyrobitumen for this example are 8.2 (mg/1 g TOC), 12.4 (mg/1 g TOC), and 18.4 (mg/1 g TOC), respectively. Thus, the assumption that the hydrocarbon liberated from organic matter by pyrolysis is dominated by aliphatic units result in an under-reporting of hydrogen that is relatively large for tar (14.2%) and even larger for pyrobitumen (34%). Moreover, assuming that the average composition of the measured hydrocarbons from pyrolysis-FID has a lower H/C ratio, would only increase the hydrogen under-reporting from pyrolysis, which would necessitate a larger correction.
In order to simplify the approach, in the method of the invention the —CnH2n— structural units are assumed to be the dominant form that is produced by pyrolysis of organic matter from petroleum reservoirs. However, it should be noted that the error involved in misrepresenting the amount of hydrogen in the sample is relatively small. For example, the 18.4 mg H/1 g TOC under-reporting for hydrogen in pyrobitumen would result in only a 1.75% error if totally disregarded. When the purpose of the calculations is to assess the volume of tar in a reservoir, and the critical amount of tar present is around 5% of the volume, it will be understood that the difference between 4.9% and 5.1% tar by volume is not significant.
Once the respective weights for hydrogen in the sample have been determined, the weight of carbon in the pyrolyzable OM and non-pyrolyzable OM can be calculated as follows:

Wt. CHCPy (mg/1 g TOC)=THI (mg/1 g TOC)×1 g TOC−Wt.HHCPy (mg/1 g TOC)  (10)

Wt. CNon-Py(mg/1 g TOC)=1000 mg TOC−Wt. CHCPy (mg/1 g TOC)  (11)

The elemental composition analysis of crude oils shows that elemental NSO typically comprise between 1-4% and asphaltene fractions will typically range from 5-12%. In the present example, the NSO values determined for the oil, tar and pyrobitumen were 2.8%, 7.3%, and 7.6%, respectively, which is within the range noted in the literature. For example, see Ancheyta, et al., Energy and Fuels, Vol. 16, pp 1121-27, 2002; Holleran, VSS Technology Library, Valley Slurry Seal Company, www.slurry.com/techpapers contrbit.shtml, 2000. Thus, from CHNOS elemental analysis, the weight-percent attributed to elemental nitrogen, sulfur and oxygen in the sample can be readily determined. Once these values have been determined, the weight of this material in the organic end-member can be calculated as follows:




  
    
      
        
          Wt
          .
          
            NSO
            OM
          
        
        =
        
          
            
              
                
                  
                    (
                    
                      %
                      ⁢
                      
                        
                      
                      ⁢
                      
                        NSO
                        OM
                      
                      ⁢
                      
                        /
                      
                      ⁢
                      100
                    
                    )
                  
                  ×
                
              
            
            
              
                
                  (
                  
                    
                      Wt
                      .
                      
                        C
                        
                          HCP
                          y
                        
                      
                    
                    +
                  
                
              
            
            
              
                
                  
                    
                      Wt
                      .
                      
                        
                      
                      ⁢
                      
                        C
                        
                          Non
                          -
                          Py
                        
                      
                    
                    +
                    
                      Wt
                      .
                      
                        
                      
                      ⁢
                      
                        H
                        
                          HCP
                          y
                        
                      
                    
                    +
                    
                      Wt
                      .
                      
                        H
                        
                          Non
                          -
                          Py
                        
                      
                    
                  
                  )
                
              
            
          
          
            (
            
              1
              -
              
                (
                
                  %
                  ⁢
                  
                    NSO
                    OM
                  
                  ⁢
                  
                    /
                  
                  ⁢
                  100
                
                )
              
            
            )
          
        
      
    
    
      
        (
        12
        )
      
    
  




The total weight of organic matter for each end-member can be calculated relative to 1 gram of TOC and the ratio of organic matter to pyrolyzable hydrocarbon (OM/HCPy) can be determined as follows:

Wt. OM (mg/1 g TOC)=Wt. HHCPy+Wt. HNon-Py+Wt. CHCPy+Wt. CNon-Py+Wt.NSOOM  (13)

OM/HCPy=Wt. OM/THI, or  (14)

OM/HCPy=Wt. OM/(Wt. HHCPy+Wt. CHCPy)  (15)

As with the weight of hydrogen in the sample, variations in the amount of elemental NSO in organic matter within a reservoir are not likely to affect the estimated volume of organic matter by more than a few percent. It is important to determine suitable end-members that produce accurate CoMod results. Since OM/HCPy ratios for the example shown were 1.13 for oil, 2.23 for tar, and 4.56 for pyrobitumen, errors in the modeling process are far more important than minor errors associated with the chemical make-up of these materials. Therefore, it is particularly important to confirm modeled results with laboratory separations that show that the relative amount of soluble versus insoluble materials is similar. This is because a 5% volume of pyrobitumen has a response that is equal to about 2.5% volume of tar.
Finally, the quantity of inert carbon (TOCinert) in the samples is determined by plotting TOC versus hydrocarbon yield in accordance with well known and commonly applied methods, such as those described by Langford, F. F. and M.-M. Blank-Valleron, 1990, Interpreting Rock-Eval Pyrolysis Data Using Graphs of Pyrolizable Hydrocarbons vs. Total Organic Carbon, Bulletin of the American Association of Petroleum Geologists, v. 74, p. 799-804; Dahl, B., J. Bojesen-Koefoed, A. Holm, H. Justwan, E. Rasmussen, and E. Thomsen, 2004, A New Approach to Interpreting Rock-Eval S2 and TOC Data for Kerogen Quality Assessment, Organic Geochemistry, v. 35, pp. 1461-1477, and others.
With these parameters and other data entered in the memory of the appropriately programmed general purpose computer, the results from applying compositional modeling as described in U.S. Pat. No. 7,363,206 for each end member (EM) are used to calculate the hydrocarbon yield that is attributable to each end member (YieldEMx) in a group of samples from a source rock or unconventional reservoir. From this data, the weight of organic matter represented by each sample's hydrocarbon yield can be calculated by utilizing the ratio of the organic matter to pyrolizable hydrocarbons (OM/HCpyx) as follows:

EMx Weight OM=YieldEMx*OM/HCpyx  (16)

The weight of elemental nitrogen, sulfur, and oxygen is calculated as follows:

EMx Weight NSOs=NOSx%*EMx Weight OM  (17)

Next, the weight of Hydrogen is calculated as follows:

EMx Weight H=Wt % HOM*EMx Weight OM  (18)

Finally, the weight of total organic carbon for the end member (TOCEMx) is calculated as follows:

TOCEMx (wt. %)=(EMx Weight OM (mg/g Rock)−EMx Weight NSOs (mg/g Rock)−EMx Weight H (mg/g Rock))×100/1000 mg/g Rock  (19)

After the TOC for each end member component is calculated, they are summed and the inert carbon present in the formation as determined from a nearby well using standard methods is added to yield the reconstructed total organic carbon (TOCRCN or TOCCoMod) as follows:

TOCRCN=TOCEM1+TOCEM2 . . . TOCEMx+TOCinert  (20)

The method for determining of the quantity of “inert carbon” in a rock sample is well established in the art of source rock analysis, and is described, for example, in Langford, F. F. and M.-M. Blank-Valleron, 1990, Interpreting Rock-Eval Pyrolysis Data Using Graphs of Pyrolizable Hydrocarbons vs. Total Organic Carbon. For the example from AMAD_2, reference is made to FIG. 7 where there is shown a plot of S2 Hydrocarbon Yield versus TOC for a nearby well, AMAD_1, which is used to assess the amount of inert carbon by noting the X-intercept of the resultant trend line.

Example 1
The following example of the method of the present invention uses the compositional modeling analysis method based on POPI data to calculate reconstructed total organic carbon (TOC) values. It has been found that reconstructed TOC values closely correspond to TOC measurements by standard instruments. The method of the invention can be used in unconventional oil reservoirs to provide relevant data in real time at the well site for use in making operational decisions while drilling. In addition to providing TOC data, the method also provides data that discriminates among the organic matter found the reservoir rock, such as free oil, residual bitumen and kerogen.
The data for this example is based upon data derived from core samples collected from an Arabian well identified as DBYT-1 that were analyzed for TOC and also by standard POPI analysis. The POPI data were processed by the GC-ROXTM commercial software available from the Saudi Arabian Oil Company and the compositional modeling method in order to assess the relative contribution of each organic matter end member to the total pyrolytic signature.
Shown in FIG. 8 is a plot of the end member pyrograms used to model the DBYT-1 well. The curve labeled “Free Oil” is representative of the free oil in the reservoir and was obtained by a cyclohexane extraction of the rock samples collected during drilling. Notably, the API gravity values for these samples as assessed by ranged between 32° and 34°. API was determined by pyrolysis methods disclosed in U.S. Pat. No. 6,823,298 B1. The “Bitumen” curve is representative of residual bitumen present in the formation that was obtained by extraction of the rock samples using methylene chloride. Its appearance is typical of what is observed for asphaltene fractions. The “Kerogen” curve is characteristic of the kerogen found in the reservoir and was determined by pyrolytic analysis of the rock samples after the above extractions by cyclohexane and methylene chloride.
Referring to FIG. 9, there are shown the results from modeling a typical pyrogram from the DBYT-1 well. The “Adjusted: (DBYT1-211)” curve is the actual sample data and the pattern-filled curve shows the results from the application of compositional modeling as a summation of the relative components of the end members. Notably, the model solution under estimates the amount of light components present in the LV (light volatile) portion of the pyrogram. This lower value reflects losses typically encountered in the light ends from extraction processes and also suggests that significant light components are in fact present the reservoir rock. The modeling results show high levels of free oil throughout the entire interval sampled.
Shown in FIG. 10 is a plot comparing the determination of TOC from modeling POPI data to measured TOC from instrumental data. The slope for the least squares fit of the line is very close to 1 with a high correlation coefficient. The results indicated a very robust correlation and provide confirmation of the utility and reliability of compositional modeling in determining the relative quantities of organic matter components in the source rocks.
In order to convert POPI data into TOC values, compositional modeling was used to determine the relative abundance of each end member in a rock sample. The chemical and physical properties of those end members were then used to calculate the weight percentage of both the pyrolyzable and non-pyrolyzable portions of each component. Significantly, the contribution of these three relatively simple components can be varied to achieve viable acceptable matches with actual rock samples. This implies that the model being used is a plausible representation of the organic matter distribution found in the rocks. The results can also form the basis to infer the level of transformation of organic matter into hydrocarbons that has occurred in the source rocks.
Shown in FIG. 11 is a composite log of five tracks from the GC-ROXTM software for a 290 foot section of the DBYT-1 well. In track one, a stack presentation of the volume of organic matter data from each end member is shown. In track two, the stack presentation of the modeled percentages of the end members from compositional modeling (CoMod %) are plotted. In track three, the value of TOC calculated from compositional modeling is plotted (circles) along with instrument measurements of TOC (triangles) for a large selection of the samples. In track four, TOC is compared with the quantity of the free oil in rock samples expressed as mg oil per gram of rock sample. In track five, CoMod TOC and light volatiles (LV) hydrocarbon yield per gram of TOC are plotted to illustrate the crossover significance as described below relative to the reference by Jarvie effect. Note that the “LV adj.” entry on plot is basically LV, but GC-ROX performs quality control during sample import and adjusts the yields slightly due to baselining procedures and smoothing algorithms.
Examination of the data presented in tracks one and two show a very high level of free oil throughout the DBYT-1 well. However, this data is difficult to compare with published accounts of unconventional oil reservoirs that use other parameters. This comparison is difficult because prior art methods rely on simple bulk integration and temperature cutoffs and do not distinguish the actual organic matter types present in the samples.
In track three, as shown in the cross plot of TOC data, there is a very close correspondence between the values of measured TOC and TOC calculated from CoMod. This plot justifies a high confidence level in the patterns that are shown in the plots from the pyrolysis data.
In tracks four and five, plots were constructed that are similar to published methods used for the assessment of unconventional oil reservoirs. From the standpoint of organic matter in hydrocarbon, the most commonly relied upon method of the prior art is to compare the TOC obtained by the Rock-Eval method and the S1 yield. As shown in Jarvie 2011, in press, AAPG Memoir 97, Shale reservoirs—Giant resources for the 21st century, J. Breyer, ed., in press, Jarvie, Daniel M., 2012, Shale resource systems for oil and gas: Part 1—Shale gas resource systems; Part 2—Shale oil resource systems, AAPG Memoir 97, p. 69-119, the preferred form of presenting the data is to plot S1 values and TOC values using the same scale. When S1 exceeds the value for TOC, these regions are said to exhibit crossover and are inferred to be favorable for unconventional oil production.
In the analysis shown in FIG. 11, two somewhat different plots are shown: one is TOC plotted against the yield of free oil and the other is TOC plotted against the LV (light volatile) yield. Neither one of these reflects precisely the determination of S1 yield. In the Rock-Eval method, S1 yield is determined by the hydrocarbons that are liberated at 300° C., whereas the LV yield is the amount of hydrocarbons liberated at 195° C. Thus, the LV yield is more conservative than S1 yield. The free oil yield as determined from compositional modeling represents the actual oil in any sample; however, it is possible that this results in a greater value than the simple S1 yield. However, the use of S1 yield at 300° C. is meant to represent the light hydrocarbons present in rock samples, and utilizing free oil values from compositional modeling is at present believed to be a viable indicator with good accuracy.
Both techniques of plotting light hydrocarbon yield versus TOC result in significant observed crossover effects in the DBYT-1 well. This suggests that the sampled section at the DBYT-1 well compares very favorably to productive unconventional oil reservoirs found as proven resources in North America.
Actual testing of the DBYT-1 well based on the method of the invention resulted in the recovery of some oil despite difficulties that occurred during well site operations and the testing procedures. In view of the fact that a great deal of the success in completing unconventional oil reservoirs in North America come only after significant investments in experimental well drilling and completion methods, the geochemical results from the DBYT-1 represent a viable method for future exploration and development drilling.
Example 2 and 3
The method described above was applied to samples from the DBYT-1 and AMAD-2 wells to yield the results set forth in Table 2 shown in FIGS. 12A-B and Table 3 shown in FIGS. 13A-B, respectively.
FIG. 14 is a cross-plot of instrument measurements of the TOC versus reconstructed TOC (TOCRCN) from the DBYT-1 well data reported in Table 2 and FIG. 6 is a cross plot of instrument measurements of the TOC versus reconstructed TOC (TOCRCN) from the AMAD-2 Well data reported in Table 3. The plots of FIGS. 14 and 15 demonstrate that TOC calculations that are reconstructed from compositional modeling provide a very close correlation to TOC values that are based on instrument measurements. The least squares fit of the line for the respective plots were as follows: DBYT-1, the R2=0.9418 and for AMAD-2 the R2=0.9493.
Example 4
An example of the improved interpretations provided through the reconstructed TOC method is illustrated by comparison of data from the DBYT-1 well to another potential unconventional oil resource well (MZLJ-29). FIG. 16 is a plot of the pyrogram for a representative sample from DBYT-1 produced during standard Rock-Eval analysis; i.e., using a starting temperature of 300° C. to capture the S1 hydrocarbon yield (in mg hydrocarbons/g Rock) and then programmed pyrolysis at 25° C./minute to 600° C. to capture the S2 hydrocarbon yield. As mentioned previously, the S1 hydrocarbon yield is commonly taken as representative of the “Free Oil” content of a sample. For example, see Jarvie, Daniel M., 2012, Shale resource systems for oil and gas: Part 1—Shale gas resource systems; Part 2—Shale oil resource systems, AAPG Memoir 97, p. 69-119.
FIG. 17 is the corresponding plot for a representative sample from MZLJ-29. The corresponding parameters commonly used for source rock/unconventional reservoir analysis are annotated on the plot. For both of these wells, the relative abundance of the “Free Oil” component as represented by S1 hydrocarbon yield is similar, with the Production Index (PI: S1/(S1+S2)) equal to 0.23 for DBYT-1 and 0.27 for MZLJ-29. This data would suggest a somewhat greater “Free Oil” content for MZLJ-29 than DBYT-1; however, the difference is not significant. Other parameters, such as the value of Tmax (445.7 for DBYT-1 vs. 451.3 for MZLJ-29) and the Hydrogen Index (HI; 373 for DBYT-1 vs. 193 for MZLJ-29) also suggest that greater organic matter transformation has occurred at the MZLJ-29 well.
FIGS. 18 and 19 are the corresponding plots of the same samples as assessed by the method of the invention. These samples were run with the standard analytical procedure as provided in the POPI method (U.S. Pat. No. 5,866,814) and subsequent application of the Compositional Modeling method (U.S. Pat. No. 7,363,206). In these plots, the original data are shown by the black solid pyrogram line and the results of Compositional Modeling are shown by the pattern-filled pyrogram, where the amount of each end-member for “Free Oil,” “Bitumen,” and “Kerogen” are reflected by the respective patterns in the legend. In this case, it can be clearly seen that the sample from the DBYT-1 well shows abundant “Free Oil,” whereas, the sample from the MZLJ-29 well shows almost no “Free Oil.” Also annotated on these plots, is the approximate division between material that would normally be comprised of S1 and S2 hydrocarbon yields under the standard Rock-Eval analysis. As can be clearly observed, a large amount of “Bitumen” can be associated with the early portion of a pyrolysis run. Significant differences are observed in the curve shape by using a lower starting temperature, which allows the differentiation between “Free Oil” and “Bitumen” end members.
FIG. 20 is a composite log of the data provided from the method of the invention that corresponds to FIG. 11 from the DBYT-1 well. As can be readily observed, the results for the well as a whole are similar to the example for single samples. Most notably, in these plots is the lack of “Free Oil” end member in the MZLJ-29 well samples. Also, Tracks 4 & 5 that show the crossover effect for DBYT-1 when plotting TOC and “Free Oil” or TOC and LV, show no crossover for the MZLJ-29 well. Lastly, well tests and stimulation of the section shown for the MZLJ-29 well did not recover any hydrocarbon, whereas, the DBYT-1 well was able to flow oil. Thus, the subject methods are able to predict well performance in cases where the standard industry method provides equivocal information.
The above description is intended to be illustrative and not exhaustive. This description will suggest many modifications, variations, and alternatives to those of ordinary skill in this art that are within the scope of the present invention. Those familiar with the art may recognize other equivalents to the specific embodiments described above. Accordingly, the scope of protection is to be determined by the claims that follow.

## Claims

1. A method of estimating the total organic carbon (TOC) content in a sample of reservoir rock from a known depth taken from a specific oil field, the sample containing a plurality of distinguishable end member (EM) components that can include oil, bitumen, soluble tar, pyrobitumen, kerogen, coal, drilling mud additives such as lubricants, mud stabilizers, lost circulation material, and other carbon-containing contaminants, the improvement comprising:
a. determining a set of local end members for compositional modeling;
b. applying compositional modeling to determine the relative proportion of each of the local end members in the sample;
c. recording the values of the percent of total yield of each end member in the sample;
d. obtaining parameters for the total hydrogen index (THIOM), the ratio of organic matter-to-pyrolizable hydrocarbons (OM/HCpy), the ratio of hydrogen to carbon in the organic matter (H/COM), the percent of Nitrogen, Oxygen, and Sulfur in the organic matter (NOSx%), and the weight of Hydrogen in the organic matter (Wt % Hom),
e. applying compositional modeling to determine the hydrocarbon yield (YieldEMx) that is attributable to each EM in the sample;
f. determining the weight of organic matter (OM) represented by the hydrocarbon yield in accordance with the following:


EMx Weight OM=YieldEMx×OM/HCPyx 



where HCPyx is a pyrolizable hydrocarbon end member;
g. determining the weight of elemental nitrogen, sulfur and oxygen in accordance with the following:


EMx WeightNSO=NSOx%×EMx Weight OM;


h. determining weight of hydrogen for each end member in a sample by the equation:


EMx WeightH=Wt%HOM×EMx Weight OM;


i. calculating the total organic carbon for the end member (TOCEMx) as follows:


TOCEMx=EMx Weight OM−EMx Weight NSOs−EMx Weight H;


j. summing the TOC for each end member and the inert carbon present in the formation to get the reconstructed total organic carbon (TOCRCN or TOCCoMod) as follows:


TOCRCN=TOCEM1+TOCEM2 . . . TOCEMx+TOCinert.

2. The method of claim 1 in which the set of local end members is determined by (i) analyzing the sample to determine the presence of each end member or (ii) selecting the end members based on pre-existing analytical results from rock samples obtained from one or more comparable zones of interest in adjacent wells in the reservoir rock formation.

3. The method of claim 1 which includes the further preliminary steps of collecting a plurality of samples from the same depth, identifying one or more samples that exhibit a uniformity in the content of at least one end member, selecting one sample having substantially uniform end member content, and subjecting that sample to steps (b) through (j).

4. The method of claim 1 which further includes:
k. selecting a plurality of samples from different intervals in the reservoir rock formation;
l. repeating steps (b) through (j) on each of the samples;
m. displaying the results in a graphic and/or tabular form for visual analysis by a user.

5. The method of claim 4 in which the results are displayed in tabular form for each of the plurality of samples and the table includes the estimated TOC and depth from which the sample was obtained, and related data selected from the group consisting of free oil, bitumen, kerogen and inert organic matter.

6. The method of claim 4 in which the TOC data are displayed graphically in stack form corresponding to the depth of the respective samples and the display includes volume of organic matter, % of kerogen, TOC vs free oil and TOC vs light volatiles (LV) as measured by the Pyrolytic Oil-Productivity Index (POPI) pyrolysis method.

7. A method for the reconstruction of total organic carbon content from compositional modeling analysis, the method comprising:
a. collecting a plurality of samples typical of unconventional oil reservoirs and source rocks;
b. performing pyrolysis by the Pyrolytic Oil-Productivity Index (POPI) pyrolysis method;
c. obtaining and reviewing preexisting pyrolysis data for comparative samples obtained from a nearby well and the same reservoir to assess the likely end-member components that appear to be present in those comparative samples;
d. selecting samples obtained in step (a) that appear to have a substantially uniform composition of an end-member, selected from the group of end members consisting of free oil, bitumen, kerogen such as appearing to be composed mostly of “free oil”, bitumen, kerogen and contaminants;
e. separating and moveable hydrocarbon component by extraction with a non-polar solvent, saving both the extract and extracted rock;
f. separating and soluble tar from extracted rock by extraction with polar solvent, saving both extract and extracted rock;
g. performing separation of residual organic matter on a portion of the polar solvent extracted rock by demineralization methods;
h. analyzing all extracts and organic matter separations by elemental analysis to determine the percent of C, H, N, O, and S in each sample;
i. analyzing all extracts, extracted rocks, and separated organic matter by pyrolysis-TOC analysis to determine THI;
j. analyzing the pyrolysis data by applying compositional modeling;
k. recording the value for the percentage of total yield for end-members in the group of samples;
l. obtaining parameters for the total hydrogen index (THIOM), the ratio of organic matter-to-pyrolizable hydrocarbons (OM/HCpy), the ratio of hydrogen-to-carbon in the organic matter (H/COM), the percent of Nitrogen, Oxygen, and Sulfur in the organic matter (NOSx%), and the weight of Hydrogen in the organic matter (Wt% HOM);
m. analyzing a plurality of samples from a well that has a source rock or unconventional reservoir where the reconstructed TOC values are needed;
n. applying compositional modeling to obtain the hydrocarbon yield that is attributable to each end member (YieldEMx) in a group of samples from a source rock or unconventional reservoir;
o. determining the weight of organic matter represented by the hydrocarbon yield by applying the equation:


EMx Weight OM=YieldEMx×OM/HCpyx;


p. determining the weight of elemental Nitrogen, Sulfur, and Oxygen by the equation:


EMx Weight NSOs=NOSx%×EMx Weight OM;


q. determining weight of Hydrogen for each end member in a sample by the equation:


EMxWeight H=Wt% HOM×EMx Weight OM;


r. calculating the total organic carbon for the end member (TOCEMx) as follows:


TOCEMx=EMx Weight OM−EMx Weight NSOs−EMx Weight H; and


s. summing the TOC for each end member and the inert carbon present in the formation to obtain the reconstructed total organic carbon (TOCRCN or TOCCoMod) as follows:


TOCRCN=TOCEM1+TOCEM2 . . . TOCEMx+TOCinert.