CN105116001A - 基于二维核磁共振确定油泥中油水含量的方法及装置 - Google Patents

Abstract

本发明实施例提供一种基于二维核磁共振确定油泥中油水含量的方法及装置。该方法包括：对油泥样品施加DEFIR脉冲序列获得T₁分布的第一幅值A_DEFIR(lgT₁)；对所述油泥样品施加FIR脉冲序列获得回波信号，对所述回波信号进行反演处理获得T₁分布的第二幅值A_FIR(lgT₁)；依据所述T₁分布的第一幅值A_DEFIR(lgT₁)和所述T₁分布的第二幅值A_FIR(lgT₁)获得T₁分布与T₁/T₂分布的函数关系；依据所述T₁分布与T₁/T₂分布的函数关系获得T₁-T₂分布，并依据所述T₁-T₂分布获得油泥中的油含量和水含量。本发明实施例通过对油泥样品先后施加DEFIR脉冲序列和FIR脉冲序列，获得油泥样品二维弛豫时间分布即T₁-T₂分布，依据T₁-T₂分布获得油泥中的油含量和水含量，相比于采用一维核磁共振方法确定油泥中油水含量的方法，提高了油泥中油水含量的测量精度。

Description

基于二维核磁共振确定油泥中油水含量的方法及装置

技术领域

本发明实施例涉及石油勘探领域，尤其涉及一种基于二维核磁共振确定油泥中油水含量的方法及装置。

背景技术

油泥是石油化学工业主要的污染物之一，它来自于包括原油生产、运输、储存和炼制在内的许多环节。在油泥处理工程中，确定油泥中油水含量至关重要。

现有技术中，确定油泥中油水含量的方法包括一维核磁共振方法，一维核磁共振方法具体采用T2谱确定油泥中油水含量，采用该方法的前提条件是油泥中的油峰和水峰必须在T2谱上区分开，否则，测量的油泥中油水含量与实际的油泥中油水含量偏差大，导致油泥中油水含量的测量精度低。

发明内容

本发明实施例提供一种基于二维核磁共振确定油泥中油水含量的方法及装置，以提高油泥中油水含量的测量精度。

本发明实施例的一个方面是提供一种基于二维核磁共振确定油泥中油水含量的方法，包括：

对油泥样品施加DEFIR脉冲序列获得T₁分布的第一幅值A_DEFIR(lgT₁)；

对所述油泥样品施加FIR脉冲序列获得回波信号，对所述回波信号进行反演处理获得T₁分布的第二幅值A_FIR(lgT₁)；

依据所述T₁分布的第一幅值A_DEFIR(lgT₁)和所述T₁分布的第二幅值A_FIR(lgT₁)获得T₁分布与T₁/T₂分布的函数关系；

依据所述T₁分布与T₁/T₂分布的函数关系获得T₁-T₂分布，并依据所述T₁-T₂分布获得油泥中的油含量和水含量。

本发明实施例的另一个方面是提供一种基于二维核磁共振确定油泥中油水含量的装置，包括：

探头，用于对油泥样品施加DEFIR脉冲序列；对所述油泥样品施加FIR脉冲序列；

电子线路，用于获得T₁分布的第一幅值A_DEFIR(lgT₁)；获得回波信号，对所述回波信号进行反演处理获得T₁分布的第二幅值A_FIR(lgT₁)；

上位机，用于依据所述T₁分布的第一幅值A_DEFIR(lgT₁)和所述T₁分布的第二幅值A_FIR(lgT₁)获得T₁分布与T₁/T₂分布的函数关系；依据所述T₁分布与T₁/T₂分布的函数关系获得T₁-T₂分布，并依据所述T₁-T₂分布获得油泥中的油含量和水含量。

本发明实施例提供的基于二维核磁共振确定油泥中油水含量的方法及装置，通过对油泥样品先后施加DEFIR脉冲序列和FIR脉冲序列，获得油泥样品二维弛豫时间分布即T₁-T₂分布，依据T₁-T₂分布获得油泥中的油含量和水含量，相比于采用一维核磁共振方法确定油泥中油水含量的方法，提高了油泥中油水含量的测量精度。

附图说明

图1为本发明实施例提供的基于二维核磁共振确定油泥中油水含量的方法流程图；

图2为本发明实施例提供的基于二维核磁共振确定油泥中油水含量的装置的结构图；

