井下稠油分子链长的测量方法和装置的制造方法
【技术领域】
[0001] 本发明实施例设及稠油分析技术领域，尤其设及一种井下稠油分子链长的测量方 法和装置。
【背景技术】
[0002] 稠油是一种组成成分极其复杂的地下储层流体，具有粘度高，密度较大开采、运输 难度大等特点。稠油的开采价值已得到越来越多的关注。了解稠油内部原油组分是实际开 采，甚至分解、提炼和运移过程中提高重质油的开采利用效率重要的一步。
[0003] 目前常用原油组分及分子链长分析方法主要在实验室内进行，主要包括气相色谱 分析方法、化学位移方法及核磁共振方法等。然而，对于重质油来说，很难将重质组分由液 相或半固相转变为气相进行分析，运也导致气相色谱分析方法无法得到重质油全部的组分 信息。即使是对于化学位移方法来说，虽然可W得到样品精准的化学组分及内部分子结构， 但由于重质油组分过于复杂，其化学谱分布高度重叠，因而难W对其进行解释。对于核磁共 振方法来说，通过编辑扩散系数D、纵向弛豫时间Τι、横向弛豫时间T2,直接采用Ξ维反拉普 拉斯变换求取Ξ维联合概率分布函数，造成了数据处理速度慢。
[0004] 另一方面在实验室内进行原油组分及分子链长的分析方法，在原油从井下取样运 送至地面实验室的过程中，由于外界环境溫度、压力的变化，样品的物性参数也会发生变 化，个别重质成分可能析出，不利于稠油组分的准确测量。
【发明内容】
[000引本发明实施例提供一种井下稠油分子链长的测量方法和装置，实现了对稠油样品 中各组分分子链长的原位无损测量，无需对稠油组分进行气相分析，提高了对稠油样品中 各组成分子链长进行测量的准确度，并采用快速反演算法计算Ξ维联合概率分布函数，加 快了数据处理速度。
[0006] 本发明实施例提供一种井下稠油分子链长的测量方法，包括：
[0007] 向填充有稠油样品的核磁共振流体分析仪的探头内发射Ξ维核磁共振脉冲序列， 所述Ξ维核磁共振脉冲序列中包含Ξ个独立的分别用于编辑纵向弛豫时间Τι的信息、扩散 系数D的信息、横向弛豫时间Τ2的信息的窗口；
[0008] 采集回波数据；
[0009] 采用快速反演算法对所述回波数据进行反演，计算关于所述纵向弛豫时间Τι、所 述扩散系数D、所述横向弛豫时间Τ2的联合概率分布函数；
[0010] 根据预先搭建的所述纵向弛豫时间Τι与稠油样品中各组分分子链长关系的模型 或所述扩散系数D与稠油样品中各组分分子链长关系的模型或所述横向弛豫时间Τ2与稠油 样品中各组分分子链长关系的模型，求解所述稠油样品中各组分分子链长的概率分布函 数，W获得所述稠油样品中各组分的分子链长。
[0011] 本发明实施例提供一种井下稠油分子链长的测量装置，包括：
[0012] 发射模块，用于向填充有稠油样品的核磁共振流体分析仪的探头内发射Ξ维核磁 共振脉冲序列，所述Ξ维核磁共振脉冲序列中包含Ξ个独立的分别用于编辑纵向弛豫时间 Τι的信息、扩散系数D的信息、横向弛豫时间Τ2的信息的窗口；
[0013] 采集模块，用于采集回波数据；
[0014] 计算模块，用于采用快速反演算法对所述回波数据进行反演，计算关于所述纵向 弛豫时间Τι、所述扩散系数D、所述横向弛豫时间Τ2的联合概率分布函数；
[0015] 求解模块，用于根据预先搭建的所述纵向弛豫时间Τι与稠油样品中各组分分子链 长关系的模型或所述扩散系数D与稠油样品中各组分分子链长关系的模型或所述横向弛豫 时间Τ2与稠油样品中各组分分子链长关系的模型，求解所述稠油样品中各组分分子链长的 概率分布函数，W获得所述稠油样品中各组分的分子链长。
[0016] 本发明实施例提供一种井下稠油分子链长的测量方法和装置，通过向填充有稠油 样品的核磁共振流体分析仪的探头内发射Ξ维核磁共振脉冲序列，Ξ维核磁共振脉冲序 列中包含Ξ个独立的分别用于编辑纵向弛豫时间Τι的信息、扩散系数D的信息、横向弛豫时 间Τ2的信息的窗口；采集回波数据;采用快速反演算法对回波数据进行反演，计算关于纵向 弛豫时间Τι、扩散系数D、横向弛豫时间Τ2的联合概率分布函数;根据预先搭建的纵向弛豫时 间Τι与稠油样品中各组分分子链长关系的模型或扩散系数D与稠油样品中各组分分子链长 关系的模型或横向弛豫时间Τ2与稠油样品中各组分分子链长关系的模型，求解稠油样品中 各组分分子链长的概率分布函数，W获得稠油样品中各组分的分子链长。实现了对稠油样 品中各组分分子链长的原位无损测量，无需对稠油组分进行气相分析，提高了对稠油样品 中各组成分子链长进行测量的准确度，并采用快速反演算法计算Ξ维联合概率分布函数， 加快了数据处理速度。
【附图说明】
[0017] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现 有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发 明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可W 根据运些附图获得其他的附图。
[0018] 图1为本发明井下稠油分子链长的测量方法实施例一的流程图；
[0019] 图2为本发明实施例中；维核磁共振脉冲序列及回波数据的示意图；
[0020] 图3为本发明井下稠油分子链长的测量方法实施例二的流程图；
[0021] 图4为本发明中求解Ξ种不同粘度的油样的联合概率分布函数的结果示意图；
[0022] 图5为采用本发明提供的方法和气相色谱方法对1号油样各组分分子链长测量的 结果示意图；
[0023] 图6为采用本发明提供的方法和气相色谱方法对2号油样各组分分子链长测量的 结果示意图；
[0024] 图7为采用本发明提供的方法和气相色谱方法对3号油样各组分分子链长测量的 结果示意图；
[0025] 图8为本发明井下稠油分子链长的测量装置实施例一的结构示意图；
[0026] 图9为本发明井下稠油分子链长的测量装置实施例二的结构示意图。
【具体实施方式】
[0027] 为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例 中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是 本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员 在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0028] 在下述实施例中，公式中的连续的多个参数或变量表示该多个参数或变量相乘。 如两个变量为X，y，"巧"则表示X X y。
[0029] 图1为本发明井下稠油分子链长的测量方法实施例一的流程图，如图1所示，该方 法可W由井下稠油分子链长的测量装置来执行，该井下稠油分子链长的测量装置可W集成 在核磁共振流体分析仪中，如图1所示，则本实施例提供的井下稠油分子链长的测量方法包 括W下步骤：
[0030] 步骤101，向填充有稠油样品的核磁共振流体分析仪的探头内发射Ξ维核磁共振 脉冲序列，Ξ维核磁共振脉冲序列中包含Ξ个独立的分别用于编辑纵向弛豫时间Τι的信 息、扩散系数D的信息、横向弛豫时间Τ2的信息的窗口。
[0031] 本实施例中，将核磁共振流体分析仪直接下入储层，并从储层中抽吸稠油样品到 探头内进行探测，探测控件封闭，受外界环境影响小。该核磁共振流体分析仪的探头包括高 导磁外壳、磁体和天线。流体分析仪抽吸储层中的稠油样品进入探头后，对稠油样品进行Ξ 维核磁共振测量。采用"outside-in"的探测模式。该磁体能够产生均匀梯度的磁场，为进行 Ξ维核磁共振测量提供必要的静磁场环境。天线采用"单发单收"模式，长天线用于发射脉 冲序列，短天线用于接收回波数据，发射天线采用螺旋线圈结构，接收天线采用马鞍形结 构，W避免两天线之间的禪合作用。
[0032] 本实施例中，图2为本发明实施例中Ξ维核磁共振脉冲序列及回波数据的示意图， 如图2所示，Ξ维核磁共振脉冲序列是一种实现静止状态下测量Ξ维T1-D-T2联合概率分布 的脉冲序列。该Ξ维核磁共振脉冲序列中包括Ξ个独立的窗口。第一窗口 I通过改变极化时 间τι来编辑纵
向弛豫时间Τι的信息，第二窗口 II为扩散编辑脉冲序列，通过改变回波间隔T2 来编辑扩散系数D的信息，第Ξ窗口 ΙΠ 为常规的CPMG脉冲序列，通过短回波间隔Τε来编辑 横向弛豫时间Τ2的信息。
[0033] 本实施例中，Τ3 = ηΤΕ，其中，η可设置为任意自然数。
[0034] 本实施例中，在Ξ维核磁共振脉冲序列中改变第一窗口中的极化时间τι来编辑纵 向弛豫时间Τι的信息时，保持第二窗口的回波间隔Τ2和第Ξ窗口的短回波间隔Τε不变。在改 变第二窗口中的回波间隔Τ2来编辑扩散系数D的信息时，保持第一窗口的极化时间τι和第Ξ 窗口中的短回波间隔Τε不变。在改变第Ξ窗口中的短回波间隔Τε来编辑横向弛豫时间Τ2的 信息时，保持第一窗口中的极化时间τι和第二窗口中的回波间隔Τ2不变。
[00巧]步骤102,采集回波数据。
[0036]本实施实例中，采集的回波数据包括各回波的幅值数据。如图2所示，在发射第一 窗口的脉冲序列后，采集到了第一回波数据201，在发射Ξ维核磁共振脉冲序列的同时，依 次采集回波数据，在图2中采集到的回波数据分别为201，202,203,204和205。运五个回波数 据的最大幅值依次减小。
[0037] 步骤103,采用快速反演算法对回波数据进行反演，计算关于纵向弛豫时间Τι、扩 散系数D、横向弛豫时间Τ2的联合概率分布函数。
[0038] 本实施例中，采集到的回波数据为关于纵向弛豫时间Τι的信息、扩散系数D的信 息、横向弛豫时间Τ2的信息的Ξ维回波数据，大量的Ξ维回波数据构成了回波数据体，对运 些Ξ维回波数据进行反演，在反演时采用快速反演算法，首先对Ξ维回波数据降维压缩，然 后计算关于纵向弛豫时间Τι、扩散系数D、横向弛豫时间Τ2的联合概率分布函数。
[0039] 步骤104,根据预先搭建的纵向弛豫时间Τι与稠油样品中各组分分子链长关系的 模型或扩散系数D与稠油样品中各组分分子链长关系的模型或横向弛豫时间Τ2与稠油样品 中各组分分子链长关系的模型，求解稠油样品中各组分分子链长的概率分布函数，W获得 稠油样品中各组分的分子链长。
[0040] 本实施例中，流体的核磁共振自旋弛豫响应由整个流体中分子的纵向弛豫时间 Τι、扩散系数D的信息、横向弛豫时间Τ2的信息表征。特定组分的分子的扩散系数与分子大小 或分子链长相关。同时整体的流体环境也影响分子的扩散系数。因此根据扩散系数与分子 链长、整体的流体环境的关系，搭建扩散系数D与稠油样品中各组分分子链长关系的模型， 该模型定义为第一模型。
[0041] 本实施例中，分子的弛豫通常由分子间核子的禪合作用决定。纵向弛豫时间Τι的 信息表征自旋系统与核子的环境的能量交换。横向弛豫时间Τ2描述了自旋系统内部的偶 极-偶极相关作用。纵向弛豫时间Τι、横向弛豫时间Τ2也分别与分子链长及整体的流体环境 有关，所W根据纵向弛豫时间Τι与分子链长及整体的流体环境的关系，搭建纵向弛豫时间Τι 与稠油样品中各组分分子链长关系的模型，该模型定义为第二模型。根据横向弛豫时间Τ2 与分子链长及整体的流体环境的关系，搭建横向弛豫时间Τ2与稠油样品中各组分分子链长 关系的模型，该模型定义为第Ξ模型。
[0042] 本实施例中，在获取纵向弛豫时间Τι、扩散系数D、横向弛豫时间Τ2的联合概率分布 函数后，可获得稠油中各组分的纵向弛豫时间Τι、扩散系数D、横向弛豫时间Τ2。根据上述中 Ξ个模型中的任意模型，求解稠油样品中各组分分子链长的概率分布函数，在获得稠油样 品中各组分分子链长的概率分布函数后，可获得稠油样品中各组分的分子链长。
[0043] 本实施例中，Ξ个模型都是根据大量实验搭建出来的，Ξ个模型中的各参数的取 值也是根据大量实验确定的。
[0044] 本实施例提供的井下稠油分子链长的测量方法，通过向填充有稠油样品的核磁共 振流体分析仪的探头内发射Ξ维核磁共振脉冲序列，Ξ维核磁共振脉冲序列中包含Ξ个独 立的分别用于编辑纵向弛豫时间Τι的信息、扩散系数D的信息、横向弛豫时间Τ2的信息的窗 口；采集回波数据;采用快速反演算法对回波数据进行反演，计算关于纵向弛豫时间Τι、扩 散系数D、横向弛豫时间Τ2的联合概率分布函数;根据预先搭建的纵向弛豫时间Τι与稠油样 品中各组分分子链长关系的模型或扩散系数D与稠油样品中各组分分子链长关系的模型或 横向弛豫时间Τ2与稠油样品中各组分分子链长关系的模型，求解稠油样品中各组分分子链 长的概率分布函数，W获得稠油样品中各组分的分子链长。由于在井下封闭环境下进行测 量，所W实现了对稠油样品中各组分分子链长的原位无损测量，无需对稠油组分进行气相 分析，提高了对稠油样品中各组成分子链长进行测量的准确度，并采用快速反演算法计算 Ξ维联合概率分布函数，加快了数据处理速度。
[0045] 图3为本发明井下稠油分子链长的测量方法实施例二的流程图，如图3所示，该方 法可W由井下稠油分子链长的测量装置来执行，该井下稠油分子链长的测量装置可W集成 在核磁共振流体分析仪中，如图3所示，则本实施例提供的显示信息的防偷看方法包括W下 步骤：
[0046] 步骤301，向填充有稠油样品的核磁共振流体分析仪的探头内发射Ξ维核磁共振 脉冲序列，Ξ维核磁共振脉冲序列中包含Ξ个独立的分别用于编辑纵向弛豫时间Τι的信 息、扩散系数D的信息、横向弛豫时间Τ2的信息的窗口。
[0047] 步骤302,采集回波数据。
[0048] 本实施例中，步骤301-步骤302和本发明井下稠油分子链长的测量方法实施例一 中的步骤101-步骤102相同，在此不再一一寶述。
[0049] 步骤303,针对采集的回波数据的第一类化e化olm方程，采用奇异值分解和张量积 的算法，对回波数据进行压缩降维处理，W获取关于回波数据的降维方程。
[0050] 本实施例中，关于回波数据的第一类Fre化olm方程可表示为式(1)所示：
[0051] Μ(τι,τ2,τ3)=Π ]·?α(χ,τι)1?2(γ，τ2)1?(ζ,τ3)?·(χ,γ,ζ)(1χ(1γ(1ζ+Ε(τι,τ2,τ3) (1)
[0052] 其中，Μ(τι，τ2，τ3)表示回波数据，回波数据为Ξ维数据，大量的回波数据构成回波 数据体，τ康示第一窗口中的极化时间，τ康示第二窗口中的回波间隔，Τ3 = ηΤΕ，ΤΕ表示第立 窗口中的短回波间隔，η为预设的自然数，X表示为纵向弛豫时间Ti，y表示扩散系数D，z表示 横向弛豫时间T2，kl(Ti，X)表示Ξ维核磁共振脉冲序列的第一反演核函数，k2(T2，y)表示Ξ 维核磁共振脉冲序列的第二反演核函数，k3(Z，T3)表示Ξ维核磁脉冲序列的第Ξ反演核函 数，f(x，y，z)表示Ξ维联合概率密度函数，ε(τ^2，Τ3)为噪声项。
[0053] 本实施例中，ε(τιτ2，Τ3)为可W为高斯白噪声。
[0054] 具体地，ki(Ti，x)为关于纵向弛豫时间Τι的反演核函数，k2(T2，y)为关于扩散系数D 的反演核函数，k3(Z，T3)为关于横向弛豫时间T2的反演核函数，kl(Ti,X)、k2(T2，y)、k3(Z，T3) 可表示为式(2)-式(4)所示。
[005引其中，G表示均匀梯度磁场的梯度，丫表示质子旋磁比。由于在式（3)中e邱(-2T2/ Z)通常情况下近似为1，所W式（3)可W写为A::0-:，J)=紛的-.1''护' r '6)，即k2(T2，y) 为关于扩散系数D的反演核函数。
[0059] 本实施例中，由于采集的大量的回波数据为Ξ维的回波数据体，所W对回波数据 进行压缩降维处理，加快数据处理速度。在进行压缩降维处理时，采用奇异值分解和张量积 的算法。把纵向弛豫时间Τι、扩散系数D两部分禪合成一个新函数，获取关于回波数据的降 维方程表示为式(5)所示：
[0060] M(TiT2,T3)=JJkl2Uy，T^2)k3(Z,T3)f2DUy，Z)ckydz + e(T：^2，T3) 巧）
[0061] 其中，Μ(τιτ2，τ3)表示压缩降维后的回波数据，ε(τιτ2，Τ3)表示降维压缩后的噪声 项，kl2(xy，τlτ2)=kl(τl,χ)k2(τ2，y)，f2D(xy
，z)表示二维联合概率密度函数。
[0062] 本实施例中，通过压缩降维处理后，将求解Ξ维联合概率分布函数的问题转变为 求解二维联合概率分布函数的问题。
[0063] 步骤304,利用零阶Ti化onov正则化方法，求解关于纵向弛豫时间Τι和扩散系数D 的二维联合概率分布函数。
[0064] 本实施例中，流体的核磁共振弛豫或扩散响应满足多指数衰减规律，由回波数据 反演得到弛豫或扩散分布是一个病态问题，使反演过程的解不唯一，并且回波数据中小的 扰动就会对求解结果造成较大影响，所W在求解过程中引入正则化处理，W对反演过程进 行稳定处理。
[006引进一步地，步骤304具体为:利用零阶Ti化onov正则化方法，求解使惩罚函数最小 的关于纵向弛豫时间Τι和扩散系数D的二维联合概率分布函数。
[0066] 其中，惩罚函数表示为式(6)所示：
[0067]
(0)
[006引本实施例中，X表示为纵向弛豫时间Τι，y表示扩散系数D，Ζ表示横向弛豫时间Τ2，τι 表示第一窗口中的极化时间，Τ2表示第二窗口中的回波间隔，Τ3 = ηΤΕ，Τε表示第Ξ窗口中的 短回波间隔，11为预设的自然数山2(巧，1仍）=41(11，义化2(12,7)山（11，义)表示；维核磁共 振脉冲序列的第一反演核函数，k2(T2，y)表示Ξ维核磁共振脉冲序列的第二反演核函数，k3 (z，T3)表示Ξ维核磁共振脉冲序列的第Ξ反演核函数，0<a<l，F2D(xy，z)表示关于纵向弛豫 时间Τι和扩散系数D的二维联合概率分布函数，务表示惩罚函数，Μ(τιτ2，τ3)表示降维压缩 后的回波数据。
[0069] 步骤305,根据关于横向弛豫时间Τ2的概率分布函数和二维联合概率分布函数，采 用Ξ维反拉普拉斯变换的方法，求解关于纵向弛豫时间Τι、扩散系数D、横向弛豫时间Τ2的联 合概率分布函数。
[0070] 本实施例中，首先根据第Ξ反演核函数k3(z，T3)求解关于横向弛豫时间Τ2的概率 分布函数，然后根据关于横向弛豫时间T2的概率分布函数和二维联合概率分布函数，采用 Ξ维反拉普拉斯变换的方法，求解关于纵向弛豫时间Τι、扩散系数D、横向弛豫时间T2的联合 概率分布函数。
[0071] 步骤306,根据预先搭建的纵向弛豫时间Τι与稠油样品中各组分分子链长关系的 模型或扩散系数D与稠油样品中各组分分子链长关系的模型或横向弛豫时间Τ2与稠油样品 中各组分分子链长关系的模型，求解稠油样品中各组分分子链长的概率分布函数，W获得 稠油样品中各组分的分子链长。
[0072] 进一步地，预先搭建的纵向弛豫时间Τι与稠油样品中各组分分子链长关系的模型 表示为式(7)所示：
[007引
巧
[0074] 预先搭建的扩散系数D与稠油样品中各组分分子链长关系的模型表示为式(8)所 示：
[0075]
(8)
[0076] 预先搭建的横向弛豫时间Τ2与稠油样品中各组分分子链长关系的模型表示为式 (9)所示：
[0077]
(9)
[007引其中，Di、Tii、T2i分别为第i种组分的扩散系数，纵向弛豫时间，横向弛豫时间，劳 为平均分子链长，Ni为第i种组分的分子链长，A为第一经验相关系数，B为第二经验相关系 数，β表示第一整体流体环境的分子尺寸信息，V表示第i种组分分子在所有组分中的第一权 重，丫 1表示第二整体流体环境的分子尺寸信息，ki表示第i种组分分子在所有组分中的第二 权重，丫 2表示第Ξ整体流体环境的分子尺寸信息，k2表示第i种组分分子在所有组分中的第 Ξ权重。
[0079] 本实施例中，Α、Β、β、ν、丫 i、ki、丫 2、k2的取值都为进行大量实验后的经验值，可预 先进行设定。
[0080] 本实施例中，可根据Ξ个模型中的任一个求解稠油样品中各组分分子链长的概 率分布函数。
[0081] 本实施例提供的井下稠油分子链长的测量方法，在采用快速反演算法对回波数据 进行反演，计算关于纵向弛豫时间Τι、扩散系数D、横向弛豫时间T2的联合概率分布函数时， 首先针对采集的回波数据的第一类Fre化olm方程，采用奇异值分解和张量积的算法，对回 波数据进行压缩降维处理，W获取关于回波数据的降维方程;然后利用零阶Ti化onov正则 化方法，求解关于纵向弛豫时间Τι和扩散系数D的二维联合概率分布函数;最后根据关于横 向弛豫时间Τ2的概率分布函数和二维联合概率分布函数，采用Ξ维反拉普拉斯变换的方 法，求解关于纵向弛豫时间Τι、扩散系数D、横向弛豫时间Τ2的联合概率分布函数，通过压缩 降维处理后，将求解Ξ维联合概率分布函数的问题转变为求解二维联合概率分布函数的问 题，加快了数据处理速度，并且在求解过程中引入正则化处理，使反演过程更稳定，由于可 选择Ξ个模型中的任意模型求解稠油样品中各组分分子链长的概率分布函数，增加了测量 稠油样品中各组分分子链长的途径。
[0082] 为了能说明本实施例提供的井下稠油分子链长的测量方法的有益效果，可通过W 下实验进行验证。
[0083] 1实验内容:采用本发明提供的井下稠油分子链长的测量方法及气相色谱方法对 Ξ种不同粘度的油样进行各组分分子链长的测量。
[0084] 2实验结果
[0085] 图4为本发明中求解Ξ种不同粘度的油样的联合概率分布函数的结果示意图。图4 中的第一行为1号油样的联合概率分布示意图，第1行的第1列为1号油样的Ξ维联合概率分 布示意图，第1行的第2-4列分别为1号油样的二维联合概率分布示意图，第2行为2号油样的 联合概率分布示意图，第2行的第1列为2号油样的Ξ维联合概率分布示意图，第2行的第2-4 列分别为第2号油样的二维联合概率分布示意图，第3行为3号油样的联合概率分布示意图， 第3行的第1列为3号油样的Ξ维联合概率分布示意图，第3行的第2-4列分别为第3号油样的 二维联合概率分布示意图。
[0086] 图5为采用本发明提供的方法和气相色谱方法对1号油样各组分分子链长测量的 结果示意图。图6为采用本发明提供的方法和气相色谱方法对2号油样各组分分子链长测量 的结果示意图。图7为采用本发明提供的方法和气相色谱方法对3号油样各组分分子链长 测量的结果示意图。
[0087] 其中，利用气相色谱方法对上述Ξ个油样各组分分子链长进行测量时采用 Agilent 7890A气相色谱仪。利用纯度为99.999%的氮气作为气相载体。系统初始溫度为40 °C，持续10分钟后W4°C/min的速度增加至70°C，然后W8°C/min的速度再持续增加至300 °C。受到气相转换的限制，碳个数大于36的分子无法转化为气相，故而无法在气相色谱分析 上得到体现。因此对于重质油样品来说，该方法只能作为一种参考测量手段。有前面的分析 可知，分子弛豫及扩散信息与其碳链长度成指数关系。因此将上述的Ξ个油样的气相色谱 分析结果转换至对数坐标系中，便于和Ξ维核磁共振结果相关联。
[0088] 如图4所示，随着油样粘度的增加，油样的Ξ维联合概率分布越来越复杂。对于1号 油样，Ξ维联合概率分布中存在两个信号。运两个信号对应的扩散系数在l(Tiim 2/s左右，具 有不同的Τι和T2值，T2值相差一个数量级而Τι值在同一数量级。对于2号油样，出现了 Ξ个信 号，其中，两个信号对应扩散系数大约在ICTiVVs范围内，另一个信号对应的扩散系数大约 为2Xl(TiV/s。粘度更大的3号油样在；维联合概率分布空间中出现了四个信号。运四个 信号可W被分为两组，分别对应的扩散系数为l(riV/s和1(TV/s^^个油样的信号中， 主要的信号对应扩散系数为l(Tiim 2/s，并且纵向弛豫时间Τι均大于横向弛豫时间T2,暗示着 Ξ个油样中均存在长链的大分子。一些比较特殊的分子团，例如渐青质，可能是造成重质油 成分复杂性的原因。而对应于扩散系数为区间的信号，其弛豫时间几乎 相等，运表明运些信号由轻质成分贡献而来的。尽管3号原油的粘度最大，但是3号原油在运 个扩散系数范围内的信号最强。表明该样品中的轻质成分最多。通过分馈的结果，可W发 现，3号重质油样的饱和烧控质量比重为50%，高于1号油样的40.8%和
2号油样的43.6%。
[0089] 如图5,图6和图7所示，在图5,图6和图7中的核磁共振测量结果即为采用本实施例 提供的测量方法的结果。碳链长度即为各组分的分子链长，由于本实施例测量方法在无需 对稠油组分进行气相分析，能够得到油样中的全部的组分信息，并且对油样中各组分分子 链长的测量在井下封闭环境，为原位无损测量，所W本发明提供井下稠油分子链长的测量 方法测量的分子链长准确度更高。
[0090] 图8为本发明井下稠油分子链长的测量装置实施例一的结构示意图，如图8所示， 本实施例提供的井下稠油分子链长的测量装置包括:发射模块801，采集模块802,计算模块 803和求解模块804。
[0091] 其中，发射模块801，用于向填充有稠油样品的核磁共振流体分析仪的探头内发射 Ξ维核磁共振脉冲序列，Ξ维核磁共振脉冲序列中包含Ξ个独立的分别用于编辑纵向弛豫 时间Τι的信息、扩散系数D的信息、横向弛豫时间T2的信息的窗口。采集模块802,用于采集回 波数据。计算模块803,用于采用快速反演算法对回波数据进行反演，计算关于纵向弛豫时 间Τι、扩散系数D、横向弛豫时间Τ2的联合概率分布函数。求解模块804,用于根据预先搭建的 纵向弛豫时间Τι与稠油样品中各组分分子链长关系的模型或扩散系数D与稠油样品中各组 分分子链长关系的模型或横向弛豫时间Τ2与稠油样品中各组分分子链长关系的模型，求解 稠油样品中各组分分子链长的概率分布函数，W获得稠油样品中各组分的分子链长。
[0092] 具体地，本实施例中的发射模块801可W为发射天线，该发射天线采用螺旋线圈结 构采集模块802可W为接收天线，接收天线采用马鞍形结构，W避免两天线之间的禪合作 用。
[0093] 本实施例中，该井下稠油分子链长的测量装置可集成在核磁共振流体分析仪中， 该核磁共振流体分析仪的探头包括高导磁外壳、磁体和天线。流体分析仪抽吸储层中的稠 油样品进入探头后，对稠油样品进行Ξ维核磁共振测量。采用"outside-in"的探测模式。该 磁体能够产生均匀梯度的磁场。为进行Ξ维核磁共振测量提供必要的静磁场环境。
[0094] 本实施例提供的井下稠油分子链长的测量装置可W执行图1所示方法实施例的技 术方案，其实现原理和技术效果类似，此处不再寶述。
[0095] 图9为本发明井下稠油分子链长的测量装置实施例二的结构示意图；如图9所示， 本实施例提供的井下稠油分子链长的测量装置包括:发射模块901，采集模块902,计算模块 903和求解模块904。
[0096] 发射模块901，用于向填充有稠油样品的核磁共振流体分析仪的探头内发射Ξ维 核磁共振脉冲序列，Ξ维核磁共振脉冲序列中包含Ξ个独立的分别用于编辑纵向弛豫时间 Τι的信息、扩散系数D的信息、横向弛豫时间T2的信息的窗口。采集模块902,用于采集回波数 据。计算模块903,用于采用快速反演算法对回波数据进行反演，计算关于纵向弛豫时间 Τι、扩散系数D、横向弛豫时间Τ2的联合概率分布函数。求解模块904,用于根据预先搭建的纵 向弛豫时间Τι与稠油样品中各组分分子链长关系的模型或扩散系数D与稠油样品中各组分 分子链长关系的模型或横向弛豫时间Τ2与稠油样品中各组分分子链长关系的模型，求解稠 油样品中各组分分子链长的概率分布函数，W获得稠油样品中各组分的分子链长。
[0097] 进一步地，计算模块903包括压缩降维单元903a和求解单元903b。
[0098] 其中，压缩降维单元903曰，用于针对采集的回波数据的第一类Fre化olm方程，采用 奇异值分解和张量积的算法，对回波数据进行压缩降维处理，W获取关于回波数据的降维 方程。求解单元903b，用于利用零阶Ti化onov正则化方法，求解关于纵向弛豫时间Τι和扩散 系数D的二维联合概率分布函数。求解单元903b，还用于根据关于横向弛豫时间Τ2的概率分 布函数和二维联合概率分布函数，采用Ξ维反拉普拉斯变换的方法，求解关于纵向弛豫时 间Τι、扩散系数D、横向弛豫时间T2的联合概率分布函数。
[0099] 进一步地，求解单元903b，具体用于：利用零阶Ti化onov正则化方法，求解使惩罚 函数最小的关于纵向弛豫时间Τι和扩散系数D的二维联合概率分布函数。
[0100] 其中，惩罚函数表示为：
X表示为纵向弛豫时间Ti，y表示扩散系数D，z表示横向弛豫时间Τ2，τι表示第一窗口中的极 化时间，τ康示第二窗口中的回波间隔，Τ3 = ηΤΕ，Τε表示第；窗口中的短回波间隔，η为预设 的自然数山2(巧，1：112)=41(1：1,义化2(12,7)山(1：1,义)表示立维核磁共振脉冲序列的第一反 演核函数，k2(T2，y)表示Ξ维核磁共振脉冲序列的第二反演核函数，k3(Z，T3)表示Ξ维核磁 共振脉冲序列的第Ξ反演核函数，0<a<l，F2D(xy，z)表示关于纵向弛豫时间Τι和扩散系数D 的二维联合概率分布函数，#表示惩罚函数，Μ(τιτ2，τ3)表示降维压缩后的回波数据。
[0101] 进一步地，第一类 Fre 化olm 方程表示为:M(Ti，T2，T3) = nAi(x，Ti)k2(y，T2)k3(z, T3)f(x，y，z)dxdydz+e(Ti,T2,T3);
[010引关于回波数据的降维方程表示为:M(T1T2，T3)=J化12(Xy，TiT2)k3(Z，T3)f2D(Xy，Z) (kydz+e(T^2，T3);
[0103] 其中，M( τι，T2，T3)表示回波数据，f (X，y，z)表示Ξ维联合概率密度函数，ε (τι，T2， T3)表示噪声项，f2D(xy，z)表示二维联合概率密度函数，ε(τιτ2，Τ3)表示压缩降维后的噪声 项。
[0104] 进一步地，纵向弛豫时间Τι与稠油样品中各组分分子链长关系的模型表示为： ?[，二 Β.矿7、·，
[0105] 扩散系数D与稠油样品中各组分分子链长关系的模型表示为：巧方-;
[0106] 横向弛豫时间T2与稠油样品中各组分分子链长关系的模型表示为： 二公 '矿&-典
[0107] 其中，Di、Tii、T2i分别为第i种组分的扩散系数，纵向弛豫时间，横向弛豫时间，茨 为平均分子链长，Ni为第i种组分的分子链长，A为第一经验相关系数，B为第二经验相关系 数，β表示第一整体流体环境的分子尺寸信息，V表示第i种组分分子在所有组分中的第一权 重，丫 1表示第二整体流体环境的分子尺寸信息，ki表示第i种组分分子在所有组分中的第二 权重，丫 2表示第Ξ整体流体环境的分子尺寸信息，k2表示第i种组分分子在所有组分中的第 二权重。
[0108] 本实施例的装置可W执行图3所示方法实施例的技术方案，其实现原理和技术效 果类似，此处不再寶述。
[0109] 最后应说明的是：W上各实施例仅用W说明本发明的技术方案，而非对其限制;尽 管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解:其依 然可W对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进 行等同替换；而运些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术 方案的范围。
【主权项】
1. 一种井下稠油分子链长的测量方法，其特征在于，包括： 向填充有稠油样品的核磁共振流体分析仪的探头内发射三维核磁共振脉冲序列，所述 三维核磁共振脉冲序列中包含三个独立的分别用于编辑纵向弛豫时间T1的信息、扩散系数 D的信息、横向弛豫时间!^的信息的窗口； 采集回波数据； 采用快速反演算法对所述回波数据进行反演，计算关于所述纵向弛豫时间!^、所述扩散 系数D、所述横向弛豫时间!^的联合概率分布函数； 根据预先搭建的所述纵向弛豫时间T1与稠油样品中各组分分子链长关系的模型或所述 扩散系数D与稠油样品中各组分分子链长关系的模型或所述横向弛豫时间T2与稠油样品中 各组分分子链长关系的模型，求解所述稠油样品中各组分分子链长的概率分布函数，以获 得所述稠油样品中各组分的分子链长。2. 根据权利要求1所述的方法，其特征在于，所述采用快速反演算法对所述回波数据进 行反演，计算关于所述纵向弛豫时间T1、所述扩散系数D、所述横向弛豫时
间T 2的联合概率分 布函数具体包括： 针对采集的所述回波数据的第一类Fredholm方程，采用奇异值分解和张量积的算法， 对所述回波数据进行压缩降维处理，以获取关于回波数据的降维方程； 利用零阶Tikhonov正则化方法，求解关于所述纵向弛豫时间TdP所述扩散系数D的二维 联合概率分布函数； 根据关于横向弛豫时间T2的概率分布函数和所述二维联合概率分布函数，采用三维反 拉普拉斯变换的方法，求解关于所述纵向弛豫时间!^、所述扩散系数D、所述横向弛豫时间T2 的联合概率分布函数。3. 根据权利要求2所述的方法，其特征在于，所述利用零阶Tikhonov正则化方法，求解 关于所述纵向弛豫时间TdP所述扩散系数D的二维联合概率分布函数具体包括： 利用零阶Tikhonov正则化方法，求解使惩罚函数最小的关于所述纵向弛豫时间T1和所 述扩散系数D的二维联合概率分布函数； 其中，所述惩罚函数表示为：X表示为纵向弛豫时间Thy表示扩散系数D，z表示横向弛豫时间Τ2，τι表示第一窗口中的极 化时间，12表示第二窗口中的回波间隔，τ3 = ηΤΕ，ΤΕ表示第三窗口中的短回波间隔，η为预设 的自然数，ki2(xy，T 1T2) = ki (T1，X)k2 (Τ2，y)，ki (T1，X)表示所述三维核磁共振脉冲序列的第 一反演核函数，k2(i2，y)表示所述三维核磁共振脉冲序列的第二反演核函数，k 3(z，T3)表示 所述三维核磁共振脉冲序列的第三反演核函数，〇〈〇〈1^邪(灯， 2)表示关于所述纵向弛豫时 间TjP所述扩散系数D的二维联合概率分布函数，#表示惩罚函数，M(T1Tm)表示降维压 缩后的回波数据。4. 根据权利要求3所述的方法，其特征在于，所述第一类Fredholm方程表示为:M( T1，T2， T3) = iiiki(x,Ti)k2(y,T2)k3(z,T3)f(x,y,z)dxdydz+e(Ti,T2,T3)； 所述关于回波数据的降维方程表示为:Μ(τιτ2，τ3) = Π?^(Χ7，τ?τ2)1?(Ζ，τ3)?· 2?(Χ7，Ζ) dxydz+e(TiT2,T3)； 其中，Μ(τι，τ2，τ3)表示回波数据，f(x，y，z)表示三维联合概率密度函数，以^仏"表 示噪声项，f2D(xy，z)表示二维联合概率密度函数，ε(τ?τ2，τ3)表示压缩降维后的噪声项。5. 根据权利要求1-4任一项所述的方法，其特征在于，所述纵向弛豫时间T1与稠油样品 中各组分分子链长关系的模型表示为：； 所述扩散系数D与稠油样品中各组分分子链长关系的模型表示为：所述横向弛豫时间Τ2与稠油样品中各组分分子链长关系的模型表示为：其中，D1J11J21分别为第i种组分的扩散系数，纵向弛豫时间，横向弛豫时间，歹为平均 分子链长，N1为第i种组分的分子链长，A为第一经验相关系数，B为第二经验相关系数，β表 示第一整体流体环境的分子尺寸信息，V表示第i种组分分子在所有组分中的第一权重， γι 表示第二整体流体环境的分子尺寸信息，kl表示第i种组分分子在所有组分中的第二权重， γ 2表示第三整体流体环境的分子尺寸信息，k2表示第i种组分分子在所有组分中的第三权 重。6. -种井下稠油分子链长的测量装置，其特征在于，包括： 发射模块，用于向填充有稠油样品的核磁共振流体分析仪的探头内发射三维核磁共振 脉冲序列，所述三维核磁共振脉冲序列中包含三个独立的分别用于编辑纵向弛豫时间1^的 信息、扩散系数D的信息、横向弛豫时间!^的信息的窗口； 采集模块，用于采集回波数据； 计算模块，用于采用快速反演算法对所述回波数据进行反演，计算关于所述纵向弛豫 时间T1、所述扩散系数D、所述横向弛豫时间!^的联合概率分布函数； 求解模块，用于根据预先搭建的所述纵向弛豫时间T1与稠油样品中各组分分子链长关 系的模型或所述扩散系数D与稠油样品中各组分分子链长关系的模型或所述横向弛豫时间 T2与稠油样品中各组分分子链长关系的模型，求解所述稠油样品中各组分分子链长的概率 分布函数，以获得所述稠油样品中各组分的分子链长。7. 根据权利要求6所述的装置，其特征在于，所述计算模块包括： 压缩降维单元，用于针对采集的所述回波数据的第一类Fredholm方程，采用奇异值分 解和张量积的算法，对所述回波数据进行压缩降维处理，以获取关于回波数据的降维方程； 求解单元，用于利用零阶Tikhonov正则化方法，求解关于所述纵向弛豫时间TdP所述扩 散系数D的二维联合概率分布函数； 所述求解单元，还用于根据关于横向弛豫时间T2的概率分布函数和所述二维联合概率 分布函数，采用三维反拉普拉斯变换的方法，求解关于所述纵向弛豫时间!^、所述扩散系数 D、所述横向弛豫时间!^的联合概率分布函数。8. 根据权利要求7所述的装置，其特征在于，所述求解单元具体用于： 利用零阶Tikhonov正则化方法，求解使惩罚函数最小的关于所述纵向弛豫时间T1和所 述扩散系数D的二维联合概率分布函数； 其中，所述惩罚函数表示为：X表示为纵向弛豫时间Thy表示扩散系数D，z表示横向弛豫时间Τ2，τι表示第一窗口中的极 化时间，12表示第二窗口中的回波间隔，τ3 = ηΤΕ，ΤΕ表示第三窗口中的短回波间隔，η为预设 的自然数，ki2(xy，T1T2) = ki (T1，X)k2 (T2，y)，ki (T1，X)表示所述三维核磁共振脉冲序列的第 一反演核函数，k2(i2，y)表示所述三维核磁共振脉冲序列的第二反演核函数，k3(z，T 3)表示 所述三维核磁共振脉冲序列的第三反演核函数，〇〈〇〈1^邪(灯，2)表示关于所述纵向弛豫时 间TjP所述扩散系数D的二维联合概率分布函数，#表示惩罚函数，Μ( τιτ2, τ3)表示降维压 缩后的回波数据。9. 根据权利要求8所述的装置，其特征在于，所述第一类Fredholm方程表示为:M( T1，τ2， T3) = iiiki(x,Ti)k2(y,T2)k3(z,T3)f(x,y,z)dxydz+e(Ti,T2,T3)； 所述关于回波数据的降维方程表示为:Μ(τ?τ2，τ3) = Π?α2(Χ7，τ?τ2)1?(Ζ，τ3)?· 2?(Χ7，Ζ) dxydz+e(TiT2,T3)； 其中，M(Hi3)表示回波数据，f(x，y，z)表示三维联合概率密度函数，"^，"表 示噪声项，f2D(xy，z)表示二维联合概率密度函数，ε(τιτ2，τ3)表示压缩降维后的噪声项。10. 根据权利要求6-9任一项所述的装置，其特征在于，所述纵向弛豫时间T1与稠油样品 中各组分分子链长关系的模型表示为：所述扩散系数D与稠油样品中各组分分子链长关系的模型表示为：所述横向弛豫时间T2与稠油样品中各组分分子链长关系的模型表示为：其中，D1J11J21分别为第i种组分的扩散系数，纵向弛豫时间，横向弛豫时间，F为平均 分子链长，N1为第i种组分的分子链长，A为第一经验相关系数，B为第二经验相关系数，β表 示第一整体流体环境的分子尺寸信息，V表示第i种组分分子在所有组分中的第一权重， γι 表示第二整体流体环境的分子尺寸信息，kl表示第i种组分分子在所有组分中的第二权重， γ 2表示第三整体流体环境的分子尺寸信息，k2表示第i种组分分子在所有组分中的第三权 重。
【专利摘要】本发明提供了一种井下稠油分子链长的测量方法和装置。该方法包括：向填充有稠油样品的核磁共振流体分析仪的探头内发射三维核磁共振脉冲序列，三维核磁共振脉冲序列中包含三个独立的分别用于编辑纵向弛豫时间T1的信息、扩散系数D的信息、横向弛豫时间T2的信息的窗口；采集回波数据；采用快速反演算法对回波数据进行反演，计算关于纵向弛豫时间T1、扩散系数D、横向弛豫时间T2的联合概率分布函数；根据预先搭建的纵向弛豫时间T1与稠油样品中各组分分子链长关系的模型或扩散系数D与稠油样品中各组分分子链长关系的模型或横向弛豫时间T2与稠油样品中各组分分子链长关系的模型，求解稠油样品中各组分分子链长的概率分布函数。
【IPC分类】G01N24/08
【公开号】CN105486709
【申请号】CN201510958700
【发明人】肖立志, 邓峰, 刘化冰, 李新, 张宗富
【申请人】中国石油大学(北京)
【公开日】2016年4月13日
【申请日】2015年12月17日
