US20100155078A1 - Method For Predicting Composition of Petroleum - Google Patents
Method For Predicting Composition of Petroleum Download PDFInfo
- Publication number
- US20100155078A1 US20100155078A1 US12/573,428 US57342809A US2010155078A1 US 20100155078 A1 US20100155078 A1 US 20100155078A1 US 57342809 A US57342809 A US 57342809A US 2010155078 A1 US2010155078 A1 US 2010155078A1
- Authority
- US
- United States
- Prior art keywords
- petroleum
- expulsion
- recited
- petroleum expulsion
- predicting
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000003208 petroleum Substances 0.000 title claims abstract description 100
- 238000000034 method Methods 0.000 title claims abstract description 68
- 239000000203 mixture Substances 0.000 title claims abstract description 38
- 239000000047 product Substances 0.000 claims abstract description 82
- 239000000126 substance Substances 0.000 claims abstract description 50
- 238000006243 chemical reaction Methods 0.000 claims abstract description 44
- 238000010438 heat treatment Methods 0.000 claims abstract description 37
- 238000005336 cracking Methods 0.000 claims abstract description 33
- 239000007795 chemical reaction product Substances 0.000 claims abstract description 27
- 150000002430 hydrocarbons Chemical class 0.000 claims description 46
- 229930195733 hydrocarbon Natural products 0.000 claims description 41
- 238000004519 manufacturing process Methods 0.000 claims description 19
- 239000007789 gas Substances 0.000 description 29
- 239000011435 rock Substances 0.000 description 27
- 230000008569 process Effects 0.000 description 21
- 238000010586 diagram Methods 0.000 description 18
- 239000000571 coke Substances 0.000 description 13
- 238000000197 pyrolysis Methods 0.000 description 13
- 238000004227 thermal cracking Methods 0.000 description 13
- 239000010426 asphalt Substances 0.000 description 11
- 238000002474 experimental method Methods 0.000 description 11
- 230000004913 activation Effects 0.000 description 10
- 230000015572 biosynthetic process Effects 0.000 description 9
- 150000001875 compounds Chemical class 0.000 description 9
- 238000005755 formation reaction Methods 0.000 description 9
- 230000005012 migration Effects 0.000 description 9
- 238000013508 migration Methods 0.000 description 9
- 238000005194 fractionation Methods 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 7
- 238000013459 approach Methods 0.000 description 7
- 238000009533 lab test Methods 0.000 description 7
- 230000037361 pathway Effects 0.000 description 7
- 238000012545 processing Methods 0.000 description 7
- 230000000717 retained effect Effects 0.000 description 7
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 6
- 230000006870 function Effects 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 6
- 238000013213 extrapolation Methods 0.000 description 5
- 230000007246 mechanism Effects 0.000 description 5
- 230000015654 memory Effects 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 5
- 239000004215 Carbon black (E152) Substances 0.000 description 4
- 230000015556 catabolic process Effects 0.000 description 4
- 230000001186 cumulative effect Effects 0.000 description 4
- 238000000354 decomposition reaction Methods 0.000 description 4
- 238000006731 degradation reaction Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 4
- 239000000463 material Substances 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 239000002243 precursor Substances 0.000 description 4
- NINIDFKCEFEMDL-UHFFFAOYSA-N Sulfur Chemical compound [S] NINIDFKCEFEMDL-UHFFFAOYSA-N 0.000 description 3
- 230000009471 action Effects 0.000 description 3
- 150000004945 aromatic hydrocarbons Chemical class 0.000 description 3
- 230000006399 behavior Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 229910052757 nitrogen Inorganic materials 0.000 description 3
- 239000012071 phase Substances 0.000 description 3
- 229920006395 saturated elastomer Polymers 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 239000007787 solid Substances 0.000 description 3
- 241000894007 species Species 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 229910052717 sulfur Inorganic materials 0.000 description 3
- 239000011593 sulfur Substances 0.000 description 3
- 208000035126 Facies Diseases 0.000 description 2
- 230000035508 accumulation Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 239000003575 carbonaceous material Substances 0.000 description 2
- 239000012530 fluid Substances 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 230000005484 gravity Effects 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 230000014759 maintenance of location Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 239000011148 porous material Substances 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 238000005979 thermal decomposition reaction Methods 0.000 description 2
- 239000004079 vitrinite Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000001644 13C nuclear magnetic resonance spectroscopy Methods 0.000 description 1
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 1
- 238000004998 X ray absorption near edge structure spectroscopy Methods 0.000 description 1
- 238000004833 X-ray photoelectron spectroscopy Methods 0.000 description 1
- 239000008186 active pharmaceutical agent Substances 0.000 description 1
- 229910052586 apatite Inorganic materials 0.000 description 1
- 238000006065 biodegradation reaction Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 239000004927 clay Substances 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000004590 computer program Methods 0.000 description 1
- 238000006482 condensation reaction Methods 0.000 description 1
- 239000000470 constituent Substances 0.000 description 1
- CVSVTCORWBXHQV-UHFFFAOYSA-N creatine Chemical compound NC(=[NH2+])N(C)CC([O-])=O CVSVTCORWBXHQV-UHFFFAOYSA-N 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004992 fission Effects 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 230000000155 isotopic effect Effects 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 239000007791 liquid phase Substances 0.000 description 1
- 230000035800 maturation Effects 0.000 description 1
- 230000000877 morphologic effect Effects 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 239000004058 oil shale Substances 0.000 description 1
- 239000005416 organic matter Substances 0.000 description 1
- 239000003960 organic solvent Substances 0.000 description 1
- VSIIXMUUUJUKCM-UHFFFAOYSA-D pentacalcium;fluoride;triphosphate Chemical compound [F-].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O VSIIXMUUUJUKCM-UHFFFAOYSA-D 0.000 description 1
- 229920005547 polycyclic aromatic hydrocarbon Polymers 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 238000000045 pyrolysis gas chromatography Methods 0.000 description 1
- 238000010791 quenching Methods 0.000 description 1
- 230000000171 quenching effect Effects 0.000 description 1
- 230000035484 reaction time Effects 0.000 description 1
- 230000003362 replicative effect Effects 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 230000008646 thermal stress Effects 0.000 description 1
- 238000001149 thermolysis Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B41/00—Equipment or details not covered by groups E21B15/00 - E21B40/00
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
- G01N33/241—Earth materials for hydrocarbon content
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
Definitions
- Exemplary embodiments of the present invention relate to a method for predicting composition of petroleum generation and expulsion from source rocks.
- An exemplary embodiment of the present invention takes into account reactions from secondary cracking.
- Basin modeling is the process of using either proprietary or commercially available software to assess charge risk by integrating diverse geological and engineering data types into a model of one or more petroleum systems active in an area being explored.
- the scale of the model may range in size from a single drilling prospect to an entire basin.
- compositional yields model based on open and closed-system pyrolysis.
- Models that predict compositional yields as defined by varying chemical lumps are expressed typically by the bulk kerogen yield kinetic model (discrete E a 's fixed or varying A) where the percentage of each chemical lump is defined for each individual Ea.
- One known basin simulator has models of petroleum generation/expulsion that are either two component (gas/oil) or three component (dry gas/wet gas/oil).
- Generation is modeled using standard discrete energy kinetics derived from open-system pyrolysis for main-stage oil/gas generation and a late-gas generation model for methane formation from reacted kerogen residue. Both processes are modeled using the first order kinetics with a distributed E a and fixed A.
- the expulsion model is based on filling available rock pore space to a defined threshold. Secondary cracking reactions are Oil ⁇ Gas in the two-component model or Oil ⁇ Wet Gas and Wet Gas ⁇ Dry Gas in the three-component model.
- a method for predicting petroleum expulsion comprises defining a chemical structure of a kerogen and identifying a plurality of reaction products of the kerogen under geologic heating rates.
- the exemplary method also comprises grouping the plurality of reaction products into a plurality of product lumps based on their chemical composition and predicting petroleum expulsion for each of the plurality of product lumps based on secondary cracking reactions.
- petroleum expulsion is computed at a plurality of varying heating rates.
- a production rate curve may be defined for each of the plurality of product lumps.
- a plurality of stoichiometric coefficients may be defined for each of the plurality of product lumps.
- Predicting petroleum expulsion may comprise employing a theoretical model in one exemplary embodiment of the present invention.
- An exemplary method may comprise comparing predicted petroleum expulsion to known data.
- the chemical structure may be redefined if the predicted petroleum expulsion is outside a specified range with respect to the known data.
- an outline of secondary cracking reactions is prepared.
- predicted petroleum expulsion comprises a quantity of expelled petroleum. Moreover, predicted petroleum expulsion may additionally comprise timing data regarding petroleum expulsion.
- One exemplary embodiment of the present invention relates to a method of producing hydrocarbons from an oil and/or gas field.
- Such an exemplary embodiment may comprise defining a chemical structure of a kerogen deposited in the oil and/or gas field and identifying a plurality of reaction products of the kerogen under geologic heating rates.
- One exemplary embodiment comprises grouping the plurality of reaction products into a plurality of product lumps based on their chemical composition and predicting petroleum expulsion for each of the plurality of product lumps based on secondary cracking reactions.
- Hydrocarbons may be extracted from the oil and/or gas field using the predicted petroleum expulsion.
- petroleum expulsion is computed at a plurality of varying heating rates.
- a production rate curve may be defined for each of the plurality of product lumps.
- a plurality of stoichiometric coefficients may be defined for each of the plurality of product lumps.
- Predicting petroleum expulsion may comprise employing a theoretical model in one exemplary embodiment of the present invention.
- An exemplary method of producing hydrocarbons may additionally comprise comparing predicted petroleum expulsion to known data. The chemical structure may be redefined if the predicted petroleum expulsion is outside a specified range with respect to the known data.
- an outline of secondary cracking reactions is prepared.
- predicted petroleum expulsion comprises a quantity of expelled petroleum. Moreover, predicted petroleum expulsion may additionally comprise timing data regarding petroleum expulsion.
- An exemplary tangible, machine-readable medium comprises code adapted to define a chemical structure of a kerogen and code adapted to identify a plurality of reaction products of the kerogen under geologic heating rates.
- the exemplary tangible, machine-readable medium may additionally comprise code adapted to group the plurality of reaction products into a plurality of product lumps based on their chemical composition and code adapted to predict petroleum expulsion for each of the plurality of product lumps based on secondary cracking reactions.
- FIG. 1 which includes FIGS. 1A and 1B , is a process flow diagram showing a method for predicting petroleum production in accordance with an exemplary embodiment of the present invention
- FIG. 2 is a diagram that shows the modeling of expelled petroleum in accordance with an exemplary embodiment of the present invention
- FIG. 3 is a plurality of graphs showing rate curves indicating rates of expulsion of petroleum at different heat rates
- FIG. 4 is a graph showing a family of curves corresponding to cumulative yields for a particular fraction as a function of temperature at varying heating rates in accordance with an exemplary embodiment of the present invention
- FIG. 5 is a graph showing a rate of expulsion expressed as a single frequency factor and a discrete spectrum of activation energies in accordance with an exemplary embodiment of the present invention
- FIG. 6 is a diagram showing an outline of secondary cracking reactions in accordance with an exemplary embodiment of the present invention.
- FIG. 7 is a graph showing an output from a calculation of the thermal cracking of a Polars (NSO) fraction in accordance with an exemplary embodiment of the present invention.
- FIG. 8 is a graph showing a comparison of overall thermal cracking expressions in accordance with an exemplary embodiment of the present invention.
- FIG. 9 is a graph showing a comparison between results produced by an open system pyrolysis approach and an exemplary embodiment of the present invention.
- FIG. 10 is a process flow diagram showing a method for producing hydrocarbons in accordance with an exemplary embodiment of the present invention.
- FIG. 11 is a diagram of a tangible, machine-readable medium in accordance with an exemplary embodiment of the present invention.
- FIG. 12 illustrates an exemplary computer network that may be used to perform the method for predicting composition of petroleum as disclosed herein, and is discussed in greater detail below.
- Basin model refers to a simplification of the earth and its processes with the intent being to track the dynamic evolution of one or more of those processes through time.
- the processes related to the generation and migration of hydrocarbons is commonly modeled with the intent to determine which of several possible structural culminations may be the most prospective for containing a commercial accumulation.
- Basin models use data from seismic, well control and knowledge of the geology of the area to construct a numerical model of the region and to track the changes in the various modeled parameters through time to reach a set of predictions that are then calibrated to the known information at the present.
- the model parameters are then adjusted within geologically reasonable bounds until a successful match and calibration is reached. Prediction can then be made at locations away from the calibration points.
- kerogen refers to a solid, carbonaceous material.
- the mixture is referred to as oil shale.
- the mineral is, in fact, technically shale, that is, a rock formed from compacted clay.
- Kerogen is subject to decomposing upon exposure to heat over a period of time. Upon heating, kerogen molecularly decomposes to produce oil, gas, and carbonaceous coke. Small amounts of water may also be generated. The oil, gas and water fluids are mobile within the rock matrix, while the carbonaceous coke remains essentially immobile. Kerogen may be classified into four distinct groups: Type I, Type II, Type III, and Type IV.
- Kerogen types used herein are as defined in Tissot and Welte (Tissot, B. P. and Welte, D. H., Petroleum Formation and Occurrence, second edition, Springer-Verlag, Berlin, 1984, p. 151).
- the maturation sequence for kerogen that typically occurs over geological time is due to burial leading to exposure to increased temperature and pressure.
- Classification of kerogen type may depend upon precursor materials of the kerogen.
- the precursor materials transform over time into macerals or amorphous masses. Macerals are microscopic structures that have distinguishing morphologies, different chemical structures and properties depending on the precursor materials from which they are derived.
- Amorphous kerogens have no distinguishing morphological features that can be used to characterize its precursor materials, but may have different chemical structures and properties.
- tangible machine-readable medium refers to a medium that participates in directly or indirectly providing signals, instructions and/or data to a processing system.
- a machine-readable medium may take forms, including, but not limited to, non-volatile media (e.g., ROM, disk) and volatile media (RAM).
- Machine-readable medium include, but are not limited to, a floppy disk, a flexible disk, a hard disk, a magnetic tape, other magnetic medium, a CD-ROM, other optical medium, punch cards, paper tape, other physical medium with patterns of holes, a RAM, a ROM, an EPROM, a FLASH-EPROM, or other memory chip or card, a memory stick, and other media from which a computer, a processor or other electronic device can read.
- NSO nitrogen, sulfur, and oxygen containing compounds
- a basin model according to an exemplary embodiment of the present invention comprises an interactive tool that geologists may use for a variety of purposes.
- the model may be used to test the impact of various geologic scenarios such as source characteristics, thermal history, structural tectonics and the like on potential plays and prospects.
- Another use of the model may include assessing the charge risk associated with each of the various elements of a given petroleum system, such as determining if there is adequate source richness, sufficient source volume, adequate source maturity, and appropriate timing of generation relative to timing of trap formation, to name just a few examples.
- Another exemplary use of model includes estimating the volume of petroleum generated from a given volume of source rock(s), the volume of petroleum expelled from the source rock, the petroleum losses during migration, and ultimately the quantity of petroleum delivered to the trap(s).
- the basin model may help the geologist address issues such as the timing of hydrocarbon generation relative to the timing a trap formation.
- the model may also be useful in predicting the volumes of hydrocarbons delivered to a trap and the hydrocarbon type (liquid vs. gas and their relative proportions) likely to be currently present in the trap.
- the model may be used to predict the physical state of the hydrocarbon charge (Gas-Liquid phase behavior at specific pressure-temperature) and the composition of the gas/liquids to assess economic value.
- the model may assist in determining the possible migration paths to the prospect.
- the model may be used to predict the possible post-charge loss of hydrocarbons from the trap (leakage through seal, tilt and spill of trap, thermal cracking, biodegradation and the like).
- a basin model in accordance with an exemplary embodiment of the present invention may be constructed by integrating various types of geological and geochemical data.
- a thermal history is defined that is based on geophysical models/measurement (for example, tectonic reconstruction, basinal heat flow, well temperatures) and or geochemical measurements (for example, vitrinite reflectance data, Rock-Eval data, apatite fission track data, fluid inclusion data, oil/gas molecular and isotopic compositions or the like).
- the generative potential of source rock units are specified by their organic richness (for example, percentage of total organic content or % TOC, Rock-Eval HI) and their spatial distribution.
- a generative model (kinetics and product yields) is assigned to individual source rock units based on either direct measurement or extrapolation of organic facies.
- the thermal history is then applied to the kinetic model and the timing and amount of oil/gas generated is calculated over the course of the basin evolution. These calculations are then compared to models of potential migration pathways and the timing of trap and seal formation, resulting in a prediction of the petroleum accumulations.
- the prediction of the amount and composition of oil and gas generated and expelled from organic rich source rocks is a key factor in basin simulations of petroleum systems. Other aspects of basin modeling, such as migration, trap formation, and seal integrity, are equally important, but are not discussed further as they are not relevant to an exemplary embodiment of the present invention.
- Determining the composition of oil and gas as it is expelled from the source rocks consists of modeling three major processes.
- the first process comprises the manner in which kerogen (the solid organic matter within source rocks that generates oil) reacts at the specified time-temperature history experienced by the source rocks over geologic time. This includes determining the rate of the reaction (kinetics) and the composition of the products (bitumen) that is formed.
- the second process comprises the mechanisms responsible for expelling the bitumen from the kerogen out of the source rock. Expulsion is known to be controlled by the amount and chemistry of the kerogen. Chemical fractionation occurs whereby saturated and aromatic hydrocarbons are selectively expelled while polar compounds and asphaltenes (collectively known as NSO's for compounds that contain nitrogen, sulfur, and/or oxygen) are selectively retained.
- the third process comprises secondary thermal cracking of the bitumen (retained within the source rock matrix) or expelled petroleum (either along migration pathways or within reservoir rocks) may occur that further alters the composition of the petroleum.
- the timing of petroleum (oil and gas) generation from source rocks is dependent on the kinetic model.
- These kinetic models may be derived from laboratory experiments where a source rock, kerogen, asphaltene, or other carbonaceous material is artificially matured and the product yields are measured as a function of varying heating rate, temperature, and/or time.
- most models of oil and gas generation are based on measurements of pyrolyzate yields from experiments conducted on a specific sample that are then fitted to a series of parallel first-order reactions or other kinetic reaction schemes.
- these laboratory studies accelerate the temperature (heating rate) to compensate for short reaction time.
- models are derived from reactions occurring at temperatures from ⁇ 250 to 650° C.
- kinetic models based on bulk kerogen thermal decomposition can only approximate actual product yields and do not model product composition.
- a separate kinetic model may be determined for each product lump that may have different frequency factors.
- product lump means a group of reaction products that are similar in chemical composition and behavior.
- a method in accordance with an exemplary embodiment of the present invention may desirably be performed such that the individual products are quantified as a function of time/temperature and may involve complex analyses of closed-system pyrolysates or customized instrumentation.
- each individual activation energy is divided into portions representative of the lumped components.
- This approach is considered easier to define as the bulk kerogen decomposition kinetics can be determined by routine open-system pyrolysis methods and the products can be independently quantified.
- the compositional model may reflect either primary bitumen or expelled petroleum.
- a bitumen model may desirably calculate the chemical fractionation that occurs upon expulsion, while a petroleum model has already adjusted the lumped component compositions to reflect this fractionation.
- an exemplary embodiment of the present invention models compositional changes in petroleum resulting from secondary cracking.
- FIG. 1 which includes FIG. 1A and 1B , is a process flow diagram showing a method for predicting differences in subsurface conditions in accordance with an exemplary embodiment of the present invention.
- the method is generally referred to by the reference number 100 .
- a chemical analysis of a kerogen is performed.
- a chemical structure referred to as CS in FIG. 1 , is populated.
- a decision is made regarding whether the chemical structure has been optimized to match known analytical data, as shown at block 106 . If the chemical structure is not optimized with respect to known analytical data, process flow returns to block 104 , where the chemical structure is re-populated.
- modeling of the amount and composition of generated and expelled oil and gas is performed on the chemical structure using a process known as chemical structure-chemical yields modeling or CY CSM.
- a CS-CYM is generally described in U.S. Pat. No. 7,344,889, entitled “Chemical Structural and Compositional Yields Model for Predicting Hydrocarbon Thermolysis Products”, which issued to Kelemen, et al. on Mar. 18, 2008.
- a CY-CSM is also described in the following publication: Freund H., Walters, C. C. Kelemen S. R., Siskin M., Gorboty, M. L., Curry D. J., Bence A. E., 2006. Organic Geochemistry, 38 288-305.
- CY-CSM relies on constructing a representative model of the molecular structure of a kerogen using constraints from numerous chemical analyses, modeling the thermal reactions that occur within this representative structure, and applying a thermodynamic-based theory of expulsion and chemical fractionation to calculate the product yields.
- This model has been shown to predict the yields and composition of oil and gas generated in high temperature experiments. When basinal heating rates are used, the compositions of both expelled oil and retained bitumen are modeled following the reaction pathways, kinetics, and expulsion mechanisms that occur under these conditions and are not extrapolations of high-temperature experimental data.
- Yield curves for individual fractions are extracted, as shown at block 110 of FIG. 1A . For example, yield curves for a group of hydrocarbons 112 may be individually extracted, as set forth below.
- each product may be considered independently within the basin simulator.
- hydrocarbons for example, C 1 , C 2 , C 3 , C 4 , and C 5 hydrocarbons
- lumped hydrocarbons for example, C 6 -C 14 Aromatics, C 6 -C 15 Saturates, C 15+ Saturates, C 15+ Aromatics, and C 15+ polars.
- a kinetic model for each of the group of hydrocarbons 112 is produced, as shown at block 114 . Because the kinetics of each expelled product are modeled using the yields calculated from a coupled generation and expulsion model, each product may be considered independently within the basin simulator.
- a kinetic model for each compositional component is defined.
- the implementation is extremely flexible as new components (for example, non-HC gases) or refinements of a lumped component (for example, C 15+ polars modeled as C 15+ resins and C 15+ asphaltenes) can be easily changed without affecting the other component models.
- the kinetics of the individual components are determined by running CS-CYM at varying heating rates typical of geologic basins (0.5 to 32° C./million years (Ma) are given as examples in block 108 ) and fitting the calculated yield curves to a discrete activation energy model.
- additional CS-CYM modeling is performed on a subset of the group of hydrocarbons 112 , as shown at block 116 .
- Yield curves for the subset of the group of hydrocarbons are extracted for the resulting individual fractions, as shown at block 118 .
- Examples of these individual fractions include a first fraction 120 , a second fraction 122 and a third fraction 124 .
- the first fraction 120 includes input from additional ones of the group of hydrocarbons 112 .
- kinetic models of secondary cracking are applied to a second subset of the group of hydrocarbons 112 .
- kinetic models and stoichiometric coefficients are derived based on input data from the first fraction 120 , the second fraction 122 and the third fraction 124 .
- modeling of secondary cracking reactions which further alter expelled oil composition via thermal reactions that occur within source rock mineral matrices, along migration pathways, and within reservoir rock, also are calculated using CS-CYM.
- generation/expulsion of a compositional lump for example, C 15+ Saturates
- the individual product is trapped.
- the trapped product is then re-run within CS-CYM under closed conditions and the product yields are calculated.
- FIG. 2 is a diagram that shows the modeling of expelled petroleum in accordance with an exemplary embodiment of the present invention.
- the diagram is generally referred to by the reference number 200 .
- CS-CYM is used to derive the data needed to determine the yields and kinetic model for the coupled generation/expulsion of individual lumped products and the yields and secondary cracking kinetics of the C 14+ lumped products.
- a first panel 202 of the diagram 200 shows a snapshot of a Type 1 kerogen substructure.
- kerogen is isolated from its rock matrix and analyzed using numerous analytical methods (for example, elemental CHNOS analysis, solid-state 13 C NMR, X-ray photoelectron spectroscopy, Sulfur-XANES, pyrolysis-GC or the like) that provide constraining input values used to construct a representation ( ⁇ 1 million atoms) of the kerogen structure.
- the chemical structure of the kerogen is then placed within a reaction network that predicts the composition of the generated and expelled products under specified temperature conditions, as shown at a second panel 204 of the diagram 200 .
- CS-CYM For simulating high temperature, open system experiments such as those typically used for determining bulk kerogen degradation kinetics, CS-CYM first calculates the products produced from kerogen cracking and then determines if the product is sufficiently volatile to be considered as an expelled product. Yield curves for expelled products may be compared to known geologic data, as shown in a third panel 206 of the diagram 200 .
- the third panel 206 includes a left graph and a right graph.
- the left graph includes a y-axis 208 that represents a maximum temperature in degrees Celsius.
- An x-axis 210 of the left graph represents a heating rate in degrees Celsius per minute.
- the right graph of the third panel 206 includes a y-axis 212 that represents the rate of expelled products from a simulation of an open-system laboratory pyrolysis in units of mg expelled product/g total carbon/0.2 minutes.
- An x-axis 214 of the right graph represents temperature in units of degrees Celsius. The calculated yields compare very well with laboratory yields when kerogens are reacted at high temperatures and heating rates, suggesting that the reaction network is replicating the basic chemical reaction and expulsion processes that occur under these conditions.
- the products replicate high temperature laboratory results.
- a CS-CYM program is used to determine the product yield under geologic heating rates.
- the generated products are calculated as several thousand individual species, then lumped into chemical classes defined as: C 1 , C 2 -C 5 , C 6 -C 14 saturates, C 6 -C 14 aromatics+polars, C 15+ saturates, C 15+ aromatics, C 15+ polars.
- the C 2 -C 5 component may be further divided into C 2 , C 3 , C 4 , C 5 hydrocarbons if the final model is to be used in an Equation of State (EOS) calculation to determine phase behavior.
- EOS Equation of State
- C 15+ polars may be split into C 15+ NSO and C 15+ asphaltenes based on a solubility test of individual components.
- the composition of the expelled products are calculated using a simplified formulation of a thermodynamic theory of expulsion and chemical fractionation that is an extended Flory-Rehner Regular Solution.
- This expulsion theory couples concepts of polymer theory to describe the amount of bitumen that can dissolve in a kerogen matrix with Regular Solution theory that explains which molecules are preferentially retained and which are expelled.
- a simplification of the theory is desirable as an exact calculation of the composition of the equilibrated retained and expelled is too currently computationally intensive to be run for all components at all time steps.
- FIG. 3 is a plurality of graphs showing rate curves indicating rates of expulsion of petroleum at different heat rates.
- the plurality of graphs is generally referred to by the reference number 300 .
- the expelled products are modeled at varying heating rates. Typical heating rates used are 0.5, 1, 2, 4, 8, 16, and 32° C./Ma, which span the conditions most commonly seen in nature.
- a family of curves is shown for heat rates of 0.5° C./Ma, 2° C./Ma and 8° C./Ma.
- a first graph 302 comprises a y-axis 304 that corresponds to a rate if expelled oil at a heating rate of 0.5° C./Ma.
- the first graph 302 includes an x-axis 306 that corresponds to temperature in degrees Celsius.
- a family of curves 308 shows the rate of expelled oil for a plurality of different reaction products. Moreover, each curve within the family of curves 308 corresponds to a different reaction product.
- a second graph 310 comprises a y-axis 312 that corresponds to a rate if expelled oil at a heating rate of 2° C./Ma.
- the second graph 310 includes an x-axis 314 that corresponds to temperature in degrees Celsius.
- a family of curves 316 shows the rate of expelled oil for a plurality of different reaction products. Moreover, each curve within the family of curves 316 corresponds to a different reaction product.
- a third graph 318 comprises a y-axis 320 that corresponds to a rate if expelled oil at a heating rate of 8° C./Ma.
- the second graph 318 includes an x-axis 322 that corresponds to temperature in degrees Celsius.
- a family of curves 324 shows the rate of expelled oil for a plurality of different reaction products. Moreover, each curve within the family of curves 324 corresponds to a different reaction product.
- FIG. 4 is a graph showing a family of curves corresponding to cumulative yields for a particular fraction as a function of temperature at varying heating rates in accordance with an exemplary embodiment of the present invention.
- each individual fraction or product lump is analyzed as a separate entity.
- the graph is generally referred to by the reference number 400 .
- the graph 400 includes a y-axis 402 that corresponds to a fraction reacted index ranging from 0.0 to 1.0.
- An x-axis 404 corresponds to temperature in degrees Celsius.
- a family of evolution curves 406 represent cumulative yields as a function of temperature at varying heating rates.
- the family of evolution curves may be fitted using fitting software. From this analysis, a kinetic model is generated whereby the rate of expulsion of a specific product lump is expressed by a single frequency factor and a discrete spectrum of activation energies.
- FIG. 5 is a graph showing a rate of expulsion expressed as a single frequency factor and a discrete spectrum of activation energies in accordance with an exemplary embodiment of the present invention.
- the graph is generally referred to by the reference number 500 .
- the graph 500 includes a y-axis 502 corresponding to a percentage rate of expulsion of a particular fraction.
- An x-axis 504 corresponds to activation energy in units of kcal/mol.
- the bars on the graph represent a percentage of overall rate of expulsion for a typical reaction product across a spectrum of activation energies.
- the model is completed by determining the stoichiometric coefficients for the finally cumulative product.
- the amount of each lumped expelled product i that is expelled from a source kerogen is calculated as: TOC ⁇ HI ⁇ f i. ⁇ a i , where a i is a fixed coefficient. Since the kinetics of expulsion of C 2 -C 5 are nearly identical, they can be calculated as a single C 2 -C 5 lump.
- the individual C 2 , C 3 , C 4 , and C 5 components are then separated from a fixed stoichiometric relationship (b i ).
- the residual kerogen termed coke, is a non-mobile product that can convert to inert char and CH 4 at higher temperatures.
- the amount of coke is calculated by grams of kerogen ⁇ f i ⁇ a coke where a coke is a fixed coefficient.
- the discrete energy kinetic models for expelled products are readily incorporated into basin simulators as the form is similar to what is currently used for kerogen degradation kinetics.
- the models are for expelled product.
- the generation and expulsion mechanisms, along with chemical fractionation, are inherently built into the model and the simulator need not calculate generation and expulsion using different models and/or principles.
- FIG. 6 is a diagram showing an outline of secondary cracking reactions in accordance with an exemplary embodiment of the present invention.
- the diagram is generally referred to by the reference number 600 .
- the diagram shows a vertical column of expelled products 602 and a corresponding series of secondary cracking reaction products 604 .
- CS-CYM may be performed in a manner that the expelled products for each chemical lump were mathematically trapped as they were expelled. These molecular pools were then subjected to a second round of heating at various heating rates under closed-system conditions and the degradation of each product lump calculated. From these rates curves, a kinetic model was fashioned in the same manner as described above. This procedure was applied to the C 15+ fractions. While the models for the thermal cracking of C 15+ saturated and aromatic hydrocarbons were similar regardless of the composition of the original source kerogen, the thermal cracking of the C 15+ polar compounds differ. Consequently, the origin of the polar compounds must be tracked within the basin simulator. The kinetics for the thermal cracking of ⁇ C 15+ fractions were derived from literature values and laboratory heating experiments.
- FIG. 7 is a graph showing an output of mobile products from a calculation of the thermal cracking of a Polars (NSO) fraction in accordance with an exemplary embodiment of the present invention.
- the graph is generally referred to by the reference number 700 .
- the Polars (NSO) fraction shown in FIG. 7 is from Green River Shale (Type I) kerogen at a 0.5° C./Ma heating rate.
- the graph 700 comprises a y-axis 702 that corresponds to NSO secondary cracking products measured in mg/gNSO.
- An x-axis 704 corresponds to temperature in degrees Celsius.
- Each of a family of curves 706 corresponds to a different component of the Polars (NSO) fraction.
- FIG. 8 is a graph showing a comparison of overall thermal cracking expressions in accordance with an exemplary embodiment of the present invention.
- the graph is generally referred to by the reference number 800 .
- the graph 800 comprises a y-axis 802 corresponding to a percentage of fractional conversion ranging from 0 to 1.
- An x-axis 804 corresponds to temperature in degrees Celsius.
- Each of a family of curves 806 corresponds to a different component of a fraction.
- FIG. 9 is a graph showing a comparison between results produced by an open system pyrolysis approach and an exemplary embodiment of the present invention.
- the graph is generally referred to by the reference number 900 .
- the graph 900 comprises a y-axis 902 corresponding to a percentage of fractional conversion ranging from 0 to 1.
- An x-axis 904 corresponds to temperature in degrees Celsius.
- a first curve 906 shows results produced by an exemplary embodiment of the present invention.
- a second curve 908 shows results produced using an open system pyrolysis approach.
- LOM 8 to 10.5 vitrinite reflectance
- the open-system pyrolysis kinetic model predicts petroleum is generated from ⁇ 0.7% to ⁇ 1.2% (LOM 9 to 11).
- the calculations were conducted using a 4° C./Ma heating rate that approximates the geothermal history of the Uinta Basin during Paleocene-Eocene subsidence.
- FIG. 10 is a process flow diagram showing a method for producing hydrocarbons in accordance with an exemplary embodiment of the present invention.
- the method is generally represented by the reference number 1000 .
- the method begins.
- a chemical structure of a kerogen deposited in the oil and/or gas field is defined, as shown at block 1004 .
- a plurality of reaction products of the kerogen under geologic heating rates is identified.
- the reaction products are grouped into a plurality of product lumps based on their chemical composition, as shown at block 1008 .
- petroleum expulsion for each of the plurality of product lumps is predicted based on secondary cracking reactions.
- Hydrocarbons are extracted from the oil and/or gas field using the predicted petroleum expulsion, as shown at block 1012 .
- the process ends.
- FIG. 11 is a diagram of a tangible, machine-readable medium in accordance with an exemplary embodiment of the present invention.
- the tangible, machine-readable medium is generally referred to by the reference number 1100 .
- the tangible, machine-readable medium 110 may comprise a disk drive such as a magnetic or optical disk or the like.
- the tangible, machine-readable medium 1100 comprises code 1102 adapted to define a chemical structure of a kerogen.
- Code 1104 adapted to identify a plurality of reaction products of the kerogen under geologic heating rates is stored on the tangible, machine-readable medium.
- code 1106 adapted to group the plurality of reaction products into a plurality of product lumps based on their chemical composition is stored on the tangible, machine-readable medium.
- the tangible, machine-readable medium stores code 1108 adapted to predict petroleum expulsion for each of the plurality of product lumps based on secondary cracking reactions.
- An exemplary embodiment of the present invention relates to a method for calculating the generation and expulsion of oil and gas and its secondary thermal cracking that differs from previously published techniques.
- determination of the quantity and timing of petroleum generation is performed according to a theoretical model as opposed to extrapolation of laboratory experiments.
- the generation and expulsion models are coupled into a single model for expelled petroleum.
- a theoretical model to describe the quantity and timing of generation of hydrocarbons from the secondary cracking of hydrocarbons and polar compounds may additionally be employed.
- An exemplary embodiment of the present invention provides an approximation of an expulsion model for chemical fractionation.
- the stoichiometry and kinetics of the secondary cracking of individual lumps are modeled.
- FIG. 12 illustrates an exemplary computer system 1200 on which software for performing processing operations of embodiments of the present invention may be implemented.
- a central processing unit (CPU) 1201 is coupled to system bus 1202 .
- the CPU 1201 may be any general-purpose CPU. The present invention is not restricted by the architecture of CPU 1201 (or other components of exemplary system 1200 ) as long as CPU 1201 (and other components of system 1200 ) supports the inventive operations as described herein.
- the CPU 1201 may execute the various logical instructions according to embodiments. For example, the CPU 1201 may execute machine-level instructions for performing processing according to the exemplary operational flow described above in conjunction with FIG. 10 . For instance, CPU 1201 may execute machine-level instructions for performing operational block 1006 of FIG. 10 , as an example.
- the computer system 1200 also preferably includes random access memory (RAM) 1203 , which may be SRAM, DRAM, SDRAM, or the like.
- the computer system 1200 preferably includes read-only memory (ROM) 1204 which may be PROM, EPROM, EEPROM, or the like.
- the RAM 1203 and the ROM 1204 hold user and system data and programs, as is well-known in the art.
- the computer system 1200 also preferably includes an input/output (I/O) adapter 1205 , a communications adapter 1211 , a user interface adapter 1208 , and a display adapter 1209 .
- the I/O adapter 1205 , the user interface adapter 1208 , and/or communications adapter 1211 may, in certain embodiments, enable a user to interact with computer system 1200 in order to input information.
- the I/O adapter 1205 preferably connects to a storage device(s) 1206 , such as one or more of hard drive, compact disc (CD) drive, floppy disk drive, tape drive, etc. to computer system 1200 .
- the storage devices may be utilized when the RAM 1203 is insufficient for the memory requirements associated with storing data for operations of embodiments of the present invention.
- the data storage of the computer system 1200 may be used for storing information and/or other data used or generated in accordance with embodiments of the present invention.
- the communications adapter 1211 is preferably adapted to couple the computer system 1200 to a network 1212 , which may enable information to be input to and/or output from system 1200 via such network 1212 (e.g., the Internet or other wide-area network, a local-area network, a public or private switched telephony network, a wireless network, any combination of the foregoing).
- the user interface adapter 1208 couples user input devices, such as a keyboard 1213 , a pointing device 1207 , and a microphone 1214 and/or output devices, such as a speaker(s) 1215 to the computer system 1200 .
- the display adapter 1209 is driven by the CPU 1201 to control the display on a display device 1210 to, for example, display information or a representation pertaining to a portion of a subsurface region under analysis, such as displaying a generated 3D representation of a target area, according to certain embodiments.
- the present invention is not limited to the architecture of system 1200 .
- any suitable processor-based device may be utilized for implementing all or a portion of embodiments of the present invention, including without limitation personal computers, laptop computers, computer workstations, and multi-processor servers.
- embodiments may be implemented on application specific integrated circuits (ASICs) or very large scale integrated (VLSI) circuits.
- ASICs application specific integrated circuits
- VLSI very large scale integrated circuits
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Chemical & Material Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Analytical Chemistry (AREA)
- Crystallography & Structural Chemistry (AREA)
- Medicinal Chemistry (AREA)
- Remote Sensing (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Food Science & Technology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Investigating Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)
Abstract
A method for predicting petroleum expulsion is provided. An exemplary embodiment of the method comprises defining a chemical structure of a kerogen and identifying a plurality of reaction products of the kerogen under geologic heating rates. The exemplary method also comprises grouping the plurality of reaction products into a plurality of product lumps based on their chemical composition and predicting petroleum expulsion for each of the plurality of product lumps based on secondary cracking reactions.
Description
- This application claims the benefit of U.S. Provisional Patent Application 61/140,249 filed Dec. 23, 2008 entitled METHOD FOR PREDICTING COMPOSITION OF PETROLEUM, the entirety of which is incorporated by reference herein.
- Exemplary embodiments of the present invention relate to a method for predicting composition of petroleum generation and expulsion from source rocks. An exemplary embodiment of the present invention takes into account reactions from secondary cracking.
- This section is intended to introduce various aspects of the art, which may be associated with exemplary embodiments of the present invention. This discussion is believed to assist in providing a framework to facilitate a better understanding of particular aspects of the present invention. Accordingly, it should be understood that this section should be read in this light, and not necessarily as admissions of prior art.
- Basin modeling is the process of using either proprietary or commercially available software to assess charge risk by integrating diverse geological and engineering data types into a model of one or more petroleum systems active in an area being explored. The scale of the model may range in size from a single drilling prospect to an entire basin.
- There are several commercial basin simulators that are widely used in academia and industry. These known simulators base the kinetics of petroleum generation on some form of a series of parallel first order reactions, usually expressed by the Arrhenius Equation as a discrete spectrum of activation energies (Ea). These models comprise a distribution of Ea's at fixed or variable spacing (typically 1 kcal/mole) using a single or varying frequency factor (A). The parameters for these models are based on various laboratory experiments. Some known models rely on simple compositional models that are rooted in a petroleum modeling computer program called PMOD developed and made available by the Lawrence Livermore National Laboratory. Additionally, some models offer compositional yield models based on forms of pyrolysis, including open-system pyrolysis. Other models use a compositional yields model based on open and closed-system pyrolysis. Models that predict compositional yields as defined by varying chemical lumps are expressed typically by the bulk kerogen yield kinetic model (discrete Ea's fixed or varying A) where the percentage of each chemical lump is defined for each individual Ea.
- One known basin simulator has models of petroleum generation/expulsion that are either two component (gas/oil) or three component (dry gas/wet gas/oil). Generation is modeled using standard discrete energy kinetics derived from open-system pyrolysis for main-stage oil/gas generation and a late-gas generation model for methane formation from reacted kerogen residue. Both processes are modeled using the first order kinetics with a distributed Ea and fixed A. The expulsion model is based on filling available rock pore space to a defined threshold. Secondary cracking reactions are Oil→Gas in the two-component model or Oil→Wet Gas and Wet Gas→Dry Gas in the three-component model. Hence, only the gas/oil ratio or GOR changes as a function of oil cracking with the oil remaining at a fixed composition and quality. In reality, the composition of petroleum changes in a systematic fashion with increasing thermal stress. Polar compounds and large polynuclear aromatic hydrocarbons or PAH are depleted, while the distribution of the surviving hydrocarbons shifts toward smaller species. The net-effect on oil quality is an increase in American Petroleum Institute gravity (referred to as API gravity herein) and GOR and a decrease in sulfur and nitrogen.
- Many known modeling techniques employ data that is derived from laboratory experiments. A potential problem relating to laboratory data involving high temperatures is that the coupled processes of generation, expulsion, and secondary cracking may be perturbed and follow mechanisms and pathways that do not occur under geologic conditions. For example, open system pyrolysis experiments measure gas-phase products that are transported from the source rock place in a heated sample holder to a flame ionization or mass spectrometer detector. Hence, expulsion is controlled mostly by volatility. In closed-system pyrolysis experiments, product yields may be defined in terms of gas-phase volatility, solubility in organic solvents, or as free-floating bitumen in hydrous systems. However, the factors that control petroleum expulsion and chemical fractionation in geologic settings, such as kerogen retention and relative solubilities are typically greatly perturbed by the high temperatures required for generation under observable time intervals (typically less than three days though experiments lasting up to five years have been conducted).
- Since expulsion and generation cannot be reliably simulated as a coupled process, laboratory experiments have great difficulty in resolving primary generation from secondary, potentially non-realistic, reactions. In closed system experiments, it is often very difficult to resolve whether the measured products evolve directly from kerogen decomposition or from the thermal cracking of confined products that are not normally retained in a natural system environment. The artificial confinement of metastable species also may promote condensation reactions, forming new organic solids. Open system experiments pose the opposite effect; that is, the rapid removal and quenching or detection of metastable products that may remain in natural systems and undergo secondary cracking reactions. Laboratory simulations typically fail to account accurately the retention that occurs in natural geologic systems at much lower temperatures and removal of products as kerogen matures under geologic conditions. Research has shown that all laboratory artifacts cannot be eliminated and that some models may produce unrealistic bitumen compositions.
- If laboratory experiments are inherently flawed, an accurate prediction of the thermal decomposition of kerogen and its product yields rests in constructing a theoretical framework based on fundamental principles. That is, the definition of mechanisms and kinetic rates of elementary reactions, which may be measured with high confidence, to molecular assemblages that reflect the complexity of kerogen. With such a model, it should be possible to predict the generation and expulsion under various laboratory and geologic conditions. While there have been numerous studies that formulated reaction schemes and associated kinetics to explain the results of specific laboratory experiments or field observations, there have been few attempts at constructing a comprehensive model that can account for petroleum generation and expulsion under a wide range of conditions. Some attempts have been made at constructing reaction models that account for product yields from experiments conducted under a variety of laboratory conditions and apply the model to a geologic setting. However, these models are phenomenological and rely on the extrapolation of laboratory-derived kinetics to geologic heating rates. An improved method of modeling basin performance, including predicting petroleum production, is desirable.
- A method for predicting petroleum expulsion is provided. An exemplary embodiment of the method comprises defining a chemical structure of a kerogen and identifying a plurality of reaction products of the kerogen under geologic heating rates. The exemplary method also comprises grouping the plurality of reaction products into a plurality of product lumps based on their chemical composition and predicting petroleum expulsion for each of the plurality of product lumps based on secondary cracking reactions.
- In one exemplary method, petroleum expulsion is computed at a plurality of varying heating rates. A production rate curve may be defined for each of the plurality of product lumps. A plurality of stoichiometric coefficients may be defined for each of the plurality of product lumps.
- Predicting petroleum expulsion may comprise employing a theoretical model in one exemplary embodiment of the present invention. An exemplary method may comprise comparing predicted petroleum expulsion to known data. The chemical structure may be redefined if the predicted petroleum expulsion is outside a specified range with respect to the known data. According to one exemplary embodiment, an outline of secondary cracking reactions is prepared.
- In one exemplary embodiment, predicted petroleum expulsion comprises a quantity of expelled petroleum. Moreover, predicted petroleum expulsion may additionally comprise timing data regarding petroleum expulsion.
- One exemplary embodiment of the present invention relates to a method of producing hydrocarbons from an oil and/or gas field. Such an exemplary embodiment may comprise defining a chemical structure of a kerogen deposited in the oil and/or gas field and identifying a plurality of reaction products of the kerogen under geologic heating rates. One exemplary embodiment comprises grouping the plurality of reaction products into a plurality of product lumps based on their chemical composition and predicting petroleum expulsion for each of the plurality of product lumps based on secondary cracking reactions.
- Hydrocarbons may be extracted from the oil and/or gas field using the predicted petroleum expulsion. In one exemplary method of producing hydrocarbons, petroleum expulsion is computed at a plurality of varying heating rates. A production rate curve may be defined for each of the plurality of product lumps. A plurality of stoichiometric coefficients may be defined for each of the plurality of product lumps.
- Predicting petroleum expulsion may comprise employing a theoretical model in one exemplary embodiment of the present invention. An exemplary method of producing hydrocarbons may additionally comprise comparing predicted petroleum expulsion to known data. The chemical structure may be redefined if the predicted petroleum expulsion is outside a specified range with respect to the known data. According to one exemplary embodiment of a method for producing hydrocarbons, an outline of secondary cracking reactions is prepared.
- In one exemplary embodiment of a method for producing hydrocarbons, predicted petroleum expulsion comprises a quantity of expelled petroleum. Moreover, predicted petroleum expulsion may additionally comprise timing data regarding petroleum expulsion.
- One exemplary embodiment of the present invention relates to a tangible, machine-readable medium. An exemplary tangible, machine-readable medium comprises code adapted to define a chemical structure of a kerogen and code adapted to identify a plurality of reaction products of the kerogen under geologic heating rates. The exemplary tangible, machine-readable medium may additionally comprise code adapted to group the plurality of reaction products into a plurality of product lumps based on their chemical composition and code adapted to predict petroleum expulsion for each of the plurality of product lumps based on secondary cracking reactions.
- The foregoing and other advantages of the present invention may become apparent upon reviewing the following detailed description and drawings of non-limiting examples of embodiments in which:
-
FIG. 1 , which includesFIGS. 1A and 1B , is a process flow diagram showing a method for predicting petroleum production in accordance with an exemplary embodiment of the present invention; -
FIG. 2 is a diagram that shows the modeling of expelled petroleum in accordance with an exemplary embodiment of the present invention; -
FIG. 3 is a plurality of graphs showing rate curves indicating rates of expulsion of petroleum at different heat rates; -
FIG. 4 is a graph showing a family of curves corresponding to cumulative yields for a particular fraction as a function of temperature at varying heating rates in accordance with an exemplary embodiment of the present invention; -
FIG. 5 is a graph showing a rate of expulsion expressed as a single frequency factor and a discrete spectrum of activation energies in accordance with an exemplary embodiment of the present invention; -
FIG. 6 is a diagram showing an outline of secondary cracking reactions in accordance with an exemplary embodiment of the present invention; -
FIG. 7 is a graph showing an output from a calculation of the thermal cracking of a Polars (NSO) fraction in accordance with an exemplary embodiment of the present invention; -
FIG. 8 is a graph showing a comparison of overall thermal cracking expressions in accordance with an exemplary embodiment of the present invention; -
FIG. 9 is a graph showing a comparison between results produced by an open system pyrolysis approach and an exemplary embodiment of the present invention; -
FIG. 10 is a process flow diagram showing a method for producing hydrocarbons in accordance with an exemplary embodiment of the present invention; -
FIG. 11 is a diagram of a tangible, machine-readable medium in accordance with an exemplary embodiment of the present invention; and -
FIG. 12 illustrates an exemplary computer network that may be used to perform the method for predicting composition of petroleum as disclosed herein, and is discussed in greater detail below. - In the following detailed description section, the specific embodiments of the present invention are described in connection with preferred embodiments. However, to the extent that the following description is specific to a particular embodiment or a particular use of the present invention, this is intended to be for exemplary purposes only and simply provides a description of the exemplary embodiments. Accordingly, the invention is not limited to the specific embodiments described below, but rather, it includes all alternatives, modifications, and equivalents falling within the true spirit and scope of the appended claims.
- At the outset, and for ease of reference, certain terms used in this application and their meanings as used in this context are set forth. To the extent a term used herein is not defined below, it should be given the broadest definition persons in the pertinent art have given that term as reflected in at least one printed publication or issued patent.
- As used herein, the term “basin model” refers to a simplification of the earth and its processes with the intent being to track the dynamic evolution of one or more of those processes through time. For example, the processes related to the generation and migration of hydrocarbons is commonly modeled with the intent to determine which of several possible structural culminations may be the most prospective for containing a commercial accumulation. Basin models use data from seismic, well control and knowledge of the geology of the area to construct a numerical model of the region and to track the changes in the various modeled parameters through time to reach a set of predictions that are then calibrated to the known information at the present. The model parameters are then adjusted within geologically reasonable bounds until a successful match and calibration is reached. Prediction can then be made at locations away from the calibration points.
- As used herein, the term “kerogen” refers to a solid, carbonaceous material. When kerogen is imbedded in rock formations, the mixture is referred to as oil shale. This is true whether or not the mineral is, in fact, technically shale, that is, a rock formed from compacted clay. Kerogen is subject to decomposing upon exposure to heat over a period of time. Upon heating, kerogen molecularly decomposes to produce oil, gas, and carbonaceous coke. Small amounts of water may also be generated. The oil, gas and water fluids are mobile within the rock matrix, while the carbonaceous coke remains essentially immobile. Kerogen may be classified into four distinct groups: Type I, Type II, Type III, and Type IV. Kerogen types used herein are as defined in Tissot and Welte (Tissot, B. P. and Welte, D. H., Petroleum Formation and Occurrence, second edition, Springer-Verlag, Berlin, 1984, p. 151). The maturation sequence for kerogen that typically occurs over geological time is due to burial leading to exposure to increased temperature and pressure. Classification of kerogen type may depend upon precursor materials of the kerogen. The precursor materials transform over time into macerals or amorphous masses. Macerals are microscopic structures that have distinguishing morphologies, different chemical structures and properties depending on the precursor materials from which they are derived. Amorphous kerogens have no distinguishing morphological features that can be used to characterize its precursor materials, but may have different chemical structures and properties.
- As used herein, “tangible machine-readable medium” refers to a medium that participates in directly or indirectly providing signals, instructions and/or data to a processing system. A machine-readable medium may take forms, including, but not limited to, non-volatile media (e.g., ROM, disk) and volatile media (RAM). Common forms of a machine-readable medium include, but are not limited to, a floppy disk, a flexible disk, a hard disk, a magnetic tape, other magnetic medium, a CD-ROM, other optical medium, punch cards, paper tape, other physical medium with patterns of holes, a RAM, a ROM, an EPROM, a FLASH-EPROM, or other memory chip or card, a memory stick, and other media from which a computer, a processor or other electronic device can read.
- As used herein, “NSO” or “NSOs” refers to nitrogen, sulfur, and oxygen containing compounds.
- Some portions of the detailed descriptions which follow are presented in terms of procedures, steps, logic blocks, processing and other symbolic representations of operations on data bits within a computer memory. These descriptions and representations are the means used by those skilled in the data processing arts to most effectively convey the substance of their work to others skilled in the art. In the present application, a procedure, step, logic block, process, or the like, is conceived to be a self-consistent sequence of steps or instructions leading to a desired result. The steps are those requiring physical manipulations of physical quantities. Usually, although not necessarily, these quantities take the form of electrical or magnetic signals capable of being stored, transferred, combined, compared, and otherwise manipulated in a computer system.
- It should be borne in mind, however, that all of these and similar terms are to be associated with the appropriate physical quantities and are merely convenient labels applied to these quantities. Unless specifically stated otherwise as apparent from the following discussions, it is appreciated that throughout the present application, discussions utilizing the terms such as “processing”, “computing”, “defining”, “redefining”, “identifying”, “grouping”, “preparing”, “predicting”, “comparing”, or the like, refer to the action and processes of a computer system, or similar electronic computing device, that transforms data represented as physical (electronic) quantities within the computer system's registers and memories into other data similarly represented as physical quantities within the computer system memories or registers or other such information storage, transmission or display devices. Example methods may be better appreciated with reference to flow diagrams.
- While for purposes of simplicity of explanation, the illustrated methodologies are shown and described as a series of blocks, it is to be appreciated that the methodologies are not limited by the order of the blocks, as some blocks can occur in different orders and/or concurrently with other blocks from that shown and described. Moreover, less than all the illustrated blocks may be required to implement an example methodology. Blocks may be combined or separated into multiple components. Furthermore, additional and/or alternative methodologies can employ additional, not illustrated blocks. While the figures illustrate various actions occurring in serial, it is to be appreciated that various actions could occur concurrently, substantially in parallel, and/or at substantially different points in time.
- A basin model according to an exemplary embodiment of the present invention comprises an interactive tool that geologists may use for a variety of purposes. For example, the model may be used to test the impact of various geologic scenarios such as source characteristics, thermal history, structural tectonics and the like on potential plays and prospects. Another use of the model may include assessing the charge risk associated with each of the various elements of a given petroleum system, such as determining if there is adequate source richness, sufficient source volume, adequate source maturity, and appropriate timing of generation relative to timing of trap formation, to name just a few examples. Another exemplary use of model includes estimating the volume of petroleum generated from a given volume of source rock(s), the volume of petroleum expelled from the source rock, the petroleum losses during migration, and ultimately the quantity of petroleum delivered to the trap(s).
- Depending on the sophistication of the model (and of the data upon which the model is based), the basin model may help the geologist address issues such as the timing of hydrocarbon generation relative to the timing a trap formation. The model may also be useful in predicting the volumes of hydrocarbons delivered to a trap and the hydrocarbon type (liquid vs. gas and their relative proportions) likely to be currently present in the trap. In addition, the model may be used to predict the physical state of the hydrocarbon charge (Gas-Liquid phase behavior at specific pressure-temperature) and the composition of the gas/liquids to assess economic value. In one exemplary embodiment of the present invention, the model may assist in determining the possible migration paths to the prospect. Also, the model may be used to predict the possible post-charge loss of hydrocarbons from the trap (leakage through seal, tilt and spill of trap, thermal cracking, biodegradation and the like).
- A basin model in accordance with an exemplary embodiment of the present invention may be constructed by integrating various types of geological and geochemical data. A thermal history is defined that is based on geophysical models/measurement (for example, tectonic reconstruction, basinal heat flow, well temperatures) and or geochemical measurements (for example, vitrinite reflectance data, Rock-Eval data, apatite fission track data, fluid inclusion data, oil/gas molecular and isotopic compositions or the like). The generative potential of source rock units are specified by their organic richness (for example, percentage of total organic content or % TOC, Rock-Eval HI) and their spatial distribution. A generative model (kinetics and product yields) is assigned to individual source rock units based on either direct measurement or extrapolation of organic facies. The thermal history is then applied to the kinetic model and the timing and amount of oil/gas generated is calculated over the course of the basin evolution. These calculations are then compared to models of potential migration pathways and the timing of trap and seal formation, resulting in a prediction of the petroleum accumulations.
- The prediction of the amount and composition of oil and gas generated and expelled from organic rich source rocks is a key factor in basin simulations of petroleum systems. Other aspects of basin modeling, such as migration, trap formation, and seal integrity, are equally important, but are not discussed further as they are not relevant to an exemplary embodiment of the present invention. Determining the composition of oil and gas as it is expelled from the source rocks consists of modeling three major processes. The first process comprises the manner in which kerogen (the solid organic matter within source rocks that generates oil) reacts at the specified time-temperature history experienced by the source rocks over geologic time. This includes determining the rate of the reaction (kinetics) and the composition of the products (bitumen) that is formed. The second process comprises the mechanisms responsible for expelling the bitumen from the kerogen out of the source rock. Expulsion is known to be controlled by the amount and chemistry of the kerogen. Chemical fractionation occurs whereby saturated and aromatic hydrocarbons are selectively expelled while polar compounds and asphaltenes (collectively known as NSO's for compounds that contain nitrogen, sulfur, and/or oxygen) are selectively retained. The third process comprises secondary thermal cracking of the bitumen (retained within the source rock matrix) or expelled petroleum (either along migration pathways or within reservoir rocks) may occur that further alters the composition of the petroleum.
- The timing of petroleum (oil and gas) generation from source rocks is dependent on the kinetic model. These kinetic models may be derived from laboratory experiments where a source rock, kerogen, asphaltene, or other carbonaceous material is artificially matured and the product yields are measured as a function of varying heating rate, temperature, and/or time. Presently, most models of oil and gas generation are based on measurements of pyrolyzate yields from experiments conducted on a specific sample that are then fitted to a series of parallel first-order reactions or other kinetic reaction schemes. As set forth above, these laboratory studies accelerate the temperature (heating rate) to compensate for short reaction time. Hence, models are derived from reactions occurring at temperatures from ˜250 to 650° C. over very short time periods (typically minutes to less than 3 days) and are then extrapolated to occur at ˜100 to 170° C. over millions of years. This extrapolation requires spanning rates that differ by about 14 orders of magnitude. An alternative approach to running individual source samples is to develop generic models that represent either global or basinal average kinetics for a broadly defined organic facies (for example, marine shale, terrestrial coal, and carbonate).
- Those of ordinary skill in the art will recognize the importance of selecting the appropriate kinetic model. However, there is generally no consensus as to how the best kinetic model should be determined. The fact that significant differences in pyrolysate yields and modeled results emerge depending on whether the reactions are conducted under open or closed and under hydrous or anhydrous conditions has lead to an extended and unresolved discussion of the accuracy of such models when applied to geologic thermal histories and to whether the parallel reaction concept is sufficient to capture the actual nature of thermal reactions of kerogen. Indeed, reactions are known to occur under high temperature laboratory conditions that are not the same as those that occur in sedimentary basins at lower temperatures over much greater time. Furthermore, some cases studies have shown empirically that laboratory-based kinetic models of kerogen decomposition must be artificially adjusted to confirm to observations.
- With respect to petroleum composition, kinetic models based on bulk kerogen thermal decomposition can only approximate actual product yields and do not model product composition. In a product generation approach, a separate kinetic model may be determined for each product lump that may have different frequency factors. As used herein, the term “product lump” means a group of reaction products that are similar in chemical composition and behavior. A method in accordance with an exemplary embodiment of the present invention may desirably be performed such that the individual products are quantified as a function of time/temperature and may involve complex analyses of closed-system pyrolysates or customized instrumentation. In the kerogen degradation approach, where a discrete distribution of activation energies (Ea) usually with a fixed frequency factor (A) is determined for the bulk kerogen, each individual activation energy is divided into portions representative of the lumped components. This approach is considered easier to define as the bulk kerogen decomposition kinetics can be determined by routine open-system pyrolysis methods and the products can be independently quantified.
- With respect to petroleum expulsion, it is typical for the majority of the products formed in the laboratory heating experiments to be polar (NSO-containing) compounds rather than saturated and aromatic hydrocarbons that dominate expelled petroleum. Research and field observation has shown that the bitumen retained in a source rock is enriched in polar compounds (NSOs and asphaltenes) while the expelled petroleum is enriched in hydrocarbons. Depending on the expulsion model, the compositional model may reflect either primary bitumen or expelled petroleum. A bitumen model may desirably calculate the chemical fractionation that occurs upon expulsion, while a petroleum model has already adjusted the lumped component compositions to reflect this fractionation.
- After petroleum has been expelled, it may undergo secondary thermal cracking reactions that further alter its composition. This secondary cracking may occur within the source rock mineral matrix, along migration pathways, or within reservoir rocks. Beyond simple Oil→Gas models, an exemplary embodiment of the present invention models compositional changes in petroleum resulting from secondary cracking.
-
FIG. 1 , which includesFIG. 1A and 1B , is a process flow diagram showing a method for predicting differences in subsurface conditions in accordance with an exemplary embodiment of the present invention. The method is generally referred to by thereference number 100. Atblock 102, a chemical analysis of a kerogen is performed. Atblock 104, a chemical structure, referred to as CS inFIG. 1 , is populated. A decision is made regarding whether the chemical structure has been optimized to match known analytical data, as shown atblock 106. If the chemical structure is not optimized with respect to known analytical data, process flow returns to block 104, where the chemical structure is re-populated. - At
block 108, modeling of the amount and composition of generated and expelled oil and gas is performed on the chemical structure using a process known as chemical structure-chemical yields modeling or CY CSM. One example of a CS-CYM is generally described in U.S. Pat. No. 7,344,889, entitled “Chemical Structural and Compositional Yields Model for Predicting Hydrocarbon Thermolysis Products”, which issued to Kelemen, et al. on Mar. 18, 2008. A CY-CSM is also described in the following publication: Freund H., Walters, C. C. Kelemen S. R., Siskin M., Gorboty, M. L., Curry D. J., Bence A. E., 2006. Organic Geochemistry, 38 288-305. Predicting Oil and Gas Compositional Yields via Chemical Structure-Chemical Yield Modeling (CS-CYM):Part 1—Concepts and Implementation. Another example of a CY-CSM is set forth in U.S. Provisional Patent Application No. 61/140,246, Attorney Docket No. 2008EM357 entitled “Method for Predicting Petroleum Expulsion”, by Walters, et al., filed on Dec. 23, 2008. - In an exemplary embodiment of the present invention, CY-CSM relies on constructing a representative model of the molecular structure of a kerogen using constraints from numerous chemical analyses, modeling the thermal reactions that occur within this representative structure, and applying a thermodynamic-based theory of expulsion and chemical fractionation to calculate the product yields. This model has been shown to predict the yields and composition of oil and gas generated in high temperature experiments. When basinal heating rates are used, the compositions of both expelled oil and retained bitumen are modeled following the reaction pathways, kinetics, and expulsion mechanisms that occur under these conditions and are not extrapolations of high-temperature experimental data. Yield curves for individual fractions are extracted, as shown at
block 110 ofFIG. 1A . For example, yield curves for a group ofhydrocarbons 112 may be individually extracted, as set forth below. - In performing CY-CSM in accordance with an exemplary embodiment of the present invention, by providing the amount and kinetic model for individual hydrocarbons (for example, C1, C2, C3, C4, and C5 hydrocarbons) or lumped hydrocarbons (for example, C6-C14 Aromatics, C6-C15 Saturates, C15+ Saturates, C15+ Aromatics, and C15+ polars). A kinetic model for each of the group of
hydrocarbons 112 is produced, as shown atblock 114. Because the kinetics of each expelled product are modeled using the yields calculated from a coupled generation and expulsion model, each product may be considered independently within the basin simulator. Hence, instead of using a kerogen decomposition kinetic model where the stoichiometric yields for each activation energy bin may vary, a kinetic model for each compositional component is defined. As such, the implementation is extremely flexible as new components (for example, non-HC gases) or refinements of a lumped component (for example, C15+ polars modeled as C15+ resins and C15+ asphaltenes) can be easily changed without affecting the other component models. The kinetics of the individual components are determined by running CS-CYM at varying heating rates typical of geologic basins (0.5 to 32° C./million years (Ma) are given as examples in block 108) and fitting the calculated yield curves to a discrete activation energy model. - In the exemplary embodiment shown in
FIG. 1A , additional CS-CYM modeling is performed on a subset of the group ofhydrocarbons 112, as shown atblock 116. Yield curves for the subset of the group of hydrocarbons are extracted for the resulting individual fractions, as shown atblock 118. Examples of these individual fractions include afirst fraction 120, asecond fraction 122 and athird fraction 124. As shown inFIG. 1B , thefirst fraction 120 includes input from additional ones of the group ofhydrocarbons 112. Atblock 126, kinetic models of secondary cracking are applied to a second subset of the group ofhydrocarbons 112. Atblock 128, kinetic models and stoichiometric coefficients are derived based on input data from thefirst fraction 120, thesecond fraction 122 and thethird fraction 124. - As shown at
block 130, modeling of secondary cracking reactions, which further alter expelled oil composition via thermal reactions that occur within source rock mineral matrices, along migration pathways, and within reservoir rock, also are calculated using CS-CYM. Here, the generation/expulsion of a compositional lump (for example, C15+ Saturates) is first modeled as described above, and the individual product is trapped. The trapped product is then re-run within CS-CYM under closed conditions and the product yields are calculated. -
FIG. 2 is a diagram that shows the modeling of expelled petroleum in accordance with an exemplary embodiment of the present invention. The diagram is generally referred to by thereference number 200. In an exemplary embodiment of the present invention, CS-CYM is used to derive the data needed to determine the yields and kinetic model for the coupled generation/expulsion of individual lumped products and the yields and secondary cracking kinetics of the C14+ lumped products. - A
first panel 202 of the diagram 200 shows a snapshot of aType 1 kerogen substructure. To create this snapshot, kerogen is isolated from its rock matrix and analyzed using numerous analytical methods (for example, elemental CHNOS analysis, solid-state 13C NMR, X-ray photoelectron spectroscopy, Sulfur-XANES, pyrolysis-GC or the like) that provide constraining input values used to construct a representation (˜1 million atoms) of the kerogen structure. The chemical structure of the kerogen is then placed within a reaction network that predicts the composition of the generated and expelled products under specified temperature conditions, as shown at asecond panel 204 of the diagram 200. For simulating high temperature, open system experiments such as those typically used for determining bulk kerogen degradation kinetics, CS-CYM first calculates the products produced from kerogen cracking and then determines if the product is sufficiently volatile to be considered as an expelled product. Yield curves for expelled products may be compared to known geologic data, as shown in athird panel 206 of the diagram 200. Thethird panel 206 includes a left graph and a right graph. The left graph includes a y-axis 208 that represents a maximum temperature in degrees Celsius. Anx-axis 210 of the left graph represents a heating rate in degrees Celsius per minute. The right graph of thethird panel 206 includes a y-axis 212 that represents the rate of expelled products from a simulation of an open-system laboratory pyrolysis in units of mg expelled product/g total carbon/0.2 minutes. Anx-axis 214 of the right graph represents temperature in units of degrees Celsius. The calculated yields compare very well with laboratory yields when kerogens are reacted at high temperatures and heating rates, suggesting that the reaction network is replicating the basic chemical reaction and expulsion processes that occur under these conditions. - With respect to identifying the constituent products as shown in the second panel of the diagram 204, the products replicate high temperature laboratory results. A CS-CYM program is used to determine the product yield under geologic heating rates. The generated products are calculated as several thousand individual species, then lumped into chemical classes defined as: C1, C2-C5, C6-C14 saturates, C6-C14 aromatics+polars, C15+ saturates, C15+ aromatics, C15+ polars. The C2-C5 component may be further divided into C2, C3, C4, C5 hydrocarbons if the final model is to be used in an Equation of State (EOS) calculation to determine phase behavior. C15+ polars may be split into C15+ NSO and C15+ asphaltenes based on a solubility test of individual components.
- As shown in the
second panel 204, the composition of the expelled products are calculated using a simplified formulation of a thermodynamic theory of expulsion and chemical fractionation that is an extended Flory-Rehner Regular Solution. This expulsion theory couples concepts of polymer theory to describe the amount of bitumen that can dissolve in a kerogen matrix with Regular Solution theory that explains which molecules are preferentially retained and which are expelled. A simplification of the theory is desirable as an exact calculation of the composition of the equilibrated retained and expelled is too currently computationally intensive to be run for all components at all time steps. -
FIG. 3 is a plurality of graphs showing rate curves indicating rates of expulsion of petroleum at different heat rates. The plurality of graphs is generally referred to by thereference number 300. For a given kerogen, the expelled products are modeled at varying heating rates. Typical heating rates used are 0.5, 1, 2, 4, 8, 16, and 32° C./Ma, which span the conditions most commonly seen in nature. In the exemplary plurality ofgraphs 300, a family of curves is shown for heat rates of 0.5° C./Ma, 2° C./Ma and 8° C./Ma. In particular, afirst graph 302 comprises a y-axis 304 that corresponds to a rate if expelled oil at a heating rate of 0.5° C./Ma. Thefirst graph 302 includes anx-axis 306 that corresponds to temperature in degrees Celsius. A family ofcurves 308 shows the rate of expelled oil for a plurality of different reaction products. Moreover, each curve within the family ofcurves 308 corresponds to a different reaction product. - A
second graph 310 comprises a y-axis 312 that corresponds to a rate if expelled oil at a heating rate of 2° C./Ma. Thesecond graph 310 includes anx-axis 314 that corresponds to temperature in degrees Celsius. A family ofcurves 316 shows the rate of expelled oil for a plurality of different reaction products. Moreover, each curve within the family ofcurves 316 corresponds to a different reaction product. - A
third graph 318 comprises a y-axis 320 that corresponds to a rate if expelled oil at a heating rate of 8° C./Ma. Thesecond graph 318 includes anx-axis 322 that corresponds to temperature in degrees Celsius. A family ofcurves 324 shows the rate of expelled oil for a plurality of different reaction products. Moreover, each curve within the family ofcurves 324 corresponds to a different reaction product. -
FIG. 4 is a graph showing a family of curves corresponding to cumulative yields for a particular fraction as a function of temperature at varying heating rates in accordance with an exemplary embodiment of the present invention. In an exemplary embodiment of the present invention, each individual fraction or product lump is analyzed as a separate entity. The graph is generally referred to by thereference number 400. Thegraph 400 includes a y-axis 402 that corresponds to a fraction reacted index ranging from 0.0 to 1.0. Anx-axis 404 corresponds to temperature in degrees Celsius. A family of evolution curves 406 represent cumulative yields as a function of temperature at varying heating rates. The family of evolution curves may be fitted using fitting software. From this analysis, a kinetic model is generated whereby the rate of expulsion of a specific product lump is expressed by a single frequency factor and a discrete spectrum of activation energies. -
FIG. 5 is a graph showing a rate of expulsion expressed as a single frequency factor and a discrete spectrum of activation energies in accordance with an exemplary embodiment of the present invention. The graph is generally referred to by thereference number 500. Thegraph 500 includes a y-axis 502 corresponding to a percentage rate of expulsion of a particular fraction. Anx-axis 504 corresponds to activation energy in units of kcal/mol. Thus, the bars on the graph represent a percentage of overall rate of expulsion for a typical reaction product across a spectrum of activation energies. - Once the kinetics of all compositional lumps are determined, the model is completed by determining the stoichiometric coefficients for the finally cumulative product. The amount of each lumped expelled product i that is expelled from a source kerogen is calculated as: TOC×HI×fi.×ai, where ai is a fixed coefficient. Since the kinetics of expulsion of C2-C5 are nearly identical, they can be calculated as a single C2-C5 lump. The individual C2, C3, C4, and C5 components are then separated from a fixed stoichiometric relationship (bi). The residual kerogen, termed coke, is a non-mobile product that can convert to inert char and CH4 at higher temperatures. The amount of coke is calculated by grams of kerogen×fi×acoke where acoke is a fixed coefficient.
- An example of a complete compositional yields model for expelled product is shown below in Table 1:
-
TABLE 1 Kinetic Models Type I (based on Green River Shale) Expelled Products and Coke Fraction C1 C2-C5 C6-C14Sat C6-C14Ar C15+Sat C15+Aro Polars Coke A, sec-1 Kcal/mol 5.44E+11 1.12E+14 3.40E+13 2.16E+13 4.49E+13 1.22E+13 3.52E+16 4.41E+13 47 0 0 0 0 0 0 0 0 48 0 0 0 0 0 0 0 0 49 0 0 0 13.9 0 81.44 0 0 50 100 0 87.93 10.02 87.93 1.34 0 75.86 51 0 46.02 0 0 0 15.05 0 0 52 0 5.06 0 0 0 0 0 0 53 0 0 10.24 53.33 10.24 2.16 0 17.53 54 0 48.91 1.83 9.86 1.83 0 50 6.43 55 0 0 0 8.14 0 0 50 0 56 0 0 0 4.76 0 0 0 0 57 0 0 0 0 0 0 0 0 58 0 0 0 0 0 0 0 0 59 0 0 0 0 0 0 0 0.18 indicates data missing or illegible when filed - Stoichiometric yield coefficients are shown below in Table 2.
- Secondary coefficients for the C2-C5 product lump are shown below in Table 3.
-
TABLE 3 % C2-C5. Fraction of individual HC's C2 C3 C4 C5 bi 0.15 0.25 0.3 0.3 - Coefficients for coke are shown below in Table 4.
-
TABLE 4 Coke C1 acoke 0.253051 - Coefficients relating to coke formation are shown below in Table 5.
- The discrete energy kinetic models for expelled products are readily incorporated into basin simulators as the form is similar to what is currently used for kerogen degradation kinetics. In accordance with an exemplary embodiment of the present invention, the models are for expelled product. Hence, the generation and expulsion mechanisms, along with chemical fractionation, are inherently built into the model and the simulator need not calculate generation and expulsion using different models and/or principles.
- Secondary cracking occurs after petroleum has migrated out of the source kerogen. These cracking reactions may occur to bitumen trapped in pore spaces within the source rock matrix, along migration pathways, and in reservoir rocks. Thermal cracking can significantly alter the composition of petroleum.
-
FIG. 6 is a diagram showing an outline of secondary cracking reactions in accordance with an exemplary embodiment of the present invention. The diagram is generally referred to by thereference number 600. The diagram shows a vertical column of expelledproducts 602 and a corresponding series of secondary crackingreaction products 604. - To model the secondary cracking reactions, CS-CYM may be performed in a manner that the expelled products for each chemical lump were mathematically trapped as they were expelled. These molecular pools were then subjected to a second round of heating at various heating rates under closed-system conditions and the degradation of each product lump calculated. From these rates curves, a kinetic model was fashioned in the same manner as described above. This procedure was applied to the C15+ fractions. While the models for the thermal cracking of C15+ saturated and aromatic hydrocarbons were similar regardless of the composition of the original source kerogen, the thermal cracking of the C15+ polar compounds differ. Consequently, the origin of the polar compounds must be tracked within the basin simulator. The kinetics for the thermal cracking of <C15+ fractions were derived from literature values and laboratory heating experiments.
-
FIG. 7 is a graph showing an output of mobile products from a calculation of the thermal cracking of a Polars (NSO) fraction in accordance with an exemplary embodiment of the present invention. The graph is generally referred to by thereference number 700. The Polars (NSO) fraction shown inFIG. 7 is from Green River Shale (Type I) kerogen at a 0.5° C./Ma heating rate. - The
graph 700 comprises a y-axis 702 that corresponds to NSO secondary cracking products measured in mg/gNSO. Anx-axis 704 corresponds to temperature in degrees Celsius. Each of a family ofcurves 706 corresponds to a different component of the Polars (NSO) fraction. -
FIG. 8 is a graph showing a comparison of overall thermal cracking expressions in accordance with an exemplary embodiment of the present invention. The graph is generally referred to by thereference number 800. Thegraph 800 comprises a y-axis 802 corresponding to a percentage of fractional conversion ranging from 0 to 1. Anx-axis 804 corresponds to temperature in degrees Celsius. Each of a family ofcurves 806 corresponds to a different component of a fraction. -
FIG. 9 is a graph showing a comparison between results produced by an open system pyrolysis approach and an exemplary embodiment of the present invention. The graph is generally referred to by thereference number 900. Thegraph 900 comprises a y-axis 902 corresponding to a percentage of fractional conversion ranging from 0 to 1. Anx-axis 904 corresponds to temperature in degrees Celsius. Afirst curve 906 shows results produced by an exemplary embodiment of the present invention. Asecond curve 908 shows results produced using an open system pyrolysis approach. - An exemplary embodiment of the present invention provides solutions that differ from known published or commercial products. For instance, in the example provided above for Type I kerogen, an exemplary embodiment of the present invention yields hydrocarbons at substantially lower temperatures than a conventional model derived from open-system pyrolysis (A=7.86×1013 sec−1, Emean=53 kcal/mol, std. dev.=1 kcal/mol) at a 4° C./Ma heating rate. Moreover, an exemplary embodiment of the present invention predicts petroleum is generated and expelled from a vitrinite reflectance of ˜0.55% to ˜0.95% (
LOM 8 to 10.5) with a 50% fractional conversion at ˜0.67%. The open-system pyrolysis kinetic model predicts petroleum is generated from ˜0.7% to ˜1.2% (LOM 9 to 11). The calculations were conducted using a 4° C./Ma heating rate that approximates the geothermal history of the Uinta Basin during Paleocene-Eocene subsidence. -
FIG. 10 is a process flow diagram showing a method for producing hydrocarbons in accordance with an exemplary embodiment of the present invention. The method is generally represented by thereference number 1000. Atblock 1002, the method begins. - A chemical structure of a kerogen deposited in the oil and/or gas field is defined, as shown at
block 1004. Atblock 1006, a plurality of reaction products of the kerogen under geologic heating rates is identified. The reaction products are grouped into a plurality of product lumps based on their chemical composition, as shown atblock 1008. Atblock 1010, petroleum expulsion for each of the plurality of product lumps is predicted based on secondary cracking reactions. Hydrocarbons are extracted from the oil and/or gas field using the predicted petroleum expulsion, as shown atblock 1012. Atblock 1014, the process ends. -
FIG. 11 is a diagram of a tangible, machine-readable medium in accordance with an exemplary embodiment of the present invention. The tangible, machine-readable medium is generally referred to by thereference number 1100. The tangible, machine-readable medium 110 may comprise a disk drive such as a magnetic or optical disk or the like. In an exemplary embodiment of the present invention, the tangible, machine-readable medium 1100 comprisescode 1102 adapted to define a chemical structure of a kerogen.Code 1104 adapted to identify a plurality of reaction products of the kerogen under geologic heating rates is stored on the tangible, machine-readable medium. In addition,code 1106 adapted to group the plurality of reaction products into a plurality of product lumps based on their chemical composition is stored on the tangible, machine-readable medium. Finally, the tangible, machine-readablemedium stores code 1108 adapted to predict petroleum expulsion for each of the plurality of product lumps based on secondary cracking reactions. - An exemplary embodiment of the present invention relates to a method for calculating the generation and expulsion of oil and gas and its secondary thermal cracking that differs from previously published techniques. According to an exemplary embodiment of the present invention, determination of the quantity and timing of petroleum generation is performed according to a theoretical model as opposed to extrapolation of laboratory experiments. In addition, the generation and expulsion models are coupled into a single model for expelled petroleum. A theoretical model to describe the quantity and timing of generation of hydrocarbons from the secondary cracking of hydrocarbons and polar compounds may additionally be employed.
- An exemplary embodiment of the present invention provides an approximation of an expulsion model for chemical fractionation. In addition, the stoichiometry and kinetics of the secondary cracking of individual lumps are modeled.
-
FIG. 12 illustrates anexemplary computer system 1200 on which software for performing processing operations of embodiments of the present invention may be implemented. A central processing unit (CPU) 1201 is coupled tosystem bus 1202. TheCPU 1201 may be any general-purpose CPU. The present invention is not restricted by the architecture of CPU 1201 (or other components of exemplary system 1200) as long as CPU 1201 (and other components of system 1200) supports the inventive operations as described herein. TheCPU 1201 may execute the various logical instructions according to embodiments. For example, theCPU 1201 may execute machine-level instructions for performing processing according to the exemplary operational flow described above in conjunction withFIG. 10 . For instance,CPU 1201 may execute machine-level instructions for performingoperational block 1006 ofFIG. 10 , as an example. - The
computer system 1200 also preferably includes random access memory (RAM) 1203, which may be SRAM, DRAM, SDRAM, or the like. Thecomputer system 1200 preferably includes read-only memory (ROM) 1204 which may be PROM, EPROM, EEPROM, or the like. TheRAM 1203 and theROM 1204 hold user and system data and programs, as is well-known in the art. Thecomputer system 1200 also preferably includes an input/output (I/O)adapter 1205, acommunications adapter 1211, auser interface adapter 1208, and adisplay adapter 1209. The I/O adapter 1205, theuser interface adapter 1208, and/orcommunications adapter 1211 may, in certain embodiments, enable a user to interact withcomputer system 1200 in order to input information. - The I/
O adapter 1205 preferably connects to a storage device(s) 1206, such as one or more of hard drive, compact disc (CD) drive, floppy disk drive, tape drive, etc. tocomputer system 1200. The storage devices may be utilized when theRAM 1203 is insufficient for the memory requirements associated with storing data for operations of embodiments of the present invention. The data storage of thecomputer system 1200 may be used for storing information and/or other data used or generated in accordance with embodiments of the present invention. Thecommunications adapter 1211 is preferably adapted to couple thecomputer system 1200 to anetwork 1212, which may enable information to be input to and/or output fromsystem 1200 via such network 1212 (e.g., the Internet or other wide-area network, a local-area network, a public or private switched telephony network, a wireless network, any combination of the foregoing). Theuser interface adapter 1208 couples user input devices, such as akeyboard 1213, apointing device 1207, and amicrophone 1214 and/or output devices, such as a speaker(s) 1215 to thecomputer system 1200. Thedisplay adapter 1209 is driven by theCPU 1201 to control the display on adisplay device 1210 to, for example, display information or a representation pertaining to a portion of a subsurface region under analysis, such as displaying a generated 3D representation of a target area, according to certain embodiments. - It shall be appreciated that the present invention is not limited to the architecture of
system 1200. For example, any suitable processor-based device may be utilized for implementing all or a portion of embodiments of the present invention, including without limitation personal computers, laptop computers, computer workstations, and multi-processor servers. Moreover, embodiments may be implemented on application specific integrated circuits (ASICs) or very large scale integrated (VLSI) circuits. In fact, persons of ordinary skill in the art may utilize any number of suitable structures capable of executing logical operations according to the embodiments. - While the present invention may be susceptible to various modifications and alternative forms, the exemplary embodiments discussed above have been shown only by way of example. However, it should again be understood that the invention is not intended to be limited to the particular embodiments disclosed herein. Indeed, the present invention includes all alternatives, modifications, and equivalents falling within the true spirit and scope of the appended claims.
Claims (20)
1. A method for predicting petroleum expulsion, the method comprising:
defining a chemical structure of a kerogen;
identifying a plurality of reaction products of the kerogen under geologic heating rates;
grouping the plurality of reaction products into a plurality of product lumps based on their chemical composition; and
predicting petroleum expulsion for each of the plurality of product lumps based on secondary cracking reactions.
2. The method for predicting petroleum expulsion recited in claim 1 , comprising computing petroleum expulsion at a plurality of varying heating rates.
3. The method for predicting petroleum expulsion recited in claim 1 , comprising defining a production rate curve for each of the plurality of product lumps.
4. The method for predicting petroleum expulsion recited in claim 1 , comprising defining a plurality of stoichiometric coefficients for each of the plurality of product lumps.
5. The method for predicting petroleum expulsion recited in claim 1 , wherein predicting petroleum expulsion comprises employing a theoretical model.
6. The method for predicting petroleum expulsion recited in claim 1 , comprising comparing predicted petroleum expulsion to known data.
7. The method for predicting petroleum expulsion recited in claim 6 , comprising redefining the chemical structure if the predicted petroleum expulsion is outside a specified range with respect to the known data.
8. The method for predicting petroleum expulsion recited in claim 1 , comprising preparing an outline of secondary cracking reactions.
9. The method for predicting petroleum expulsion recited in claim 1 , wherein predicted petroleum expulsion comprises a quantity of expelled petroleum.
10. The method for predicting petroleum expulsion recited in claim 1 , wherein predicted petroleum expulsion comprises timing data regarding petroleum expulsion.
11. A method for producing hydrocarbons from an oil and/or gas field, the method comprising:
defining a chemical structure of a kerogen deposited in the oil and/or gas field;
identifying a plurality of reaction products of the kerogen under geologic heating rates;
grouping the plurality of reaction products into a plurality of product lumps based on their chemical composition;
predicting petroleum expulsion for each of the plurality of product lumps based on secondary cracking reactions; and
extracting hydrocarbons from the oil and/or gas field using the predicted petroleum expulsion.
12. The method for producing hydrocarbons recited in claim 11 , comprising computing petroleum expulsion at a plurality of varying heating rates.
13. The method for producing hydrocarbons recited in claim 11 , comprising defining a production rate curve for each of the plurality of product lumps.
14. The method for producing hydrocarbons recited in claim 11 , wherein predicting petroleum expulsion comprises employing a theoretical model.
15. The method for producing hydrocarbons recited in claim 11 , comprising comparing predicted petroleum expulsion to known data.
16. The method for producing hydrocarbons recited in claim 15 , comprising redefining the chemical structure if the predicted petroleum expulsion is outside a specified range with respect to the known data.
17. The method for producing hydrocarbons recited in claim 11 , comprising preparing an outline of secondary cracking reactions.
18. The method for producing hydrocarbons recited in claim 11 , wherein predicted petroleum expulsion comprises a quantity of expelled petroleum.
19. The method for producing hydrocarbons recited in claim 11 , wherein predicted petroleum expulsion comprises timing data regarding petroleum expulsion.
20. A tangible, machine-readable medium, comprising:
code adapted to define a chemical structure of a kerogen;
code adapted to identify a plurality of reaction products of the kerogen under geologic heating rates;
code adapted to group the plurality of reaction products into a plurality of product lumps based on their chemical composition; and
code adapted to predict petroleum expulsion for each of the plurality of product lumps based on secondary cracking reactions.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/573,428 US9552462B2 (en) | 2008-12-23 | 2009-10-05 | Method for predicting composition of petroleum |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US14024908P | 2008-12-23 | 2008-12-23 | |
US12/573,428 US9552462B2 (en) | 2008-12-23 | 2009-10-05 | Method for predicting composition of petroleum |
Publications (2)
Publication Number | Publication Date |
---|---|
US20100155078A1 true US20100155078A1 (en) | 2010-06-24 |
US9552462B2 US9552462B2 (en) | 2017-01-24 |
Family
ID=42264389
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US12/573,428 Active 2031-11-15 US9552462B2 (en) | 2008-12-23 | 2009-10-05 | Method for predicting composition of petroleum |
Country Status (2)
Country | Link |
---|---|
US (1) | US9552462B2 (en) |
EP (1) | EP2202514A3 (en) |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100250468A1 (en) * | 2009-03-24 | 2010-09-30 | Spl, Inc. | System to form an actual sales or delivery value for all components of a commingled hydrocarbon fluid stream |
US20100252270A1 (en) * | 2007-12-18 | 2010-10-07 | Chul-Sung Kim | Determining Connectivity Architecture In 2-D and 3-D Heterogeneous Data |
US7895052B1 (en) * | 2009-03-24 | 2011-02-22 | Spl, Inc. | Computer instructions to form an actual sales or delivery value for all components of a commingled hydrocarbon fluid stream |
US7895051B1 (en) * | 2009-03-24 | 2011-02-22 | Spl, Inc. | Method to form an actual sales or delivery value for all components of a commingled hydrocarbon fluid stream |
US8370122B2 (en) | 2007-12-21 | 2013-02-05 | Exxonmobil Upstream Research Company | Method of predicting connectivity between parts of a potential hydrocarbon reservoir and analyzing 3D data in a subsurface region |
WO2014040622A1 (en) * | 2012-09-12 | 2014-03-20 | Bp Exploration Operating Company Limited | System and method for determining retained hydrocarbon fluid |
WO2015031088A1 (en) * | 2013-08-27 | 2015-03-05 | Schlumberger Canada Limited | Determining phase behavior of a reservoir fluid |
US20150063650A1 (en) * | 2013-08-30 | 2015-03-05 | Halliburton Energy Services, Inc. | Identifying Potential Fracture Treatment Locations In a Formation Based on Production Potential |
US9128076B2 (en) | 2010-04-30 | 2015-09-08 | Exxonmobil Upstream Research Company | Measurement of isotope ratios in complex matrices |
US9146225B2 (en) | 2011-11-11 | 2015-09-29 | Exxonmobil Upstream Research Company | Exploration method and system for detection of hydrocarbons with an underwater vehicle |
WO2017105875A1 (en) | 2015-12-15 | 2017-06-22 | Exxonmobil Upstream Research Company | Methods for the determination of biogenic gas |
US9726012B2 (en) | 2013-03-08 | 2017-08-08 | Halliburton Energy Services, Inc. | Systems and methods for optical fluid identification approximation and calibration |
US9891331B2 (en) | 2014-03-07 | 2018-02-13 | Scott C. Hornbostel | Exploration method and system for detection of hydrocarbons from the water column |
CN108609789A (en) * | 2014-09-30 | 2018-10-02 | 北京合众有益科技有限责任公司 | A kind of gas field sewage method for innocent treatment |
US10309217B2 (en) | 2011-11-11 | 2019-06-04 | Exxonmobil Upstream Research Company | Method and system for reservoir surveillance utilizing a clumped isotope and/or noble gas data |
US10415367B2 (en) | 2012-12-27 | 2019-09-17 | Halliburton Energy Services, Inc. | System and methods for estimation of intra-kerogen porosity of downhole formation samples from pyrolysis tests and basin modeling data |
US10515715B1 (en) | 2019-06-25 | 2019-12-24 | Colgate-Palmolive Company | Systems and methods for evaluating compositions |
CN110869762A (en) * | 2017-06-29 | 2020-03-06 | 沙特阿拉伯石油公司 | Pyrolysis for determining hydrocarbon expulsion efficiency of source rock |
CN112211620A (en) * | 2019-06-24 | 2021-01-12 | 中国石油化工股份有限公司 | Method and system for generating ancient oil reservoir storage capacity evaluation chart |
US20210035659A1 (en) * | 2019-07-30 | 2021-02-04 | Chengdu University Of Technology | Method for forming history of natural gas accumulation by using carbon isotopes by pyrolysis experiment |
US11093842B2 (en) | 2018-02-13 | 2021-08-17 | International Business Machines Corporation | Combining chemical structure data with unstructured data for predictive analytics in a cognitive system |
CN113406307A (en) * | 2021-07-21 | 2021-09-17 | 吉林大学 | Conversion method of resistivity index and relative permeability of argillaceous sandstone reservoir |
CN114034838A (en) * | 2020-10-21 | 2022-02-11 | 核工业北京地质研究院 | Oil-gas dissipation and sandstone-type uranium ore body space positioning method in multi-energy basin |
US11518941B2 (en) | 2015-05-20 | 2022-12-06 | Saudi Arabian Oil Company | Pyrolysis to determine hydrocarbon expulsion efficiency of hydrocarbon source rock |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11073504B2 (en) | 2015-07-27 | 2021-07-27 | Oil India Limited | Method to determine transformation of kerogen to oil/gas at different maturities using step-wise pyrolysis-gas chromatography |
US11099292B1 (en) * | 2019-04-10 | 2021-08-24 | Vinegar Technologies LLC | Method for determining the composition of natural gas liquids, mean pore-size and tortuosity in a subsurface formation using NMR |
US11846184B2 (en) | 2020-08-05 | 2023-12-19 | ExxonMobil Technology and Engineering Company | Systems and methods for predicting the composition of petroleum hydrocarbons |
US11584889B2 (en) | 2021-01-04 | 2023-02-21 | Saudi Arabian Oil Company | Synthetic source rock with tea |
US11885790B2 (en) | 2021-12-13 | 2024-01-30 | Saudi Arabian Oil Company | Source productivity assay integrating pyrolysis data and X-ray diffraction data |
Citations (72)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5040414A (en) * | 1989-06-29 | 1991-08-20 | Peter Graebner | Analyzing a hydrocarbon reservoir by determining the response of that reservoir to tidal forces |
US5586082A (en) * | 1995-03-02 | 1996-12-17 | The Trustees Of Columbia University In The City Of New York | Method for identifying subsurface fluid migration and drainage pathways in and among oil and gas reservoirs using 3-D and 4-D seismic imaging |
US5757663A (en) * | 1995-09-26 | 1998-05-26 | Atlantic Richfield Company | Hydrocarbon reservoir connectivity tool using cells and pay indicators |
US5774381A (en) * | 1992-03-04 | 1998-06-30 | Meier; Paul F. | Modeling and simulation of catalytic cracking |
US5798982A (en) * | 1996-04-29 | 1998-08-25 | The Trustees Of Columbia University In The City Of New York | Method for inverting reflection trace data from 3-D and 4-D seismic surveys and identifying subsurface fluid and pathways in and among hydrocarbon reservoirs based on impedance models |
US5835882A (en) * | 1997-01-31 | 1998-11-10 | Phillips Petroleum Company | Method for determining barriers to reservoir flow |
US5835982A (en) * | 1995-10-17 | 1998-11-10 | Lawer S.P.A. | Device for withdrawing and weighing products such as, for example, dye products in textile installations |
US5905657A (en) * | 1996-12-19 | 1999-05-18 | Schlumberger Technology Corporation | Performing geoscience interpretation with simulated data |
US6013172A (en) * | 1997-11-13 | 2000-01-11 | The University Of Chicago | Methodology for extracting local constants from petroleum cracking flows |
US6246963B1 (en) * | 1999-01-29 | 2001-06-12 | Timothy A. Cross | Method for predicting stratigraphy |
US20020013687A1 (en) * | 2000-03-27 | 2002-01-31 | Ortoleva Peter J. | Methods and systems for simulation-enhanced fracture detections in sedimentary basins |
US20020049575A1 (en) * | 2000-09-28 | 2002-04-25 | Younes Jalali | Well planning and design |
US6393906B1 (en) * | 2001-01-31 | 2002-05-28 | Exxonmobil Upstream Research Company | Method to evaluate the hydrocarbon potential of sedimentary basins from fluid inclusions |
US20020067373A1 (en) * | 2000-06-29 | 2002-06-06 | Eric Roe | System and method for defining and displaying a reservoir model |
US20020099504A1 (en) * | 1999-01-29 | 2002-07-25 | Cross Timothy A. | Method of predicting three-dimensional stratigraphy using inverse optimization techniques |
US20020120429A1 (en) * | 2000-12-08 | 2002-08-29 | Peter Ortoleva | Methods for modeling multi-dimensional domains using information theory to resolve gaps in data and in theories |
US6514915B1 (en) * | 1999-09-29 | 2003-02-04 | Baker Hughes Incorporated | Synthetic base fluid for enhancing the results of crude oil characterization analyses |
US6661000B2 (en) * | 2001-12-12 | 2003-12-09 | Exxonmobil Upstream Research Company | Method for measuring absorbed and interstitial fluids |
US20040148147A1 (en) * | 2003-01-24 | 2004-07-29 | Martin Gregory D. | Modeling in-situ reservoirs with derivative constraints |
US20040210547A1 (en) * | 2002-07-12 | 2004-10-21 | Chroma Energy, Inc. | Pattern recognition template application applied to oil exploration and production |
US6810332B2 (en) * | 2003-01-31 | 2004-10-26 | Chevron U.S.A. Inc. | Method for computing complexity, confidence and technical maturity indices for reservoir evaluations |
US20040220790A1 (en) * | 2003-04-30 | 2004-11-04 | Cullick Alvin Stanley | Method and system for scenario and case decision management |
US6826483B1 (en) * | 1999-10-13 | 2004-11-30 | The Trustees Of Columbia University In The City Of New York | Petroleum reservoir simulation and characterization system and method |
US20040254734A1 (en) * | 2003-06-02 | 2004-12-16 | Isabelle Zabalza-Mezghani | Method for optimizing production of an oil reservoir in the presence of uncertainties |
US20050096893A1 (en) * | 2003-06-02 | 2005-05-05 | Mathieu Feraille | Decision support method for oil reservoir management in the presence of uncertain technical and economic parameters |
US20050149307A1 (en) * | 2000-02-22 | 2005-07-07 | Schlumberger Technology Corporation | Integrated reservoir optimization |
US20050199391A1 (en) * | 2004-02-03 | 2005-09-15 | Cudmore Julian R. | System and method for optimizing production in an artificially lifted well |
US20050209912A1 (en) * | 2004-03-17 | 2005-09-22 | Schlumberger Technology Corporation | Method system and program storage device for automatically calculating and displaying time and cost data in a well planning system using a Monte Carlo simulation software |
US20050209866A1 (en) * | 2004-03-17 | 2005-09-22 | Schlumberger Technology Corporation | Method and apparatus and program storage device adapted for visualization of qualitative and quantitative risk assessment based on technical wellbore design and earth properties |
US6950751B2 (en) * | 2003-03-31 | 2005-09-27 | Conocophillps Company | Method and apparatus for the assimilation and visualization of information from 3D data volumes |
US20050234690A1 (en) * | 2004-04-14 | 2005-10-20 | Marc Mainguy | Method of constructing a geomechanical model of an underground zone intended to be coupled with a reservoir model |
US20050256647A1 (en) * | 2004-05-14 | 2005-11-17 | Leroy Ellis | Mud gas isotope logging interpretative process utilizing mixing lines in oil and gas drilling operations |
US6985841B2 (en) * | 2000-07-10 | 2006-01-10 | Institut Francais Du Petrole | Modelling method allowing to predict as a function of time the detailed composition of fluids produced by an underground reservoir under production |
US20060014647A1 (en) * | 2002-11-13 | 2006-01-19 | Baker Hughes Incorporated | Synthetic base fluid for enhancing the results of crude oil characterization analyses |
US20060041409A1 (en) * | 2004-08-20 | 2006-02-23 | Chevron U.S.A. Inc. | Method for making a reservoir facies model utilizing a training image and a geologically interpreted facies probability cube |
US20060047489A1 (en) * | 2004-08-30 | 2006-03-02 | Celine Scheidt | Method of modelling the production of an oil reservoir |
US20060052938A1 (en) * | 2004-08-20 | 2006-03-09 | Chevron U.S.A. Inc. | Method for creating facies probability cubes based upon geologic interpretation |
US20060092766A1 (en) * | 2004-11-03 | 2006-05-04 | Halliburton Energy Services, Inc. | Method and system for predicting production of a well |
US7124030B2 (en) * | 2004-05-14 | 2006-10-17 | Leroy Ellis | Mud gas isotope logging interpretive method in oil and gas drilling operations |
US20060235668A1 (en) * | 2005-04-14 | 2006-10-19 | Swanson Donald C | Modeling clastic reservoirs |
US20060235667A1 (en) * | 2005-04-14 | 2006-10-19 | Fung Larry S | Solution method and apparatus for large-scale simulation of layered formations |
US20060241867A1 (en) * | 2005-04-26 | 2006-10-26 | Fikri Kuchuk | System and methods of characterizing a hydrocarbon reservoir |
US20060265204A1 (en) * | 2005-04-26 | 2006-11-23 | Schlumberger Technology Corporation | Apparatus, method and system for improved reservoir simulation using a multiplicative overlapping Schwarz preconditioning for adaptive implicit linear systems |
US20060277012A1 (en) * | 2005-06-02 | 2006-12-07 | Ludovic Ricard | Method of upscaling absolute permeabilities to construct a flow simulation model |
US20060277013A1 (en) * | 2005-06-02 | 2006-12-07 | Chakib Bennis | Method for simulating fluid flows within a reservoir by means of a chimera type discretization |
US20060282243A1 (en) * | 2004-11-29 | 2006-12-14 | Chevron U.S.A. Inc. Schlumberger Technology Company | Method, system and program storage device for simulating fluid flow in a physical system using a dynamic composition based extensible object-oriented architecture |
US20060287201A1 (en) * | 2005-06-03 | 2006-12-21 | Baker Hughes Incorporated | Pore-scale geometric models for interpretation of downhole formation evaluation data |
US20060293872A1 (en) * | 2004-03-16 | 2006-12-28 | M-I L.L.C. | Three-dimensional wellbore analysis and visualization |
US20070005253A1 (en) * | 2005-06-03 | 2007-01-04 | Alexandre Fornel | Method for updating a geologic model by seismic and production data |
US20070011646A1 (en) * | 2005-06-24 | 2007-01-11 | College Of William And Mary | Parallel Decoupled Mesh Generation |
US20070013690A1 (en) * | 2005-07-18 | 2007-01-18 | Jean-Jacques Grimaud | Method, system and software for visualizing 3D models |
US20070016389A1 (en) * | 2005-06-24 | 2007-01-18 | Cetin Ozgen | Method and system for accelerating and improving the history matching of a reservoir simulation model |
US7210342B1 (en) * | 2001-06-02 | 2007-05-01 | Fluid Inclusion Technologies, Inc. | Method and apparatus for determining gas content of subsurface fluids for oil and gas exploration |
US20070143024A1 (en) * | 2005-12-12 | 2007-06-21 | Anthony Michel | Method for Determining the Acid Gas Storage Capacities of a Geologic Medium by Means of a Multiphase Reactive Transport Model |
US7249009B2 (en) * | 2002-03-19 | 2007-07-24 | Baker Geomark Llc | Method and apparatus for simulating PVT parameters |
US20070219724A1 (en) * | 2004-07-01 | 2007-09-20 | Dachang Li | Method for Geologic Modeling Through Hydrodynamics-Based Gridding (Hydro-Grids) |
US20070219725A1 (en) * | 2004-09-10 | 2007-09-20 | Tao Sun | Method For Evaluating Sedimentary Basin Properties By Numerical Modeling Of Sedimentation Processes |
US20070242564A1 (en) * | 2006-04-11 | 2007-10-18 | Devi Kattelmalvadi R S | Method of processing seismic data and method of producing a mineral hydrocarbon fluid and a computer program product |
US20070265778A1 (en) * | 2000-10-26 | 2007-11-15 | Suter James R | Method for automated management of hydrocarbon gathering systems |
US7297661B2 (en) * | 1999-09-29 | 2007-11-20 | Baker Hughes Incorporated | Synthetic base fluid for enhancing the results of crude oil characterization analyses |
US20080040086A1 (en) * | 2006-08-09 | 2008-02-14 | Schlumberger Technology Corporation | Facilitating oilfield development with downhole fluid analysis |
US7337660B2 (en) * | 2004-05-12 | 2008-03-04 | Halliburton Energy Services, Inc. | Method and system for reservoir characterization in connection with drilling operations |
US20080059140A1 (en) * | 2006-08-04 | 2008-03-06 | Elodie Salmon | Method of quantifying hydrocarbon formation and retention in a mother rock |
US7344889B2 (en) * | 2002-05-01 | 2008-03-18 | Exxonmobil Upstream Research Company | Chemical structural and compositional yields model for predicting hydrocarbon thermolysis products |
US7387021B2 (en) * | 2005-05-24 | 2008-06-17 | Baker Hughes Incorporated | Method and apparatus for reservoir characterization using photoacoustic spectroscopy |
US20080147326A1 (en) * | 2004-05-14 | 2008-06-19 | Leroy Ellis | Method and system of processing information derived from gas isotope measurements in association with geophysical and other logs from oil and gas drilling operations |
US20080173804A1 (en) * | 2007-01-19 | 2008-07-24 | Schlumberger Technology Corporation | Methods and apparatus for multi dimension fluorescence spectrum measurement and correlations downhole |
US20090071239A1 (en) * | 2007-09-13 | 2009-03-19 | Schlumberger Technology Corp. | Methods for optimizing petroleum reservoir analysis |
US7520158B2 (en) * | 2005-05-24 | 2009-04-21 | Baker Hughes Incorporated | Method and apparatus for reservoir characterization using photoacoustic spectroscopy |
US7526418B2 (en) * | 2004-08-12 | 2009-04-28 | Saudi Arabian Oil Company | Highly-parallel, implicit compositional reservoir simulator for multi-million-cell models |
US7529626B1 (en) * | 2004-05-14 | 2009-05-05 | Leroy Ellis | Method of integration and displaying of information derived from a mud gas isotope logging interpretative process in association with geophysical and other logs from oil and gas drilling operations |
US8352228B2 (en) * | 2008-12-23 | 2013-01-08 | Exxonmobil Upstream Research Company | Method for predicting petroleum expulsion |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1983003676A1 (en) | 1982-04-20 | 1983-10-27 | Bather, John, Maxwell | Analysis of kerogens |
GB8321073D0 (en) | 1983-08-04 | 1983-09-07 | Boc Group Plc | Refrigeration method |
WO2004102156A2 (en) | 2003-05-07 | 2004-11-25 | Saudi Arabian Oil Company | Compositional modeling and pyrolysis data analysis methods |
US20100214293A1 (en) | 2005-12-02 | 2010-08-26 | Koninklijke Philips Electronics, N.V. | System and method for user interation in data-driven mesh generation for parameter reconstruction from imaging data |
CN101689102B (en) | 2007-02-16 | 2014-01-29 | 沙特阿拉伯石油公司 | Method for determining volume of organic matter in reservoir rock |
WO2009094064A1 (en) | 2008-01-22 | 2009-07-30 | Exxonmobil Upstream Research Company | Dynamic connectivity analysis |
CA2726491C (en) | 2008-07-14 | 2016-09-27 | Exxonmobil Upstream Research Company | Systems and methods for determining geologic properties using acoustic analysis |
-
2009
- 2009-10-05 US US12/573,428 patent/US9552462B2/en active Active
- 2009-11-18 EP EP09176276.5A patent/EP2202514A3/en not_active Withdrawn
Patent Citations (79)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5159833A (en) * | 1989-06-29 | 1992-11-03 | Chevron Research And Technology Company | Analyzing a hydrocarbon reservoir by determining the response of that reservoir to tidal forces |
US5040414A (en) * | 1989-06-29 | 1991-08-20 | Peter Graebner | Analyzing a hydrocarbon reservoir by determining the response of that reservoir to tidal forces |
US5774381A (en) * | 1992-03-04 | 1998-06-30 | Meier; Paul F. | Modeling and simulation of catalytic cracking |
US5586082A (en) * | 1995-03-02 | 1996-12-17 | The Trustees Of Columbia University In The City Of New York | Method for identifying subsurface fluid migration and drainage pathways in and among oil and gas reservoirs using 3-D and 4-D seismic imaging |
US5757663A (en) * | 1995-09-26 | 1998-05-26 | Atlantic Richfield Company | Hydrocarbon reservoir connectivity tool using cells and pay indicators |
US5835982A (en) * | 1995-10-17 | 1998-11-10 | Lawer S.P.A. | Device for withdrawing and weighing products such as, for example, dye products in textile installations |
US5798982A (en) * | 1996-04-29 | 1998-08-25 | The Trustees Of Columbia University In The City Of New York | Method for inverting reflection trace data from 3-D and 4-D seismic surveys and identifying subsurface fluid and pathways in and among hydrocarbon reservoirs based on impedance models |
US5905657A (en) * | 1996-12-19 | 1999-05-18 | Schlumberger Technology Corporation | Performing geoscience interpretation with simulated data |
US5835882A (en) * | 1997-01-31 | 1998-11-10 | Phillips Petroleum Company | Method for determining barriers to reservoir flow |
US6013172A (en) * | 1997-11-13 | 2000-01-11 | The University Of Chicago | Methodology for extracting local constants from petroleum cracking flows |
US6246963B1 (en) * | 1999-01-29 | 2001-06-12 | Timothy A. Cross | Method for predicting stratigraphy |
US20020099504A1 (en) * | 1999-01-29 | 2002-07-25 | Cross Timothy A. | Method of predicting three-dimensional stratigraphy using inverse optimization techniques |
US6754588B2 (en) * | 1999-01-29 | 2004-06-22 | Platte River Associates, Inc. | Method of predicting three-dimensional stratigraphy using inverse optimization techniques |
US7297661B2 (en) * | 1999-09-29 | 2007-11-20 | Baker Hughes Incorporated | Synthetic base fluid for enhancing the results of crude oil characterization analyses |
US6514915B1 (en) * | 1999-09-29 | 2003-02-04 | Baker Hughes Incorporated | Synthetic base fluid for enhancing the results of crude oil characterization analyses |
US6826483B1 (en) * | 1999-10-13 | 2004-11-30 | The Trustees Of Columbia University In The City Of New York | Petroleum reservoir simulation and characterization system and method |
US20050149307A1 (en) * | 2000-02-22 | 2005-07-07 | Schlumberger Technology Corporation | Integrated reservoir optimization |
US20070156377A1 (en) * | 2000-02-22 | 2007-07-05 | Gurpinar Omer M | Integrated reservoir optimization |
US20020013687A1 (en) * | 2000-03-27 | 2002-01-31 | Ortoleva Peter J. | Methods and systems for simulation-enhanced fracture detections in sedimentary basins |
US20020067373A1 (en) * | 2000-06-29 | 2002-06-06 | Eric Roe | System and method for defining and displaying a reservoir model |
US6985841B2 (en) * | 2000-07-10 | 2006-01-10 | Institut Francais Du Petrole | Modelling method allowing to predict as a function of time the detailed composition of fluids produced by an underground reservoir under production |
US20020049575A1 (en) * | 2000-09-28 | 2002-04-25 | Younes Jalali | Well planning and design |
US20070265778A1 (en) * | 2000-10-26 | 2007-11-15 | Suter James R | Method for automated management of hydrocarbon gathering systems |
US20020120429A1 (en) * | 2000-12-08 | 2002-08-29 | Peter Ortoleva | Methods for modeling multi-dimensional domains using information theory to resolve gaps in data and in theories |
US6393906B1 (en) * | 2001-01-31 | 2002-05-28 | Exxonmobil Upstream Research Company | Method to evaluate the hydrocarbon potential of sedimentary basins from fluid inclusions |
US7395691B2 (en) * | 2001-06-02 | 2008-07-08 | Fluid Inclusion Technologies, Inc. | Method and apparatus for determining gas content of subsurface fluids for oil and gas exploration |
US7210342B1 (en) * | 2001-06-02 | 2007-05-01 | Fluid Inclusion Technologies, Inc. | Method and apparatus for determining gas content of subsurface fluids for oil and gas exploration |
US6661000B2 (en) * | 2001-12-12 | 2003-12-09 | Exxonmobil Upstream Research Company | Method for measuring absorbed and interstitial fluids |
US7249009B2 (en) * | 2002-03-19 | 2007-07-24 | Baker Geomark Llc | Method and apparatus for simulating PVT parameters |
US7344889B2 (en) * | 2002-05-01 | 2008-03-18 | Exxonmobil Upstream Research Company | Chemical structural and compositional yields model for predicting hydrocarbon thermolysis products |
US20040210547A1 (en) * | 2002-07-12 | 2004-10-21 | Chroma Energy, Inc. | Pattern recognition template application applied to oil exploration and production |
US20060014647A1 (en) * | 2002-11-13 | 2006-01-19 | Baker Hughes Incorporated | Synthetic base fluid for enhancing the results of crude oil characterization analyses |
US20040148147A1 (en) * | 2003-01-24 | 2004-07-29 | Martin Gregory D. | Modeling in-situ reservoirs with derivative constraints |
US6810332B2 (en) * | 2003-01-31 | 2004-10-26 | Chevron U.S.A. Inc. | Method for computing complexity, confidence and technical maturity indices for reservoir evaluations |
US6950751B2 (en) * | 2003-03-31 | 2005-09-27 | Conocophillps Company | Method and apparatus for the assimilation and visualization of information from 3D data volumes |
US20040220790A1 (en) * | 2003-04-30 | 2004-11-04 | Cullick Alvin Stanley | Method and system for scenario and case decision management |
US20050096893A1 (en) * | 2003-06-02 | 2005-05-05 | Mathieu Feraille | Decision support method for oil reservoir management in the presence of uncertain technical and economic parameters |
US20040254734A1 (en) * | 2003-06-02 | 2004-12-16 | Isabelle Zabalza-Mezghani | Method for optimizing production of an oil reservoir in the presence of uncertainties |
US20050199391A1 (en) * | 2004-02-03 | 2005-09-15 | Cudmore Julian R. | System and method for optimizing production in an artificially lifted well |
US20060293872A1 (en) * | 2004-03-16 | 2006-12-28 | M-I L.L.C. | Three-dimensional wellbore analysis and visualization |
US20050209912A1 (en) * | 2004-03-17 | 2005-09-22 | Schlumberger Technology Corporation | Method system and program storage device for automatically calculating and displaying time and cost data in a well planning system using a Monte Carlo simulation software |
US20050209866A1 (en) * | 2004-03-17 | 2005-09-22 | Schlumberger Technology Corporation | Method and apparatus and program storage device adapted for visualization of qualitative and quantitative risk assessment based on technical wellbore design and earth properties |
US20050234690A1 (en) * | 2004-04-14 | 2005-10-20 | Marc Mainguy | Method of constructing a geomechanical model of an underground zone intended to be coupled with a reservoir model |
US20080099241A1 (en) * | 2004-05-12 | 2008-05-01 | Halliburton Energy Services, Inc., A Delaware Corporation | Characterizing a reservoir in connection with drilling operations |
US7337660B2 (en) * | 2004-05-12 | 2008-03-04 | Halliburton Energy Services, Inc. | Method and system for reservoir characterization in connection with drilling operations |
US20080097735A1 (en) * | 2004-05-12 | 2008-04-24 | Halliburton Energy Services, Inc., A Delaware Corporation | System for predicting changes in a drilling event during wellbore drilling prior to the occurrence of the event |
US7529626B1 (en) * | 2004-05-14 | 2009-05-05 | Leroy Ellis | Method of integration and displaying of information derived from a mud gas isotope logging interpretative process in association with geophysical and other logs from oil and gas drilling operations |
US7124030B2 (en) * | 2004-05-14 | 2006-10-17 | Leroy Ellis | Mud gas isotope logging interpretive method in oil and gas drilling operations |
US20080147326A1 (en) * | 2004-05-14 | 2008-06-19 | Leroy Ellis | Method and system of processing information derived from gas isotope measurements in association with geophysical and other logs from oil and gas drilling operations |
US7174254B2 (en) * | 2004-05-14 | 2007-02-06 | Leroy Ellis | Mud gas isotope logging interpretative process utilizing mixing lines in oil and gas drilling operations |
US20050256647A1 (en) * | 2004-05-14 | 2005-11-17 | Leroy Ellis | Mud gas isotope logging interpretative process utilizing mixing lines in oil and gas drilling operations |
US20070219724A1 (en) * | 2004-07-01 | 2007-09-20 | Dachang Li | Method for Geologic Modeling Through Hydrodynamics-Based Gridding (Hydro-Grids) |
US7526418B2 (en) * | 2004-08-12 | 2009-04-28 | Saudi Arabian Oil Company | Highly-parallel, implicit compositional reservoir simulator for multi-million-cell models |
US20060041409A1 (en) * | 2004-08-20 | 2006-02-23 | Chevron U.S.A. Inc. | Method for making a reservoir facies model utilizing a training image and a geologically interpreted facies probability cube |
US20060052938A1 (en) * | 2004-08-20 | 2006-03-09 | Chevron U.S.A. Inc. | Method for creating facies probability cubes based upon geologic interpretation |
US20060047489A1 (en) * | 2004-08-30 | 2006-03-02 | Celine Scheidt | Method of modelling the production of an oil reservoir |
US20070219725A1 (en) * | 2004-09-10 | 2007-09-20 | Tao Sun | Method For Evaluating Sedimentary Basin Properties By Numerical Modeling Of Sedimentation Processes |
US20060092766A1 (en) * | 2004-11-03 | 2006-05-04 | Halliburton Energy Services, Inc. | Method and system for predicting production of a well |
US20060282243A1 (en) * | 2004-11-29 | 2006-12-14 | Chevron U.S.A. Inc. Schlumberger Technology Company | Method, system and program storage device for simulating fluid flow in a physical system using a dynamic composition based extensible object-oriented architecture |
US20060235667A1 (en) * | 2005-04-14 | 2006-10-19 | Fung Larry S | Solution method and apparatus for large-scale simulation of layered formations |
US20060235668A1 (en) * | 2005-04-14 | 2006-10-19 | Swanson Donald C | Modeling clastic reservoirs |
US20060265204A1 (en) * | 2005-04-26 | 2006-11-23 | Schlumberger Technology Corporation | Apparatus, method and system for improved reservoir simulation using a multiplicative overlapping Schwarz preconditioning for adaptive implicit linear systems |
US20060241867A1 (en) * | 2005-04-26 | 2006-10-26 | Fikri Kuchuk | System and methods of characterizing a hydrocarbon reservoir |
US7520158B2 (en) * | 2005-05-24 | 2009-04-21 | Baker Hughes Incorporated | Method and apparatus for reservoir characterization using photoacoustic spectroscopy |
US7387021B2 (en) * | 2005-05-24 | 2008-06-17 | Baker Hughes Incorporated | Method and apparatus for reservoir characterization using photoacoustic spectroscopy |
US20060277012A1 (en) * | 2005-06-02 | 2006-12-07 | Ludovic Ricard | Method of upscaling absolute permeabilities to construct a flow simulation model |
US20060277013A1 (en) * | 2005-06-02 | 2006-12-07 | Chakib Bennis | Method for simulating fluid flows within a reservoir by means of a chimera type discretization |
US20060287201A1 (en) * | 2005-06-03 | 2006-12-21 | Baker Hughes Incorporated | Pore-scale geometric models for interpretation of downhole formation evaluation data |
US20070005253A1 (en) * | 2005-06-03 | 2007-01-04 | Alexandre Fornel | Method for updating a geologic model by seismic and production data |
US20070016389A1 (en) * | 2005-06-24 | 2007-01-18 | Cetin Ozgen | Method and system for accelerating and improving the history matching of a reservoir simulation model |
US20070011646A1 (en) * | 2005-06-24 | 2007-01-11 | College Of William And Mary | Parallel Decoupled Mesh Generation |
US20070013690A1 (en) * | 2005-07-18 | 2007-01-18 | Jean-Jacques Grimaud | Method, system and software for visualizing 3D models |
US20070143024A1 (en) * | 2005-12-12 | 2007-06-21 | Anthony Michel | Method for Determining the Acid Gas Storage Capacities of a Geologic Medium by Means of a Multiphase Reactive Transport Model |
US20070242564A1 (en) * | 2006-04-11 | 2007-10-18 | Devi Kattelmalvadi R S | Method of processing seismic data and method of producing a mineral hydrocarbon fluid and a computer program product |
US20080059140A1 (en) * | 2006-08-04 | 2008-03-06 | Elodie Salmon | Method of quantifying hydrocarbon formation and retention in a mother rock |
US20080040086A1 (en) * | 2006-08-09 | 2008-02-14 | Schlumberger Technology Corporation | Facilitating oilfield development with downhole fluid analysis |
US20080173804A1 (en) * | 2007-01-19 | 2008-07-24 | Schlumberger Technology Corporation | Methods and apparatus for multi dimension fluorescence spectrum measurement and correlations downhole |
US20090071239A1 (en) * | 2007-09-13 | 2009-03-19 | Schlumberger Technology Corp. | Methods for optimizing petroleum reservoir analysis |
US8352228B2 (en) * | 2008-12-23 | 2013-01-08 | Exxonmobil Upstream Research Company | Method for predicting petroleum expulsion |
Cited By (40)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100252270A1 (en) * | 2007-12-18 | 2010-10-07 | Chul-Sung Kim | Determining Connectivity Architecture In 2-D and 3-D Heterogeneous Data |
US8365831B2 (en) | 2007-12-18 | 2013-02-05 | Exxonmobil Upstream Research Company | Determining connectivity architecture in 2-D and 3-D heterogeneous data |
US8370122B2 (en) | 2007-12-21 | 2013-02-05 | Exxonmobil Upstream Research Company | Method of predicting connectivity between parts of a potential hydrocarbon reservoir and analyzing 3D data in a subsurface region |
US7895052B1 (en) * | 2009-03-24 | 2011-02-22 | Spl, Inc. | Computer instructions to form an actual sales or delivery value for all components of a commingled hydrocarbon fluid stream |
US7895134B2 (en) * | 2009-03-24 | 2011-02-22 | Spl, Inc. | System to form an actual sales or delivery value for all components of a commingled hydrocarbon fluid stream |
US7895051B1 (en) * | 2009-03-24 | 2011-02-22 | Spl, Inc. | Method to form an actual sales or delivery value for all components of a commingled hydrocarbon fluid stream |
US20100250468A1 (en) * | 2009-03-24 | 2010-09-30 | Spl, Inc. | System to form an actual sales or delivery value for all components of a commingled hydrocarbon fluid stream |
US9128076B2 (en) | 2010-04-30 | 2015-09-08 | Exxonmobil Upstream Research Company | Measurement of isotope ratios in complex matrices |
US10527601B2 (en) | 2011-11-11 | 2020-01-07 | Exxonmobil Upstream Research Company | Method for determining the location, size, and fluid composition of a subsurface hydrocarbon accumulation |
US9146225B2 (en) | 2011-11-11 | 2015-09-29 | Exxonmobil Upstream Research Company | Exploration method and system for detection of hydrocarbons with an underwater vehicle |
US10330659B2 (en) | 2011-11-11 | 2019-06-25 | Exxonmobil Upstream Research Company | Method for determining the location, size, and fluid composition of a subsurface hydrocarbon accumulation |
US9612231B2 (en) | 2011-11-11 | 2017-04-04 | Exxonmobil Upstream Research Company | Exploration method and system for detection of hydrocarbons |
US10309217B2 (en) | 2011-11-11 | 2019-06-04 | Exxonmobil Upstream Research Company | Method and system for reservoir surveillance utilizing a clumped isotope and/or noble gas data |
US20150212235A1 (en) * | 2012-09-12 | 2015-07-30 | Bp Exploration Operating Company Limited | System and method for determining retained hydrocarbon fluid |
US10012764B2 (en) * | 2012-09-12 | 2018-07-03 | Bp Exploration Operating Company Limited | System and method for determining retained hydrocarbon fluid |
RU2604565C2 (en) * | 2012-09-12 | 2016-12-10 | Бп Эксплорейшн Оперейтинг Компани Лимитед | System and method for determining amount of retained hydrocarbon fluid |
WO2014040622A1 (en) * | 2012-09-12 | 2014-03-20 | Bp Exploration Operating Company Limited | System and method for determining retained hydrocarbon fluid |
US10415367B2 (en) | 2012-12-27 | 2019-09-17 | Halliburton Energy Services, Inc. | System and methods for estimation of intra-kerogen porosity of downhole formation samples from pyrolysis tests and basin modeling data |
US9726012B2 (en) | 2013-03-08 | 2017-08-08 | Halliburton Energy Services, Inc. | Systems and methods for optical fluid identification approximation and calibration |
US10088597B2 (en) | 2013-08-27 | 2018-10-02 | Schlumberger Technology Corporation | Determining phase behavior of a reservoir fluid |
WO2015031088A1 (en) * | 2013-08-27 | 2015-03-05 | Schlumberger Canada Limited | Determining phase behavior of a reservoir fluid |
US9250173B2 (en) * | 2013-08-30 | 2016-02-02 | Halliburton Energy Services, Inc. | Identifying potential fracture treatment locations in a formation based on production potential |
US20150063650A1 (en) * | 2013-08-30 | 2015-03-05 | Halliburton Energy Services, Inc. | Identifying Potential Fracture Treatment Locations In a Formation Based on Production Potential |
US9891331B2 (en) | 2014-03-07 | 2018-02-13 | Scott C. Hornbostel | Exploration method and system for detection of hydrocarbons from the water column |
CN108609789A (en) * | 2014-09-30 | 2018-10-02 | 北京合众有益科技有限责任公司 | A kind of gas field sewage method for innocent treatment |
CN108622961A (en) * | 2014-09-30 | 2018-10-09 | 北京合众有益科技有限责任公司 | A kind of gas field sewage method for innocent treatment |
CN108622962A (en) * | 2014-09-30 | 2018-10-09 | 北京合众有益科技有限责任公司 | A kind of gas field sewage method for innocent treatment |
CN108773875A (en) * | 2014-09-30 | 2018-11-09 | 北京合众有益科技有限责任公司 | A kind of gas field sewage method for innocent treatment |
US11518941B2 (en) | 2015-05-20 | 2022-12-06 | Saudi Arabian Oil Company | Pyrolysis to determine hydrocarbon expulsion efficiency of hydrocarbon source rock |
US10544442B2 (en) | 2015-12-15 | 2020-01-28 | Exxonmobil Upstream Research Company | Methods for the determination of biogenic gas |
WO2017105875A1 (en) | 2015-12-15 | 2017-06-22 | Exxonmobil Upstream Research Company | Methods for the determination of biogenic gas |
CN110869762A (en) * | 2017-06-29 | 2020-03-06 | 沙特阿拉伯石油公司 | Pyrolysis for determining hydrocarbon expulsion efficiency of source rock |
US11093842B2 (en) | 2018-02-13 | 2021-08-17 | International Business Machines Corporation | Combining chemical structure data with unstructured data for predictive analytics in a cognitive system |
US11100413B2 (en) * | 2018-02-13 | 2021-08-24 | International Business Machines Corporation | Combining chemical structure data with unstructured data for predictive analytics in a cognitive system |
CN112211620A (en) * | 2019-06-24 | 2021-01-12 | 中国石油化工股份有限公司 | Method and system for generating ancient oil reservoir storage capacity evaluation chart |
US10515715B1 (en) | 2019-06-25 | 2019-12-24 | Colgate-Palmolive Company | Systems and methods for evaluating compositions |
US20210035659A1 (en) * | 2019-07-30 | 2021-02-04 | Chengdu University Of Technology | Method for forming history of natural gas accumulation by using carbon isotopes by pyrolysis experiment |
US11705221B2 (en) * | 2019-07-30 | 2023-07-18 | Chengdu University Of Technology | Method for forming history of natural gas accumulation by using carbon isotopes by pyrolysis experiment |
CN114034838A (en) * | 2020-10-21 | 2022-02-11 | 核工业北京地质研究院 | Oil-gas dissipation and sandstone-type uranium ore body space positioning method in multi-energy basin |
CN113406307A (en) * | 2021-07-21 | 2021-09-17 | 吉林大学 | Conversion method of resistivity index and relative permeability of argillaceous sandstone reservoir |
Also Published As
Publication number | Publication date |
---|---|
EP2202514A3 (en) | 2017-03-22 |
US9552462B2 (en) | 2017-01-24 |
EP2202514A2 (en) | 2010-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9552462B2 (en) | Method for predicting composition of petroleum | |
US8352228B2 (en) | Method for predicting petroleum expulsion | |
US10360994B2 (en) | Method of quantifying hydrocarbon formation and retention in a source rock | |
CA2853295C (en) | Method for determining the presence and location of a subsurface hydrocarbon accumulation and the origin of the associated hydrocarbons | |
US11237147B2 (en) | Methods for using isotopic signatures to determine characteristics of hydrocarbon sources | |
Chen et al. | Quick evaluation of source rock kerogen kinetics using hydrocarbon pyrograms from regular rock-eval analysis | |
US10494923B2 (en) | Method to perform hydrocarbon system analysis for exploration, production and development of hydrocarbons | |
US9890617B2 (en) | Method to determine the presence of source rocks and the timing and extent of hydrocarbon generation for exploration, production and development of hydrocarbons | |
US20160084817A1 (en) | Method to Enhance Exploration, Development and Production of Hydrocarbons Using Multiply Substituted Isotopologue Geochemistry, Basin Modeling and Molecular Kinetics | |
US10416349B2 (en) | Phase predictions using geochemical data | |
EP3195022B1 (en) | Method to enhance exploration, development and production of hydrocarbons using multiply substituted isotopologue geochemistry, basin modeling and molecular kinetics | |
Wood | Thermal maturity and burial history modelling of shale is enhanced by use of Arrhenius time-temperature index and memetic optimizer | |
Walters et al. | Predicting oil and gas compositional yields via chemical structure–chemical yield modeling (CS–CYM): Part 2–Application under laboratory and geologic conditions | |
Yu et al. | Analysis of factors of productivity of tight conglomerate reservoirs based on random forest algorithm | |
US11846184B2 (en) | Systems and methods for predicting the composition of petroleum hydrocarbons | |
Hamdi et al. | Importance of multiple-contact and swelling tests for huff-n-puff simulations: A montney shale example | |
Liao et al. | Experimental study on fractional compositions of residual oil from shale and coal of china using grain-based MSSV pyrolysis | |
Acunzo et al. | NC92Soil: A computer code for deterministic and stochastic 1D equivalent linear seismic site response analyses | |
Xu et al. | Hydrocarbon generation kinetics and expulsion models of the Triassic Chang 9 lacustrine shales, Ordos Basin, China: Implications for the tight sandstone oil occurrence | |
Allouache et al. | Developing Equations of State for Unconventional Source Rock Fluids | |
Hazra et al. | Kerogen’s Potential to Be Converted into Petroleum: Reaction Kinetics and Modelling Thermal Maturity Plus Petroleum Transformation Processes | |
Rodríguez et al. | Estimation Coalbed Methane Content from Elemental Analysis for" Los Cuervos" Formations High Volatile Bituminous C Coals Applying Nonparametric Regression Techniques |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1551); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 4 |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 8TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1552); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 8 |