图3为本发明实施例提供的DEFIR脉冲序列的示意图；

图4为本发明实施例提供的油泥样品的磁化矢量的幅值的变化示意图；

图5为本发明实施例提供的回波信号幅值随响应时间的变化示意图；

图6为本发明实施例提供的T₁-T₂分布示意图；

图7为本发明实施例提供的油泥样品的含水率拟合示意图；

图8为本发明实施例提供的油泥样品的含油率拟合示意图。

具体实施方式

图1为本发明实施例提供的基于二维核磁共振确定油泥中油水含量的方法流程图；图2为本发明实施例提供的基于二维核磁共振确定油泥中油水含量的装置的结构图；图3为本发明实施例提供的DEFIR脉冲序列的示意图；图4为本发明实施例提供的油泥样品的磁化矢量的幅值的变化示意图。本发明实施例针对现有技术测量油泥中油水含量精度低的问题，提供了基于二维核磁共振确定油泥中油水含量的方法，该方法具体步骤如下：

步骤S101、对油泥样品施加DEFIR脉冲序列获得T₁分布的第一幅值A_DEFIR(lgT₁)；

执行本发明实施例方法的设备为基于二维核磁共振确定油泥中油水含量的装置，如图2所示，基于二维核磁共振确定油泥中油水含量的装置包括上位机10、探头12和电子线路13，油泥样品11放置在探头12上，上位机10控制电子线路13产生如图3所示的DEFIR脉冲序列，通过DEFIR脉冲序列对油泥样品11进行测量。

所述DEFIR脉冲序列包括DE脉冲序列、180度X轴正向脉冲和FIR脉冲序列；所述对油泥样品施加DEFIR脉冲序列，包括：对所述油泥样品依次施加x个所述DE脉冲序列、一个所述180度X轴正向脉冲和y个所述FIR脉冲序列，x≥2,y≥2。

如图3所示，DEFIR脉冲序列包括DE脉冲序列、180度X轴正向脉冲和FIR脉冲序列，其中，DE脉冲序列循环x次，FIR脉冲序列循环y次。

所述对所述油泥样品依次施加x个所述DE脉冲序列后，所述油泥样品的磁化矢量达到平衡值M_eq，其中，M₀表示所述油泥样品被完全极化时的磁化矢量，T₁表示T₁分布，T₂表示T₂分布，τ₁表示x个所述DE脉冲序列中最后一个所述DE脉冲序列的90度X轴负向脉冲与所述180度X轴正向脉冲之间的时间间隔，τ₂＝4τ_DE，τ_DE表示一个所述DE脉冲序列中相邻两个脉冲之间的时间间隔。

通过可知，油泥样品的磁化矢量对应的平衡值M_eq只与τ_DE、τ₁、T₁/T₂有关，与初始磁化矢量无关，因此，无论油泥样品在开始DE脉冲序列测量时初始磁化矢量是多少，在DE脉冲序列测量过程中或结束时，油泥样品的磁化矢量将趋于一个平衡值M_eq，且平衡值M_eq是固定值。如图4所示，对于不同的检测点其初始的磁化矢量的幅值不同，但随着施加DE脉冲序列时间的延长，油泥样品的磁化矢量的幅值均趋近于平衡值M_eq。

其中，x>>(τ₁/T₁+τ₂/T₂)^-1。

对所述油泥样品施加DEFIR脉冲序列获得所述油泥样品的磁化矢量M_DEFIR随极化时间TW_FIR的变化量M_DEFIR(TW_FIR)，依据M_DEFIR(TW_FIR)获得T₁分布的第一幅值A_DEFIR(lgT₁)。

步骤S102、对所述油泥样品施加FIR脉冲序列获得回波信号，对所述回波信号进行反演处理获得T₁分布的第二幅值A_FIR(lgT₁)；

所述对所述油泥样品施加FIR脉冲序列，包括：

对所述油泥样品依次施加y个所述FIR脉冲序列。

本发明实施例先执行步骤S101后执行S102，即先对油泥样品施加DEFIR脉冲序列，后对油泥样品施加FIR脉冲序列，此处，FIR脉冲序列重复y次。另外，对油泥样品先后施加的DEFIR脉冲序列与FIR脉冲序列之间可以有较长的时间间隔，也可以没有时间间隔。

