Summary
同步快速断层扫描是用来以CO 2 -饱和盐水油藏条件下存在的石灰石动态图像溶解。 100次扫描以6.1微米分辨率在一段2小时服用。
Abstract
地下存储持久性是碳捕获和储存的主要问题。抽的 CO 2成碳酸盐岩储层有权解散地质密封，并允许二氧化碳逃逸的可能性。然而，在油藏条件下溶解过程知之甚少。因此，需要时间分辨实验在孔隙尺度观察和预测溶出的性质和速度。同步加速器快速层析是采取复杂的孔隙结构的高分辨率时间分辨图像的速度远远超过传统的μ-CT的方法。钻石Lightsource粉红色光束被用来以CO 2 -饱和盐水油藏条件下存在的石灰石动态图像溶解。 100次扫描以6.1微米的分辨率超过2小时内服用。的图像进行分割，并使用图像分析和网络萃取测定的孔隙度和渗透性。孔隙沿lengt均匀地提高样品的小时然而，无论是孔隙度和渗透率的增加速率减慢在稍后的时间。
Introduction
碳捕获和储存（CCS）的一个主要关注的是长期贮存的安全1，2。二氧化碳，CO 2注入到地下将在宿主盐水溶解，并形成碳酸3，4，5。该酸性盐水必须与反应并溶解围岩，尤其是当主岩是石灰石6的电位。溶解可以有利，并允许继续地层渗透率7和更大的存储持久8。然而，地质密封完整性可通过该溶解受到损害，并允许CO 2迁移到表面9。存储永久性的准确预测建模因此依赖于盐水岩系统充分理解溶解和分布和流体的运动在地下10，11，12的速度。
然而，性质和溶解在碳酸盐的速率依赖于盐水13，14，15，16和主岩17的两个属性。溶解率也强烈地依赖于盐水的温度和压力6，使实验技术的发展在重要代表性油藏条件下测量复杂的时间依赖的过程。
之前的实验观察到，现场规模的反应率通常要比实验间歇反应器测量18，19低几个数量级。风化，矿物heterogene性，并且在非均相流场不完全混合是这种现象的可能的解释。然而，这是不可能的评估最显著因素而不反应期间不断变化的孔隙空间的直接观察。因此，需要动态孔隙规模的实验，以提供两者见解运输和反应之间的相互影响，并验证的预测模型。
在碳存储应用研究孔隙尺度过程建立的实验方法是X射线微断层（μCT）20，21。 μ-CT有几个好处：它达到下的高空间分辨率到约1微米，其是非侵入性的，并提供了三维图像。石灰石溶解已经研究在芯（〜厘米）比例22和发现岩盐水洗涤反应增加物理异质性。要提前多不同的T的理解ransport和反应条件改变复杂固体和孔结构，有必要测量在孔隙空间几何，拓扑反应引起的变化，并在贮存器的温度和压力在地下岩石系统流，并以更高的分辨率，以详细调查孔隙规模的过程。本文介绍了在复杂的孔隙结构的岩石学活性溶解过程，专注于测量油藏条件下二氧化碳 -acidified盐水和石灰岩之间的时间和空间依赖性反应速率的方法。
已经有多项研究已在复杂碳酸盐23，24，25，26，27看了看反应，但由于实验或成像约束他们已经不是仅限于前，后的反应图像或没有完成在代表性的地下条件。 Menke等 28已经在孔隙尺度历时几个小时，在有代表性的一个含水层在约1公里的深度的温度和压力下进行反应的原位成像动态二氧化碳 -acidified盐水和Ketton石灰石之间。然而，Ketton是大颗粒比较均匀的岩石在很短的时间容易图像（〜17分钟），并与一些预测（〜400）。大多数碳酸盐岩都需要许多预测准确地解决它可以使用传统的μ-CT一个非常耗时的过程孔隙结构复杂 - 无论是与在同步辐射光源的单色光束或台式X射线扫描仪。因此，需要断层的一个快速的方法来查看动态异构碳酸盐反应引起的变化。
的时间量花费图像样本由第通量控制E X射线源。扫描的一种方法很快是使用同步加速器源20的多色光束。这种所谓的'粉红梁“，提供了幅度较多的光级比台式来源，因此，可采取在数万的秒，而不是小时的时间尺度上的图像。一个由偶极磁体的周期结构的波荡产生粉红梁。电子束被迫接受振荡，因为它穿过磁体并因此辐射能量。产生的能量是浓缩缩小波长频带和是非常激烈的。反射镜和过滤器然后用于缩小光的光谱，以适应实验需要。镜吸收高能量光谱而过滤器吸收低能量。因此，可以对光谱缩小到仅使用这些工具的辐射的期望的频带。
但是，使用该强X射线通量是并非没有挑战。该粉红梁频谱的低能量X射线被示例作为热吸收。这可以用的原位装置中的温度控制干扰而引起的 CO 2从溶液中20至exsolve。 的 CO 2 -饱和食盐水是热和压力，因此，在热平衡的变化小可以显著改变的原位流体5的pH值非常敏感。因此，对于X射线光谱仔细设计和控制元件必须纳入之前成像的光束线设备。
快速断层也以高速率产生数据的大量。从照相机和随后的存储读出的数据的限制提供大幅的技术挑战。有些人通过采取几个连续的扫描和读取它们对外部数据服务器之前，相机内存存储他们克服这一点。然而，这需要该实验是relatively短相机内存只能容纳数据的有限体积。分箱相机上的数据也减少，因为它减少了需要被传输的数据量的传送时间，但它具有降低图像的质量的潜力。或者，数据可以在每次扫描后开始下一，这将增加扫描之间的总时间之前传送关闭照相机。该研究中使用与每个图像采集服用〜45秒，并读出取额外〜30秒的数据后一种方法。
当以高速率服用扫描，样品台必须旋转比传统扫描要快得多，因此，在芯保持器的电位角应力是很大的。碳纤维，而X射线透明，强调当是灵活的。如果图像获取图像模糊在采样动作可以发生。芯保持套筒被设计成尽可能的短，以减少这些潜在的应力。此外，灵活的polyeth呃醚酮（PEEK）油管在靠近舞台的实验装置中的所有元素使用，以至于阶段是自由转动。使用的PEEK管的一个缺点是，它是可渗透CO 2上扩散的时间尺度。居住在长时间的线流体将逐渐在一段约24小时成为去饱和。所没有的芯保持器附近的所有系由不锈钢制成的并且流体是预平衡在剧烈混合哈氏反应器加热加压，实验条件23，29，30。
实验设备如图1所示。储温度由包裹套管的外部中的X射线透明加热胶带，并通过电池的径向端口并进入限制流体插入热电偶保持在芯保持器。比例积分ðerivative（PID）控制器然后调节温度到1内 C。压力和流量条件，使用三个高压注射泵是精确到0.001毫升/分钟的流速维持。两种盐被用于实验中，将高吸收25％（重量）的KI不反应盐水和低的吸收1重量％的KCl，5％（重量）的NaCl的反应盐水。在衰减的差异使得它很容易看到反应盐水中的制芯死体积计算不必要的到来。
Protocol
1.成像策略设计
- 计算在最高粉束能量和磁通束线的X射线谱，以便利用实验调谐曲线和测量滤波器的传输进行预测的成像性能。钻石Lightsource I13-2'粉红梁“的X射线谱的一个例子示于图2。
- 作为低能量的X射线引起的样品的加热和不添加到成像对比，滤出的X射线频谱的较低部分，使得只有最高能量的X射线在成像样品使用。选择的那个通过计算在现有的光波长31,32理论滤波器透射吸收在适合于光源的光的所希望的低波长的材料线滤波器。这里，可以使用铝和金在该光源的光束线。
- 使用由combina的带通滤波器的X射线滤波器的高通和一个X射线反射镜接近临界角的低通滤波操作和灰。在这种情况下，使用一组为0.2mm热解碳0.2毫米铝过滤器和用于反射镜的铂涂覆的带材1.15毫弧度的入射角下。反射镜反射低于30千电子伏特和附加在线过滤器被安装为2mm的Al和0.1微米的Au其中在13和22千电子伏有吸收峰，分别筛选更多的低能量X射线的唯一的光。 图3描绘了束线成像装置。
- 选择在现有的光的频率和流量的光束线丰富scintillates闪烁体。在这种情况下，闪烁屏幕是由250微米厚的镉tungstenate（CdWO 4），其上叠加了750微米厚的铅tungstenate（钨酸铅4）。然后选择一个客观的镜头和相机的视适当的领域，抢购时间分辨率为experimental要求。在这种情况下，夫妇与0.04与PCO EDGE 5.5 CMOS摄像的孔并用一1.25X物镜捕获的视4毫米场与0.001秒的帧速率。
- 选择图像采集的'flyscan“技术为舞台旋转这种方法减少了样品的震动。传统采集要求阶段停留在每一角增量，取投影，然后移动到下一个角度。在这些动态tomographies图像获取用'flyscan'，它接受作为舞台移动断层扫描，并假定一个角增量，使得每个相继投影之间的差小完成。该'flyscan'方法消除了起动和停止运动的小振动的效果和更快地提供更高质量的图像。
2.设备和电池的组装
- 装载在用于岩心驱替制备芯进入细胞。
- 首先，包芯在铝箔的一个层和插入一个套筒（ 例如 ，氟橡胶）（ 图4）。
- 切断套筒尺寸，使得它比纤芯和内部端部装配件的组合长度短2毫米。端部装配件是1/16“国家管螺纹（NPT），对已加工成5mm的外径联合接头，而套筒是4毫米内径。
- 伸过来5 mm页端接头套筒来创建一个密封。确保没有端部配件和芯之间的任何空间，以确保围压不超过压缩套管和夹断流动。
- 包裹在铝两个附加层中的配件和套筒既能防止气态的 CO 2扩散进入限制流体并保持套筒就位在接头和防止液压通路从连接围和孔隙流体。
- 把核心HOL滑动管子和密封件放回原处明镜回到一起并通过更换螺栓密封端盖和端部配件。
- 扫描在大约4毫米的宽度和长度的重叠部分中的核心。校准扫描曝光时间为15000左右的平均计数值，以确保高的信号噪声比，而不过度饱和闪烁器。每次取干扫描至少有2400的预测保持相衬和边缘清晰度。
- 取闪烁器，使得任何损害和外部噪声可以重建过程中占的平面和暗图像。由芯保持器移出的视场取单位并考虑只是闪烁体的图像与在横梁上。采取利用与梁关的相同方法深色衣物。
3.系统加压
- 通过反应器容器的顶部倾倒流体在加载在1％（重量）氯化钾（KCl）中氯化物的5％（重量）氯化钠（NaCl）食盐水进入拆卸反应器中。
- 加入粉状碳酸盐岩以实现所需盐水酸度。在这种情况下不加入碳酸酯。
- 通过拧紧螺栓和加热胶带重新打包，并插入温度探头进入前重新组装反应器。
- 负载CO 2的进入喷射泵由（ 图1中V 1）打开阀1。
- 关闭阀1和加压喷射泵100巴。
- 开阀2 CO 2淹没反应堆。加热该反应器50℃，用在组合一个PID控制的加热卷绕温度探头和反对对具有夹带搅拌器由外部电动机驱动tinuously搅拌。平衡用CO 2的盐水在10MPa和50℃2和6小时之间，以确保该盐水是完全用CO 2饱和的和碳酸酯完全溶解。
- 前连接核心保持器，彻底吹扫空气和可能沉淀的系统在先前实验的行。要做到这一点，连接线的上方和芯保持器下面绕过芯保持器（U1和U2）。
- 负载去离子（DI）水成通过设置在接收泵笔芯通过阀11接收泵。
- 开阀7,4和3，并使用在恒压模式的接收泵向后通过该系统，并从反应器下面阀3的驱动DI水。使用大约十系统卷，以确保线条都清晰的空气，并冲洗干净。
- 清空接收泵次，然后通过阀11和负载DI水进入围泵通过阀10装入25％（重量）的KI盐水到接收泵。
- 关闭阀10并打开阀8和6，使用限制泵在2MPa限制的核心。
- 关闭阀门11和加压接收泵10栏。
- 打开阀9，7,4，和3，并使用所产生的压力降来驱动KI通过芯掺杂盐水洗涤。
- 逐步加强围和孔隙压力，直到一个合理的流量建立。通过核心驱动盐水的约两个完整系统卷，并通过在反应器下面阀3排出的流体。在这种方式，所有空气被从系统中排出和芯充斥着高对比度盐水洗涤，使未掺杂的反应盐水的到来易于观察。
- 关闭阀3和逐步增加围和孔隙压力，直到芯12兆帕局限和孔隙压力为10MPa。开关PID控制器上以使芯50 C。
- 停止接收泵，关闭阀3，并打开阀5在反应器的底部到反应器系统连接到的核心。
4.流体流动和图像采集
- 中心芯部的中间在视场并为核心充斥跟踪岩心驱替进展连续服用的2-D投影。通过打开粉光束，并使用相机拍摄的图像不旋转阶段采取的2-D投影。反应液注射前开始的2-D的预测给出了前盐水图像清晰，稍后将相对于后续的盐水填充的图像。
- 设置接收泵以从而通过磁芯在所需流动条件拉动流体从反应器中的期望的流速再填充而离开喷射泵从前端调节压力。
- 监控的2-D凸起用于衰减该信号功盐水的到来的变化。当反应盐水到达时，芯体的传输将增加，并且随着更多的光照射到闪烁器和掺杂盐水是由高度X-射线透明反应流体位移的二维凸起将大大变亮。如果存在反应性和非反应性盐水之间没有衰减差异然后，取决于束线谱，重新从步骤2.1具有较高盐浓度的KI的实验或使用不同的高度吸收的盐可以是必需的。
- 停止的2-D扫描和3-D取tomographies依次一样快的成像装置允许。使用每次扫描大约1000预测。扫描只使用180°的旋转（相对于传统的360度）的核心。同时使用旋转更少度降低了信噪比，这是更快的，并有助于避免拉伸和缠结的流动和电气线路。以3-D扫描，直到eith器的时间限制达到或芯看起来充分溶解有内部结构倒塌的迫在眉睫的危险（从而导致这两个围压和未来整个芯干扫描数据的损失）。
- 上次扫描拍摄后，有效地将系统减压，以避免任何核心的进一步反应。
- 首先停止接收泵。然后关闭阀5的反应器连接到系统的其余部分。
- 步骤的系统压力下使用约束和接收泵周围兆帕的压力保持在限制液体。
- 一旦内达到大气压力为1MPa，打开围和接收泵用阀10和11以及在恒流模式运行排出任何剩余流体。
- 关闭PID控制器和在芯保持器以释放任何残留的系统压力的顶部打开工会4路（U2）。
- 慢慢松开围线而醒目过量的二围水用吸水纸。关闭阀6和7和断开工会1和电力线。
- 松开阶段夹，从舞台上拆下核心持有人。
- 小心地从芯保持核心组件，然后断开内部端接头套管。不要从袖子去掉内核，因为这样做可能会破坏脆弱的反应的核心。放置在一个完整的去离子水的烧杯中的套筒覆盖的核心，以稀释的任何潜在的反应盐水洗涤，并停止所有反应。
- 干燥在60℃的烘箱整个芯为至少12小时。然后用传统的样品安装重新安装在舞台上的核心，并在相同的分辨率和预测作为初始干扫描扫描一遍。
5.图像处理
- 通过假设任何影响更正使用多色束关联的任何束硬化的重建图像是放射状符号公制高斯函数33。
- 筛选使用边缘保留滤波器，如非本地化装置以增加信号噪声34，35中的图像（见参考文件 ）。
- 段使用的一个分水岭干扫描图像分割 算法 36 和界定种子岩石和无效的 （见 补充文件 ）。
- 采取具有反应盐水芯的第一图像和每个后续图像注册到第一图像，并使用与第一图像作为参考兰克泽斯37重采样方法重新取样它。作为正在进行的反应趋于模糊的边缘，对图像分水岭分割不足够精确分割。
- 从第一图像中减去每个反应芯图像得到的差图象。段差图像转换成变化和没有变化。寄存器分割干扫描到第一反应扫描，然后减去从分割干扫描分段变化以实现分段反应图像38。
6.造型
- 使用二值图像作为输入到一个无论是直接的纳维-斯托克斯流求解器39，40，或网络抽取模型41（ 图8）来表征通透性的改变，并提供物理见解溶解的动态。
Representative Results
将反应物在毫米直径的4 1.2厘米长的波特兰碳酸芯42方解石和无缓冲SCCO 2饱和盐水之间进行成像。波特兰碳酸盐是相对纯的（<99％）方解石鲕具有复杂异质孔结构43。低能量的X射线通过使光束穿过的Al 2毫米和Au为0.1μm过滤。一个CdWO 4的闪烁体用1.25X物镜和PCO EDGE相机在检测器组件被使用。干扫描，以4000预测收购而动态扫描相依为命1000预测。总采集时间为约1分钟，每次扫描15秒〜100次扫描所历时2小时。
重建和工件的去除采用钻石Lightsource专有软件完成。每个图像由2000年第3体素，WHI沟道然后分级以6.1微米（ 图5）的分辨率，以增加信号到所得的1000 3体素的图像中的噪声。然后，将图像用在Avizo 8.1和ImageJ的节目图像处理模块处理（见补充文件）。每个图像需要大约12个小时CPU和处理3小时GPU采用了3.0 GHz的CPU和特斯拉K20C GPU的计算机上。
分割图像进行分析，作为通过计算孔隙和岩石的体素的数目的孔隙率变化的时间系列。期间随时间溶解孔隙率的增加（ 图6）。划分图像的目视检查（ 图7）示出了一个信道的在流动方向上的存在。当孔隙率被绘制为的时间和从样品的距离的函数入口显然一个信道在第一小时形成，然后加宽作为实验继续（F igure 8）。
划分图像然后作为输入到网络提取模型来分析通透性改变（ 图9）。据发现，有在初始小时期间急剧增加渗透性，但随后的渗透性稳定在稍后的时间。
图1. 原位实验装置。 CO 2由注射泵加压并用于在反应器中平衡盐水洗涤。反应盐水由接收泵通过芯组件拉出。细胞通过在围泵DI水并通过在限制流体的热电偶控制的加热用加热胶带密闭。实验系统被连接在一起用管路和流体流动用阀（V）和联盟（U）的定向。PS：//www.jove.com/files/ftp_upload/53763/53763fig1large.jpg“目标=”_空白“>点击此处查看该图的放大版本。
图2.使用两个实验调谐曲线和理论镜面反射率和过滤器传输所计算的钻石LightsourceⅠ-13I粉束的X射线谱。镜子吸收超过30千电子伏的能量;铝和Au分别过滤器吸收低于13和22千电子伏的能量。 请点击此处查看该图的放大版本。
图3.束线成像装置。在Al和Au片过滤粉红光束，而其余的X射线击中第E磁芯装配。 X射线的一部分被样品吸收，而其余穿过样品并击中该荧光在可见光谱的闪烁体。此可见光然后由物镜聚焦到CCD，它转换该光进入其中像素强度值是由该闪烁体所吸收的X射线的数量的函数的像素化的数字图像。 请点击此处查看该图的放大版本。
图4.芯保持器内的芯组件。 PEEK管附连到内部端部配件并通过钢端盖螺纹。芯包裹在铝箔和插入套管。套管然后，在端部配件拉伸以创建一个水密本身人与铝箔的两个附加层被添加到保持在适当位置的一切，并防止气体扩散。热电偶被固定到芯组件的外部与粘合剂的铝箔的外层。图从Menke 等改性。 42。 请点击此处查看该图的放大版本。
（a）之前和（b）后溶解重构图像的图5 2-D切片。较轻的地区是粮食和暗区是孔隙。在晶/孔边界的边缘模糊所用的孔隙空间（B）的反应的部分可以看到。 请点击此处查看该图的放大版本。
图6.孔隙绘出时间。线性溶出的第二个小时的斜率略有减少孔隙增大。图从Menke 等改性。 42。 请点击此处查看该图的放大版本。
在60分钟时在孔隙率变化到实验，其中绿色代表孔隙率变化最大和红最少的图7. 3-D呈现。由流体 - 固体化学反应产生一个清晰的多孔通道被看作在溶解是最大的芯的中心。图从Menke 等改性。“外部参照”> 42。 请点击此处查看该图的放大版本。
图8概要孔隙率作为距离从样品入口的功能。孔隙率是沿溶出的轴线均匀，但溶解速率的变化作为时间的函数。图从Menke 等改性。 42。 请点击此处查看该图的放大版本。
图9（A）上划分图像进行网络提取显示在60分钟，显示大孔隙（球）和它们的连接（管）。 （B）中所计算的渗透性被示出为增加与时间用40至60分钟之间的急剧上升作为宽溶解信道被建立。图从Menke 等改性。 42。 请点击此处查看该图的放大版本。
Discussion
在储层条件异构孔结构反应的动态影像中最关键的步骤是：1）的粉红色光束内的细胞的精确的温度控制; 2）快速移动阶段成功的核心持有人的稳定性; 3）有效的数据处理和存储技术; 4）时间分辨图像的有效分割。
温度控制是用粉红色光束油藏条件下成像的关键。如果温度上升到高于在反应器温度，CO 2将在孔隙空间exsolve和既改变盐水的pH值，并在可能改变溶出44的性质孔隙空间创建超临界的 CO 2的神经节。使用过滤器来吸收低能量的X射线是用于除去此额外的温度应力允许热电偶和加热缠绕，有效地控制温度在外部的关键。然而，过滤器降低光束的总能量可以通过，因此必须谨慎使用，以便不显著增加总采集时间。此外，滤波器类型和厚度必须进行调整以符合特定能量的波长和吞吐量的束线的。
芯保持器经受断层采集期间旋转和振动应力，可能会导致该碳纤维套筒阶段转动期间摇动和模糊的突起。为了尽量减少这种潜力，芯保持器被设计成一个短6厘米套筒在同步加速器使用。这个套筒不会与台式扫描仪使用有利，因为钢端部装配件将抑制源 - 样品距离和几何放大的最小化。然而，在一个平行光源这些都不顾虑。
在一系列的拍摄各断层扫描可具有尺寸超过20 GB的意思是一个系列100的扫描将在大小2 TB。当在连续服用多种扫描非常归仁ckly两个仪器的带宽和存储选项提供了大量的数据管理挑战。实验成像设备的设计必须考虑到这些约束条件，以充分实现快速断层的动态显像潜力。数据传输瓶颈，必须在开始实验和技术基础设施适应这样如相机问题，读出速度快，传输带宽和存储写入速度不抑制采集速度的潜力之前予以确定。
溶解的时间分辨图像的有效分割提供了一个挑战。当一个断层扫描是在不断变化的系统中的固 - 液边界的边缘可以变得模糊。这使得模糊传统的分割技术，如分水岭，这是假设的界限将是最高的衰减梯度的区域，更成功的作品。为了规避这一点，unreac的差分图像泰德并反应图像被计算，提供变化的唯一的区域的图像。此方法允许连续地改变孔结构的成功分割。
同步加速器加上一个水库规模装置快速层析是一个功能强大的实验方法，该方法可以适用于探索一个范围的应用，包括多相流过程中，对流 - 扩散，并在化学上异构介质传输。然而，目前的装置限制为秒，单相的实验，和小样本大小量级的时间分辨率。未来设计升级可以包括用于三相功能的其他的泵，增加光通量，以便能够穿透较大介质，更好的重建技术，允许更少突起每次扫描进行拍摄，并多元办法图像采集以及可以进一步改善信息分割深度，宽度，和准确性。
Materials
| Name | Company | Catalog Number | Comments |
| NaCl salt | Sigma Aldrich | S7653-1KG | |
| KCl salt | Sigma Aldrich | P9333-1KG | |
| KI salt | Sigma Aldrich | 30315-1KG | |
| Core holder | Airbourne Composites | 110 mm Core holder | Constructed in conjunction with Imperial College |
| PEEK tubing | Kinesis | 1560xL | |
| Thermocouple | Omega Engineering | KMTSS-IM300U-150 | |
| Flexible Heating Tape | Omega Engineering | KH-112/10-P | |
| 1/16" Needle Valve | Hydrasun Ltd | MVE1002 | |
| High Pressure Syringe Pump | Teledyne ISCO | 1000D | |
| 600 mL Parr Reactor | Parr Instrument Company | 4547A - hastelloy | |
| CO2 Cylinder | BOC | CO2 - size E | |
| Viton | Fisher Scientific | 11572583 | |
| Aluminium Foil | Coroplast | 1510AWX | |
| ImageJ - image processing | NIH | ImageJ | |
| Matlab | Mathworks | Matlab | Used for data analysis |
| Avizo | FEI | Avizo | |
| Snoop Leak Detector | Swagelok | MS-SNOOP-8OZ |
References
- Herzog, H., Caldeira, K., Reilly, J. An issue of permanence: Assessing the effectiveness of temporary carbon storage. Clim. Change. 59, 293-310 (2003).
- Metz, B., Davidson, O., de Coninck, H., Loos, M., Meyer, L. Carbon Dioxide Capture and Storage. IPCC. , Cambridge University Press. 431 (2005).
- Langmuir, D., Hall, P., Drever, J. Aqueous Environmental Geochemistry. , Prentice Hall. (1997).
- Morse, J. W., Mackenzie, F. T. Geochemistry of sedimentary carbonates. , Elsevier. (1990).
- Peng, C., Crawshaw, J. P., Maitland, G. C., Trusler, J. M., Vega-Maza, D. The pH of CO 2-saturated water at temperatures between 308K and 423K at pressures up to 15MPa. J Supercrit Fluid. 82, 129-137 (2013).
- Peng, C., Crawshaw, J. P., Maitland, G. C., Trusler, J. M. Kinetics of calcite dissolution in CO 2-saturated water at temperatures between (323 and 373) K and pressures up to 13.8 MPa. Chem. Geol. 403, 74-85 (2015).
- Bachu, S., Nordbotten, J. M., Celia, M. A. Evaluation of the spread of acid gas plumes injected in deep saline aquifers in western Canada as an analogue to CO2 injection in continental sedimentary basins. Proceedings of 7th International Conference on Greenhouse Gas Control Technologies. 1, (2004).
- Bachu, S. Review of CO 2 storage efficiency in deep saline aquifers. Int J Greenh Gas Con. , (2015).
- Marland, G., Fruit, K., Sedjo, R. Accounting for sequestered carbon: the question of permanence. Environ Sci Policy. 4, 259-268 (2001).
- Daccord, G., Lenormand, R., Lietard, O. Chemical Dissolution of a Porous-Medium by a Reactive Fluid .1. Model for the Wormholing Phenomenon. Chem. Eng. Sci. 48, 169-178 (1993).
- Daccord, G., Lietard, O., Lenormand, R. Chemical Dissolution of a Porous-medium by a Reactive Fluid .2. Convection vs Reaction, Behavior Diagram. Chem. Eng. Sci. 48, 179-186 (1993).
- Maheshwari, P., Ratnakar, R., Kalia, N., Balakotaiah, V. 3-D simulation and analysis of reactive dissolution and wormhole formation in carbonate rocks. Chem. Eng. Sci. 90, 258-274 (2013).
- El-Maghraby, R., Pentland, C., Iglauer, S., Blunt, M. A fast method to equilibrate carbon dioxide with brine at high pressure and elevated temperature including solubility measurements. J Supercrit Fluid. 62, 55-59 (2012).
- Fredd, C., Fogler, S. Influence of Transport and Reaction on Wormhole Formations in Porous Media. AIChE. 44, (1998).
- Gharbi, O., Toth, A., Bijeljic, B., Boek, E., Blunt, M. PGE Seminar Series. , Imperial College London. (2013).
- Luquot, L., Gouze, P. Experimental determination of porosity and permeability changes induced by injection of CO2 into carbonate rocks. Chem. Geol. 265, 148-159 (2009).
- Cohen, C. E., Ding, D., Quintard, M., Bazin, B. From pore scale to wellbore scale: Impact of geometry on wormhole growth in carbonate acidization. Chem. Eng. Sci. 63, 3088-3099 (2008).
- Li, L., Peters, C. A., Celia, M. A. Upscaling geochemical reaction rates using pore-scale network modeling. Adv Water Resour. 29, 1351-1370 (2006).
- Swoboda-Colberg, N. G., Drever, J. I. Mineral dissolution rates in plot-scale field and laboratory experiments. Chem. Geol. 105, 51-69 (1993).
- Berg, S., et al. Real-time 3D imaging of Haines jumps in porous media flow. Proc. Natl. Acad. Sci. U.S.A. 110, 3755-3759 (2013).
- Blunt, M. J., et al.
- Ott, H., et al. Core-flood experiment for transport of reactive fluids in rocks. Rev. Sci. Instrum. 83, 084501 (2012).
- Gharbi, O. Fluid-Rock Interactions in Carbonates: Applications to CO2 storage. , Imperial College London. (2014).
- Noiriel, C., Gouze, P., Made, B. 3D analysis of geometry and flow changes in a limestone fracture during dissolution. J Hydrol. 486, 211-223 (2013).
- Hao, Y., Smith, M., Sholokhova, Y., Carroll, S. CO2-induced dissolution of low permeability carbonates. Part II: Numerical modeling of experiments. Adv Water Resour. 62, 388-408 (2013).
- Smith, M. M., Sholokhova, Y., Hao, Y., Carroll, S. A. CO2-induced dissolution of low permeability carbonates. Part I: Characterization and experiments. Adv Water Resour. 62, 370-387 (2013).
- Gouze, P., Luquot, L. X-ray microtomography characterization of porosity, permeability and reactive surface changes during dissolution. J. Contam. Hydrol. 120-121, 45-55 (2011).
- Menke, H. P., Bijeljic, B., Andrew, M. G., Blunt, M. J. Dynamic Three-Dimensional Pore-Scale Imaging of Reaction in a Carbonate at Reservoir Conditions. Environ. Sci. Technol. 49, 4407-4414 (2015).
- Andrew, M., Bijeljic, B., Blunt, M. Reservoir Condition Pore-scale Imaging of Multiple Fluid Phases Using X-ray Microtomography. J Vis Exp. , (2015).
- Andrew, M., Bijeljic, B., Blunt, M. J. Pore-scale imaging of trapped supercritical carbon dioxide in sandstones and carbonates. Int J Greenh Gas Con. 22, 1-14 (2014).
- Henke, B. L. Filter Transmission. , Available from: http://henke.lbl.gov/optical_constants/filter2.html (2015).
- Henke, B. L., Gullikson, E. M., Davis, J. C. X-ray interactions: photoabsorption, scattering, transmission, and reflection at E= 50-30,000 eV, Z= 1-92. At. Data Nucl. Data Tables. 54, 181-342 (1993).
- Schlüter, S., Sheppard, A., Brown, K., Wildenschild, D. Image processing of multiphase images obtained via X-ray microtomography: A review. Water Resour. Res. 50, 3615-3639 (2014).
- Buades, A., Coll, B., Morel, J. -M. A non-local algorithm for image denoising. Computer Vision and Pattern Recognition, 2005. CVPR 2005. IEEE Computer Society Conference on. 2, 60-65 (2005).
- Buades, A., Coll, B., Morel, J. -M. Nonlocal image and movie denoising. Int J Comput Vision. 76, 123-139 (2008).
- Sheppard, A. P., Sok, R. M., Averdunk, H. Techniques for image enhancement and segmentation of tomographic images of porous materials. Physica A. 339, 145-151 (2004).
- Lanczos, C. An iteration method for the solution of the eigenvalue problem of linear differential and integral operators. , United States Governm. Press Office. (1950).
- Andrew, M., Menke, H., Blunt, M. J., Bijeljic, B. The Imaging of Dynamic Multiphase Fluid Flow Using Synchrotron-Based X-ray Microtomography at Reservoir Conditions. Transport Porous Med. , 1-24 (2015).
- Raeini, A. Q., Blunt, M. J., Bijeljic, B. Modelling two-phase flow in porous media at the pore scale using the volume-of-fluid method. J. Comput. Phys. 231, 5653-5668 (2012).
- Bijeljic, B., Raeini, A., Mostaghimi, P., Blunt, M. J. Predictions of non-Fickian solute transport in different classes of porous media using direct simulation on pore-scale images. Phys Rev E. 87, 013011 (2013).
- Dong, H., Blunt, M. J. Pore-network extraction from micro-computerized-tomography images. Phys Rev E. 80, 036307 (2009).
- Menke, H. P., Bijeljic, B., Andrew, M. G., Blunt, M. J. Dynamic Pore-scale Imaging of Reactive Transport in Heterogeneous Carbonates at Reservoir Conditions. Energy Procedia. 63, 5503-5511 (2014).
- Bijeljic, B., Mostaghimi, P., Blunt, M. J. Signature of non-Fickian solute transport in complex heterogeneous porous media. Phys. Rev. Lett. 107, 204502 (2011).
- Ott, H., Oedai, S. Wormhole formation and compact dissolution in single-and two-phase CO2-brine injections. Geophys. Res. Lett. 42, 2270-2276 (2015).
