CN107966465B - 一种基于三维脉冲序列的岩心核磁信号采集及反演方法 - Google Patents

一种基于三维脉冲序列的岩心核磁信号采集及反演方法 Download PDF

Info

Publication number
CN107966465B
CN107966465B CN201711177772.2A CN201711177772A CN107966465B CN 107966465 B CN107966465 B CN 107966465B CN 201711177772 A CN201711177772 A CN 201711177772A CN 107966465 B CN107966465 B CN 107966465B
Authority
CN
China
Prior art keywords
echo
amplitude
window
string signal
pulse train
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.)
Active
Application number
CN201711177772.2A
Other languages
English (en)
Other versions
CN107966465A (zh
Inventor
范宜仁
刘建宇
葛新民
邓少贵
邢东辉
吴飞
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.)
China University of Petroleum East China
Original Assignee
China University of Petroleum East China
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 China University of Petroleum East China filed Critical China University of Petroleum East China
Priority to CN201711177772.2A priority Critical patent/CN107966465B/zh
Publication of CN107966465A publication Critical patent/CN107966465A/zh
Application granted granted Critical
Publication of CN107966465B publication Critical patent/CN107966465B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/08—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using nuclear magnetic resonance
    • G01N24/081—Making measurements of geologic samples, e.g. measurements of moisture, pH, porosity, permeability, tortuosity or viscosity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/448—Relaxometry, i.e. quantification of relaxation times or spin density
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/46—NMR spectroscopy
    • G01R33/4625—Processing of acquired signals, e.g. elimination of phase errors, baseline fitting, chemometric analysis

Landscapes

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

Abstract

本发明公开一种基于三维脉冲序列的岩心核磁信号采集及反演方法，步骤如下：将岩心样品放入恒定磁场强度的磁体中，设置三维脉冲序列参数并施加于岩心样品，三维脉冲序列在时间轴上分为三个窗口，第一个窗口使用反转恢复脉冲序列，其中极化时间可调；第二个窗口使用梯度线圈给测试区域施加脉冲梯度磁场并采集第一回波串信号的幅值，其中脉冲梯度大小、窗口时间、回波间隔可调；第三个窗口使用仪器最短回波间隔的CPMG脉冲序列采集第二回波串信号的幅值。对第一回波串信号的幅值和第二回波串信号的幅值进行联合反演获得岩心中各弛豫组分的含量。本发明可同时获取岩心样品中快弛豫信号和慢弛豫信号，避免有效信号的缺失，降低了反演结果的误差。

Description