步骤S103、依据所述T₁分布的第一幅值A_DEFIR(lgT₁)和所述T₁分布的第二幅值A_FIR(lgT₁)获得T₁分布与T₁/T₂分布的函数关系；

在本发明实施例中T₁分布与T₁/T₂分布的函数关系为 $<\frac{T_1}{T_2}{>}_T=\frac{{\mathrm{\&tau;}}_1}{{\mathrm{\&tau;}}_2}\&lsqb;\frac{2A_{FIR}\left({\mathrm{logT}}_1\right)-2A_{DEFIR}\left({\mathrm{logT}}_1\right)}{2A_{DEFIR}\left({\mathrm{logT}}_1\right)-A_{FIR}\left(1{\mathrm{ogT}}_1\right)}\&rsqb;.$

步骤S104、依据所述T₁分布与T₁/T₂分布的函数关系获得T₁-T₂分布，并依据所述T₁-T₂分布获得油泥中的油含量和水含量。

根据步骤S103获取到T₁分布与T₁/T₂分布的函数关系，另外结合T₁分布，采用现有技术中的任意一种可实现算法获得T₁-T₂分布，并依据所述T₁-T₂分布获得油泥中的油含量和水含量。

本发明实施例通过对油泥样品先后施加DEFIR脉冲序列和FIR脉冲序列，获得油泥样品二维弛豫时间分布即T₁-T₂分布，依据T₁-T₂分布获得油泥中的油含量和水含量，相比于采用一维核磁共振方法确定油泥中油水含量的方法，提高了油泥中油水含量的测量精度。

在上述实施例的基础上，所述对油泥样品施加DEFIR脉冲序列获得T₁分布的第一幅值A_DEFIR(lgT₁)，包括：

对所述油泥样品施加DEFIR脉冲序列获得所述油泥样品的磁化矢量M_DEFIR随极化时间TW_FIR的变化量M_DEFIR(TW_FIR)，M_DEFIR(TW_FIR)表示为公式(1)和公式(2)：

$M_{DEFIR}\left({\mathrm{TW}}_{FIR}\right)=\&Integral;d\lg T_1{A}_{DEFIR}\left(\lg T_1\right)\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0(1-2\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\left\}\right)---\left(1\right)$

$M_{DEFIR}\left({\mathrm{TW}}_{FIR}\right)=\&Integral;\&Integral;d\lg T_1{f}_1\left(\lg T_1\right)\frac{M_{eq}}{M_0}\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0---\left(2\right)$

依据所述公式(1)和所述公式(2)获得

其中，M_eq是对所述油泥样品依次施加x个所述DE脉冲序列后，所述油泥样品的磁化矢量达到的平衡值，M₀表示所述油泥样品被完全极化时的磁化矢量，TW_FIR表示所述DE脉冲序列中所述180度X轴正向脉冲与第一个所述FIR脉冲序列的90度X轴负向脉冲之间的时间间隔，f₁(lgT₁)表示一维T₁分布函数。

对所述油泥样品施加DEFIR脉冲序列获得所述油泥样品的磁化矢量M_DEFIR随极化时间TW_FIR的变化量为M_DEFIR(TW_FIR)， $M_{DEFIR}\left({\mathrm{TW}}_{FIR}\right)=(M_0+M_{y=0})\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0(1-2\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\left\}\right),$ 其中，M_y＝0＝-M_eq。

对所述油泥样品施加FIR脉冲序列获得所述油泥样品的磁化矢量M_FIR随极化时间TW_FIR的变化量为M_FIR(TW_FIR)， $M_{FIR}\left({\mathrm{TW}}_{FIR}\right)=(M_0+M_{y=0})\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0(1-2\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\left\}\right),$ 其中，M_y＝0＝-M₀。

将T₁分布的第一幅值A_DEFIR(lgT₁)引入公式 $M_{DEFIR}\left({\mathrm{TW}}_{FIR}\right)=(M_0+M_{y=0})\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0(1-2\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\left\}\right)$ 中得到公式(1)：

$M_{DEFIR}\left({\mathrm{TW}}_{FIR}\right)=\&Integral;d\lg T_1{A}_{DEFIR}\left(\lg T_1\right)\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0(1-2\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\left\}\right)---\left(1\right).$

