CN104428662B - 用于确定地质地层中的分子结构的系统和方法 - Google Patents

用于确定地质地层中的分子结构的系统和方法 Download PDF

Info

Publication number
CN104428662B
CN104428662B CN201380032782.0A CN201380032782A CN104428662B CN 104428662 B CN104428662 B CN 104428662B CN 201380032782 A CN201380032782 A CN 201380032782A CN 104428662 B CN104428662 B CN 104428662B
Authority
CN
China
Prior art keywords
chemical shift
intensity
magnetic resonance
nuclear magnetic
conversion
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.)
Expired - Fee Related
Application number
CN201380032782.0A
Other languages
English (en)
Other versions
CN104428662A (zh
Inventor
M·O·维甘德
孙伯勤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Chevron USA Inc
Original Assignee
Chevron USA Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Chevron USA Inc filed Critical Chevron USA Inc
Publication of CN104428662A publication Critical patent/CN104428662A/zh
Application granted granted Critical
Publication of CN104428662B publication Critical patent/CN104428662B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01RMEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00Arrangements or instruments for measuring magnetic variables
    • G01R33/20Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/46NMR spectroscopy
    • G01R33/4633Sequences for multi-dimensional NMR
    • AHUMAN NECESSITIES
    • A45HAND OR TRAVELLING ARTICLES
    • A45CPURSES; LUGGAGE; HAND CARRIED BAGS
    • A45C11/00Receptacles for purposes not provided for in groups A45C1/00-A45C9/00
    • A45C11/20Lunch or picnic boxes or the like
    • AHUMAN NECESSITIES
    • A47FURNITURE; DOMESTIC ARTICLES OR APPLIANCES; COFFEE MILLS; SPICE MILLS; SUCTION CLEANERS IN GENERAL
    • A47JKITCHEN EQUIPMENT; COFFEE MILLS; SPICE MILLS; APPARATUS FOR MAKING BEVERAGES
    • A47J47/00Kitchen containers, stands or the like, not provided for in other groups of this subclass; Cutting-boards, e.g. for bread
    • A47J47/14Carriers for prepared human food
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G01N24/081Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity

Landscapes

  • Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Health & Medical Sciences (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Geology (AREA)
  • Food Science & Technology (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

确定地质地层中的有机分子的分子结构。该有机分子可以包括:油母岩质、煤炭和/或其它有机分子。具体来说,所实现的技术可以操作以将核磁共振数据转换成多维空间,其准许通过比较跨该多维空间的强度信息与该空间的截断图来识别分子结构。这不仅可以简化有机分子的分子结构的识别,而且还使用针对混合物样本的精确数学模型,以导出结构和动态参数两者以及它们的变化。

Description

用于确定地质地层中的分子结构的系统和方法
技术领域
本公开涉及通过核磁共振测量分析来确定地质地层中的有机分子的分子结构。
背景技术
已知确定诸如油母岩质或煤炭的地质地层中的有机分子的分子结构。这些常规确定中的一些基于在该地质地层中采取的核磁共振测量。
常规技术中的第一步骤,获取给定样本的高分辨率13C CP/MAS谱,接着在离散的处理中将多个不同截断(cutoff)应用至该谱,以便确定不同的分子结构和动力学的组的摩尔分数(例如,芳香碳和脂肪碳;脂肪族状态下的甲基、乙烯/次甲基、以及甲氧基;羟基碳和芳香碳;醛/酮和酸/酯酰胺基;质子化和非质子化芳环;以及桥头酯、烷基连接芳环酯,以及苯酚/酚酯;和/或其它结构)。这些离散的处理在计算和/或人力资源方面可能耗时且昂贵。常规技术可以假定NMR谱来自所有分子都具有相同结构的纯样本。常规技术通常不测量和/或使用静态和动态参数分布形状,该形状表征分子结构的异质性和分子迁移率的变化。
发明内容
本公开的一方面涉及确定地质地层中的有机分子的分子结构的计算机实现方法。该方法在包括被配置成执行计算机程序模块的一个或多个物理处理器的计算机系统中实现。所述方法包括以下步骤:针对混合时间的多个不同长度,获得从在关注地质体积处采取的核磁共振测量结果所导出的多个化学位移谱,其中,已经针对对应于一给定化学位移谱的混合时间,从核磁共振测量结果导出了所述给定化学位移谱;向所述多个化学位移谱应用变换,所述变换产生作为化学位移和所述核磁共振测量结果的一个或多个动态参数的函数的强度的多维分布,其中,所述变换被单独地应用至具有共同化学位移值的化学位移测量结果的分离的集合,以使所述变换针对一给定化学位移值包括:从跨所述多个化学位移谱获取在所述给定化学位移值和获取所述测量结果的混合时间值下的强度测量结果;并且反演所述给定化学位移处的所述强度测量结果,以提供作为所述化学位移的所述一个或多个动态参数的函数的强度,其中,这种反演包括执行拉普拉斯变换;以及基于强度的多维分布及其形状的分析,确定所述地质地层中的有机分子的分子结构。
一种系统,被配置成确定地质地层中的有机分子的分子结构。该系统包括:一个或多个处理器,该一个或多个处理器被配置成执行计算机程序模块。该计算机程序模块包括:谱模块、变换模块、以及分析模块。该谱模块被配置成,针对混合时间的多个不同长度,获取从在关注地质体积处采取的核磁共振测量结果所导出的多个化学位移谱,其中,已经从针对对应于一给定化学位移谱的混合时间的核磁共振测量结果导出了所述给定化学位移谱。该变换模块被配置成,向产生作为化学位移和所述核磁共振测量结果的一个或多个动态参数的函数的强度的多维分布的所述多个化学位移谱应用变换。所述变换模块被配置成,使得所述变换被分别应用至分离的具有共同化学位移值的化学位移测量结果集,以使所述变换针对一给定化学位移值包括:跨所述多个化学位移谱获取所述给定化学位移值下的强度测量结果、和获取所述测量结果的混合时间的值;并且反演所述给定化学位移处的所述强度测量结果,以提供针对所述给定化学位移的作为一个或多个动态参数的函数的强度,其中,这种反演包括执行拉普拉斯变换。该分析模块被配置成基于强度的多维分布的分析,确定所述地质地层中的所述有机分子的分子结构。
一种存储有处理器可执行指令的非暂时电子存储介质,该处理器可执行指令被配置成使一个或多个处理器执行确定地质地层中的有机分子的分子结构的方法。该方法包括以下步骤:针对混合时间的多个不同长度,获取从在关注地质体积处采取的核磁共振测量结果所导出的多个化学位移谱,其中,已经针对对应于一给定化学位移谱的混合时间,从核磁共振测量结果导出了所述给定化学位移谱;向产生强度的多维分布的所述多个化学位移谱应用变换,所述强度的多维分布是化学位移和所述核磁共振测量结果的一个或多个动态参数的函数,其中,所述变换被分别地应用至分离的具有共同化学位移值的化学位移测量结果集,以使所述变换针对一给定化学位移值包括:跨所述多个化学位移谱获取所述给定化学位移值下的强度测量结果和获取所述测量结果的混合时间的值;并且反演所述给定化学位移下的所述强度测量结果,以提供作为所述给定化学位移的一个或多个动态参数的函数的强度,其中,这种反演包括执行拉普拉斯变换;以及基于强度的多维分布分析,确定所述地质地层中的所述有机分子的分子结构。
当参照附图考虑下面的描述和所附权利要求书时,在此公开的系统和/或方法的这些和其它目的、特征以及特性,连同操作方法和结构与组合部分的相关部件的功能以及制造的经济性将变得更清楚,附图形成了本说明书的一部分,其中,相同标号指定各个图中的对应部分。然而,应当明白,附图仅仅是出于例示和描述的目的,而非旨在作为对本发明的限制的解说。如在本说明书和权利要求书中使用的,单数形式“一(a)”、“一个(an)”,以及“该/所述(the)”包括多个指示物,除非上下文另外清楚地规定。
附图说明
图1例示了一种确定地质地层中的有机分子的分子结构的方法。
图2例示了化学位移谱的多个标绘图,并且每一个谱对应于混合时间的不同长度。
图3例示了化学位移谱的多个标绘图,并且每一个谱对应于混合时间的不同长度。
图4例示了强度的二维分布。
图5例示了强度的三维分布。
图6例示了针对通过化学位移和偶极去相时间常数参数化的强度的二维分布的截断图。
图7例示了强度的二维分布。
图8例示了将化学位移信息变换成通过化学位移和一个或多个动态参数所参数化的强度的多维分布的方法。
图9例示了被配置成确定地质地层中的有机分子的分子结构的系统。
具体实施方式
本技术可以按一系统和要通过一计算机执行的计算机方法的一般背景来描述和实现。这种计算机可执行指令可以包括程序、例程、对象,组件、数据结构、以及可以被用于执行特定任务和处理抽象数据类型的计算机软件技术。本技术的软件实现可以按用于多种计算平台和环境中的应用的不同语言来编码。应当清楚,本技术的范围和基本原理不限于任何特定计算机软件技术。
而且,本领域技术人员应当清楚,本技术可以利用硬件和软件配置中的任一个或组合来实践,包括但不限于,具有单一和/或多处理器计算机处理器系统的系统、手持式装置、可编程消费类电子设备、迷你计算机、大型计算机等。本技术还可以在其中通过经由一个或多个数据通信网络链接的服务器或其它处理装置执行任务的分布式计算环境中实践。在分布式计算环境中,程序模块可以位于包括存储器存储装置的本地和远程计算机存储介质两者中。
而且,用于与计算机处理器一起使用的制造品(如CD、预记录盘或其它等同装置)可以包括计算机程序存储介质和记录在其上的、用于引导计算机处理器以易于实现和实践本技术的程序装置。这种装置和制造品也落入本技术的精神和范围内。
下面,参照附图,对本技术的实施例进行描述。本技术可以按许多方式来实现,例如包括系统(包括计算机处理系统)、方法(包括计算机实现方法)、装置、计算机可读介质、计算机程序产品、图形用户接口、门户网站、或者有形地固定在计算机可读存储器中的数据结构。下面,对本技术的几个实施例进行讨论。附图仅例示了本技术的典型实施例,并由此,不应被视为对其范围和宽度的限制。
图1例示了确定地质地层中的有机分子的分子结构的方法10。该有机分子可以包括:油母岩质、煤炭、和/或其它有机分子。具体来说,方法10可以操作以将核磁共振数据转换成多维空间,其准许通过比较跨该多维空间的强度和形状信息与该空间的截断图来识别分子结构。这不仅可以简化有机分子的分子结构的识别,而且还使用针对混合物样本的精确数学模型,以导出结构和动态参数两者加上它们的变化。然而,现有技术仅使用针对纯样本的数学方程,而不考虑任何变化。
在操作12,获取地质地层的核磁共振数据。获取核磁共振数据可以包括:执行核磁共振测量,从核磁共振测量结果导出信息,访问先前存储的核磁共振数据,接收核磁共振数据的用户输入和/或用于获取核磁共振数据的其它技术。
核磁共振数据通过在地质地层上执行核磁共振测量来生成。核磁共振测量通常涉及两个阶段,极化和获取。在极化期间,将一个或多个磁场应用至该地质地层,以定向该地质地层内的氢原子。这可以包括施加静态磁场,以向氢原子提供初始取向,跟着施加一个或多个动态磁场(例如,振荡磁场)。该动态磁场使氢原子产生进动,其产生一信号。该信号的衰变通过核磁共振测量来测量。该测量可以包括确定一个或多个动态参数,其特征化分子中的原子或原子核的连接性和/或迁移率。该动态参数例如可以包括:偶极去相时间常数、旋转坐标系纵向驰豫时间、交叉极化传递时间常数、和/或其它动态参数中的一个或多个。它们的值反映了分子的迁移率,并且它们的形状指示因分子结构的变化而造成的迁移率变化。操作12处获取的核磁共振数据可以包括针对混合时间的多个不同长度采取的测量结果、和/或由其导出的信息。例如,测量可以利用混合期期间的不同接触时间、混合时段期间的不同偶极去相时间、和/或其它不同混合时间来采取。核磁共振测量结果可以在基于魔角旋转(magic angle spinning)原理的偶极去相实验中、基于魔角旋转原理的常规交叉极化实验中、和/或其它核磁共振实验和/或方法中采取。
在操作14,获取多个化学位移谱。从在多个不同的混合时间在地质地层处采取的核磁共振测量结果导出各个化学位移谱。该化学位移谱指示作为化学位移的函数的强度。该化学位移对应于相对于一标准的核子的谐振频率,如通过核磁共振数据指示的。化学位移值因其唯一性而在结构确定中用作分子的指纹识别。并且,其单峰值的变化指示类似分子结构的变化。
通过例示的方式,图2示出了化学位移谱的多个标绘图,并且每一个谱对应于不同的混合时间。在一些实现中,图2所示的标绘图对应于从在不同偶极去相时间基于魔角旋转原理的偶极去相实验中采取的核磁共振测量结果所导出的化学位移谱。该各个谱分别对应于各个偶极去相时间。
作为另一例示,图3示出了根据不同组的核磁共振数据的化学位移谱的多个标绘图。图3所示的标绘图已经根据利用不同接触时间在魔角旋转方法下的常规交叉极化实验中采取的核磁共振测量结果导出。该各个谱分别对应于各个接触时间。
返回至图1,在操作16,向产生强度的多维分布的所述多个化学位移谱应用变换,所述强度的多维分布作为化学位移和一个或多个动态参数的函数,所述参数特征化分子中的原子或核子的连接性和/或迁移率。在此,参照图6对该变换的一些实现进行描述。
通过例示的方式,图4描绘了强度的二维分布(其中强度用阴影指示)。该二维分布根据化学位移和一动态参数来参数化。例示例中的动态参数是偶极去相时间常数。图4所示强度的分布可以是在图2描绘的化学位移谱上执行变换的结果。
作为另一例示,图5描绘了强度的三维分布(并且强度用彩色指示)。该三维分布根据化学位移和两个动态参数来参数化。这两个动态参数是旋转坐标系纵向驰豫时间和交叉极化传递时间常数。图5所示强度的分布可以是在图3描绘的化学位移谱上执行变换的结果。
返回至图1,在操作18,可以获取针对在操作16生成的强度的多维分布的截断图。该截断图可以关联该多维分布内的特定区域与特定分子结构。给定区域内的升高的强度指示在地质地层中的有机分子中存在与给定区域相关联的分子结构。获取截断图可以包括基于分布信息确定截断图(例如,后期核磁共振测量和/或对应分析、校准、和/或基于其它信息)、存取先前生成的截断图(例如,从电子存储部、通过网络、和/或从其它源)、和/或获取来自其它源的截断图。
通过例示的方式,图6示出了针对通过化学位移和偶极去相时间常数参数化的强度的二维分布的截断图20。截断图20包括多个区域22,其中各个区域22分别对应于不同分子结构(例如,如图6标注的)。
图7示出了通过化学位移和偶极去相时间常数参数化的强度的二维分布(例如,通过图1所示操作16生成的和在此描述的),并且强度用阴影指示。该标绘被截断图20交叠。截断图20的各个区域20中的升高强度指示地质地层中存在与强度升高的区域22相对应的分子结构。
返回至图1,在操作24,地质地层中存在的分子结构基于在操作16生成的多维分布和在操作18获取的截断图来识别。该识别通过确定截断图中哪些区域的强度升高来执行。检测升高强度可以自动地和/或通过可视检查人工地进行。作为一非限制例,升高强度可以根据超过阈值强度的单个强度读数、违背阈值等级的给定区域内的强度聚集(例如,平均值、加权平均值,和/或其它聚集)、和/或用于识别给定区域内的升高强度的其它分析工具来识别。
图8例示了将化学位移信息变换成通过化学位移和一个或多个动态参数所参数化的强度的多维分布的方法30。在一些实现中,方法30可以被实现为方法10的操作16(图1所示并且在此描述)。应当清楚,这不是旨在进行限制,因为方法30可以在多种其它背景下应用。针对方法30的输入包括表示按不同混合时间采取的核磁共振测量结果的多个化学位移谱,以使各个谱分别对应于不同混合时间。这例如可以包括在方法10的操作14处获取的一组化学位移谱(图1所示并且在此描述)。
在操作32,对于当前化学位移值,跨所述多个化学位移谱获取强度和对应混合时间。这包括:针对给定化学位移谱,获取当前化学位移值下该给定化学位移谱中的强度、和与该给定化学位移谱相对应的混合时间的长度。
在操作34,当前化学位移值下的强度和混合时间的长度被反演,以提供作为特征化分子中的原子和/或核子的连接性和迁移率的一个或多个动态参数的函数的强度。所述一个或多个动态参数可以特征化分子中的特定类型的核子。该反演可以包括对强度和混合时间的长度执行拉普拉斯变换。该动态参数例如可以包括:偶极去相时间常数、旋转坐标系纵向驰豫时间、交叉极化传递时间常数、和/或其它动态参数中的一个或多个。
在操作36,进行有关是否存在尚未对其执行操作32和34的另外的化学位移值的确定。响应于确定存在要处理的另外的化学位移值,方法300针对下一化学位移值返回至操作32。响应于确定不存在要处理的另外的化学位移值,方法30进行至操作38。
在操作38,被确定为所述一个或多个动态参数的函数的、针对各个化学位移值的强度被用于生成根据化学位移和所述一个或多个动态参数所参数化的强度的多维分布。如上讨论的,图4-6示出了这种分布的例示例。
在此呈现的方法10和30的操作旨在例示。在一些实施方式中,方法10和/或30可以利用未描述的一个或多个附加操作和/或不利用所讨论操作中的一个或多个来完成。另外,其中图1和8例示和在此描述的方法10和30的操作的次序不旨在进行限制。
在一些实施方式中,方法10和/或30可以在一个或多个处理装置(例如,数字处理器、模拟处理器、被设计成处理信息的数字电路、被设计成处理信息的模拟电路、状态机、和/或用于电子地处理信息的其它机构)中实现。该一个或多个处理装置可以包括响应于电子地存储在电子存储介质上的指令来执行方法10和/或30的一些或全部操作的一个或多个装置。该一个或多个处理装置可以包括通过要具体设计用于执行方法10和/或30的一个或多个操作的硬件、固件和/或软件配置的一个或多个装置。
图9例示了被配置成确定地质地层中的有机分子的分子结构的系统40。在一些实现中,系统40可以包括至少一个处理器42、用户接口44、电子存储部46、和/或其它组件中的一个或多个。
处理器42被配置成执行计算机程序模块。计算机程序模块可以包括以下中的一个或多个:测量模块48、谱模块50、变换模块52、截断图模块54、结构模块56、和/或其它模块。
测量模块48被配置成获取针对地质地层的核磁共振数据。在一些实现中,测量模块48被配置成提供在此与操作12相关联的一些或全部功能(如图1所示和在此所述的)。
谱模块50被配置成获取根据核磁共振测量结果导出的多个化学位移谱。针对混合时间的多个不同长度,从在地质地层处采取的核磁共振测量结果导出一组化学位移谱。在一些实现中,谱模块50被配置成提供在此与操作14相关联的一些或全部功能(如图1所示和在此所述的)。
变换模块52被配置成向所述多个化学位移谱应用变换。该变换产生强度的多维分布,作为化学位移和一个或多个动态参数的函数,所述动态参数特征化分子中的原子或核子的连接性和/或迁移率。变换模块52被配置成使得将变换分别地应用于分离的具有共同化学位移值的化学位移测量结果集。在一些实现中,变换模块52被配置成提供在此与操作16(如图1所示和在此所述的)和/或方法30(如图8所示和在此所述的)相关联的一些或全部功能。
该截断图模块54被配置成获取针对通过变换模块52生成的强度的多维分布的截断图。该截断图可以关联该多维分布内的特定区域与特定分子结构。给定区域内的升高强度和形状指示在地质地层中的有机分子中存在与给定区域相关联的分子结构。在一些实现中,该截断图模块54被配置成提供与操作18相关联的一些或全部功能(如图1所示和在此所述的)。
结构模块56被配置成,基于通过变换模块52生成的强度和形状的多维分布的分析和/或通过截断图模块54获取的截断图来确定地质地层中的有机分子的分子结构。在一些实现中,结构模块56被配置成提供与操作24相关联的一些或全部功能(如图1所示和在此所述的)。
处理器42被配置成在系统40中提供信息处理能力。同样地,处理器42可以包括数字处理器、模拟处理器、被设计成处理信息的数字电路、被设计成处理信息的模拟电路、状态机、和/或用于电子地处理信息的其它机构中的一个或多个。尽管处理器42在图9中被示出为单一实体,但这仅仅是出于例示性目的。在一些实现中,处理器42可以包括多个处理单元。这些处理单元可以物理地位于同一装置内,或者处理器42可以表示协作操作的多个装置的处理功能。处理器42可以被配置成通过软件,硬件,固件,软件、硬件和/或固件的某一组合,和/或用于配置处理器42上的处理能力的其它机制来执行模块48、50、52、54、和/或56。
应当清楚,尽管模块48、50、52、54和/或56在图9中被例示为共同位于单一处理单元内,但在其中处理器42包括多个处理单元的实现中,模块48、50、52、54和/或56中的一个或多个可以相对于其它模块远程定位。下面描述的对由不同模块48、50、52、54和/或56提供的功能的描述出于例示性目的,而非旨在进行限制,如模块48、50、52、54和/或56中的任一个可以提供比所述的更多或更少的功能。例如,模块48、50、52、54和/或56中的一个或多个可以消除,并且其一些或全部功能可以由模块48、50、52、54和/或56中的其它模块提供。作为另一实施例,处理器42可以被配置成执行可以执行下面归因于模块48、50、52、54和/或56中的一个的一些或全部功能的一个或多个附加模块。
电子存储部46可以包括电子地存储信息的电子存储介质。电子存储部46的电子存储介质可以包括与系统40集成地(即,基本上不可去除)设置的系统存储部和/或例如经由端口(例如,USB端口、固件端口等)或驱动器(例如,盘驱动器等)可去除地连接至系统40的可去除存储部。电子存储部46可以包括一个或多个光学可读存储介质(例如,光盘等)、磁可读存储介质(例如,磁带、磁硬盘驱动器、软盘驱动器等)、基于电荷的存储介质(例如、EEPROM、RAM等)、固态存储介质(例如,闪存驱动器等)、和/或其它电子可读存储介质。该电子存储部46可以包括一个或多个虚拟存储资源(例如,云存储部、虚拟专用网络、和/或其它虚拟存储资源)。电子存储部46可以存储软件算法、由处理器42确定的信息、经由用户接口44接收的信息、和/或使得系统40能够如在此所述地起作用的其它信息。
尽管基于当前被认为是最有用且优选的实现,出于例示的目的对本公开的系统和/或方法进行了详细描述,但要明白的是,这种细节仅用于该目的,并且本公开不限于所公开的实现,而是正相反,其旨在覆盖在所附权利要求书的精神和范围内的修改例和等同布置。例如,要明白的是,本公开设想,在尽可能的情况下,可以将任何实现的一个或多个特征与任何其它实现的一个或多个特征相组合。

Claims (12)

1.一种确定地质地层中的有机分子的分子结构的计算机实现的方法,该方法在计算机系统中实现,所述计算机系统包括被配置成执行计算机程序模块的一个或多个物理处理器,所述方法包括以下步骤:
针对多个不同长度的混合时间,获取从在关注地质体积处采取的核磁共振测量结果所导出的多个化学位移谱,其中,已针对对应于一给定化学位移谱的混合时间,从核磁共振测量结果导出了所述给定化学位移谱;
向所述多个化学位移谱应用变换,所述变换产生作为化学位移和所述核磁共振测量结果的一个或多个动态参数的函数的强度的多维分布,其中,所述变换被分别地应用至分离的具有共同化学位移值的化学位移测量结果集,以使所述变换针对一给定化学位移值包括:
跨所述多个化学位移谱获取所述给定化学位移值下的强度测量结果和获取所述测量结果的混合时间的值;并且
反演所述给定化学位移值下的所述强度测量结果,以提供作为针对所述给定化学位移值的所述一个或多个动态参数的函数的强度,其中,这种反演包括执行拉普拉斯变换;以及
基于对所述强度的多维分布及其形状的分析,确定所述地质地层中的所述有机分子的分子结构,
其中,所述多维分布是三维的,并且其中,所述一个或多个动态参数包括旋转坐标系纵向驰豫时间和交叉极化传递时间常数。
2.根据权利要求1所述的方法,其中,所述一个或多个动态参数包括偶极去相时间常数。
3.根据权利要求1所述的方法,其中,所述核磁共振测量结果是利用混合期期间的多个去相时间基于魔角旋转原理在偶极去相实验中获取的。
4.根据权利要求1所述的方法,其中,所述核磁共振测量结果是利用混合期期间的多个不同接触时间基于魔角旋转在常规交叉极化实验中获取的。
5.一种被配置成确定地质地层中的有机分子的分子结构的系统,该系统包括:
一个或多个处理器,被配置成执行计算机程序模块,所述计算机程序模块包括:
谱模块,被配置成针对多个不同长度的混合时间,获取从在关注地质体积处采取的核磁共振测量结果所导出的多个化学位移谱,其中,已经针对对应于一给定化学位移谱的混合时间,从核磁共振测量结果导出了所述给定化学位移谱;
变换模块,被配置成向所述多个化学位移谱应用变换,所述变换产生作为化学位移和所述核磁共振测量结果的一个或多个动态参数的函数的强度的多维分布,其中,所述变换模块被配置成使得所述变换被分别地应用至分离的具有共同化学位移值的化学位移测量结果集,以使所述变换针对一给定化学位移值包括:
跨所述多个化学位移谱获取所述给定化学位移值下的强度测量结果和获取所述测量结果的混合时间的值;并且
反演所述给定化学位移值下的所述强度测量结果,以提供作为针对所述给定化学位移值的一个或多个动态参数的函数的强度,其中,这种反演包括执行拉普拉斯变换;以及
分析模块,被配置成基于对所述强度的多维分布的分析,确定所述地质地层中的所述有机分子的分子结构,
其中,所述变换模块被配置成使得所述多维分布是三维的,并且使得所述一个或多个动态参数包括旋转坐标系纵向驰豫时间和交叉极化传递时间常数。
6.根据权利要求5所述的系统,其中,所述变换模块被配置成使得所述一个或多个动态参数包括偶极去相时间常数。
7.根据权利要求5所述的系统,其中,所述谱模块被配置成,使得利用混合期期间的多个去相时间,基于魔角旋转原理在偶极去相实验中获取所述核磁共振测量结果。
8.根据权利要求5所述的系统,其中,所述谱模块被配置成,使得利用混合期期间的多个不同接触时间基于魔角旋转在常规交叉极化实验中获取所述核磁共振测量结果。
9.一种确定地质地层中的有机分子的分子结构的系统,所述系统包括:
用于针对多个不同长度的混合时间,获取从在关注地质体积处采取的核磁共振测量结果所导出的多个化学位移谱的装置,其中,已针对对应于一给定化学位移谱的混合时间,从核磁共振测量结果导出了所述给定化学位移谱;
用于向所述多个化学位移谱应用变换的装置,所述变换产生作为化学位移和所述核磁共振测量结果的一个或多个动态参数的函数的强度的多维分布,其中,所述变换被分别地应用至分离的具有共同化学位移值的化学位移测量结果集,以使所述变换针对一给定化学位移值包括:
跨所述多个化学位移谱获取所述给定化学位移值下的强度测量结果和获取所述测量结果的混合时间的值;并且
反演所述给定化学位移值下的所述强度测量结果,以提供作为针对所述给定化学位移值的一个或多个动态参数的函数的强度,其中,这种反演包括执行拉普拉斯变换;以及
用于基于对所述强度的多维分布的分析,确定所述地质地层中的所述有机分子的分子结构的装置,
其中,所述多维分布是三维的,并且其中,所述一个或多个动态参数包括旋转坐标系纵向驰豫时间和交叉极化传递时间常数。
10.根据权利要求9所述的系统,其中,所述一个或多个动态参数包括偶极去相时间常数。
11.根据权利要求9所述的系统,其中,所述核磁共振测量结果是利用混合期期间的多个去相时间基于魔角旋转原理在偶极去相实验中获取的。
12.根据权利要求9所述的系统,其中,所述核磁共振测量结果是利用混合期期间的多个不同接触时间基于魔角旋转在常规交叉极化实验中获取的。
CN201380032782.0A 2012-06-22 2013-06-14 用于确定地质地层中的分子结构的系统和方法 Expired - Fee Related CN104428662B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US13/530,766 US9291690B2 (en) 2012-06-22 2012-06-22 System and method for determining molecular structures in geological formations
US13/530,766 2012-06-22
PCT/US2013/045924 WO2013192039A1 (en) 2012-06-22 2013-06-14 System and method for determining molecular structures in geological formations

Publications (2)

Publication Number Publication Date
CN104428662A CN104428662A (zh) 2015-03-18
CN104428662B true CN104428662B (zh) 2017-03-22

Family

ID=49769240

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201380032782.0A Expired - Fee Related CN104428662B (zh) 2012-06-22 2013-06-14 用于确定地质地层中的分子结构的系统和方法

Country Status (8)

Country Link
US (1) US9291690B2 (zh)
EP (1) EP2864766B1 (zh)
CN (1) CN104428662B (zh)
AU (1) AU2013277466B2 (zh)
BR (1) BR112014029422A2 (zh)
CA (1) CA2875665C (zh)
RU (1) RU2015101709A (zh)
WO (1) WO2013192039A1 (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104091504A (zh) * 2014-07-29 2014-10-08 神华集团有限责任公司 煤结构模型的构建方法
EP3329263B1 (en) * 2015-07-28 2023-10-11 ConocoPhillips Company Nmr in kinetics of hydrocarbon generation
CN105486709B (zh) * 2015-12-17 2017-11-24 中国石油大学(北京) 井下稠油分子链长的测量方法和装置
CN113848156A (zh) * 2021-09-23 2021-12-28 滨州学院 一种褐煤分子结构模型的构建与优化方法
CN114550850A (zh) * 2022-01-28 2022-05-27 厦门大学 基于机器学习模拟动态核磁共振谱的方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6027941A (en) * 1996-05-15 2000-02-22 Curagen Corporation Method for distance measurements with solid-state NMR
CA2651108A1 (en) * 2006-05-03 2007-11-15 Schlumberger Canada Limited Downhole micro magnetic resonance analyzer
CN101479619A (zh) * 2006-06-24 2009-07-08 于利希研究中心有限公司 在超低磁场中通过高分辨率nmr谱法测量同核j耦合
CN101896834A (zh) * 2007-12-12 2010-11-24 雪佛龙美国公司 从核磁共振数据获得质子密度分布
WO2011034656A2 (en) * 2009-09-16 2011-03-24 Chevron U.S.A. Inc. Method for creating a 3d rock representation using petrophysical data
CN102232198A (zh) * 2008-12-01 2011-11-02 雪佛龙美国公司 处理井眼核磁共振测井记录以增强t2分布的连续性的方法
CN103098062A (zh) * 2010-09-13 2013-05-08 雪佛龙美国公司 用于页岩气储集层中的低硫区标识的系统和方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2002088662A2 (en) * 2001-04-25 2002-11-07 Bristol-Myers Squibb Company Method of molecular structure recognition
US6833698B2 (en) 2002-05-15 2004-12-21 Chevrontexaco U.S.A. Inc. Methods of decoupling diffusion effects from relaxation times to determine properties of porous media containing fluids

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6027941A (en) * 1996-05-15 2000-02-22 Curagen Corporation Method for distance measurements with solid-state NMR
CA2651108A1 (en) * 2006-05-03 2007-11-15 Schlumberger Canada Limited Downhole micro magnetic resonance analyzer
CN101479619A (zh) * 2006-06-24 2009-07-08 于利希研究中心有限公司 在超低磁场中通过高分辨率nmr谱法测量同核j耦合
CN101896834A (zh) * 2007-12-12 2010-11-24 雪佛龙美国公司 从核磁共振数据获得质子密度分布
CN102232198A (zh) * 2008-12-01 2011-11-02 雪佛龙美国公司 处理井眼核磁共振测井记录以增强t2分布的连续性的方法
WO2011034656A2 (en) * 2009-09-16 2011-03-24 Chevron U.S.A. Inc. Method for creating a 3d rock representation using petrophysical data
CN103098062A (zh) * 2010-09-13 2013-05-08 雪佛龙美国公司 用于页岩气储集层中的低硫区标识的系统和方法

Also Published As

Publication number Publication date
CA2875665C (en) 2019-06-11
AU2013277466B2 (en) 2017-01-05
US9291690B2 (en) 2016-03-22
CN104428662A (zh) 2015-03-18
EP2864766A1 (en) 2015-04-29
WO2013192039A1 (en) 2013-12-27
RU2015101709A (ru) 2016-08-10
CA2875665A1 (en) 2013-12-27
BR112014029422A2 (pt) 2017-06-27
AU2013277466A1 (en) 2014-12-04
EP2864766B1 (en) 2016-04-27
US20130342201A1 (en) 2013-12-26

Similar Documents

Publication Publication Date Title
CN104428662B (zh) 用于确定地质地层中的分子结构的系统和方法
CN105378714A (zh) 时间序列的快速分组
Zhihong et al. Evaluation of granular particle roundness using digital image processing and computational geometry
CN109541172A (zh) 土壤属性值的计算方法及装置
Larsen Global poles of the two-loop six-point N<? format?>=<? format?> 4 SYM integrand
Nikolaou et al. Detection of early warning signals in paleoclimate data using a genetic time series segmentation algorithm
Żaba et al. Solving fractional Schrödinger-type spectral problems: Cauchy oscillator and Cauchy well
Krone et al. Interactive extraction and tracking of biomolecular surface features
Maiorino et al. Multifractal characterization of protein contact networks
Fofonov et al. Projected Field Similarity for Comparative Visualization of Multi‐Run Multi‐Field Time‐Varying Spatial Data
Madan Shape-related characteristics of age-related differences in subcortical structures
Carrea et al. MATLAB Virtual Toolbox for Retrospective Rockfall Source Detection and Volume Estimation Using 3D Point Clouds: A Case Study of a Subalpine Molasse Cliff
Dickinson et al. HydroClimATe: hydrologic and climatic analysis toolkit
Wang et al. Forecast verification and visualization based on Gaussian mixture model co‐estimation
CN109100165B (zh) 桥梁运营模态分析方法、装置、终端与计算机可读存储介质
Lamanna et al. Territorial origin of olive oil: representing georeferenced maps of olive oils by NMR profiling
Xie et al. Variable dimension fractal-based conversion method between the nuclear magnetic resonance T2 spectrum and capillary pressure curve
Bosler et al. Stride Search: a general algorithm for storm detection in high-resolution climate data
Zacharatou et al. RUBIK: efficient threshold queries on massive time series
Hanson et al. XSophe–Sophe–XeprView and Molecular Sophe: computer simulation software suites for the analysis of continuous wave and pulsed EPR and ENDOR spectra
Rahimov et al. Use of local binary pattern in texture classification of carbonate rock micro-CT images
Hou et al. Weighted iterative inversion method for T2 spectrum in nuclear magnetic resonance
Audibert et al. Reconstruction of averaging indicators for highly heterogeneous media
McKnight et al. Estimating loop length from CryoEM images at medium resolutions
Fricke et al. Magnetic resonance insights into the heterogeneous, fractal-like kinetics of chemically recyclable polymers

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170322

Termination date: 20180614

CF01 Termination of patent right due to non-payment of annual fee