本公开大体上涉及地层评估。更具体地说，本公开涉及地层评估技术，如地下地层和其中的流体的成像。勘探可以用于定位有价值的烃，如石油和天然气。钻机安置于井场以钻井眼并且在井下部署工具来定位地下储层。井下工具，如钻探工具被推进到井眼中。井下工具，如电缆工具通过电缆部署到井眼中以收集地下样品用于评估。井下工具通常具备测井装置，如核磁共振装置，以对地下地层进行成像。井下工具还具备采样工具，如探针，以将地下流体的样品抽取到井下工具中，以及取芯工具，如轴向和侧壁取芯装置，以切割地下地层的样品。井下工具的实例在美国专利第6047239号和第6897652号中提供。在井下工具中捕获收集的样品并且将其取回到地面。样品被放入实验室用于测试。对样品进行测试以确定烃的存在。在一些情况下，可以使用核磁共振来测试岩心样品。测试的实例在美国专利第9133709号、第8499856号、第6220371号和第4769602号中提供尽管在地层测试和取样方面取得了进展，但仍需要能够精确评估地下地层的技术和工具。技术实现要素：在至少一个方面中，本发明涉及一种用于使地下地层的流体成像的成像器。所述成像器包含壳体，所述壳体具有限定接收穿过其中的所述地下地层的岩心样品的通道的侧壁。所述壳体安置于井下工具中并且具有流体入口以从所述地下地层接收流体并且进入所述通道。所述成像器还包含安置于所述壳体的所述侧壁中，定向为引导磁场通过所述通道的永磁体；安置于所述壳体的所述侧壁中，在所述永磁体与所述通道之间，定向为引导射频场通过所述通道的射频线圈；定位于所述壳体的所述侧壁中，在所述永磁体与所述射频线圈之间的磁场梯度，以选择性地引导梯度场通过所述通道；以及以操作方式连接到所述射频线圈的化学选择性成像器，以根据脉冲序列选择性地脉冲频率，由此产生所述岩心样品的个别流体测量值。在另一方面中，本发明涉及一种使定位于地下地层中的流体成像的方法。所述方法涉及将所述地下地层的岩心样品放置在井下工具中的成像器的通道中，通过使流体样品从所述地层传递到所述通道中来驱替岩心样品，并且使所述岩心样品成像。所述成像涉及在沿所述通道的纵向轴线的方向上引导磁场通过所述通道，选择性地引导梯度场通过所述通道，通过引导射频场以与所述磁场和所述通道的所述纵向轴线正交的方向通过所述通道来选择性地脉冲，以及在所述脉冲期间产生所述岩心样品中的所述流体的图像。最后，在另一方面中，本发明涉及一种使定位于地下地层中的流体成像的方法。所述方法涉及将所述地下地层的岩心样品放置在成像器的流体填充的通道中，引导磁场在沿所述通道的纵向轴线的方向上通过通道，选择性地引导梯度场通过所述通道，引导射频场沿所述通道的所述纵向轴线方向通过所述通道，根据预定的k空间采样图通过选择性地将所述射频场的频率脉冲到所述岩心样品并且将所述磁场梯度施加到所述岩心样品，选择性地获取所述岩心样品中的所述流体的核磁共振测量值，以及通过对所述获取的核磁共振测量进行压缩感测来产生所述岩心样品中的所述流体的图像。附图说明通过参考在附图中示出的本公开的实施例，可以获得上文简要概括的本公开的更具体的描述。然而，应注意，附图说明本公开的示例实施例并且因此不认为限制其范围。图不一定按比例绘制，并且为了清楚和简明起见，图的某些特征和某些视图可能在比例上放大地示出或示意性地示出。图1a是描绘井场的示意图，所述井场具有部署到穿透其中具有流体的地下地层的井眼中的井下工具，所述井下工具在其中具有成像器。图1b是地下地层中的孔隙的放大图。图2是描绘用于使从地下地层取得的岩心样品成像的成像器的示意图。图3a和图3b是描绘用于使岩心样品中的流体成像的成像方法的示意图。图4a到图4c是分别描绘在成像期间产生的1d、2d和3d脉冲序列的示意图。图5a到图5c是岩心样品中的流体的图像。图6是描绘使在地下地层中的流体成像的方法的流程图。具体实施方式在以下描述中，阐述大量细节以提供本发明的理解。然而，所属领域的技术人员应理解，本发明可以在没有这些细节的情况下实施并且所描述的实施例的大量变化或修改是可能的。公开了用于对地下地层进行化学选择性成像的技术。这些技术涉及对地下地层的岩心样品进行磁共振成像(mri)(或核磁共振(nmr)成像)。可以使用能够在岩心样品处选择性地施加以射频振荡的磁场脉冲的装置在原位和/或在表面进行成像。使用单调谐或多调谐探针，化学选择性成像可以使任何核磁共振活性物质成像(例如，1h或23na)。还可以进行对比成像(例如，弛豫和/或弥散)用于比较。可以进行成像以选择性地测量岩心样品中的各种流体，如烃(例如原油或十二烷)和水性流体(例如，水、盐水等)。可以分别或组合地使用这类技术以使地层中的各种流体成像。具体来说，成像可以用于区分岩心样品中的水性流体和烃。这些图像可以用于例如对流体参数，如流动速率和产生的烃的类型特征化。从这种成像收集的信息可以用于例如识别特定流体、分别地使流体成像、评估含有流体的地层、确定井下参数、检测有价值的烃、提供用于规划油田作业的信息等等。可以使用例如1d、2d或3d脉冲序列选择性地使成像流体成像。为了促进成像(例如，减少采集时间)，可以使用各种成像序列，如用于收集数据的经减小的样品尺寸的快速成像(利用弛豫增强(rare)脉冲排序的快速采集)和用于从经减小的样品尺寸重构图像的压缩感测(cs)。快速成像技术可以与压缩感测组合使用以减少图像采集时间，其可以用于例如使工具在井下进行成像所花费的时间减到最少。图1a是描绘用于进行地下操作的井场100的示意图。井场100置于地下地层102周围，地下地层102具有其中具有流体的储层104。地层102具有一个或多个孔隙108，其中具有流体。如图1b孔隙108的放大图所示，流体的部分106a可以位于孔隙108内的中心，而部分106b可以位于孔隙108的凹部110中。中心部分106a可以表示可获取流体并且部分106b可以表示孔隙108内的截留流体。如所展示，井场100包含钻机112、井下工具114和地面单元116。井下工具114被部署到井眼118中以测量地下参数。如所展示的井下工具114是经由电缆线120部署到井眼118中的电缆工具，但可以使用任何井下工具(例如，钻、盘管、生产和/或其它工具)。电缆线120与地面单元116连通，用于在其间传递信号。地面单元116可以用于从井下工具114收集数据和/或将信号(例如，动力、指令等)发送到井下工具114。图1a的井下工具114可以包含用于进行各种操作的各种部件。如所展示，井下工具114包含电子件122a、流体取样器122b、岩心取样器122c、成像器122d和收集器122e。电子件可以包含各种器件，如动力、控制、处理、通信(例如遥测)和/或在井下操作中使用的其它器件。流体取样器122b可以是传统的取样工具，其能够将流体样品从地下地层抽吸到井下工具114中。岩心取样器122c可以是传统的侧壁取芯工具，其能够从围绕地层的井眼的壁切割岩心样品125b。尽管图1a示出了侧壁取芯实例，但可以在井下工具114中设置轴向取芯工具。收集器122e可以是用于存储收集的样品的容器。在美国专利第6047239号中提供了具有采样能力的电缆工具的实例。如所示的成像器122d包含壳体124，壳体124具有侧壁126，侧壁126限定穿过其中的通道128。流体取样器122b和核心取样器122c可以安置在成像器122d的周围(例如，在其上方)，以使流体样品125a和岩心样品125b分别通过入口128a并且进入通道128以进行测量，如由箭头示意性地指示。成像器122d可以包含例如磁共振成像器(mri)，其经配置以接收流体和岩心样品125a、125b，如本文中进一步描述。可以使用成像器122d使岩心样品125b成像。岩心样品125b可以传送通过位于壳体中的通道入口128a(例如，门)，通道入口128a设置成将岩心样品接收到通道中，因此岩心样品可以安置在成像器122d的通道128中。如由箭头指示，使用岩心驱替，岩心样品125b可以用流体样品(或其它流体)125a饱和。流体和岩心样品125a、125b可以从通道128通过出口128b(例如，门)释放(例如，下落)并且进入收集器122e以用于存储和/或取回。井下工具114可以具备各种装置以促进和/或控制采样和/或成像。举例来说，流体样品125a可以自由地传送通过成像器122d，或可以使用流体控制装置，如流线、阀门等来控制流体样品125a。成像器122d可能能够在井下工具114内的岩心样品上进行mri。在成像期间，岩心样品125b可以饱和有采样流体。由成像器122d取得的测量值可以由地面单元116收集，并且可以由其产生输出129，如曲线，如本文进一步所描述。可以使用传感器和/或其它装置进行额外测量以确定各种地下参数，如井下条件、地层参数、流体参数等。举例来说，井下工具114可以具备可以收集附加测量值，如组成、温度、压力等的光学流体分析仪、仪表、光谱仪、变换器等。可以评估所收集的图像和/或测量值以确定各种地下参数。图2是对比成像器222的横截面图，其可以安置在井下(参见例如图1a的成像器122d)或地面位置(例如，单机、实验室设备等)。如此视图中所示，对比成像器222包含壳体224，壳体224具有限定通道228的侧壁226，其可以类似于图1a的成像器122d。成像器可以在任何方向上定向以便于操作。通道228成形为接收岩心样品225b(例如，图1a的岩心样品125b)和流体225a(例如，图1a的流体样品125a)。岩心样品225b可以是圆柱形样品，其可以如所展示安置于壳体224的线性通道228中，或可以具有其它形状和/或尺寸。可以提供装置以将一个或多个样品自动插入和/或移出通道228。如轴向箭头所指示，在测试期间，流体225a可以传送通过通道228。可以任选地提供流线和/或其它流量控制装置以选择性地将流体225a传送到通道228中以提供所需的岩心驱替。流体225a可以以所需的流动速率流动通过通道228，或者封闭在其中以在测试期间保持静止。根据需要，可以从通道228倾倒流体225a。流体225a可以是在测试期间传送通过岩心样品225b的任何流体。在一个实例中，流体225a是来自地层的原位流体，其用于复制地下条件。如由弯曲箭头指示，成像器和/或岩心样品225b可以任选地旋转(例如，通过由马达驱动的旋转轴)，以在成像期间将岩心样品225b的取向从角α1改变为角α2。岩心样品225b相对于壳体224的旋转允许沿基于场b0和b1的方向的取向的坐标信号选择到岩心样品225b旋转时的取向。壳体224具有传感器阵列，其包含永磁体230a、射频线圈230b和施加的磁场梯度230c。永磁体230a安置于侧壁226中并且围绕通道228径向设置，以将岩心样品225b环绕在其中。永磁体230a可以是任何永磁体，如海尔贝克(halbach)磁体，其经布置以产生定向到通道228的z轴的磁场(b0)。射频线圈230b安置于侧壁226中并且围绕通道228径向设置，以将岩心样品225b环绕在其中。射频线圈230b安置于通道228与永磁体230a之间。射频线圈230b可以是布置成产生沿通道228的x或y轴以射频振荡的磁场b1的线圈。磁场梯度230c定位于永磁体230a与射频线圈230b之间。壳体224可以具备成像单元216(例如，图1a的地面单元116)或与其耦合以用于提供动力、收集数据和/或向成像器222发送指令。磁性线圈230a、射频线圈230b和磁场梯度230c可以与成像单元216耦合以向其提供测量值。成像器222和/或成像单元216可以具备通信装置，如有线和/或无线耦合，以在其间定义通信链路。成像单元216可以具有能够将测量值转换成图像以供显示的传统显示器。成像单元216可以包含例如处理器、数据库、遥测单元、电力单元和/或用于与成像器222一起操作的其它电子器件。可将成像单元216并入到井下工具(例如，图1a的122a)和/或地面单元116(图1a)的电子件中。收集的测量值可用于产生输出，如曲线229。任选地，可以提供一个或多个探针234和/或传感器s以收集测量值。举例来说，射频线圈230b的探针234可以选择性地具备单一和/或多个谐振频率，例如以允许检测多个核。可以使用成像单元216从成像器222收集图像参数(例如，流体的分布、残油饱和度等)和来自传感器(s)的地下参数(例如，组成、温度、压力等)。成像和/或收集的测量值可以用于进行各种地层评估，如成像、流体分析、流出物分析、压缩感测等。举例来说，收集的数据可以用于导出地下参数，如电阻率和磁导率。对比成像可以使用各种对比成像技术，如弛豫和弥散成像来进行评估。这类技术可以涉及例如对所产生的图像的弛豫时间t1和t2的分析。可以使用例如依赖于nmr可测量量的差异的技术产产生像参数，如弛豫时间和扩散系数(d)，以提供烃与水性流体之间的对比。nmr测量值可以包含在岩石物理工作中使用的基线测量值，如弛豫时间t2。t2可以是通过施加射频激励在系统中产生的体磁化衰减的量度。磁相干的衰减可能是由核自旋与由静电场不均匀性以及分子间和分子内运动产生的变化磁场的相互作用引起的。在多孔岩石环境中，烃和水性流体可能具有相似的t2。因此，用于提供流体相之间对比的技术可以探测第二流体特性，如弛豫时间t1和各种流体的扩散系数。t1可以是流体分子与环境交换能量的程度的量度。长t1可以指示弱耦合，而短t1可以指示强耦合。因此，t1弛豫时间可能取决于分子特性，如尺寸，并且较大的烃分子可能展现较长的t1。流体，如水性流体、液态烃和气态烃的自扩散系数可能是相当不同的，并且可以用于区分岩石样品中存在的流体相。在这些情况下，可以使用绘制t1对t2或d对t2的多维弛豫测量来提供所需的对比度。可以使用nmr岩心分析和/或光谱方法进行图像的评估。这类评估可以用于在块状样品上提供所需的流体相区分。可以对例如标准自旋回波成像序列，如自旋翘曲(spin-warp)进行这些评估。自旋翘曲的实例描述于edelstein,w.a.,hutchison,j.m.s.,johnson,g.和redpath,t.，《自旋翘曲nmr成像和人体全身成像上的应用(spinwarpnmrimagingandapplicationstohumanwhole-bodyimaging)》，《医学和生物学中的物理学(physicsinmedicineandbiology)》25，751(1980)]。为了提供关于这些流体的空间分布的信息，在岩心驱替期间一维空间分辨的t2分布可用于提供流体鉴别。为了在流体之间的t2对比度低时提供期望的分离，可以进行多维弛豫测量，如d-t2和t1-t2。这些可以提供容积测量。关于超过相对体积的烃和水性流体的分布的进一步信息可以如本文中进一步描述的那样进行。可以通过在注入的水性流体中使用的化学掺杂剂以提供弛豫对比度来确定单一岩心塞中的相的空间分布。在一个实例中，可以使用含有如cu2+、mn2+或gd3+的物质的化学掺杂剂。这些物质可用于减少水性流体的弛豫时间，并且提供各种流体，如水性流体与烃之间的区别。在另一个实例中，在展现不同t1值的系统中，可以使用t1归零来抑制来自样品中存在的t1环境之一的信号。rf激励脉冲的时序可以被设置为使得来自一个t1环境的磁化强度和所得mri信号被信号抑制。岩心样品可以饱和有具有样品流体的多个流体相，如具有单一明确定义的t1的流体。在又一实例中，nmr测量的化学选择性可以用于区分地层流体。给定物质的nmr响应取决于自旋的旋磁比，这是对于每种nmr活性物质唯一的量。在第一种情况下，可以在注入的盐水中使用d2o代替h2o以去除水性流体对图像的贡献。在另一情况下，可以对存在于地层流体中的氢(或其它nmr活性核，如钠23na或碳13c)进行成像。对比成像的实例提供于mitchell,j.,chandrasekera,t.c.,holland,d.j.,gladden,l.f.和fordham,e.j.，《岩石物理岩心分析中的磁共振成像(magneticimaginginpetrophysicalcoreanalysis)》，《物理学报告(physicsreports)》，526，第165-225页(2013)。其它现有技术可用于评估，如描述于美国专利第9133709号、第8499856号、第6220371号、和第4769602号中的那些现有技术。化学选择性成像还可以在地层内使用成像器(例如，分别为图1a和图2的成像器122d和222)使用化学选择性成像技术进行评估，以产生流体，如烃和水性(盐水)的独立图像。化学选择性成像技术利用nmr波谱中化学位移的差异来区分岩心样品中的流体。化学选择性成像包括：1)基于核磁共振(nmr)波谱中化学位移的差异将烃图像与水性(盐水)图像对比，和2)在代表性的储层流动速率下(例如，在vi＝1英尺天-1(0.304米天-1)下)的排水和渗吸实验期间，在减少连续油-水图像之间的像素模糊的时间尺度上获取1、2或3d图像。为了实现这一点，mri脉冲序列(例如，具有弛豫增强的快速采集(rare))可以与压缩感测(cs)组合使用。图3a和图3b是描绘化学选择性成像方法300a、300b的流程图，其可以分别使用图1a和/或图2的成像器122d和222来进行以使岩心样品中的流体成像。图3a的方法300a包含测量阶段336a、信号选择阶段336b、图像采集阶段336c和图像显示阶段336d。可以对岩心样品中的一种或多种流体进行方法300a。如由虚线和平滑箭头337a、337b指示，可以对岩心样品中的一种或多种流体(例如，337a水性流体、337b烃)选择性地重复部分或全部方法300a。测量阶段336a涉及收集测量值，如由成像单元216产生的图2的曲线229。如所展示的曲线229绘制由成像器222产生的信号强度(a.u.)(y轴)对频率(hz)(x轴)。所得线展示对应于流体组成的峰340a、340b。在所示的实例中，峰340a、340b对应于水性流体(w)和烃(d-十二烷)，其中峰分别在0hz、300hz。信号选择阶段336b涉及对水相或烃相磁化的选择性激发。可以进行选择以指示要使哪种流体成像。举例来说，当根据线337a进行水性流体的方法300a时，可以使用框342a选择水峰340a。在另一个实例中，当根据线337b进行烃的方法300a时，可以使用框342b选择烃峰340b。图像采集阶段336c涉及使用mri脉冲序列获取与贯穿岩心样品的流体分布对应的原始k空间数据344。样品图案344是kp1m-1(y轴)对kp2m-1(x轴)的曲线，其指示在图像采集期间必须测量的数据点。曲线344上的数据点指示待获取的点的位置。亮区的强度指示数据采样的位置；而暗区指示未采样。可以使用快速数据采集来捕获此数据，并且使用压缩感测来重构数据以产生如在本文中进一步描述的图像。图像显示阶段336d涉及产生岩心样品内的地层流体的图像346。取决于可用时间，可以使用标准成像技术获取图像346，或者通过利用所获取数据的压缩感测重构的快速采集来产生图像。尽管示出了3d图像346，但图像可以是1d或2d图像。可以显示一种或多种流体的一个或多个图像，如本文中进一步所描述。当产生图像时，图像采集阶段336c可以任选地以各种角度进行。如图3b所示，方法300b可以涉及脉冲排序400。方法300b涉及相同的测量阶段336a，组合信号选择和图像采集阶段336b、336c和图像显示阶段336d。对于各种流体，如由箭头337a、337b指示，重复方法300b的一部分。因为脉冲序列400可以选择性地捕获特定流体的数据，所以脉冲序列400可以用于确定哪些流体正在成像。可以以不同的脉冲频率重复方法300b来激励所选择的流体，如所展示的水346a和烃346b。图4a到图4c是描绘各种脉冲序列400a、400b、400c的曲线图，其可以用于在图像采集阶段336c期间获取图3的烃和水相分布的原始数据。图4a示出了1d序列。图4b示出了2d序列。图4c示出了3d序列。每个脉冲序列400a到400c包含来自rf线圈的射频脉冲(rf)和一个或多个所施加的磁场梯度(gr,gp/gp1,gs/gp2)(例如，图2的230b、230c)。可以应用脉冲序列以选择性地激发样品中存在的特定化学物质(例如油或盐水)。磁场包括读取梯度(gr)和相位梯度(gp/gp1,gs/gp2)，其被施加一段时间以实现核自旋的空间编码。每个脉冲序列还包含代表化学选择性预处理的部分449a和代表脉冲激发部分的部分449b。以各种形状并且在不同角度下进行脉冲序列400a到400c，以产生被成像的样品的不同视角。每个脉冲序列包含：psel-选择性激励脉冲、pex-非选择性激励脉冲、pref-重聚焦脉冲。举例来说，对于在快速采集中使用的脉冲序列，180°重聚焦rf脉冲可以重复nrf次，以从系统的单次激励中对多行k空间采样。如图4a到图4c中所示，每个脉冲序列具有不同的形状。射频(rf)脉冲产生由极化h原子诱导的自旋回波450。可以使用不同的相位梯度(gp/gp1,gs/gp2)来实现自旋的空间编码。如图4a的1d形式所示，脉冲序列400a包含用于激励的rf脉冲和使得仅在读取梯度的方向上实现空间分辨信息的读取梯度(gr)。可以以微秒发射的脉冲将rf场添加到磁场中。射频线r.f.上的脉冲的形状包含具有宽带的方形脉冲，其影响整个脉冲。读取梯度gr仅重复一次以产生1d图像。图4b示出了包含rf脉冲和gr脉冲的脉冲序列400b的2d形式，其具有附加的gp和gs脉冲。射频线r.f.上的脉冲的形状是高斯形，以仅影响样品的特定区域。这种形式还描绘了沿梯度线gp梯度迭代s1、s2。这些梯度迭代指示重复gp以便产生2d图像。对于迭代s1、s2，相位梯度(gp)迭代nrf次。180°重聚焦脉冲可以重复nrf次。可以按需要通过改变各种梯度的强度来获取信息。每次读取梯度(gr)和切片梯度(gs)相同时，可以改变相位梯度(gp)的幅度，并且随后通过各种梯度值(s1和s2)迭代以产生2d图像。图4c示出了包含rf脉冲和读取梯度gr的脉冲序列400c的3d形式，其具有额外的第一和第二相位编码梯度gp1和gp2。在这种形式中，gp1和gp2脉冲各自包含梯度迭代s1、s2，指示这些脉冲重复nrf次以便产生期望的3d图像。对于3d图像，可以使用相同的读取梯度迭代gp1与gp2的所有组合。脉冲序列400a、400b、400c描绘了示例性化学选择性rare脉冲序列。脉冲序列400a到400c具有在读取方向(kr)上的编码的k空间频率和在kp1和kp2上的相位编码，如图3a的图像344中所描绘。在施加读取梯度(gr)和相位梯度(gp/gp1,gs/gp2)时，产生曲线344的点。对于3d脉冲序列(图4c)，相位编码梯度的幅度gp1和gp2确定对图344上的哪些数据点进行采样。图像显示阶段336d使用压缩感测来从岩心图像重构流体图像。可以改变梯度的持续时间和幅度以控制所获取的空间分辨信息。通过改变这些，可以调节视场和成像区域以及分辨率以提供关于样品的所需聚焦级别。mri采集技术可用于欠采样数据以减少数据收集时间，并且压缩感测可用于从欠采样数据重构完整图像。图像采集技术可以包含例如具有弛豫增强(rare)的快速采集、回波平面成像(epi)、自旋翘曲和/或其它采集技术，并且通过使用欠采样和压缩感测重构可以进一步提高时间分辨率。rare的实例描述于hennig,j.,nauerth,a.和friedburg,h，《rare成像：用于临床mr的快速成像方法(rareimaging:afastimagingmethodforclinicalmr)》，《医学中的磁共振：医学磁共振学会官方期刊/医学磁共振学会(magneticresonanceinofficialjournalofthesocietyofmagneticresonanceinmedicine/societyofmagneticresonanceinmedicine)》3，823-833(1986)；epi的实例描述于mansfield,p.《使用nmr自旋回波的多平面图像形成(multi-planarimageformationusingnmrspinechoes)》，《物理学期刊c：固态物理学(journalofphysicsc:solidstatephysics)》10，l55-l58(1977)；自旋翘曲的实例描述于edelstein,w.a.,hutchison,j.m.s.,johnson,g.和redpath,t.，《自旋翘曲nmr成像和人体全身成像上的应用》，《医学和生物学中的物理学》25，751(1980)。mri中的压缩感测基于以下要求，例如，(1)线性重构中的混叠伪像(例如，样品数据)必须是非相干的并且类似噪声；(2)所需图像展现变换稀疏性；(3)使用非线性算法重构图像，所述算法强制稀疏性和与获取的k空间数据的一致性。假设重构的图像由x给出，其通过以下等式(1)与所获取的k空间测量值相关：sfx+v＝y，(1)其中s是子采样模式，f是将图像映射到k空间的傅里叶变换，v是正态分布的噪声(标准差σ和零均值)，并且y是含有所获取的k空间的测量值的向量。由于欠采样和噪声的存在，等式1可能是一个不适定的问题，并且因此，可用于完全采样k空间数据集的线性图像重构方法，如傅里叶逆变换，可能导致由于违反奈奎斯特准则而含有混叠伪像的图像。因此，可以通过使用变分正则化方法来寻求x的近似解，所述方法根据由以下给出的正则化函数j平衡模型(等式1)和x的先验假设：不等式约束的作用强制与所获取的k空间数据的一致性和正则项j并入xσ的重构的先验信息，其可能对抵消问题的任何不适定性是所需的。在cs重构的情况下，先验信息是可以隐式地或在适当的变换域中稀疏地表示图像。用于将图像映射到变换域的正则化函数(j)的选择可以取决于待重构的图像的性质。举例来说，非平滑正则化器(如全变差(tv))可能更适合于具有锐边的图像，而平滑正则化器(如多贝西(daubechies)小波变换)很适合于像素强度改变较平缓的图像。在本研究中，tv已被用作正则化函数j(x)，现将对此进行讨论。总变差惩罚图像的梯度的有限不同近似的2-范数的1-范数，如下所示：在此，为cs重构设定纽曼(neumann)边界条件。可以使用其它正则化函数，如小波变换。一般来说，用于对x近似的吉洪诺夫正则化(tikhonov-regularization)方案写成如下所示：并且正则化参数α(总是正的)对等式4中的保真度和正则化项的影响进行加权。在本研究中，如等式(5a)、(5b)所述，已经实现了包含布莱格曼(bregman)迭代的等式5的修改：使用布莱格曼方法，在每次迭代(等式5b)之后利用添加到k空间数据y的残差来解决一系列k问题(等式5a)。参见例如m.benning,l.f.gladden,d.j.holland,c.-b.schonlieb,t.valkonen，《速度编码mri测量的相位重构-稀疏促进变分方法的调查(phasereconstructionfromvelocity-encodedmrimeasurements-asurveyofsparsity-promotingvariationalapproaches)》，《磁共振期刊(journalofmagneticresonance)》，238(2014)26-43。压缩感测可以与mri脉冲序列组合使用，例如cs-rare，其中k空间欠采样。压缩感测(cs)的实例描述于lustig,m.,donoho,d.l.,santos,j.m.和pauly,j.m，《压缩感测mri(compressedsensingmri)》，《ieee信号处理杂志(ieeesignalprocessingmagazine)》25，72-82(2008)；和lustig,m.,donoho,d.和pauly,j.m.，《稀疏mri：压缩感测在快速mr成像中的应用(sparsemri:theapplicationofcompressedsensingforrapidmrimaging)》。《医学中的磁共振：医学磁共振学会官方期刊/医学磁共振学会》58，1182-95(2007)。图像采集和压缩感测技术可以与地下材料的成像组合使用。地下材料成像的实例包含chang,c.t.,edwards,c.m.,1993，《多孔介质中的流体的质子mr二组分化学位移成像(protonmrtwo-componentchemicalshiftimagingoffluidsinporousmedia)》，《测井分析者(theloganalyst)》，34，第20-28页；dereppe,j.m.,moreaux,c.，《流体填充的多孔岩石的化学位移成像(chemicalshiftimagingoffluidfilledporousrocks)》，《磁共振成像(magneticresonanceimaging)》，9，第809-813页(1991)；dereppe,j.m.，moreaux,c.，《用于分析石灰石中的油-水替换的2d自旋回波和3d化学位移成像技术(2dspin-echoand3dchemical-shift-imagingtechniquesforanalysisofoil-waterreplacementinlimestone)》，《磁共振期刊》，91，第596-603页(1991)；maudsley,a.a.,hila,l,s.k.,perman,w.h.,simon,h.e.，《利用“四维”nmr的空间分辨高分辨率波谱学(spatiallyresolvedhighresolutionspectroscopyby“four-dimensional”nmr)》，《磁共振期刊》，51，第147-152页(1983)；和dechter,jamesj.,komoroski,richarda.,ramaprasad,s.，《使用预饱和对沙岩和碳酸盐芯部中的个别流体的化学位移选择性成像(useofpresaturationforchemical-shiftselectiveimagingofindividualfluidsinsandstoneandcarbonatecores)》，《磁共振期刊》，93，第142-150页(1991)。化学选择性预处理部分449a和rare成像脉冲序列部分449b可用于促进采集。通过使用化学选择性射频(r.f.)脉冲和乱场(homospoil)梯度，可以在脉冲序列的成像部分之前有效地抑制来自各种流体(例如，烃和水性流体)的信号。图4b和4c分别示出了用于2d和3d采集的rare脉冲序列。对于2d应用，可以进行相位编码方向(p1)上的k空间的欠采样，而对于后者，可以在两个相位编码方向(p1和p2)上进行k空间的欠采样。在这两种情况下，k空间在读取方向(r)上被完全采样。使用压缩感测重构欠采样的k空间数据。在2d脉冲测序实例中，图4b的化学选择性部分449b包含用于选择性激发其中一相的高斯形rf脉冲(r.f.)，接着是乱场梯度(gr,gp,gs)来破坏磁化。为了展示化学选择性成像，已获取二维(2d)切片图像。第一图像可以包含未经化学选择性预处理的水性流体和烃的参考图像。第二图像可以独立地包含水性流体和烃。表1描述了用于产生2d图像的实验参数：表1-2d表1表明，通过抑制来自烃的信号，仅检测到水性流体。相反，通过抑制来自水性流体的信号，仅检测到烃。每次采集施加的180°重聚焦脉冲的数量可以由rare因子(nrf)确定，并且采样k空间所需的r.f.激励总数(nex)由在二相编码方向上的点np1,2除以nrf得出。可以通过使用快速成像序列(如rare、epi等)使用快速测序产生图像，并且可以通过采用欠采样并且随后使用压缩感测进行图像重构来进一步提高时间分辨率。还可以通过使用压缩感测，如cs-rare来提高标准成像序列的时间分辨率。rare和epi都可以使用从个别激发获取的多行k空间。可以从每次激发获得的数据线的数量的实际限制并且相应地采集时间加速可以通过所研究的样品的弛豫时间来确定。作为本发明主题的流体饱和岩心样品的横向弛豫时间预期在数十到数百毫秒的范围内。考虑具有t2＝150ms，回波时间te＝4ms的水饱和岩心的rare采集，可以合理地建议每次激发可以获取64行k空间。通过使用压缩感测(cs)，可以从低于奈奎斯特率采样的多个测量值中恢复具有稀疏表示的信号，如图像。因此，将cs应用于超快速mri采集，欠采样k空间可能引起采集图像时间的进一步减少，从而能够研究如实验室岩心驱替的时间分辨率更大的动态方法。在如由图4c展示的3d实例中，例如当希望获得关于岩心内的流体分布的信息时可以提供各种脉冲配置。应用于mri方案以监测实验室岩心驱替实验中的流体分布的自旋翘曲、rare和cs-rare技术的比较如下所示：表2-3d(1)自旋翘曲(2)rare(3)cs-rare读取方向上的像素，nr/-256256256第一阶段编码方向上的像素，np1/-128128128第一阶段编码方向上的像素，np2/-128128128每次激励采样的k空间的行数，nex/-16464回收延迟，tr/s1.61.61.6扫描次数，ns/-888k空间采样比，sf/％10010025总图像采集时间(tacq)可以使用等式6计算：使用等式6，表3显示了表2中每种技术的总采集时间：表3(1)自旋翘曲(2)rare(3)cs-rare采集时间，tacq/min34955514使用cs-rare可以使总采集时间减到最少。为此，可以针对理论岩心驱替的情况计算在三个采集时间的过程中注入的流体量，如表3所示。表4列出了相关的样品特性和实验条件：表4塞直径，d/mm38塞长度，l/mm76塞体积，v/ml87岩石孔隙度，φ/％23孔隙体积，p.v./ml20填隙速度，vi/英尺天-11注射剂流动速率，q/mlmin-10.06对于所考虑的三种情况中的每一种，在图像采集过程中注入的流体的孔隙体积(n.p.v)的总量可以使用等式7计算并且在表5中汇总。对于所考虑的三种情况中的每一种，在图像采集过程中注入的流体的孔隙体积(n.p.v)的总量使用等式8计算并且在表5中汇总：表5(1)自旋翘曲(2)rare(3)cs-rare注入的p.v.的量，n.p.v/ml10.500.160.04表5表明，使用压缩感测的快速成像，如cs-rare可以实现样品体积的显著减少，并且这类图像可以更好代表岩心驱替中特定时间点的流体饱和度。图5a到图5c更详细地示出了3d图像346。这些图像可用于分别地描绘在地层孔隙内驻留的每种流体。如这些图中所示，实施化学选择性3dcs-raremri方法以在代表性储层条件下在实验室岩心驱替实验中使用图2的成像器222独立地使岩心样品内的烃和水性流体成像。图4c的3d脉冲序列用于根据图3a的曲线344产生mri测量。344的白色像素确定gp1、gp2的值，其又确定需要获取哪些数据点。随后通过压缩感测处理数据以产生图5a到图5c的图像。图5c示出了由这种技术产生的组合流体图像。将原始岩心切成两块，其中在环境条件下一半在水性流体中饱和，另一半在烃中饱和。图5a示出了对烃相在沿线337b重复的方法300的施加期间产生的岩心样品中的水相的3d流体图像346a。图5b示出了对烃相在沿线337b重复的方法300a的施加期间产生的岩心样品中的烃相的3d流体图像346b的第一半。图5c示出了岩心样品的第二半，水相和烃相346a、346b在同一图上。在图5a到图5c中描绘的实例中，进行化学选择性3dcs-rare采集以用于石灰石岩心塞中吸收的a)十二烷和b)水性流体的选择性成像。a)和b)的采集时间为每个图像14分钟。已将这两个图像组合以提供水性流体和烃的整体图像。下表6示出了用于产生图5a到图5c的图像的岩心样品的实验参数：表6下表7示出了用于产生图5a到图5c的图像的岩心样品的岩石特性。表7水性流体烃岩石类型estaillades石灰石estaillades石灰石塞直径，d/mm3838塞长度，l/mm约35约35塞体积，v/ml4040吸收体积/ml9.710.6可以将使用图3a、图3b的化学选择性成像方法300a、300b产生的图像与弛豫成像方法进行比较用于验证。可以进行一个或多个成像方法。举例来说，在可能不希望使用可能对岩石样品的状况产生负面影响的d2o、任何离子掺杂剂或其它流体的情况下，在可能存在低信号(例如，由于低天然丰度或低旋磁比)的情况下，在特定rf探针可能是优选的情况下等可以使用化学选择性成像方法，以避免来自标准单回波自旋回波成像序列等的长采集时间。可能影响成像选择的其它变化可以包含通过t1归零和标准自旋回波成像序列抑制来自表现出不同t1值的物质的信号，其中可以为每个初始rf激发获取2dk空间的单行。可以进行方法300a、300b的变化。举例来说，可以在其它nmr活性核(如氢、钠等)上进行方法。这可以用于替代基于化学位移分离选择烃或盐水。可以对各种脉冲序列、测量值、图像和/或其它数据进行比较以验证结果。举例来说，可以将弛豫成像结果与化学选择性成像结果进行比较。成像可以与用传感器s，如流体分析仪感测的测量值组合使用，使得可以确定烃和盐水孔隙流体的分布。还可以使用来自光学流体分析仪的测量值来进行流出物分析。地层评价可以用于规划油田作业，如设计提高石油采收率(eor)(例如注入)以促进生产。图6是描绘进行烃操作的示例方法600的流程图。可以使用例如图1a和/或图2的成像器122d、222来进行方法600。方法包括654-将地层的岩心样品放置在成像器中(参见例如图1a和图2)。方法还包括655-通过以下来使岩心样品成像：引导磁场在沿通道的纵向轴线的磁方向上通过通道，选择性地引导梯度场通过通道，引导射频场在与通道的纵向轴线正交的方向上通过通道中的样品，以及选择性地脉冲射频场的频率(图2)。方法还包括656通过选择性地将射频场的频率脉冲到岩心样品和根据预定的k空间采样图将磁场梯度施加到岩心样品来选择性地获取岩心样品中的流体的核磁共振测量值。选择性获取可以包括进行快速采集，如rare。方法还可以包括657-通过对所获取的核磁共振测量值进行压缩感测来产生岩心样品中的地层流体的图像。选择性获取656和/或产生657可以通过以下来进行：从成像获得流体的波谱图像，选择用于化学选择性成像的相，通过从由选择性脉冲产生的曲线中选择样品点来获取图像，以及使用所选样品点的压缩感测来进行岩心样品内的地层流体的重构。方法还可以包括658-进行化学选择性成像以从流体图像中分离如烃的流体，和660-通过使用对比成像产生参考图像并且将产生的图像与参考图像进行比较来验证图像，以及基于验证的图像进行油田作业(例如，eor)。方法可以以任何顺序进行并且可以视需要重复。可以进行部分或全部方法。可以进行其它任选步骤，如还可以包括656-旋转样品并且在各种角度下重复成像(图2)。虽然参考各种实施方案和开发形式描述了实施例，但是应理解，这些实施例是说明性的，并且本发明主题的范围不限于此。许多变化形式、修改、添加和改进是可能的。举例来说，可以使用本文中的一种或多种技术来进行一种或多种成像。可以使用本文中所提供的技术的各种组合。针对本文中描述为单个实例的组件、操作或结构可以提供多个实例。总的来说，在示范性配置中呈现为独立组件的结构和功能性可以实施为组合结构或组件。类似地，呈现为单个组件的结构和功能性可实施为独立组件。这些和其它变化形式、修改、添加和改进可以属于本发明主题的范围。可以对得益于本文中的教示的所属领域的技术人员来说显而易见的不同但等效的方式来修改和实践本发明。此外，不希望对本文中所展示的构造或设计的细节进行限制。尽管根据“包括”、“含有”或“包含”各种组件或步骤来描述系统和方法，但系统和方法也可以“基本上由各种组件和步骤组成”或“由其组成”。每当公开具有下限和上限的数值范围，则特定地公开落入所述范围内的任何数字和任何包含的范围。具体来说，本文公开的每个值范围(具有形式“a到b”或等效地，“a-b”)应理解为阐述包含在更宽的值范围内的每个数字和范围。每当公开仅具有特定下限、仅具有特定上限或具有特定上限和特定下限的数值范围时，所述范围还包含在指定下限和/或指定上限“周围”的任何数值。当前第1页12