将二维分布函数f₂(lgT₁,T₁/T₂)引入公式 $M_{DEFIR}\left({\mathrm{TW}}_{FIR}\right)=(M_0+M_{y=0})\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0(1-2\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\left\}\right)$ 中得到公式(4)：

$M_{DEFIR}\left({\mathrm{TW}}_{FIR}\right)=\&Integral;\&Integral;d\lg T_{1}d\frac{T_1}{T_2}{f}_2(\lg T,\frac{T_1}{T_2})(1+\frac{M_{eq}}{M_0})\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0---\left(4\right)$

在发明实施例中由于油泥样品的T1/T2分布很窄，公式(4)可近似为公式(2)；

$M_{DEFIR}\left({\mathrm{TW}}_{FIR}\right)=\&Integral;\&Integral;d\lg T_1{f}_1\left(\lg T_1\right)\frac{M_{eq}}{M_0}\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0---\left(2\right).$

依据所述公式(1)和所述公式(2)获得

A_FIR(lgT₁)＝2*f₁(lgT₁)，T₁分布与T₁/T₂分布的函数关系系为 $<\frac{T_1}{T_2}{>}_{T_1}=\frac{{\mathrm{\&tau;}}_1}{{\mathrm{\&tau;}}_2}\&lsqb;\frac{2A_{FIR}\left({\mathrm{logT}}_1\right)-2A_{DEFIR}\left({\mathrm{logT}}_1\right)}{2A_{DEFIR}\left({\mathrm{logT}}_1\right)-A_{FIR}\left(1{\mathrm{ogT}}_1\right)}\&rsqb;.$

根据和A_FIR(lgT₁)＝2*f₁(lgT₁)获得T₁分布与T₁/T₂分布的函数关系为 $<\frac{T_1}{T_2}{>}_{T_1}=\frac{{\mathrm{\&tau;}}_1}{{\mathrm{\&tau;}}_2}\&lsqb;\frac{2A_{FIR}\left({\mathrm{logT}}_1\right)-2A_{DEFIR}\left({\mathrm{logT}}_1\right)}{2A_{DEFIR}\left({\mathrm{logT}}_1\right)-A_{FIR}\left(1{\mathrm{ogT}}_1\right)}\&rsqb;.$

图5为本发明实施例提供的回波信号幅值随响应时间的变化示意图。如图5所示，水为长弛豫组分，油为短弛豫组分，长弛豫组分越多的油泥样品产生的回波信号的幅值随响应时间衰减速度慢，短弛豫组分越多的油泥样品产生的回波信号的幅值随响应时间衰减速度快。

图6为本发明实施例提供的T₁-T₂分布示意图。如图6所示，随着油泥中含水率的增加，二维T₁-T₂谱中水峰的信号越来越强，油峰的信号会越来越弱。

图7为本发明实施例提供的油泥样品的含水率拟合示意图。图8为本发明实施例提供的油泥样品的含油率拟合示意图。根据油泥中水峰和油峰在二维核磁共振T₁-T₂分布位置不同，可以由二维核磁共振T₁-T₂得到油泥中油和水的含量。采用最小二乘方法对不同含水率的油泥信号进行拟合，相关系数达到99％以上。由最小二乘得到的含水率和含油率，同油泥实际含水率和含油率具有一致性，结果如图7和图8所示。因此，采用最小二乘拟合得到的结果可以应用于未知油水含量的油泥样品的含水率和含油率的预测。

本发明实施例通过对油泥样品先后施加DEFIR脉冲序列和FIR脉冲序列，获得油泥样品二维弛豫时间分布即T₁-T₂分布，依据T₁-T₂分布获得油泥中的油含量和水含量，相比于采用一维核磁共振方法确定油泥中油水含量的方法，提高了油泥中油水含量的测量精度。

图2为本发明实施例提供的基于二维核磁共振确定油泥中油水含量的装置的结构图。本发明实施例提供的基于二维核磁共振确定油泥中油水含量的装置可以执行基于二维核磁共振确定油泥中油水含量的方法实施例提供的处理流程，如图2所示，基于二维核磁共振确定油泥中油水含量的装置包括探头12、电子线路13和上位机10，其中，探头12用于对油泥样品施加DEFIR脉冲序列；对所述油泥样品施加FIR脉冲序列；电子线路13用于获得T₁分布的第一幅值A_DEFIR(lgT₁)；获得回波信号，对所述回波信号进行反演处理获得T₁分布的第二幅值A_FIR(lgT₁)；上位机10用于依据所述T₁分布的第一幅值A_DEFIR(lgT₁)和所述T₁分布的第二幅值A_FIR(lgT₁)获得T₁分布与T₁/T₂分布的函数关系；依据所述T₁分布与T₁/T₂分布的函数关系获得T₁-T₂分布，并依据所述T₁-T₂分布获得油泥中的油含量和水含量。