一种基于三维脉冲序列的岩心核磁信号采集及反演方法
技术领域
本发明属于核磁共振及岩石物理领域，具体地说是涉及一种基于三维脉冲序列的岩心核磁信号采集及反演方法，可应用于岩石核磁共振实验及分析。
背景技术
核磁共振技术凭借快速、无损、无侵入、无毒、只对含氢孔隙流体敏感等优点，成为岩石物理领域中的一种重要的分析测试方法。将核磁共振仪器产生的脉冲序列施加于岩心样品中，以使岩心样品对该脉冲序列产生回波串信号，依据回波串信号的幅值可分析获得岩心样品中各弛豫组分的含量。
随着能源勘探的进一步深入，有机页岩、致密砂岩等复杂及非常规储层成了研究的重点。这些储层物性差，非均质性强，孔隙结构复杂，弛豫信号弱，对传统的核磁共振技术带来巨大挑战。三维核磁共振岩心测量及分析技术以其能够同时测量并分析纵向弛豫时间、横向弛豫时间、扩散系数的优点受到广泛的重视。其中，三维脉冲序列是三维核磁共振数据采集的核心技术，其功能直接影响三维核磁共振技术的应用效果。现有的三维脉冲序列，是现有二维脉冲序列的简单组合，包括基于饱和恢复序列、多回波间隔CPMG序列的三维脉冲序列，基于饱和恢复序列、扩散编辑序列的三维脉冲序列，基于饱和恢复序列、改进CPMG序列的三维脉冲序列，基于饱和恢复序列、脉冲梯度序列、CPMG序列的三维脉冲序列。
然而，以上三维脉冲序列缺少对快弛豫信号的采集，无法兼顾岩心中快弛豫组分和慢弛豫组分的协同测量，造成有效信号的缺失。同时，现有三维核磁共振信号反演方法中，缺少对快弛豫信号的处理，无法兼顾获取岩心中快弛豫组分和慢弛豫组分的信息，造成反演结果的误差。
发明内容
针对以上问题，本发明提出一种基于三维脉冲序列的岩心核磁信号采集及反演方法。
本发明所采用的技术解决方案是：
一种基于三维脉冲序列的岩心核磁信号采集及反演方法，包括以下步骤：
S1、将岩心样品放入恒定磁场强度的磁体中；
S2、设置三维脉冲序列参数并施加于岩心样品，三维脉冲序列在时间轴上分为三个窗口，第一个窗口使用反转恢复脉冲序列，其中极化时间可调；第二个窗口使用梯度线圈给测试区域施加脉冲梯度磁场并采集第一回波串信号的幅值，其中脉冲梯度大小、窗口时间、回波间隔可调；第三个窗口使用仪器最短回波间隔的CPMG脉冲序列采集第二回波串信号的幅值；
S3、对第一回波串信号的幅值和第二回波串信号的幅值进行联合反演获得岩心中各弛豫组分的含量。
上述方法中，所述三维脉冲序列划分出的三个窗口中，第一窗口脉冲序列编辑所述岩心样品的纵向弛豫时间，第二窗口脉冲序列编辑所述岩心样品的横向弛豫时间和扩散系数，第三窗口脉冲序列编辑所述岩心样品的横向弛豫时间。
上述方法中，根据多孔介质核磁共振弛豫理论，推得所述三维脉冲序列的回波幅度，如式(1)、(2)所示，
式(1)代表第一回波串信号的幅值，式(2)代表第二回波串信号的幅值；式(1)、(2)中bis代表第s组回波串中第i个自旋回波的幅值；第s组回波串的三维脉冲序列参数中，极化时间为TWs，第二窗口的回波间隔为NE1s，脉冲梯度为Gs，回波个数为TE1s；f(T1r,T2j,Dp)为纵向弛豫时间T1r、横向弛豫时间T2j、扩散系数Dp对应的孔隙度分量；γ为氢核的旋磁比；t0为梯度脉冲持续时间；TE2为第三窗口中CPMG序列的回波间隔，是仪器可设置的最小值；NE2为第三窗口中CPMG序列的回波个数；l为纵向弛豫时间的布点数；m为横向弛豫时间的布点数；n为扩散系数的布点数。
上述方法中所提到的联合反演方法是将所述第一回波串信号的幅值和所述第二回波串信号的幅值利用式(3)联合反演：
其中，k1表示与所述第一回波串信号对应的核函数矩阵；k2表示与所述第二回波串信号对应的核函数矩阵；b1表示所述第一回波串信号的幅值；b2表示所述第二回波串信号的幅值；依据式(3)获得所述岩心样品中各弛豫组分的含量f。
在一些实施方式中，所述联合反演方法可以转换为约束最优化问题的求解，即式(4)：
其中，σ是平衡所述第一回波串信号和所述第二回波串信号噪声水平的因子；α是正则化因子；依据获得所述岩心样品中各弛豫组分的含量f；其中I代表与及维度相同的单位矩阵，上标符号T表示矩阵的转置。
本发明的有益技术效果是：
(1)本发明考虑到现有三维脉冲序列缺少对快慢弛豫组分的协同采集，提供了新型的三窗口的三维脉冲序列，在施加单次脉冲序列的过程中，在第二窗口中采集第一回波串信号的幅值，在第三窗口中采集第二回波串信号的幅值，同时获取岩心样品中快弛豫信号和慢弛豫信号，避免有效信号的缺失。
(2)本发明应用基于三维脉冲序列的联合反演方法，对第一回波串信号的幅值和第二回波串信号的幅值进行联合反演，兼顾获取岩心中快弛豫组分和慢弛豫组分的信息，降低了反演结果的误差。本发明在页岩、致密砂岩等具有弱弛豫信号、孔隙结构复杂的岩心中具有较大的应用前景。
附图说明
图1是本发明提供的一种基于三维脉冲序列的岩心核磁信号采集及反演方法的流程图；
图2是本发明提供的三维脉冲序列的示意图；
图3是某具有三弛豫组分岩心的正演模型；
图4是某具有三弛豫组分岩心的正演回波串的示意图；
图5是某具有三弛豫组分岩心的反演结果。
具体实施方式
本发明提出了一种基于三维脉冲序列的岩心核磁信号采集及反演方法，将岩心样品放入恒定磁场强度的磁体中，设置三维脉冲序列参数并施加于岩心样品，三维脉冲序列在时间轴上分为三个窗口，第一个窗口使用反转恢复脉冲序列，其中极化时间可调；第二个窗口使用梯度线圈给测试区域施加脉冲梯度磁场并采集第一回波串信号的幅值，其中脉冲梯度大小、窗口时间、回波间隔可调；第三个窗口使用仪器最短回波间隔的CPMG脉冲序列采集第二回波串信号的幅值。对第一回波串信号的幅值和第二回波串信号的幅值进行联合反演获得岩心中各弛豫组分的含量。该方法在页岩、致密砂岩等具有弱弛豫信号、孔隙结构复杂的岩心中具有较大的应用前景。
上述方法中，所述三维脉冲序列划分出的三个窗口中，第一窗口脉冲序列编辑所述岩心样品的纵向弛豫时间，第二窗口脉冲序列编辑所述岩心样品的横向弛豫时间和扩散系数，第三窗口脉冲序列编辑所述岩心样品的横向弛豫时间。
上述方法中，根据多孔介质核磁共振弛豫理论，可推得所述三维脉冲序列的回波幅度，如式(1)、(2)所示，
式(1)代表第一回波串信号的幅值，式(2)代表第二回波串信号的幅值；式(1)、(2)中bis代表第s组回波串中第i个自旋回波的幅值；第s组回波串的三维脉冲序列参数中，极化时间为TWs，第二窗口的回波间隔为NE1s，脉冲梯度为Gs，回波个数为TE1s；f(T1r,T2j,Dp)为纵向弛豫时间T1r、横向弛豫时间T2j、扩散系数Dp对应的孔隙度分量；γ为氢核的旋磁比；t0为梯度脉冲持续时间；TE2为第三窗口中CPMG序列的回波间隔，是仪器可设置的最小值；NE2为第三窗口中CPMG序列的回波个数；l为纵向弛豫时间的布点数；m为横向弛豫时间的布点数；n为扩散系数的布点数。
上述方法中所提到的联合反演方法是将所述第一回波串信号的幅值和所述第二回波串信号的幅值利用式(3)联合反演：
其中，k1表示与所述第一回波串信号对应的核函数矩阵；k2表示与所述第二回波串信号对应的核函数矩阵；b1表示所述第一回波串信号的幅值；b2表示所述第二回波串信号的幅值；依据式(3)获得所述岩心样品中各弛豫组分的含量f。
在一些实施方式中，联合反演方法可以转换为约束最优化问题的求解，即式(4)：
其中，σ是平衡所述第一回波串信号和所述第二回波串信号噪声水平的因子；α是正则化因子；依据获得所述岩心样品中各弛豫组分的含量f。
下面结合附图与具体实施方式对本发明进行说明。
如图1所示，一种基于三维脉冲序列的岩心核磁信号采集及反演方法，主要包括岩心样品核磁信号采集、岩心样品核磁信号反演两部分，该方法具体按以下步骤依次进行：
S1、将岩心样品放入恒定磁场强度的磁体中；
S2、设置三维脉冲序列参数并施加于岩心样品，三维脉冲序列在时间轴上分为三个窗口，第一个窗口使用反转恢复脉冲序列，其中极化时间可调；第二个窗口使用梯度线圈给测试区域施加脉冲梯度磁场并采集第一回波串信号的幅值，其中脉冲梯度大小、窗口时间、回波间隔可调；第三个窗口使用仪器最短回波间隔的CPMG脉冲序列采集第二回波串信号的幅值；
S3、对第一回波串信号的幅值和第二回波串信号的幅值进行联合反演获得岩心中各弛豫组分的含量。
图2是三维脉冲序列的示意图。三维脉冲序列在时间轴上分为三个窗口，分别对岩心样品有具体的编辑功能：其中第一窗口脉冲序列编辑岩心样品的纵向弛豫时间，第二窗口脉冲序列编辑岩心样品的横向弛豫时间和扩散系数，第三窗口脉冲序列编辑岩心样品的横向弛豫时间。具体地，第一个窗口使用反转恢复脉冲序列，其中极化时间可调；第二个窗口使用梯度线圈给测试区域施加脉冲梯度磁场并采集第一回波串信号的幅值，其中脉冲梯度大小、窗口时间、回波间隔可调；第三个窗口使用仪器最短回波间隔的CPMG脉冲序列采集第二回波串信号的幅值。
图3是某具有三弛豫组分岩心的正演模型。其中岩心类型为页岩，弛豫组分1代表干酪根，弛豫组分2代表有机质孔中的油，弛豫组分3代表无机孔中的油，弛豫组分1及弛豫组分2代表快弛豫组分，弛豫组分3代表慢弛豫组分。图3(a)是某具有三弛豫组分岩心正演模型的三维核磁共振图，弛豫组分1的纵向弛豫时间为2ms，横向弛豫时间为0.5ms，扩散系数为4×10-7cm2/s，占岩石总弛豫信号的30％；弛豫组分2的纵向弛豫时间为15ms，横向弛豫时间为3ms，扩散系数为4×10-6cm2/s，占岩石总弛豫信号的40％；弛豫组分3的纵向弛豫时间为150ms，横向弛豫时间为100ms，扩散系数为4×10-6cm2/s，占岩石总弛豫信号的30％。图3(b)是三弛豫组分岩心正演模型在T1-T2平面的投影图；图3(c)是三弛豫组分岩心正演模型在T1-D平面的投影图；图3(d)是三弛豫组分岩心正演模型在T2-D平面的投影图。
图4是某具有三弛豫组分岩心的正演回波串的示意图。共设置15组三维脉冲序列，其中第一窗口极化时间TW为0.10ms、0.19ms、0.37ms、0.72ms、1.39ms、2.68ms、5.18ms、10ms、19.31ms、37.28ms、71.97ms、138.95ms、268.27ms、517.95ms、1000ms；第二窗口梯度脉冲持续时间t0为10ms；回波间隔TE1分别为0.2ms、0.33ms、0.5ms、1ms、1ms、0.2ms、0.33ms、0.5ms、1ms、1ms、0.2ms、0.33ms、0.5ms、1ms、1ms；回波个数NE1分别为50、30、20、10、10、50、30、20、10、10、50、30、20、10、10；脉冲梯度为G分别为5Gs/cm、22.5Gs/cm、40Gs/cm、57.5Gs/cm、75Gs/cm、5Gs/cm、22.5Gs/cm、40Gs/cm、57.5Gs/cm、75Gs/cm、5Gs/cm、22.5Gs/cm、40Gs/cm、57.5Gs/cm、75Gs/cm；第三窗口回波间隔TE2为0.2ms；回波个数NE2为3000；氢核的旋磁比γ为2.675×104rad/s/Gs。图4为根据以上15组三维脉冲序列及图3中某具有三弛豫组分岩心的正演模型获取的回波串，其中10ms前为第一回波串信号的幅值，10ms后为第二回波串信号的幅值。第一回波串信号的幅值表征快弛豫组分的信号，如弛豫组分1及弛豫组分2；第二回波串信号的幅值表征慢弛豫组分的信号，如弛豫组分3，可见基于三维脉冲序列的岩心核磁信号采集方法，能够实现对快慢弛豫组分的协同采集。
图5是某具有三弛豫组分岩心的反演结果。根据图4中的某具有三弛豫组分岩心的正演回波串，使用联合反演方法对第一回波串信号的幅值和第二回波串信号的幅值进行联合反演获得岩心中各弛豫组分的含量，其中纵向弛豫时间的布点数l为30，横向弛豫时间的布点数m为30，扩散系数的布点数n为30。图5(a)是某具有三弛豫组分岩心反演结果的三维核磁共振图，图5(b)是三弛豫组分岩心反演结果在T1-T2平面的投影图；图5(c)是三弛豫组分岩心反演结果在T1-D平面的投影图；图5(d)是三弛豫组分岩心反演结果在T2-D平面的投影图。图中可见，弛豫组分1、弛豫组分2及弛豫组分3在三维核磁共振图及三平面的投影中清晰可辨，其相对误差为5.2％，满足解释分析精度的要求，可见基于三维脉冲序列的岩心核磁信号采集反演方法可以实现对快弛豫组分与慢弛豫组分的协同处理，同时获取岩心中快弛豫组分和慢弛豫组分的信息。
上述方式中未述及的部分采取或借鉴已有技术即可实现。
需要说明的是，在本说明书的教导下，本领域技术人员所作出的任何等同替代方式，或明显变型方式，均应在本发明的保护范围之内。

