Abstract
A tomografia de raios X foi usado para a imagem, a uma resolução de 6,6 mm, o arranjo de poros escala de gânglios residual de dióxido de carbono no espaço de poros de uma rocha de carbonato a pressões e temperaturas representante de formações mais utilizados para a armazenagem de CO 2. Equilíbrio químico entre as fases de CO 2, e salmoura rocha foi mantida utilizando uma alta temperatura do reactor de alta pressão, reproduzindo as condições de longe do local de injecção. O fluxo de fluido foi controlado usando bombas de seringa alta temperatura alta pressão. Para manter representativos condições in situ dentro do scanner de micro-TC uma fibra de carbono de alta pressão de micro-CT coreholder foi usado. Diffusive CO 2 troca em toda a manga confinando do pore-espaço do rock para o fluido de confinamento foi impedido por em torno do núcleo com um envoltório triplo de papel alumínio. Reconstruído contraste salmoura foi modelado utilizando uma fonte de raios-x policromática, e composição da salmoura wescolhido para maximizar o contraste de três fases entre os dois fluidos e a rocha. Linhas de fluxo flexíveis foram usadas para reduzir as forças sobre a amostra durante a aquisição de imagem, potencialmente causando movimento indesejável da amostra, uma grande falha em técnicas anteriores. Um termopar interno, colocado directamente adjacente ao núcleo rocha, juntamente com um envoltório flexível de aquecimento externo e um controlador PID foi utilizado para manter uma temperatura constante no interior da célula de fluxo. Quantidades substanciais de CO 2 foram presos, com uma saturação residual de 0,203 ± 0,013, e os tamanhos dos gânglios maior volume obedecer distribuições lei de potência, de forma compatível com a teoria de percolação.
Introduction
Captura e armazenamento de carbono é o processo em que o CO 2 é capturado a partir de fontes pontuais de grande e armazenado em rochas porosas, deslocando brines residentes para que ele permaneça na subsuperfície para centenas de milhares de anos 1. O CO 2 reside no subsolo como uma fase super-crítico densa (SCCO 2), com propriedades radicalmente diferentes ao CO 2 em condições ambientais. Existem quatro mecanismos principais pelos quais SCCO 2 pode ser imobilizadas na subsuperfície: estratigráfica, solubilidade, mineral e aprisionamento residual. Trapping estratigráfica é onde CO 2 é realizada sob rochas vedação impermeável; solubilidade trapping é onde CO 2 se dissolve na água salgada residente em torno do CO 2 injectado 2-4; trapping mineral é onde fases minerais de carbonato são precipitados na rocha 5; e aprisionamento residual ou capilar é onde CO 2 é detido por forças de superfíciecomo pequenas gotas (gânglios) no pore-espaço da rocha 6. Isto pode ocorrer tanto naturalmente, pela migração do CO 2 pluma 7-9, ou pode ser induzida pela injecção de perseguição Brines 10. A fim de compreender os processos que regem o fluxo e armadilhas deste CO 2 no subsolo de um novo conjunto de experimentos devem ser realizados, aproveitando os novos avanços em tecnologia para melhor compreender a física fundamental associado com o fluxo de multi-fase.
Microtomografia de raios-X foi desenvolvido como uma técnica ao longo dos últimos 25 anos das primeiras tentativas de visualizar ambas as amostras geológicas secas 11 e múltiplas fases fluidas 12 para o principal método para a imagem não invasivo de núcleos de rock, tanto para fins de modelagem e para experimental implementação 13-15. Porque microtomografia é não-invasiva, que tem a capacidade para estudar os sistemas em condições representativas, que é particularmente attractivde e para o sistema de CO 2 -brine-rocha, como o comportamento do fluxo multifásico de SCCO 2 é altamente dependente das propriedades térmicas e físicas, tais como a tensão interfacial e ângulo de contacto, que por sua vez são uma forte função das condições do sistema, tais como temperatura, pressão e salinidade 16-18. Em um sistema tão complexo, com um tal conjunto extenso e pouco compreendido de variáveis inter-dependentes, experiências usando estruturas de poros idealizadas 19 ou fluidos analógicos 20,21 pode não ser aplicável a processos de fluxo de subsuperfície. Imagiologia múltiplos fluidos em condições representativas de uma formação injecção prospectivo CO 2 tem, no entanto, manteve-se um desafio 22. Neste estudo, descrevemos uma metodologia para a análise do comportamento multi-fluido em condições de reservatório, com foco no exame de capilar prendendo 23,24. Isto incluirá a concepção de uma estratégia de imagem, a montagem da célula de fluido, a injeção strategy e processamento de imagem subsequente.
O exame experimental do comportamento de fluxo multifásico pore escala em sistemas de rochas reais centra-se na imagem de núcleos de rocha parcialmente saturados após a injeção não molhamento fase (drenagem) e molhando injeção fase (embebição). Estes fluidos são injetados, ligando os núcleos para bombas de injeção de fluidos utilizando linhas de fluxo flexíveis, enquanto confinar o núcleo usando um design coreholder Hassler-type 25. Para o arranjo de imagem com êxito in-situ de SCCO 2 e salmoura, um processo novo e altamente sensível montagem experimental foi usada, com ênfase no uso de uma alta resolução de raios-x microscópio 23,24,26. Os requisitos para a realização de experimentos em temperaturas elevadas e pressões são muito rigorosos e exigem os recentes desenvolvimentos em tecnologia de materiais e micro-CT instalações. Os principais requisitos que têm de ser cumpridos são de que qualquer detentor core / amostra precisa ser able para suportar a alta pressão e alta temperatura (HPHT) condições, permanecendo suficientemente raios-x transparente para permitir a imagiologia efectiva. Os instrumentos baseados no Laboratório impor uma restrição adicional, como titular do núcleo deve ser pequeno o suficiente de modo que a fonte de raios-x pode ser colocado perto da amostra e que suficientemente grande geométrica ampliação de raio-x pode ser alcançado de modo a que o espaço poroso é eficazmente resolvidos. Embora esta restrição foi relaxado um pouco com a introdução da óptica secundários em máquinas mais recentes de micro-CT baseados em laboratório, não tem sido completamente removido, especialmente se os tempos rápidos de aquisição são desejados, como ampliações ópticas mais elevadas tendem a aumentar o tempo necessário para a aquisição imagens.
Experiências com fluidos solúveis proporcionar um desafio adicional quando se utiliza longos tempos de aquisição, na forma de CO 2 irá difundir através das porções poliméricas da montagem experimental, a redução da saturação do fluido in situ. All estas questões com que o tempo de verificação mais do que em torno de 2 horas eram impraticáveis. A fim de manter o tempo de verificação abaixo deste requisito, particularmente rigorosos para fontes de laboratório, o titular do núcleo deve ser de cerca de 1 cm de diâmetro. Um tamanho maior coreholder teria exigido o detector a ser muito mais longe da fonte para atingir o mesmo aumento geométrico, reduzindo o fluxo incidente de raios-x sobre o detector e, por conseguinte, aumentando os tempos de exposição de projecção necessária. A célula de fluxo utilizado nestas experiências foi baseada numa concepção da célula Hassler tradicional, construído em torno de uma manga de fibra de carbono, com um desenho de manga semelhante ao utilizado por Iglauer et al 27, mas com duas alterações significativas:. 1) O compósito de fibra de carbono utilizado no fabrico da luva foi alterada a partir de fibras T700, com uma rigidez de 230 GPa, ao fibras M55, com uma rigidez de 550 GPa. Isto não só reduziu a quantidade de movimento da amostra durante a aquisição de tomografia, mas também aumentou o máximo working pressão da célula de 20 MPa a 50 MPa. 2) A manga ter sido alongado a partir de 212 milímetros a 262 mm, para permitir que a fonte e o detector para estar tão próximo quanto possível da amostra.
Um grande deficiência experimental no primeiro estudo a utilizar micro-CT para examinar CO 2 em condições de reservatório era a utilização de linhas de metal para controlar o fluxo de e para o detentor do núcleo 27. À medida que a amostra é rodado em relação às bombas, as linhas de fluxo também precisa de ser rodado. Stiff linhas de fluxo pode fazer com que a amostra se mover, reduzir a resolução da imagem efectiva ou fazendo parte ou a totalidade do conjunto de dados inutilizável. Para evitar que isso substituímos todas as linhas de fluxo perto do palco rotação com flexível éter poliéter cetona (PEEK) tubulação. Estas linhas de fluxo eram flexíveis, proporcionando forças laterais muito pequenas (de carga) para o titular do núcleo durante a aquisição. Também ligado as linhas de fluxo de válvulas associadas à fase de amostra, em vez de ligar as linhas de fluxo para o coreholder. Isto significa que qualquer linha de carga de fluxo existente foi transmitida diretamente para a etapa, em vez de para a amostra, reduzindo a probabilidade de movimento da amostra. Uma grande desvantagem de utilizar o tubo de PEEK era que o CO 2 foi capaz de se difundir lentamente através dela, através de uma escala de tempo de cerca de 24 horas. Isto significa que o CO 2 salmoura saturada deixado nas linhas de fluxo se gradualmente dessaturar.
Outra grande deficiência experimental de estudos anteriores era o controle impreciso da temperatura. Isto pode ter impacto no resultados de um número de maneiras. Em primeiro lugar, a temperatura é um forte controle sobre ambos tensão e ângulo de contato interfacial 16-18. Além disso, a solubilidade de ambos SCCO 2 e carbonato de rocha em salmoura é também altamente dependente da temperatura 28. De controlo da solubilidade é crítico, como quando SCCO 2 é injectada num aquífero carbonato salina, se dissolverá na solução salina de residente, formando um ácido carbónico altamente reactivo, o qual, por turn começar a dissolver qualquer calcite presente. Qualquer imprecisão no controle solubilidade pode, portanto, levar a SCCO 2 dissolução / exsolution ou sólido de dissolução / precipitação.
Estudos anteriores 27 utilizado um fluido confinando aquecida para aquecer a coreholder; no entanto, este foi problemática. Tem as desvantagens associadas com a dificuldade de manter rigorosamente uma pressão confinante constante utiliza um sistema de recirculação de água, exigindo banhos de aquecimento extra para que a oferta. Além disso, este sistema só mantém um controlo preciso de temperatura no ponto de banho de aquecimento (não no ponto de suporte do núcleo, e confinar o fluido iria arrefecer entre o banho de água e o suporte de núcleo). Isto também exige uma entrada e uma porta de saída para o fluido de confinamento, aumentando o número de linhas de fluido associadas à coreholder e assim aumentando a carga de linha de fluxo.
Em vez de utilizar um fluido aquecido confinante, uma aq flexívelng revestimento foi usada para envolver o suporte do núcleo. Este método de aquecimento muito simples resultou em muito pouca carga coreholder, e permitiu o aquecimento preciso e exato. Utilizou-se um filme de poliimida aquecimento extremamente fino, a fim de minimizar o tamanho da amostra. A construção deste filme é composto por um elemento cobre gravado folha de 0,0127 milímetros de espessura, encapsulado entre duas camadas de filme de poliimida 0,0508 milímetro. Os elementos de cobre presentes no revestimento não afectar a qualidade de imagem perceptível. A temperatura foi medida utilizando um termopar sentado no anel confinando da célula. Foi colocado no exterior da manga de confinamento, tão perto quanto possível do núcleo, garantindo uma leitura precisa e fiável e estável a temperatura de poros de fluido. O filme termopares e aquecimento foram conectados a um controlador de costume construído Proporcional Integral Derivativo (PID), e as temperaturas foram controladas dentro de ± 1 ° C.
Para manter o controle completo over solubilidade inter-fase, e são condições presentes no aquífero longe do local da injecção, antes da injecção da salmoura foi equilibrada com 2 SCCO por vigorosamente a mistura dos dois fluidos em conjunto com partículas pequenas (1-2 mm) da rocha hospedeira num reactor agitado e aquecido. Todos os componentes de contacto com o fluido dentro deste reactor é feito de Hastelloy para minimizar a corrosão. O reactor contém um tubo de imersão para permitir o filtrado para o fluido mais denso para ser extraído a partir da base do reactor (salmoura) e fluido menos denso para ser extraído a partir da parte superior do reactor (SCCO 2). Bombas de seringa de alta pressão foram usadas para manter o fluxo de pressão e controle no pore-espaço do rock e no reator, com uma precisão de deslocamento de 25,4 nl. O aparelho experimental utilizado neste estudo é apresentado na Figura 1. O sal iónico utilizado para a experiência de que os resultados foram representativos desenhada foi iodeto de potássio (KI), como tem um elevado peso atómico e assim umelevado coeficiente de atenuação de raios-x, tornando-se um agente de contraste eficaz. Menos atenuantes (sais tais como NaCl) ou misturas podem ser utilizados, no entanto salinidades seriam necessárias para conseguir o mesmo atenuação de raios-x.
Protocol
1. Criação de Imagens Projeto Estratégia
- A fim de prever o desempenho de imagens de diferentes opções de soluto para a salmoura, calcular o espectro de raios-x do incidente raios-x 29-31. Incluir o impacto do porta-core, conjunto de núcleo e fluidos confinantes no espectro de raios-x. Um exemplo incidente espectro de raios-x usando uma voltagem de aceleração de 80 kV e corrente de electrões de 87 uA é mostrado na Figura 1.
- Comparar este espectro para os factores de transmissão da amostra contendo diferentes fluidos de poros. Simule mudanças no fator de transmissão, devido a mudanças no pore-fluido usando a lei de Beer-Lambert, assumindo um comprimento efectivo óptica das espécies dentro da amostra, e coeficientes de atenuação de raios-x calculados (Figura 3) 32. Determine o fator de transmissão global, integrando sobre todas as energias de raios-x de incidentes. Um exemplo dos factores de transmissão resultantes eficazes para a matri rochax e material de poros de espaço, e as alterações nos factores de transmissão em relação ao caso em que o espaço dos poros é preenchida com vácuo pode ser visto na Tabela 1.
- Escolha uma concentração de soluto e água salgada de modo a que a alteração no coeficiente de transmissão associado com salmoura é aproximadamente metade da alteração no coeficiente de transmissão associado com o sólido. Isto irá maximizar três contraste de fase na imagem reconstruída. Pesar a quantidade necessária de sal (7% (w / w) de KI foi utilizado para a experiência a partir do qual foram retirados os resultados representativos) e misturar vigorosamente com água desionizada.
- Em alternativa, se uma composição de solução salina específico for desejado, alterar a voltagem de aceleração fonte de raios-x para alterar o espectro dos raios X incidentes.
2. Montagem de Equipamentos e celular
- Montagem do equipamento conforme mostrado abaixo na Figura 2. Utilização PEEK linhas de fluxo para reduzir a carga lateral sobre a amostra da célula de fluxo.Teste cada conexão cuidadosamente se há algum vazamento de fluido.
- Coloque a salmoura, da composição durante os passos decidiu 1,1-1,3, na base do reactor. Envolver o aquecedor flexível em torno da célula de fluxo.
- Construa acessórios terminais metálicos. Remova a linha do "fim de um 1/16" 1/8 a 1/8 "redutor de montagem. Em seguida, corte pequenos sulcos no rosto do "fim do encaixe 1/8 para distribuir o CO 2 injectado ao longo de toda a face do núcleo.
- Passe o termopar alta pressão através das peças metálicas nas extremidades da célula de micro-fluxo e vedação com ¼ "virolas e porca, de modo que a junção quente do termopar fica adjacente à face de entrada do núcleo, no interior do anel de confinar a célula.
- Broca de amostra desejada em um núcleo 6,5 mm de diâmetro e 30 mm a 50 mm de comprimento. Triturar as extremidades do plano core, para garantir uma boa conexão com as peças finais metal. Enrole este núcleo em papel alumínio e coloque withind flúor-polímero de elastómero de manga.
- Ligue as extremidades da manga elastômero ao arranjo finais metal. Adicionar outro invólucro de folha de alumínio para o exterior da manga de elastómero, antes de colocar a junção quente do termopar ao lado do anel tubular de confinamento da célula e a adição de um envoltório final da folha de alumínio. Isto forma o conjunto de núcleo (Figura 4).
- Montagem da célula de micro-fluxo com o conjunto de núcleo selado dentro dele e ligar a célula para a fase dentro do compartimento de micro-CT (Figura 5), utilizando um grampo montado no topo da fase rotacional CT.
3. sistema de pressurização
- Fechar todas as válvulas para além da válvula 1, 2 e 3, tal como definido na Figura 3. Carga de CO 2 a partir do cilindro na bomba 1 e o reactor, em seguida, fechar a válvula 1. elevar lentamente a temperatura e pressão no interior do reactor para que desejado para a pore fluido durante o experimento.
- Abra a válvula 14 e carregar o fluido confinando em bomba 3. Feche a válvula 14. Abra as válvulas 12 e 13. pressurizar o anel confinando da célula para pelo menos 10% maior do que a pressão de poros fluido proposto.
- Abra a válvula 11. Carga salmoura para a bomba 2. Feche a válvula 11 e válvulas abertas 9, 8 e 6.
- Lentamente pressurizar o espaço dos poros da rocha até que a pressão é a poro de líquido desejado, enchendo o espaço de poros da amostra com uma solução salina que não tenha sido equilibrada com SCCO 2.
- Abrir a válvula 4. Lavar mais de 1000 volumes de poro da solução salina equilibrada através do núcleo de reenchimento da bomba 2 a uma taxa de fluxo constante. O volume de poro é encontrada multiplicando o volume do núcleo da porosidade encontrado usando porosimetria de hélio.
NOTA: Este vai miscibly deslocar a salmoura equilibrada-un, garantindo 100% de saturação inicial salmoura e creacondições ting no núcleo semelhante às condições do subsolo em um aquífero num ponto ligeiramente à frente da frente de uma pluma SCCO 2.
4. Fluid Flow e Aquisição de Imagem
- Passar com 10 volumes de poro (cerca de 1 ml) de SCCO 2 através do núcleo a velocidades de fluxo muito baixas (1,67 x 10 -9 m 3 / seg), assegurando um baixo número de capilares de cerca de 10 -6. Continuamente tomar projeções 2D, a fim de medir com precisão o volume total injetado, observando o ponto em que SCCO 2 desloca salmoura no espaço poroso.
- Passar com 10 volumes de poro (cerca de 1 ml) de solução salina equilibrada através do núcleo à mesma taxa de fluxo baixo, causando SCCO 2 para ficar presa como uma fase residual no espaço dos poros.
- Depois dos passos 4.1 ou 4.2, tomar scans da amostra a drenagem imagem ou embebição respectivamente. Usar um tamanho de voxel de tal modo que todo o diâmetro do núcleo se encaixa dentro do campo de view.
- Reconstruir as digitalizações usando um programa de reconstrução tomográfica. Para analisar todo o comprimento do núcleo enquanto reter uma pequena tamanho do voxel, reconstruir volumes compostas por unir várias seções sobrepostas, adquiridos sequencialmente.
NOTA: Cada seção necessária cerca de 400 projeções, tendo 15-20 min para adquirir, de modo a digitalização de um volume composto inteiro levou cerca de 90 min.
5. Processamento de Imagem e Segmentação
- Aplicar uma borda de meios não-locais preservando filtro 33,34 para o conjunto de dados e corrigir as imagens para qualquer endurecimento do feixe ou amolecimento artefatos criados durante a reconstrução da imagem modelando esses artefatos funções de Gauss como radialmente simétricas 35.
- Segmento de dados (transformar a informação de escala de cinzentos em uma representação binária do CO2 no interior da imagem) pela utilização de um algoritmo de bacia com uma semente gerado utilizando um histograma 2D 36, tratando o CO
- Analise esta imagem segmentada para encontrar tanto o número total de CO 2 voxels e também os tamanhos de cada grupo ligado de CO 2 residual.
Representative Results
Os resultados para um único carbonato, Ketton calcário, um oólito do membro superior Lincolnshire pedra calcária foram analisados em 3D, a fim de identificar e medir o volume de cada gânglio único desconectado, que foi então marcado (Figura 6). Todo o processamento foi realizado no âmbito dos programas Avizo incêndio 8.0 e ImageJ 37.
As imagens parcialmente saturados segmentados foram analisados pela contagem do número de voxels de residualmente preso SCCO 2 para encontrar a percentagem do volume ocupado pela rocha preso SCCO 2 - a capacidade de aprisionamento capilar. Este pode, então, ser convertido a uma saturação residual (S r), dividindo este valor pela porosidade obtida utilizando hélio como porosimetria. Uma parte significativa das SCCO 2 foram presos como uma saturação residual, com uma saturação residual de 0,203 ± 0,013. Isso está de acordo com os resultados encontrados em estudos anteriores usando micro-CT 23. Estudos maiores núcleo escala de aprisionamento residual neste tipo de rocha mostrou uma saturação residual inferior de 0,137 ± 0,012 38.
A entrada de salmoura em um núcleo saturado SCCO 2 é um processo em que a embebição de um fluido de molhagem (salmoura) invade cada poro, deslocando fluido não molhante (SCCO 2). Em uma rocha fortemente wet-água esperamos que a água para encher as áreas do espaço poroso por ordem de tamanho 39,40, prendendo gânglios desconectado no processo chamado de snap-off. Este processo deve ser percolação como 41 também é possível fazer previsões sobre a distribuição de tamanho dos clusters isolados. O número n de aglomerados de volume de s (medido em voxels) deve ser dimensionada como, onde τ é o expoente Fisher 42. Modelagem de rede mostrou que em redes regulares cúbicos tridimensionais o valor deste expoente é de cerca de τ = 2.189 43. Uma maneira natural de extrair este expoente de dados reaisé traçar a quantidade binned, tal como definido por Dias 41 e Wilkinson.
que deve escalar como:
Este é então plotados em um gráfico log-log em função do s (Figura 7), mostrando comportamento em lei de potência para grandes gânglios, mas uma sub-representação de menor gânglios em comparação com o modelo da potência. O expoente foi calculado por exclusão gânglios menor do que 10 5 voxels (aproximadamente o início do comportamento em lei de potência) e a realização de Levenberg-Marquardt 44,45 regressão utilizando um algoritmo robusto residual mínimo absoluto encaixe 46,47. Este foi realizada utilizando um pacote de software comercial. O expoente Fisher para este sistema foi 2,287 ± 0,009, close para o valor teórico de 2,189, indicando que neste sistema de embebição é de facto como percolação. De modo mais geral estes resultados confirmam as conclusões em maiores experimentos core-inundação 38,48,49 que SCCO 2 actua como a fase não molhar em carbonatos.
Figura 1. aparelho Experimental, mostrando as bombas, válvulas e reator utilizado para controlar o fluxo e o assentamento do coreholder dentro do gabinete micro-CT. Por favor, clique aqui para ver uma versão maior desta figura.
Figura 2. O espectro de energia normalizado para raios-x incidente sobre o núcleo, fi filtradas através do coreholder, confinando manga e confinando fluido. Calculado utilizando SpekCALC 29-31.
Figura 3. Os coeficientes de atenuação linear de diferentes fluidos e materiais rochosos como uma função de energia de fotões.
Figura 4. Detalhe do conjunto de núcleo, mostrando um envoltório de alumínio triple em torno do núcleo, impedindo diffusive CO 2 troca em toda a manga flouro-elastômero. Por favor, clique aqui para ver uma versão maior desta figura.
0 / 52440fig5.jpg "/>
Figura 5. Detalhe da célula de fluxo, aparelhos de aquecimento e implantação do conjunto de núcleo dentro da célula de fluxo. O termopar deve ser colocado o mais próximo possível da face de entrada do núcleo. Por favor, clique aqui para ver uma versão maior desta figura .
Figura 6. Imagem do carbonato após a drenagem e embebição. (A) A renderização 3D do núcleo após a drenagem, onde cada cluster não molhar fase é dada uma cor diferente. (BF) A renderização em 3D do núcleo após cinco experimentos de embebição, colorido, como descrito por (A). A grande variedade de cores, indicando uma fase residual mal ligado. (G) Uma secção transversal do núcleo após a drenagem. A fase mais escura é a SCCO 2, a fase intermediária é salmoura e a fase mais leve é grão rock. A seção transversal do núcleo após a embebição (H). Por favor, clique aqui para ver uma versão maior desta figura.
Figura 7. A distribuição do tamanho dos gânglios residual mostrado na Figura 6.
| A porosidade Material de enchimento | Vácuo | CO2 | H2O | H2O - 7% em peso de NaCl | H2O - 7% em peso de KI | Sólido (CaCO 3) |
| Fator de Transmissão | 0,25 | 0,247 | 0,243 | 0,242 | 0,224 | 0,202 |
| Mudança na transmissão fator relativo ao vácuo | N / D | -0.003 | -0,007 | -0,008 | -0,026 | -0,048 |
Tabela 1. Resumo dos resultados (factor de transmissão e as alterações no factor de transmissão em relação ao caso de vácuo cheio) de simulação das propriedades ópticas de raios-x do material de rocha e de poros espaço enchimento trabalhada durante este estudo. Cada coluna representa um material diferente enchendo o poro-espaço da rocha dentro da coreholder.
Discussion
Os passos mais críticos para uma imagem com qualidade de fluidos multifásicos a pressões e temperaturas elevadas são: 1) O isolamento bem sucedido do fluido dos poros de confinar o fluido circundante; 2) o equilíbrio eficaz dos fluidos e rocha antes da injecção; 3) Controle de temperatura efetiva durante todo o experimento; e 4) a segmentação eficaz das imagens resultantes.
A utilização dos envoltórios de alumínio é crítica para o isolamento bem sucedido do poro de líquido a partir do fluido envolvente, confinando na sua troca difusiva ausência do outro lado da manga é rápido, e saturação dentro do núcleo não permanece constante durante a duração do varrimento. Este problema também pode ser evidente quando o líquido permanece nas linhas de fluxo PEEK por períodos de tempo prolongados (> 2 h) antes da injecção no núcleo do passo 4.1 e 4.2. Mais uma vez, CO 2 Mobilidade difusivamente em todo o plástico, fazendo com que a salmoura para remover a saturação. Se este desalmoura saturada é injectado para dentro do núcleo, a saturação do núcleo irá diminuir como aglomerados residuais são dissolvidos, a solução salina injectada.
Outros métodos para o equilíbrio de fluidos e rochas, incluindo fluido de recirculação 50, foram propostos na literatura. Esses métodos aumentam a complexidade da configuração experimental, que por sua vez teria aumentado a quantidade de tempo para cada experiência, o que tem por sua vez aumentou a probabilidade de que a salmoura em linhas de fluxo teria difusivamente dessaturado.
A eficácia do controlo de temperatura é essencial, e a presença de um termopar no interior do anel de confinamento da célula de fluxo é crítico para esta. Temperatura apenas é medida num ponto único, o que significa que pode haver alguns gradiente através da amostra, levando a um desequilíbrio de solubilidade e dissolução e exsolution. Isto pode ser minimizado pela localização da junção quente do termopar tão perto quanto possível a tele inlet face do núcleo rock.
A segmentação eficaz das imagens resultantes podem ser um verdadeiro desafio com estes sistemas, como a segmentação de imagens contendo uma saturação parcial de vários fluidos é significativamente mais desafiador que a segmentação de imagens secos, por isso o uso de simples limiar universal escala de cinza é insuficiente 51. O uso de segmentação divisor de águas não só dá os resultados mais fiáveis, em comparação com outros algoritmos da literatura, mas também é o mais eficaz para lidar com campainha e volume parcial artefatos 35.
Uma das limitações mais significativas desta técnica é que ele só pode acessar o espaço macro-poros de uma rocha. A microporosidade (em escalas menores do que a resolução da imagem) permanece inacessível, e pode ser importante para o escoamento multifásico. Resoluções mais elevadas revelam uma maior proporção destas partes do espaço dos poros, mas também correspondem a uma diminuição no the campo de visão. A aplicabilidade da técnica a um tipo de rocha específica poderiam ser resolvidos através da comparação da resolução da digitalização para a distribuição de tamanho dos poros garganta obtido usando o método independente, como pressão capilar injeção de mercúrio.
Este método é uma técnica de imagiologia que conduz para a poro escala de vários fluidos em condições de reservatório em sistemas realistas, com aplicações existentes, incluindo um estudo comparativo transversal de aprisionamento capilar 24 e a medição do ângulo de contacto 26, e o método é facilmente aplicável a uma vasta gama de sistemas porosos. O trabalho futuro poderia estudar, na escala dos poros, uma grande variedade de fase única e fluxo multifásico em problemas de meios porosos em condições representativas dos aquíferos do subsolo, os campos de petróleo e gás e outros sistemas geológicas profundas.
Acknowledgments
Agradecemos o financiamento do Centro de Investigação carbonatos e Armazenamento de Carbono Qatar (QCCSRC), desde conjuntamente pela Qatar Petroleum, Shell e Qatar Science & Technology Park. Reconhecemos também o financiamento do Colégio Imperial de Consortium on Pore-Scale Modelling.
Materials
| Name | Company | Catalog Number | Comments |
| High Pressure Syringe Pump | Teledyne ISCO | 1000D | |
| Parr Reactor | Parr Instrument Company | 4547A - hastelloy | |
| PEEK Tubing | Kinesis | 1560xL | |
| Potassium Iodide Salt | Sigma Aldrich | 30315-1KG | |
| Carbon Dioxide | BOC | CO2 - size E | |
| Thermocouple | Omega Engineering | KMTSS-IM300U-150 | |
| Kapton Flexible Heater | Omega Engineering | KH-112/10-P | |
| X-Ray Microscope | Zeiss | Versa XRM 500 | |
| Snoop Leak Detector | Swagelok | MS-SNOOP-8OZ | |
| Flouro-Elastomer Polymer (Viton) Sleeve | Fisher Scientific | 11572583 | |
| Micro-CT Coreholder | Airborne Composites | 262mm Coreholder | Constructed in conjunction with Imperial College |
| Tomographic program | Zeiss | XM-Reconstructor | |
| ImageJ - image processing | NIH | ImageJ | |
| Matlab | Mathworks | Matlab | Used for regression analysis |
| Avizo | FEI | Avizo |
References
- Special Report on Carbon Dioxide Capture and Storage. , Cambridge, U.K. (2005).
- Ennis-King, J., Paterson, L. Engineering aspects of geological sequestration of carbon dioxide. SPE 77809, Proceedings of the Asia Pacific Oil and Gas Conference and Exhibition. , (2002).
- Weir, G. J., White, S. P., Kissling, W. M. Reservoir storage and containment of greenhouse gases. Transport in Porous Media. 23 (1), 61-82 (1996).
- Lindeberg, E., Wessel-Berg, D. Vertical convection in an aquifer column under a gas cap of CO2. Energy Conversion & Management. 38, 229-234 (1997).
- Lin, H., Fujii, T., Takisawa, R., Takahashi, T., Hashida, T. Experimental evaluation of interactions in supercritical CO2/water/rock minerals system under geologic CO2 sequestration conditions. Journal of Materials Science. 43 (7), 2307-2315 (2007).
- Juanes, R., Spiteri, E. J., Orr, F. M. Jr, Blunt, M. J. Impact of relative permeability hysteresis on geological CO2 storage. Water Resources Research. 42 (12), (2006).
- Hesse, M., Orr, F. M. Jr, Tchelepi, H. A. Gravity Currents with Residual trapping. Journal of Fluid Mechanics. 661, 35-60 (2008).
- MacMinn, C., Szulczewski, M. L., Juanes, R. CO2 migration in saline aquifers. Part 1. Capillary trapping under slope and groundwater flow. Journal of Fluid Mechanics. 662, 329-351 (2010).
- MacMinn, C., Szulczewski, M. L., Juanes, R. CO2 migration in saline aquifers. Part 2. Capillary and solubility trapping. Journal of Fluid Mechanics. 688, 321-351 (2011).
- Qi, R., LaForce, T. C., Blunt, M. J. Design of carbon dioxide storage in aquifers. International Journal of Greenhouse Gas Control. 3 (2), 195-205 (2009).
- Flannery, B. P., Deckman, H. W., Roberge, W. G., D'Amico, K. L. Three-Dimensional X-Ray Microtomography. Science. 237 (4821), 1439-1444 (1987).
- Jatsi, J. K., Jesion, G., Feldkamp, L. Microscopic Imaging of Porous Media with X-Ray Computer Tomography. SPE 20495. SPE Formation Evaluation. , 189-193 (1993).
- Blunt, M. J., et al. Pore-scale imaging and modelling. Advances in Water Resources. 51, 197-216 (2013).
- Cnudde, V., Boone, M. N. High-resolution X-ray computed tomography in geosciences: A review of the current technology and applications. Earth-Science Reviews. 123, 1-17 (2013).
- Wildenschild, D., Sheppard, A. P. X-ray imaging and analysis techniques for quantifying pore-scale structure and processes in subsurface porous medium systems. Advances in Water Resources. 51, 217-246 (2013).
- Espinoza, D. N., Santamarina, J. C. Water-CO2-mineral systems: Interfacial tension, contact angle, and diffusion-Implications to CO2 geological storage. Water Resources Research. 46 (7), (2010).
- Li, X., Boek, E., Maitland, G., Trusler, J. P. M. Interfacial Tension of (Brines + CO2): (0.864 NaCl+0.136 KCl) at Temperatures between (298 and 448) K, Pressures between (2 and 50) MPa, and Total Molalities of (1 to 5) mol/kg. Journal of Chemical and Engineering Data. 57 (4), 1078-1088 (2012).
- Spiteri, E. J., Juanes, R., Blunt, M. J., Orr, F. M. Jr A New Model of Trapping and Relative Permeability Hysteresis for All Wettability Characteristics. SPE Journal. 13 (3), 277-288 (2008).
- Chaudhary, K., et al. Pore-scale trapping of supercritical CO2 and the role of grain wettability and shape. Geophysical Research Letters. 40, 1-5 (2013).
- Karpyn, Z. T., Piri, M., Singh, G. Experimental investigation of trapped oil clusters in a water-wet bead pack using X-ray microtomography. Water Resources Research. 46 (4), (2010).
- Wildenschild, D., Armstrong, R. T., Herring, A. L., Young, I. M., Carey, J. W. Exploring capillary trapping efficiency as a function of interfacial tension, viscosity, and flow rate. Energy Procedia. 4, 4945-4952 (2011).
- Silin, D., Tomutsa, L., Benson, S. M., Patzek, T. W. Microtomography and Pore-Scale Modeling of Two-Phase Fluid Distribution. Transport in Porous Media. 86, 495-515 (2011).
- Andrew, M. G., Bijeljic, B., Blunt, M. J. Pore-scale imaging of geological carbon dioxide storage under in situ conditions. Geophysical Research Letters. 40 (15), 3915-3918 (2013).
- Andrew, M. G., Bijeljic, B., Blunt, M. J. Pore-scale imaging of trapped supercritical carbon dioxide in sandstones and carbonates. International Journal of Greenhouse Gas Control. 22, 1-14 (2014).
- Method and Apparatus for Permeability Measurements. US Patent. Hassler, G. L. , 2,345,935 (1944).
- Andrew, M. G., Bijeljic, B., Blunt, M. J. Pore-scale contact angle measurements at reservoir conditions using X-Ray microtomography. Advances in Water Resources. 68, 24-31 (2014).
- Iglauer, S., Paluszny, A., Pentland, C. H., Blunt, M. J. Residual CO2 imaged with x-ray micro-tomography. Geophysical Research Letters. 38, (2011).
- Carroll, J. J., Slupsky, J. D., Mather, A. E. The Solubility of Carbon Dioxide in Water at Low Pressure. Journal of Physical and Chemical Reference Data. 20 (6), (1991).
- Poludniowski, G. Calculation of x-ray spectra emerging from an x-ray tube. Part II. X-ray production and filtration in x-ray targets. Medical Physics. 34 (2175), (2007).
- Poludniowski, G., Evans, P. M. Calculation of x-ray spectra emerging from an x-ray tube. Part I. Electron penetration characteristics in x-ray targets. Medical Physics. 34 (2161), (2007).
- Poludniowski, G., Landry, G., DeBlois, F., Evans, P. M., Verhaegen, F. SpekCalc: a program to calculate photon spectra from tungsten anode x-ray tubes. Physics in Medicine and Biology. 57 (19), (2009).
- Berger, M. J., et al. XCOM: Photon Cross Section Database (version 1.5). National Institute of Standards and Technology. , Gaithersburg, MD. (2007).
- Buades, A., Coll, B., Morel, J. -M. A non-local algorithm for image denoising. Computer Vision and Pattern Recognition. , (2005).
- Buades, A., Coll, B., Morel, J. -M. Nonlocal Image and Movie Denoising. International Journal of Computer Vision. 76 (2), 123-139 (2008).
- Schluter, S., Sheppard, A. P., Brown, K., Wildenschild, D. Image processing of multiphase images obtained via X-ray microtomography: A review. Water Resources Research. 50 (4), 3615-3639 (2014).
- Jones, A. C., et al. Assessment of bone ingrowth into porous biomaterials using MICRO-CT. Imaging Techniques for Biomaterials Characterization. 28 (15), 2491-2504 (2007).
- Schneider, C. A., Rasband, W. S., Eliceiri, K. W. NIH Image to ImageJ: 25 years of image analysis. Nature Methods. 9, 671-675 (2012).
- El-Maghraby, R. Measurements of CO2 trapping in Carbonate and Sandstone Rocks. , Imperial College, London. (2013).
- Roof, J. G. Snap-Off of Oil Droplets in Water-Wet Pores. SPE Journal. 10 (1), 85-90 (1970).
- Lenormand, R., Zarcone, C., Sarr, A. Mechanisms of the displacement of one fluid by another in a network of capillary ducts. Journal of Fluid Mechanics. 135, 337-353 (1983).
- Dias, M. M., Wilkinson, D. Percolation with trapping. Journal of Physics A: Mathematical and General. 19, 3131-3146 (1986).
- Fisher, M. E. The theory of equilibrium critical phenomena. Reports on Progress in Physics. 30 (2), 615 (1967).
- Lorenz, C. D., Ziff, R. M. Precise determination of the bond percolation thresholds and finite-size scaling corrections for the sc, fcc and bcc lattices. Physical Review E. 57 (1), 230-236 (1998).
- Levenberg, K. A method for the solution of certain non-linear problems in least squares. Quarterly Journal of Applied Mathmatics. 2, 164-168 (1944).
- Marquardt, D. W. An Algorithm for Least-Squares Estimation of Nonlinear Parameters. Journal of the Society for Industrial and Applied Mathematics. 11 (2), 431-441 (1963).
- Holland, P. W., Welsch, R. E. Robust regression using iteratively reweighted least-squares. Communications in Statistics - Theory and Methods. 6 (9), 813-827 (1977).
- Huber, P. J. Robust Statistics. , John Wiley & Sons, Inc. New York, NY. (1981).
- Akbarabadi, M., Piri, M. Relative permeability hysteresis and capillary trapping characteristics of supercritical CO2/brine systems: An experimental study at reservoir conditions. Advances in Water Resources. 52, 190-206 (2013).
- Pentland, C. H., El-Maghraby, R., Iglauer, S., Blunt, M. J. Measurements of the capillary trapping of super-critical carbon dioxide in Berea sandstone. Geophysical Research Letters. 38, 4 (2011).
- El-Maghraby, R. M., Blunt, M. J. Residual CO2 Trapping in Indiana Limestone. Environmental Science and Technology. 47 (1), 227-233 (2013).
- Sheppard, A. P., Sok, R. M., Averdunk, H. Techniques for image enhancement and segmentation of tomographic images of porous materials. Physica A. 339 (1-2), 145-151 (2004).