本发明实施例通过对油泥样品先后施加DEFIR脉冲序列和FIR脉冲序列，获得油泥样品二维弛豫时间分布即T₁-T₂分布，依据T₁-T₂分布获得油泥中的油含量和水含量，相比于采用一维核磁共振方法确定油泥中油水含量的方法，提高了油泥中油水含量的测量精度。

在上述实施例的基础上，所述DEFIR脉冲序列包括DE脉冲序列、180度X轴正向脉冲和FIR脉冲序列；探头12具体用于对所述油泥样品依次施加x个所述DE脉冲序列、一个所述180度X轴正向脉冲和y个所述FIR脉冲序列，x≥2,y≥2；对所述油泥样品依次施加y个所述FIR脉冲序列。

电子线路13具体用于获得所述油泥样品的磁化矢量M_DEFIR随极化时间TW_FIR的变化量M_DEFIR(TW_FIR)，M_DEFIR(TW_FIR)表示为公式(1)和公式(2)：

$M_{DEFIR}\left({\mathrm{TW}}_{FIR}\right)=\&Integral;d\lg T_1{A}_{DEFIR}\left(\lg T_1\right)\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0(1-2\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\left\}\right)---\left(1\right)$

$M_{DEFIR}\left({\mathrm{TW}}_{FIR}\right)=\&Integral;\&Integral;d\lg T_1{f}_1\left(\lg T_1\right)\frac{M_{eq}}{M_0}\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0---\left(2\right)$

上位机10还用于依据所述公式(1)和所述公式(2)获得其中，M_eq是对所述油泥样品依次施加x个所述DE脉冲序列后，所述油泥样品的磁化矢量达到的平衡值，M₀表示所述油泥样品被完全极化时的磁化矢量，TW_FIR表示所述DE脉冲序列中所述180度X轴正向脉冲与第一个所述FIR脉冲序列的90度X轴负向脉冲之间的时间间隔，f₁(lgT₁)表示一维T₁分布函数。

其中，T₁表示T₁分布，T₂表示T₂分布，τ₁表示x个所述DE脉冲序列中最后一个所述DE脉冲序列的90度X轴负向脉冲与所述180度X轴正向脉冲之间的时间间隔，τ₂＝4τ_DE，τ_DE表示一个所述DE脉冲序列中相邻两个脉冲之间的时间间隔。

A_FIR(lgT₁)＝2*f₁(lgT₁)，所述T₁分布与T₁/T₂分布的函数关系为 $<\frac{T_1}{T_2}{>}_{T_1}=\frac{{\mathrm{\&tau;}}_1}{{\mathrm{\&tau;}}_2}\&lsqb;\frac{2A_{FIR}\left({\mathrm{logT}}_1\right)-2A_{DEFIR}\left({\mathrm{logT}}_1\right)}{2A_{DEFIR}\left({\mathrm{logT}}_1\right)-A_{FIR}\left(1{\mathrm{ogT}}_1\right)}\&rsqb;.$

本发明实施例通过对油泥样品先后施加DEFIR脉冲序列和FIR脉冲序列，获得油泥样品二维弛豫时间分布即T₁-T₂分布，依据T₁-T₂分布获得油泥中的油含量和水含量，相比于采用一维核磁共振方法确定油泥中油水含量的方法，提高了油泥中油水含量的测量精度。

综上所述，本发明实施例通过对油泥样品先后施加DEFIR脉冲序列和FIR脉冲序列，获得油泥样品二维弛豫时间分布即T₁-T₂分布，依据T₁-T₂分布获得油泥中的油含量和水含量，相比于采用一维核磁共振方法确定油泥中油水含量的方法，提高了油泥中油水含量的测量精度。

在本发明所提供的几个实施例中，应该理解到，所揭露的装置和方法，可以通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性，机械或其它的形式。

所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。

另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。

上述以软件功能单元的形式实现的集成的单元，可以存储在一个计算机可读取存储介质中。上述软件功能单元存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)或处理器(processor)执行本发明各个实施例所述方法的部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(Read-OnlyMemory，ROM)、随机存取存储器(RandomAccessMemory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。

本领域技术人员可以清楚地了解到，为描述的方便和简洁，仅以上述各功能模块的划分进行举例说明，实际应用中，可以根据需要而将上述功能分配由不同的功能模块完成，即将装置的内部结构划分成不同的功能模块，以完成以上描述的全部或者部分功能。上述描述的装置的具体工作过程，可以参考前述方法实施例中的对应过程，在此不再赘述。

最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (10)

1.一种基于二维核磁共振确定油泥中油水含量的方法，其特征在于，包括：

对油泥样品施加DEFIR脉冲序列获得T₁分布的第一幅值A_DEFIR(lgT₁)；

对所述油泥样品施加FIR脉冲序列获得回波信号，对所述回波信号进行反演处理获得T₁分布的第二幅值A_FIR(lgT₁)；

依据所述T₁分布的第一幅值A_DEFIR(lgT₁)和所述T₁分布的第二幅值A_FIR(lgT₁)获得T₁分布与T₁/T₂分布的函数关系；

依据所述T₁分布与T₁/T₂分布的函数关系获得T₁-T₂分布，并依据所述T₁-T₂分布获得油泥中的油含量和水含量。

2.根据权利要求1所述的方法，其特征在于，所述DEFIR脉冲序列包括DE脉冲序列、180度X轴正向脉冲和FIR脉冲序列；

所述对油泥样品施加DEFIR脉冲序列，包括：

对所述油泥样品依次施加x个所述DE脉冲序列、一个所述180度X轴正向脉冲和y个所述FIR脉冲序列，x≥2,y≥2；

所述对所述油泥样品施加FIR脉冲序列，包括：

对所述油泥样品依次施加y个所述FIR脉冲序列。

3.根据权利要求2所述的方法，其特征在于，所述对油泥样品施加DEFIR脉冲序列获得T₁分布的第一幅值A_DEFIR(lgT₁)，包括：

对所述油泥样品施加DEFIR脉冲序列获得所述油泥样品的磁化矢量M_DEFIR随极化时间TW_FIR的变化量M_DEFIR(TW_FIR)，M_DEFIR(TW_FIR)表示为公式(1)和公式(2)：

$M_{DEFIR}\left({\mathrm{TW}}_{FIR}\right)=\&Integral;{\mathrm{dlgT}}_1{A}_{DEFIR}\left({\mathrm{lgT}}_1\right)\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0(1-2\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\left\}\right)---\left(1\right)$

$M_{DEFIR}\left({\mathrm{TW}}_{FIR}\right)=\&Integral;\&Integral;{\mathrm{dlgT}}_1{f}_1\left({\mathrm{lgT}}_1\right)\frac{M_{eq}}{M_0}\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0---\left(2\right)$

依据所述公式(1)和所述公式(2)获得

其中，M_eq是对所述油泥样品依次施加x个所述DE脉冲序列后，所述油泥样品的磁化矢量达到的平衡值，M₀表示所述油泥样品被完全极化时的磁化矢量，TW_FIR表示所述DE脉冲序列中所述180度X轴正向脉冲与第一个所述FIR脉冲序列的90度X轴负向脉冲之间的时间间隔，f₁(lgT₁)表示一维T₁分布函数。

4.根据权利要求3所述的方法，其特征在于，其中，T₁表示T₁分布，T₂表示T₂分布，τ₁表示x个所述DE脉冲序列中最后一个所述DE脉冲序列的90度X轴负向脉冲与所述180度X轴正向脉冲之间的时间间隔，τ₂＝4τ_DE，τ_DE表示一个所述DE脉冲序列中相邻两个脉冲之间的时间间隔。

5.根据权利要求4所述的方法，其特征在于，A_FIR(lgT₁)＝2*f₁(lgT₁)，所述T₁分布与T₁/T₂分布的函数关系为 $<\frac{T_1}{T_2}{>}_{T_1}=\frac{{\mathrm{\&tau;}}_1}{{\mathrm{\&tau;}}_2}\&lsqb;\frac{2A_{FIR}\left({\mathrm{logT}}_1\right)-2A_{DEFIR}\left({\mathrm{logT}}_1\right)}{2A_{DEFIR}\left({\mathrm{logT}}_1\right)-A_{FIR}\left({\mathrm{logT}}_1\right)}\&rsqb;.$