Claims (1)

1.一种基于三维脉冲序列的岩心核磁信号采集及反演方法，其特征在于包括以下步骤：
S1、将岩心样品放入恒定磁场强度的磁体中；
S2、设置三维脉冲序列参数并施加于岩心样品，三维脉冲序列在时间轴上分为三个窗口，第一个窗口使用反转恢复脉冲序列，其中极化时间可调；第二个窗口使用梯度线圈给测试区域施加脉冲梯度磁场并采集第一回波串信号的幅值，其中脉冲梯度大小、窗口时间、回波间隔可调；第三个窗口使用仪器最短回波间隔的CPMG脉冲序列采集第二回波串信号的幅值；
S3、对第一回波串信号的幅值和第二回波串信号的幅值进行联合反演获得岩心中各弛豫组分的含量；
所述三维脉冲序列划分出的三个窗口中，第一窗口脉冲序列编辑所述岩心样品的纵向弛豫时间，第二窗口脉冲序列编辑所述岩心样品的横向弛豫时间和扩散系数，第三窗口脉冲序列编辑所述岩心样品的横向弛豫时间；
根据多孔介质核磁共振弛豫理论，推得所述三维脉冲序列的回波幅度，如式(1)、(2)所示，
式(1)代表第一回波串信号的幅值，式(2)代表第二回波串信号的幅值；式(1)、(2)中bis代表第s组回波串中第i个自旋回波的幅值；第s组回波串的三维脉冲序列参数中，极化时间为TWs，第二窗口的回波间隔为NE1s，脉冲梯度为Gs，回波个数为TE1s；f(T1r,T2j,Dp)为纵向弛豫时间T1r、横向弛豫时间T2j、扩散系数Dp对应的孔隙度分量；γ为氢核的旋磁比；t0为梯度脉冲持续时间；TE2为第三窗口中CPMG序列的回波间隔，是仪器可设置的最小值；NE2为第三窗口中CPMG序列的回波个数；l为纵向弛豫时间的布点数；m为横向弛豫时间的布点数；n为扩散系数的布点数；
所述的联合反演方法是将所述第一回波串信号的幅值和所述第二回波串信号的幅值利用式(3)联合反演：
其中，k1表示与所述第一回波串信号对应的核函数矩阵；k2表示与所述第二回波串信号对应的核函数矩阵；b1表示所述第一回波串信号的幅值；b2表示所述第二回波串信号的幅值；依据式(3)获得所述岩心样品中各弛豫组分的含量f；
所述联合反演方法转换为约束最优化问题的求解，即式(4)：
其中，σ是平衡所述第一回波串信号和所述第二回波串信号噪声水平的因子；α是正则化因子；依据获得所述岩心样品中各弛豫组分的含量f，其中I代表与及维度相同的单位矩阵，上标符号T表示矩阵的转置。
CN201711177772.2A 2017-11-23 2017-11-23 一种基于三维脉冲序列的岩心核磁信号采集及反演方法 Active CN107966465B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711177772.2A CN107966465B (zh) 2017-11-23 2017-11-23 一种基于三维脉冲序列的岩心核磁信号采集及反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711177772.2A CN107966465B (zh) 2017-11-23 2017-11-23 一种基于三维脉冲序列的岩心核磁信号采集及反演方法

