The disclosure relates to the determination of molecular structures of organic molecules in a geological formation through analysis of nuclear magnetic resonance measurements.
The determination of molecular structures of organic molecules in a geological formation, such as kerogen or coal, is known. Some of these conventional determinations are based on nuclear magnetic resonance measurements taken at the geological formation.
A first step in a convention technique a high-resolution 13C CP/MAS spectrum of a given sample is obtained, and then a plurality of different cutoffs are applied to the spectrum in discrete processes in order to determine the mole fraction of different sets of molecular structures and dynamics (e.g., aromatic and aliphatic carbons; methyl group, ethylene/methine group, and methoxy group in aliphatic regime; carbonyl carbons and aromatic carbons; aldehyde/ketone and acid/ester amide groups; protonated and non-protonated aromatic rings; and bridge heads, alkyl-attached aromatic rings, and phenols/phenolic esters; and/or other structures). These discrete processes can time-consuming and costly with respect to computing and/or human resources. Conventional techniques may assume NMR spectra comes from a pure sample that all molecules have same structure. Conventional techniques typically do not measure and/or use the shapes of static and dynamic parameter distributions, which characterize the heterogeneity of molecular structures and variation of molecular mobility.
One aspect of the disclosure relates to a computer-implemented method of determining molecular structures of organic molecules in a geological formation. The method is implemented in a computer system comprising one or more physical processors configured to execute computer program modules. The method comprises obtaining a plurality of chemical shift spectra derived from nuclear magnetic resonance measurements taken at the geological volume of interest for a plurality of different lengths of a mixing time, wherein a given chemical shift spectrum has been derived from nuclear magnetic resonance measurements for a mixing time that corresponds to the given chemical shift spectrum; applying a transformation to the plurality of chemical shift spectra that yields a multi-dimensional distribution of intensity as a function of chemical shift and one or more dynamic parameters of the nuclear magnetic resonance measurements, wherein the transformation is applied individually to separate sets of chemical shift measurements having common chemical shift values such that the transformation comprises for a given chemical shift value: obtaining, from across the plurality of chemical shift spectra, the measurements of intensity at the given chemical shift value and the values of the mixing time at which the measurements were obtained; and inverting the measurements of intensity at the given chemical shift to provide intensity as a function of the one or more dynamic parameters for the given chemical shift, wherein such inversion comprises performance of a Laplace transform; and determining molecular structures of the organic molecules in the geological formation based on an analysis of the multi-dimensional distribution of intensity and its shape.
A system configured to determine molecular structures of organic molecules in a geological formation. The system comprising one or more processors configured to execute computer program modules. The computer program modules comprising a spectrum module, a transformation module, and an analysis module. The spectrum module is configured to obtain a plurality of chemical shift spectra derived from nuclear magnetic resonance measurements taken at the geological volume of interest for a plurality of different lengths of a mixing time, wherein a given chemical shift spectrum has been derived from nuclear magnetic resonance measurements for a mixing time that corresponds to the given chemical shift spectrum. The transformation module is configured to apply a transformation to the plurality of chemical shift spectra that yields a multi-dimensional distribution of intensity as a function of chemical shift and one or more dynamic parameters of the nuclear magnetic resonance measurements. The transformation module is configured such that the transformation is applied individually to separate sets of chemical shift measurements having common chemical shift values such that the transformation comprises for a given chemical shift value: obtaining, from across the plurality of chemical shift spectra, the measurements of intensity at the given chemical shift value and the values of the mixing time at which the measurements were obtained; and inverting the measurements of intensity at the given chemical shift to provide intensity as a function of the one or more dynamic parameters for the given chemical shift, wherein such inversion comprises performance of a Laplace transform. The analysis module is configured to determine molecular structures of the organic molecules in the geological formation based on an analysis of the multi-dimensional distribution of intensity.
Non-transitory electronic storage media storing processor executable instructions configured to cause one or more processors to execute a method of determining molecular structures of organic molecules in a geological formation. The method comprising obtaining a plurality of chemical shift spectra derived from nuclear magnetic resonance measurements taken at the geological volume of interest for a plurality of different lengths of a mixing time, wherein a given chemical shift spectrum has been derived from nuclear magnetic resonance measurements for a mixing time that corresponds to the given chemical shift spectrum; applying a transformation to the plurality of chemical shift spectra that yields a multi-dimensional distribution of intensity as a function of chemical shift and one or more dynamic parameters of the nuclear magnetic resonance measurements, wherein the transformation is applied individually to separate sets of chemical shift measurements having common chemical shift values such that the transformation comprises for a given chemical shift value: obtaining, from across the plurality of chemical shift spectra, the measurements of intensity at the given chemical shift value and the values of the mixing time at which the measurements were obtained; and inverting the measurements of intensity at the given chemical shift to provide intensity as a function of the one or more dynamic parameters for the given chemical shift, wherein such inversion comprises performance of a Laplace transform; and determining molecular structures of the organic molecules in the geological formation based on an analysis of the multi-dimensional distribution of intensity.
These and other objects, features, and characteristics of the system and/or method disclosed herein, as well as the methods of operation and functions of the related elements of structure and the combination of parts and economies of manufacture, will become more apparent upon consideration of the following description and the appended claims with reference to the accompanying drawings, all of which form a part of this specification, wherein like reference numerals designate corresponding parts in the various figures. It is to be expressly understood, however, that the drawings are for the purpose of illustration and description only and are not intended as a definition of the limits of the invention. As used in the specification and in the claims, the singular form of “a”, “an”, and “the” include plural referents unless the context clearly dictates otherwise.
The present technology may be described and implemented in the general context of a system and computer methods to be executed by a computer. Such computer-executable instructions may include programs, routines, objects, components, data structures, and computer software technologies that can be used to perform particular tasks and process abstract data types. Software implementations of the present technology may be coded in different languages for application in a variety of computing platforms and environments. It will be appreciated that the scope and underlying principles of the present technology are not limited to any particular computer software technology.
Moreover, those skilled in the art will appreciate that the present technology may be practiced using any one or combination of hardware and software configurations, including but not limited to a system having single and/or multi-processer computer processors system, hand-held devices, programmable consumer electronics, mini-computers, mainframe computers, and the like. The technology may also be practiced in distributed computing environments where tasks are performed by servers or other processing devices that are linked through one or more data communications networks. In a distributed computing environment, program modules may be located in both local and remote computer storage media including memory storage devices.
Also, an article of manufacture for use with a computer processor, such as a CD, pre-recorded disk or other equivalent devices, may include a computer program storage medium and program means recorded thereon for directing the computer processor to facilitate the implementation and practice of the present technology. Such devices and articles of manufacture also fall within the spirit and scope of the present technology.
Referring now to the drawings, embodiments of the present technology will be described. The technology can be implemented in numerous ways, including for example as a system (including a computer processing system), a method (including a computer implemented method), an apparatus, a computer readable medium, a computer program product, a graphical user interface, a web portal, or a data structure tangibly fixed in a computer readable memory. Several embodiments of the present technology are discussed below. The appended drawings illustrate only typical embodiments of the present technology and therefore are not to be considered limiting of its scope and breadth.
At an operation 12, nuclear magnetic resonance data for the geological formation is obtained. Obtaining the nuclear magnetic resonance data may include performing a nuclear magnetic resonance measurement, deriving information from a nuclear magnetic resonance measurement, accessing previously stored nuclear magnetic resonance data, receiving user input of nuclear magnetic resonance data, and/or other techniques for obtaining nuclear magnetic resonance data.
Nuclear magnetic resonance data is generated by performing nuclear magnetic resonance measurement on a geological formation. Nuclear magnetic resonance measurement typically involves two phases, polarization and acquisition. During polarization, one or more magnetic fields are applied to the geological formation to orient the hydrogen atoms within the geological formation. This may include applying a static magnetic field to provide an initial orientation to the hydrogen atoms, followed by the application of one or more dynamic magnetic fields (e.g., an oscillating magnetic field). The dynamic magnetic field(s) cause the hydrogen atoms to precess, which produces a signal. The decay of this signal is what is measured by the nuclear magnetic resonance measurement. This measurement may include the determination of one or more dynamic parameters that characterize the connectivity and/or mobility of atoms or nuclei in a molecule. The dynamic parameters may include, for example, one or more of dipolar dephasing time constant, a rotating-frame longitudinal relaxation time, a cross polarization transfer time constant, and/or other dynamic parameters. Their values reflect the mobility of molecules and their shapes denote the variation of the mobility due to variation of molecular structures. The nuclear magnetic resonance data obtained at operation 12 may include measurements, and/or information derived therefrom, taken for a plurality of different lengths of mixing times. For example, measurements may be taken with different contacting times during the mixing period, different dipolar dephasing times during the mixing period, and/or other different mixing times. The nuclear magnetic resonance measurements may be taken in a dipolar dephasing experiment under magic angle spinning framework, in a regular cross polarization experiment under magic angle spinning framework, and/or in other nuclear magnetic resonance experiments and/or frameworks.
At an operation 14, a plurality of chemical shift spectra are obtained. Individual ones of the chemical shift spectra are derived from nuclear magnetic resonance measurements taken at the geological formation for a plurality of different of mixing times. The chemical shift spectra indicates intensity as a function of chemical shift. The chemical shift corresponds to the resonant frequency of a nucleus relative to a standard, as indicated by the nuclear magnetic resonance data. Chemical shift values are used as a fingerprinting of a molecule in the structural determination due to its uniqueness. And its variation of single peak indicates the variation of similar molecular structures.
By way of illustration,
As another illustration,
Returning to
By way of illustration,
As another illustration,
Returning to
By way of illustration,
Returning to
At an operation 32, for a current chemical shift value, intensities and corresponding mixing times are obtained across the plurality of chemical shift spectra. This includes obtaining, for a given chemical shift spectrum, the intensity in the given chemical shift spectrum at the current chemical shift value and the length of the mixing time corresponding to the given chemical shift spectrum.
At an operation 34, the intensities and lengths of mixing time at the current chemical shift value are inverted to provide intensity as a function of one or more dynamic parameters that characterize the connectivity and mobility of atoms and/or nuclei in a molecule. The one or more dynamic parameters may characterize a specific type of nucleus in a molecule. The inversion may include performance of a Laplace transform on the intensities and lengths of mixing time. The dynamic parameters may include, for example, one or more of dipolar dephasing time constant, a rotating-frame longitudinal relaxation time, a cross polarization transfer time constant, and/or other dynamic parameters.
At an operation 36, a determination is made as to whether there are additional chemical shift values for which operations 32 and 34 have not been performed. Responsive to a determination that there are additional chemical shift values to be processed, method 30 returns to operation 32 for a next chemical shift value. Responsive to a determination that there are no additional chemical shift values to be processed, method 30 proceeds to an operation 38.
At operation 38, the intensities for the individual chemical shift values that are determined as a function of the one or more dynamic parameters are used to generate a multi-dimensional distribution of the intensities parameterized by chemical shift and the one or more dynamic parameters. As was discussed above, illustrative examples of such distributions are shown in
The operations of methods 10 and 30 presented herein are intended to be illustrative. In some embodiments, method 10 and/or 30 may be accomplished with one or more additional operations not described, and/or without one or more of the operations discussed. Additionally, the order in which the operations of method 10 and 30 are illustrated in
In some embodiments, method 10 and/or 30 may be implemented in one or more processing devices (e.g., a digital processor, an analog processor, a digital circuit designed to process information, an analog circuit designed to process information, a state machine, and/or other mechanisms for electronically processing information). The one or more processing devices may include one or more devices executing some or all of the operations of method 10 and/or 30 in response to instructions stored electronically on an electronic storage medium. The one or more processing devices may include one or more devices configured through hardware, firmware, and/or software to be specifically designed for execution of one or more of the operations of method 10 and/or 30.
Processor 42 is configured to execute computer program modules. The computer program modules may include one or more of a measurement module 48, a spectrum module 50, a transform module 52, a cutoff map module 54, a structure module 56, and/or other modules.
Measurement module 48 is configured to obtain nuclear magnetic resonance data for the geological formation. In some implementations, measurement module 48 is configured to provide some or all of the functionality associated herein with operation 12 (shown in
Spectrum module 50 is configured to obtain a plurality of chemical shift spectra derived from nuclear magnetic resonance measurements. A set of chemical shift spectra are derived from nuclear magnetic resonance measurements taken at the geological formation for a plurality of different lengths of a mixing time. In some implementations, spectrum module 50 is configured to provide some or all of the functionality associated herein with operation 14 (shown in
Transform module 52 is configured to apply a transformation to the plurality of chemical shift spectra. The transformation yields a multi-dimensional distribution of intensities as a function of chemical shift and one or more dynamic parameters that characterize the connectivity and/or mobility of atoms or nuclei in a molecule. Transform module 52 is configured such that the transformation is applied individually to separate sets of chemical shift measurements having common chemical shift values. In some implementations, transform module 52 is configured to provide some or all of the functionality associated herein with operation 16 (shown in
Cutoff map module 54 is configured to obtain a cutoff map for the multi-dimensional distribution of intensities generated by transform module 52. The cutoff map may associate specific regions within the multi-dimensional distribution with specific molecular structures. An elevated intensity and shape within a given region indicates the presence of molecular structures associated with the given region in the organic molecules in the geological formation. In some implementations, cutoff map module 54 is configured to provide some or all of the functionality associated with operation 18 (shown in
Structure module 56 is configured to determine molecular structures of the organic molecules in the geological formation based on an analysis of the multi-dimensional distribution of intensity and shape generated by transform module 52 and/or the cutoff map obtained by cutoff map module 54. In some implementations, structure module 56 is configured to provide some or all of the functionality associated with operation 24 (shown in
Processor 42 is configured to provide information processing capabilities in system 40. As such, processor 42 may include one or more of a digital processor, an analog processor, a digital circuit designed to process information, an analog circuit designed to process information, a state machine, and/or other mechanisms for electronically processing information. Although processor 42 is shown in
It should be appreciated that although modules 48, 50, 52, 54, and/or 56 are illustrated in
Electronic storage 46 may comprise electronic storage media that electronically stores information. The electronic storage media of electronic storage 46 may include one or both of system storage that is provided integrally (i.e., substantially non-removable) with system 40 and/or removable storage that is removably connectable to system 40 via, for example, a port (e.g., a USB port, a firewire port, etc.) or a drive (e.g., a disk drive, etc.). Electronic storage 46 may include one or more of optically readable storage media (e.g., optical disks, etc.), magnetically readable storage media (e.g., magnetic tape, magnetic hard drive, floppy drive, etc.), electrical charge-based storage media (e.g., EEPROM, RAM, etc.), solid-state storage media (e.g., flash drive, etc.), and/or other electronically readable storage media. The electronic storage 46 may include one or more virtual storage resources (e.g., cloud storage, a virtual private network, and/or other virtual storage resources). Electronic storage 46 may store software algorithms, information determined by processor 42, information received through user interface 44, and/or other information that enables system 40 to function as described herein.
Although the system(s) and/or method(s) of this disclosure have been described in detail for the purpose of illustration based on what is currently considered to be the most practical and preferred implementations, it is to be understood that such detail is solely for that purpose and that the disclosure is not limited to the disclosed implementations, but, on the contrary, is intended to cover modifications and equivalent arrangements that are within the spirit and scope of the appended claims. For example, it is to be understood that the present disclosure contemplates that, to the extent possible, one or more features of any implementation can be combined with one or more features of any other implementation.