6.一种基于二维核磁共振确定油泥中油水含量的装置，其特征在于，包括：

探头，用于对油泥样品施加DEFIR脉冲序列；对所述油泥样品施加FIR脉冲序列；

电子线路，用于获得T₁分布的第一幅值A_DEFIR(lgT₁)；获得回波信号，对所述回波信号进行反演处理获得T₁分布的第二幅值A_FIR(lgT₁)；

上位机，用于依据所述T₁分布的第一幅值A_DEFIR(lgT₁)和所述T₁分布的第二幅值A_FIR(lgT₁)获得T₁分布与T₁/T₂分布的函数关系；依据所述T₁分布与T₁/T₂分布的函数关系获得T₁-T₂分布，并依据所述T₁-T₂分布获得油泥中的油含量和水含量。

7.根据权利要求6所述的基于二维核磁共振确定油泥中油水含量的装置，其特征在于，所述DEFIR脉冲序列包括DE脉冲序列、180度X轴正向脉冲和FIR脉冲序列；

所述探头具体用于对所述油泥样品依次施加x个所述DE脉冲序列、一个所述180度X轴正向脉冲和y个所述FIR脉冲序列，x≥2,y≥2；对所述油泥样品依次施加y个所述FIR脉冲序列。

8.根据权利要求7所述的基于二维核磁共振确定油泥中油水含量的装置，其特征在于，所述电子线路具体用于获得所述油泥样品的磁化矢量M_DEFIR随极化时间TW_FIR的变化量M_DEFIR(TW_FIR)，M_DEFIR(TW_FIR)表示为公式(1)和公式(2)：

$M_{DEFIR}\left({\mathrm{TW}}_{FIR}\right)=\&Integral;{\mathrm{dlgT}}_1{A}_{DEFIR}\left({\mathrm{lgT}}_1\right)\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0(1-2\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\left\}\right)---\left(1\right)$

$M_{DEFIR}\left({\mathrm{TW}}_{FIR}\right)=\&Integral;\&Integral;{\mathrm{dlgT}}_1{f}_1\left({\mathrm{lgT}}_1\right)\frac{M_{eq}}{M_0}\exp \{-\frac{{\mathrm{TW}}_{FIR}}{T_1}\}+M_0---\left(2\right)$

所述上位机还用于依据所述公式(1)和所述公式(2)获得 $A_{DEFIR}\left({\mathrm{lgT}}_1\right)=f_1\left({\mathrm{lgT}}_1\right)(1+\frac{M_{eq}}{M_0});$

其中，M_eq是对所述油泥样品依次施加x个所述DE脉冲序列后，所述油泥样品的磁化矢量达到的平衡值，M₀表示所述油泥样品被完全极化时的磁化矢量，TW_FIR表示所述DE脉冲序列中所述180度X轴正向脉冲与第一个所述FIR脉冲序列的90度X轴负向脉冲之间的时间间隔，f₁(lgT₁)表示一维T₁分布函数。

9.根据权利要求8所述的基于二维核磁共振确定油泥中油水含量的装置，其中，T₁表示T₁分布，T₂表示T₂分布，τ₁表示x个所述DE脉冲序列中最后一个所述DE脉冲序列的90度X轴负向脉冲与所述180度X轴正向脉冲之间的时间间隔，τ₂＝4τ_DE，τ_DE表示一个所述DE脉冲序列中相邻两个脉冲之间的时间间隔。

10.根据权利要求9所述的基于二维核磁共振确定油泥中油水含量的装置，其特征在于，A_FIR(lgT₁)＝2*f₁(lgT₁)，所述T₁分布与T₁/T₂分布的函数关系为 $<\frac{T_1}{T_2}{>}_{T_1}=\frac{{\mathrm{\&tau;}}_1}{{\mathrm{\&tau;}}_2}\&lsqb;\frac{2A_{FIR}\left({\mathrm{logT}}_1\right)-2A_{DEFIR}\left({\mathrm{logT}}_1\right)}{2A_{DEFIR}\left({\mathrm{logT}}_1\right)-A_{FIR}\left({\mathrm{logT}}_1\right)}\&rsqb;.$