Publications (2)

Publication Number Publication Date
CN107966465A CN107966465A (zh) 2018-04-27
CN107966465B true CN107966465B (zh) 2019-12-03

Family

ID=61999749

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711177772.2A Active CN107966465B (zh) 2017-11-23 2017-11-23 一种基于三维脉冲序列的岩心核磁信号采集及反演方法

Country Status (1)

Country Link
CN (1) CN107966465B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112710688B (zh) * 2019-10-24 2023-08-22 中国石油天然气股份有限公司 核磁共振纵向弛豫获取方法及系统
CN110780249B (zh) * 2019-11-21 2020-08-11 中国科学院武汉物理与数学研究所 使用绝热射频脉冲测量射频b1场分布的磁共振成像方法
CN112462311B (zh) * 2020-11-09 2021-06-15 无锡鸣石峻致医疗科技有限公司 一种横向磁化矢量衰减时间常数的校正测量方法、装置、计算机设备及非均匀场磁共振系统
CN112798635B (zh) * 2020-12-31 2024-02-20 东南大学 一种补偿射频磁场不均匀性的核磁共振信号脉冲方法
CN112903567B (zh) * 2021-04-15 2022-06-17 中国石油大学(北京) 孔径测量方法、装置、设备及存储介质

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1301810B1 (en) * 2000-07-21 2008-09-10 Services Petroliers Schlumberger Nuclear magnetic resonance methods for extracting information about a fluid in a rock
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
US6937014B2 (en) * 2003-03-24 2005-08-30 Chevron U.S.A. Inc. Method for obtaining multi-dimensional proton density distributions from a system of nuclear spins
CN104280703B (zh) * 2014-09-28 2017-01-18 中国石油大学(华东) 一种基于脉冲梯度磁场的d‑t2脉冲施加方法
CN105334239B (zh) * 2015-11-11 2018-02-02 中国石油大学(北京) 多维核磁共振流体组分含量测量方法及装置
CN106093100B (zh) * 2016-06-16 2017-09-01 中国石油大学(华东) 一种me‑cpmg序列的岩心核磁信号采集及反演方法
CN109270107B (zh) * 2017-03-27 2021-10-08 北京青檬艾柯科技有限公司 多维核磁共振测量方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Multidimensional NMR Inverse Laplace Spectroscopy in Petrophysics;Arns et al.;《PETROPHYSICS》;20071031;第48卷(第5期);380-392 *

Also Published As

Publication number Publication date
CN107966465A (zh) 2018-04-27

Similar Documents

Publication Publication Date Title
CN107966465B (zh) 一种基于三维脉冲序列的岩心核磁信号采集及反演方法
US10551520B2 (en) Surface relaxivity calculation using nuclear magnetic resonance (NMR) measurement, three dimensional (3D) rock model and NMR response simulation
Kenyon et al. A three-part study of NMR longitudinal relaxation properties of water-saturated sandstones
US6069477A (en) Method for improving the accuracy of NMR relaxation distribution analysis with two echo trains
CN104215652B (zh) 确定油气饱和度的方法和装置
Grunewald et al. Multiecho scheme advances surface NMR for aquifer characterization
US20130060474A1 (en) Estimation of petrophysical and fluid properties using integral transforms in nuclear magnetic resonance
CN106501144A (zh) 一种基于核磁共振双截止值的致密砂岩渗透率计算方法
RU2432587C2 (ru) Способ анализа данных с общими и различающимися свойствами
CN107748126B (zh) 一种获取岩石孔隙尺寸和孔隙表面弛豫率的核磁共振方法
NO336430B1 (no) Bestemmelse av T1 relaksasjonstid fra multippel ventetid NMR logger frembragt i den samme eller forskjellige passeringer
MX2008008999A (es) Pelicula de capas multiples resistente a la intemperie.
CN105334239A (zh) 多维核磁共振流体组分含量测量方法及装置
CN106526684B (zh) 一种核磁共振测井t2谱采集方法及装置
US20170123098A1 (en) Robust multi-dimensional inversion from wellbore nmr measurements
Luo et al. A machine learning framework for low-field NMR data processing
Guo et al. Numerical simulation of multi-dimensional NMR response in tight sandstone
CN107907911A (zh) 基于核磁共振的致密储层含油量测定方法
CN102913240B (zh) 一种储层流体识别方法
CN104330433A (zh) 一种获得目的储层横向弛豫时间分布的方法及装置
CN110454153A (zh) 一种核磁共振测井弛豫反演方法
CN106644879B (zh) 一种确定岩心不同孔隙组分渗透率贡献值的方法及装置
CN106093100B (zh) 一种me‑cpmg序列的岩心核磁信号采集及反演方法
NO20120654A1 (no) Fremgangsmåte og anordning for inkorporering av intern gradient og begrenset diffusjon i NMR-inversjon ved brønnlogging
CN104280703B (zh) 一种基于脉冲梯度磁场的d‑t2脉冲施加方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant