Infiltração, profundidade da frente usando o lapso de tempo offset multi reúne de rastreamento coletados com matriz chão de antena Radar de penetração

Summary

Aqui nós apresentamos um sistema de Radar penetrante solo (GPR) com base em uma matriz de antena chão flexíveis, densamente povoadas para monitorar o processo dinâmico de infiltração de água subterrânea. Uma imagem de radar de lapso de tempo do processo de infiltração permitiu estimar a profundidade da frente de molhamento no decurso do processo de infiltração.

Abstract

Um sistema de Radar penetrante solo (GPR) com base em uma matriz de antena chão flexíveis, densamente povoadas foi usado para coletar dados durante um experimento de infiltração, conduzido em um site de teste perto de dunas de areia de Tottori, Japão. A antena utilizada neste estudo consiste de 10 antenas de transmissão (Tx) e 11 receber antenas (Rx). Para este experimento, o sistema foi configurado para usar todos os possíveis pares de Tx-Rx, resultando em um deslocamento multi reunir (MOG) constituídos por 110 Tx-Rx combinações. A matriz foi deixada parada em uma posição diretamente acima da área de infiltração e os dados foram coletados a cada 1,5 segundos usando um disparador de baseados em tempo. Cubos de dados comum-deslocamento reunir (COG) e ponto médio comum (CMP) foram reconstruídos a partir dos dados MOG durante o pós-processamento. Há poucos estudos que usaram dados CMP lapso de tempo para estimar mudanças na velocidade de propagação. Neste estudo, a velocidade da onda eletromagnética (EM) estimou heuristicamente em intervalos de 1 minuto, a partir dos dados CMP reconstruídos através da curva de encaixe, usando a equação da hipérbole. Em seguida procedeu-se calcular a profundidade da frente de molhamento. A evolução da molhadela frente ao longo do tempo obter por este método é consistente com as observações de um sensor de umidade do solo que foi colocado em uma profundidade abaixo de 20 cm. Os resultados obtidos neste estudo demonstram a capacidade de tal matriz sistema GPR para monitorar um processo dinâmico subsuperficial como infiltração de água com precisão e quantitativamente.

Introduction

Processos de transporte em massa e energia compreensão na zona do vadose são importantes para muitas aplicações em disciplinas agrícolas e ambientais. Entre esses processos, fluxo de água variável saturado é o processo essencial, como muitos dos outros processos, tais como processos físicos, geoquímicos, biológicos e até mesmo mecânicos, geralmente são acoplados com fluxo de água. Desenvolvimentos recentes em técnicas de geofísicas permitiram um para monitorar processos hidrológicos na zona do vadose não invasiva. Entre muitas técnicas geofísicas, radar de Prospecção Geotécnica (GPR) é uma das técnicas mais amplamente utilizadas para monitorar e caracterizar a dinâmica de água do solo, porque é a propagação de ondas eletromagnéticas de (EM) emitido e recebido por antenas GPR caracterizada por solo umidade conteúdo^1,2,3,4. Entre os sistemas disponíveis, superfície no chão GPR (referida como superfície GPR no restante do manuscrito) é o mais comum para usar em um campo. Sistemas tradicionais de superfície GPR com um transmissor e um receptor (sistemas de radar bistatico) são comumente usados para digitalizar o subsolo com uma separação de emissão e recepção constante (deslocamento). Conjuntos de dados coletados nesta configuração são também conhecidos como deslocamento comum reúne (COG). Dados de radar são exibidos como séries temporais, com base no tempo de viagem total entre o transmissor, eventuais refletores e volta para o receptor. Para converter o tempo de viagem com informações de profundidade, velocidade da onda na subsuperfície precisa de ser estimado. Por exemplo, isso pode ser feito através da análise de conjuntos de dados de multi deslocamento se reúnem (MOG)⁵.

Embora tenha havido um número de estudos usando GPR para monitorar infiltração subsuperficial processos^6,7,8,9, nenhum deles diretamente determinado o local da frente de molhamento ou a onda EM estrutura de velocidade que está mudando com o tempo durante a infiltração. A abordagem comum é usar objetos enterrados em profundidades conhecidas como refletores de referência para determinar a velocidade média da onda EM e profundidade da frente de umedecimento. Desde a frente de molhamento altera dinamicamente durante a infiltração, MOG lapso de tempo deve ser coletado em intervalos de tempo curto para explorar as mudanças na estrutura de velocidade de onda EM sem o uso de objetos de referência. Com comum bistatico antenas GPR superficiais, a coleção de MOG lapso de tempo em intervalos curtos entre si é difícil ou impossível como ele necessita para mover manualmente as antenas para configurar as várias configurações de deslocamento. Recentemente, uma família de antenas GPR (referido como matriz GPR, daqui em diante) tem sido amplamente utilizada para o subsolo rápido e com precisão¹⁰de imagem. O conceito básico da matriz de GPR é fornecer áreas densas com um esforço mínimo comutando eletronicamente várias antenas montadas em um único quadro. Sistemas GPR de matriz tem sido usados principalmente para gerar imagens 3D de subsuperfície de amplas áreas rapidamente. Alguns exemplos de aplicações típicas para estes sistemas são estrada e ponte de inspeção¹¹, prospecção arqueológica¹² e engenhos explosivos e minas terrestres deteção^13,14. Para tais fins, a matriz GPR é usada principalmente para fazer a varredura do subsolo com uma configuração de separação constante antena para coletar COG. Embora tem sido demonstrado que MOG coletados com uma matriz que GPR poderia ser utilizado para estimativa de velocidade¹⁵, a aplicação prática desta metodologia tem sido limitada a poucos casos. Ao colocar a antena num local fixo, tempo-lapsed MOG pode ser facilmente coletado. Como demonstrado em nossa recente publicação¹⁶, lapso de tempo radargrams coletados com a matriz do sistema GPR preferia claramente retratar os sinais de reflexão da parte dianteira de umectação enquanto se move gradualmente para baixo durante um experimento de infiltração vertical realizadas em uma duna de areia. O principal objetivo do livro foi demonstrar como usar a matriz GPR para coletar MOG lapso de tempo durante o teste de infiltração e como analisar esses dados para controlar a profundidade da frente de molhamento.

Neste estudo, utilizamos uma matriz de antena, consistindo de 10 a transmitir (Tx0 - Tx9) e 11 (Rx0 - Rx10) a receber antenas monopolo de gravata-borboleta. O deslocamento dos elementos dentro da matriz de antena é mostrado na Figura 1 (ver Tabela de materiais). A antena é controlada por uma unidade de radar de onda contínua passo-frequência (SFCW) operando na frequência de 100 MHz a 3.000 MHz. Os interruptores GPR matriz através de uma sequência definida pelo usuário de pares de Rx-Tx usando multiplexadores de radiofrequência (RF) da matriz de antena¹⁰. O número máximo de combinações de Tx-Rx para este sistema particular é 110. Para este experimento, configuramos a matriz GPR usar 110 todas as combinações, a sequência de verificação de programação, para que cada transmissor, de Tx0 a Tx9, foi pareado sequencialmente todos os 11 receptores de Rx0 para Rx10. O tempo necessário para executar uma varredura através de todas as combinações de 110 é menos de 1,5 segundos. O deslocamento entre o emissor e o receptor foi calculado com base na distância entre os pontos de entrada dos elementos da antena, onde o deslocamento vertical é 85 milímetros, conforme mostrado na Figura 1.

Protocol

O principal objetivo deste estudo foi usar uma matriz sistema GPR para visualizar processos dinâmicos subsuperficiais, tais como a infiltração de água, também é importante descrever como o teste de infiltração foi conduzido.

1. teste de infiltração de campo

1. Selecione um site com uma superfície nua para o ensaio de infiltração.
   Nota: Para este estudo, o teste de infiltração foi conduzido em uma facilidade grande estufa da Universidade árida terra Research Center de Tottori, Japão. Consulte a publicação recente para mais detalhes,¹⁶.
2. Posição sobre o paralelo de tubos muito porosa superfície seis 2,5-m a uma distância de 15 cm entre si, a fim de garantir a infiltração uniforme (Figura 2).
3. Conecte uma extremidade do tubo para uma entrada de válvula regulamentado, que é conectado a um reservatório de água ou de uma torneira, e a outra extremidade a uma tomada, a fim de remover o excesso de água.
4. Coloque uma fina 1802 x 910 mm² painel de madeira (um pouco maior que a antena), para cobrir os tubos porosos.
5. Verificar o efeito do painel de madeira do sinal de radar antes do ensaio de infiltração.
6. Instale sensores de umidade do solo ao lado do painel de madeira para monitorar alterações no teor de umidade do solo em diferentes profundidades.
   Nota: O objetivo destes sensores é fornecer dados de referência usados para verificar os resultados obtidos a partir da análise dos dados de GPR. Neste estudo, um sensor de umidade do solo de haste-tipo (ver Tabela de materiais) foi instalado ao lado da antena para monitorar a constante dielétrica de solo em 10, 20, 30, 40, 60 e 100 cm de profundidade.

2. GPR medição

1. Coloque a antena GPR de matriz no painel de madeira.
   Nota: Neste estudo, a antena ficou estacionária para garantir a consistência em sinais gravados.
2. Ligue a antena para o controlador usando cabos coaxiais.
3. Conecte o controlador para um laptop PC com um cabo Ethernet.
4. Decidir as combinações de antena e sua sequência para garantir ambos deslocamento comum se reúnem (COG) e reunir vários deslocamento (MOG) são recolhidas perfeitamente.
   Nota: Neste estudo, todas as combinações possíveis para a matriz GPR, à nossa disposição foi usado (110). Demorou menos de 1,5 s para executar uma varredura completa através de todas as combinações.
5. Inicie a gravação dos dados de GPR utilizando o software de coleta de dados.
   Nota: É importante iniciar a aquisição de dados antes de água é injectada os tubos porosos.
6. Abra a válvula para aplicar a água (ou seja, iniciar o processo de infiltração).
   Nota: A taxa de injeção de água no subsolo pode ser determinada a partir das propriedades hidráulicas (por exemplo, condutividade hidráulica) do solo alvo. A taxa de fluxo utilizada neste estudo foi 7.000 cm3/min, que foi menor que a condutividade hidráulica saturada da areia das dunas.
7. Depois de uma quantidade pré-determinada de água é injectada, feche a válvula.
   Nota: Água foi injectada continuamente por 4 h, e um total de 1.680 L de água foram usados para o teste.
8. Passar o scanner com a matriz GPR por um período adicional de tempo, a fim de monitorar o processo de redistribuição para posterior análise.
   Nota: A matriz que GPR tem sido mantido coletando dados para 4h após a válvula tinha sido fechada. No entanto, esses dados não foram analisados para o estudo aqui apresentado.

3. os dados análise para a estimativa da velocidade

1. Reconstrua dados COG e CMP extraindo as combinações de Tx-Rx relativas do cubo de dados global. Isto foi conseguido através de um código personalizado desenvolvido pelos autores que rearranjou os perfis de radar então para agrupar pares Tx-Rx no deslocamento idêntico (COG) e os que têm o mesmo ponto médio (CMP).
   Nota: Porque uma unidade SFCW GPR foi usada para este estudo, os dados brutos são no domínio da frequência; a operação de frequência para conversão de tempo foi executada usando o software fornecido pelo fabricante do sistema de GPR.
2. Estime a velocidade da onda EM ideal na zona de umectação em determinado momento, decorrido encaixando o tempo de viagem em dois sentidos computada, t,_c, da onda EM refletida na parte dianteira de infiltração para os sinais observados. A velocidade com a curva de melhor ajuste foi determinada heuristicamente.
   (1)
   onde d₀ é a profundidade do ponto de reflexão, x é a distância entre Tx e Rx, v_r é a velocidade de root-mean-square ao plano de reflexão e t₀ é o tempo de viagem em dois sentidos em zero deslocamentos determinada a partir do radargram da CMP.
   Nota: O mais comum método de estimativa da velocidade usando MOG é semelhança análise^5,15,17. Na análise de semelhança, a velocidade de raiz quadrada é escolhida usando um espectro de velocidade. Análise de semelhança não foi utilizado para a estimativa de velocidade neste estudo porque a baixa relação sinal-ruído observada na matriz de dados GPR não nos permitiu desenvolver um espectro de velocidade confiável. Encontrar um método de filtragem adequado para gerar um espectro de velocidade confiável a partir destes dados está sendo investigada em outras investigações em curso. Para o escopo deste estudo, a método de encaixe de curva foi usada para ajustar o Eq (1) para o radargrams da CMP.
3. Usando a EQ. (1), calcule a molhadela profundidade frontal d₀ em determinado momento, decorrido usando o estimado EM velocidade v_r.
   Nota: Para esta parte da análise dos dados, utilizou-se o conjunto de dados COG em um deslocamento x igual a 113 mm.

Results

A Figura 3 mostra um lapso de tempo painel diagrama de 110 combinações de Tx-Rx adquiriu cada 1,5 s durante o primeiro 60 min do experimento infiltração, depois de converter os dados de frequência para o domínio do tempo. Para reforçar a reflexões em profundidade, um filtro passa-banda foi aplicado seguido por compensação de ganho. O diagrama do painel pode ser dividido em 10 seções, cada seção correspondente a um determinado Tx. A localização do Tx é indicada por um triângulo branco, e cada fatia corresponde a um lapso de tempo sinal gravado para um Rx. O eixo vertical mostra o tempo de viagem em dois sentidos ou TWT, ou seja, o tempo necessário para a onda EM viajar a partir do transmissor, ao eventual refletor e volta para o receptor. Os perfis de exibem a amplitude do sinal em um mapa de cores em tons de cinza. Um grande contraste na cor indica alta amplitude do sinal de radar gravado. Reflexões de onda EM que são produzidos na interface entre camadas diferentes de constante dielétrica ou por objetos com diferentes características elétricas do que no meio circundante; um alto contraste dielétrico irá determinar uma reflexão de grande amplitude. Durante o processo de infiltração de água, há uma zona referida como a zona de transição onde o conteúdo de água aumenta gradualmente da borda da zona umectante, que é onde a água penetra o solo inicialmente seco. A onda EM é susceptível de ser refletido não no limite, mas dentro da zona de transição, como observado em estudos de deteção de água tabela¹⁸. No restante do manuscrito, esta área de reflexão é referida como a frente de molhamento. Na Figura 3, um sinal de alta amplitude aparece e move-se continuamente para baixo conforme o tempo avança durante o experimento. Esta reflexão é realmente produzida pela frente de molhamento como a água gradualmente para baixo penetra no subsolo. Este diagrama, COG e CMP podem ser reconstruído como mostrado na Figura 2 de Iwasaki et al ¹⁶

Análise de velocidade foi a CMP dados obtidos a cada 1 min. Para cada conjunto de dados da CMP, o tempo de viagem em dois sentidos, dado pela EQ. (1) foi montado para a reflexão da parte dianteira de umectação ajustando t₀ e v_r, assumindo uma camada uniforme na zona molhada. Zero hora foi corrigida por encaixe a onda de ar com a velocidade de 0,3 m/ns. A Figura 4 mostra dados CMP em intervalos de 5 min de tempo transcorrido t_e = 5 min para t_e = 50 min juntamente com as curvas de ajuste exibido como linhas brancas (linhas sólidas para a onda refletida) e a linha tracejada para a onda de ar. Desde que a onda de ar não é uma onda refletida, mas o sinal direto entre o emissor e o receptor, o tempo de viagem aumenta linearmente com o deslocamento. Todas as curvas foram equipadas para os picos positivos (na cor branca) das ondas refletidas. Todas as curvas bem montadas as curvas refletidas observadas mostradas no radargrams CMP em todos os deslocamentos, o que significa que os valores estimados para t₀ e v_r são boas. Para o experimento de infiltração, um painel de madeira seco foi colocado entre a antena e os tubos porosos. Porque o painel tem uma muito baixa constante dielétrica do que a de solo úmido, seus efeitos sobre a propagação de ondas EM podem não ser insignificantes, mesmo que seja fino. Um modelo de duas camadas considerou-se, então, juntamente com o modelo de camada uniforme acima mencionados, assumindo um valor de 3 para a constante dielétrica do top 5 cm. Também para este segundo modelo, a de velocidade de onda EM v_r foi estimada por encaixe para o reflexo produzido pela frente de molhamento de curva.

Em Figura 5, estimado umectante frente profundidades são plotadas em função de t_e para ambos os modelos simples e duplo. Este pode ser apreciado que a frente de molhamento se move para baixo quase linearmente com o tempo para ambos os modelos, com exceção de um abrandamento entre t_e = 10 min e t_e = 20 min. diferenças entre os dois modelos são inicialmente não significativo, mas como tempo decorrido a estimativa para o modelo uniforme se move um pouco mais rápido em comparação com o modelo de duas camadas. Na Figura 5, símbolos de diamante são usados para marcar os tempos, quando as leituras dos sensores da umidade começaram a aumentar, e quando eles mais tarde se tornou constantes; Estes são conectados com uma linha sólida para cada profundidade de sensor. Como mencionado acima, a reflexão da onda EM não necessariamente ocorre no limite da zona húmida; em outras palavras, considerando uma certa profundidade, não se pode esperar desta reflexão para ser correspondente ao ponto no tempo quando as leituras de um sensor começam a aumentar. Neste sentido, a reflexão é atribuível a um nível de profundidade na parte da frente de infiltração, onde atingiu uma certa saturação de água, em comparação com a área imediatamente abaixo. Considerando-se os sensores em 30, 40 e 60 cm de profundidade, a estimativa da profundidade da frente umectante obtida a partir dos dados GPR cai bem no intervalo mostrado pelas linhas sólidas sobre o cronograma. O tempo quando o GPR estimado umectantes frente chegaram a 20 cm de profundidade corresponde ao tempo quando foi observado o aumento súbito da leitura do sensor, enquanto a estimativa GPR alcançado a profundidade de 10 cm muito mais rápido do que o que foi produzido pelo sensor de umidade , embora o sinal de reflexão da parte dianteira de umectante é claramente observado após t_e = 5 min (Figura 4). Além disso, tem de ser mencionado que a extrapolação de GPR estimada curva não passa através da origem. Embora não seja claro o que causou esta discrepância em profundidades mais rasas, pode haver algumas explicações possíveis. Isso pode ser atribuído à heterogeneidade nas propriedades do solo, ou pode ser devido a não uniformidade na aplicação de água. Se for o caso, isto teria um efeito maior mais cedo durante o processo de infiltração do que numa fase posterior. Outra explicação seria que a rugosidade da superfície afeta a determinação do tempo zero. Além do efeito do painel de madeira e os tubos porosos, o efeito da rugosidade da superfície deve ter em conta.

Figura 1 : Esquema de matriz à configuração de antena de radar penetrante utilizada neste estudo terra. Estruturas em forma de V são antenas monopolo de gravata borboleta. Existem 10 antenas de transmissão (Tx) e 11 receber antenas (Rx) alinhadas horizontalmente. Clique aqui para ver uma versão maior desta figura.

Figura 2: esquemas do experimento infiltração. (A) uma vista superior e (B) uma vista lateral, no qual a antena matriz foi colocada em cima de seis tubos porosos de 250 cm alinhados 15 centímetros distante. Clique aqui para ver uma versão maior desta figura.

Figura 3 : Radargram de lapso de tempo obtidos durante o primeiro 60 min do experimento infiltração. Dados consistem em sinais gravados para 110 combinações de antena. Uma fatia corresponde ao lapso de tempo dados coletados com uma única combinação de Tx-Rx. Uma cor diferente é usada para a amplitude do sinal. Linhas verticais pretas de dados separada para cada triângulos Tx. White indicam os locais das Tx. clique aqui para ver uma versão maior desta figura.

Figura 4 : Radargrams da CMP no t_e = 5 min para t_e = 55 min em intervalos de 5 min. O linhas sólidas representam manualmente branco equipado tempo de viagem em dois sentidos da reflexão da frente umectante, enquanto as linhas tracejadas brancas representam o tempo de viagem da onda de ar. Clique aqui para ver uma versão maior desta figura.

Figura 5 : Umectante frente profundidades. Umectantes frente profundidades estimado a partir de matriz GPR em Time-Lapse offset multi se reúnem em função do tempo decorrido para o uniforme (triângulos) e os modelos de duas camadas (quadrados). Linhas pretas com diamantes em ambas as extremidades mostram o tempo entre o aumento inicial de leituras e quando estas atingiu um nível constante para cada profundidade de sensor (ou seja, a duração da zona de transição). Clique aqui para ver uma versão maior desta figura.

Discussion

Neste estudo, chão de matriz radar penetrante (GPR) foi usado para controlar a profundidade da frente de molhamento durante um experimento de infiltração realizada em um campo experimental perto de dunas de areia de Tottori, Japão. A matriz sistema GPR utilizado neste estudo é composto por 10 antenas de transmissão (Tx) e 11 receber antenas (Rx). O sistema pode ser configurado para usar até 110 diferentes combinações de Tx-Rx. Durante o experimento de infiltração, todas as combinações de 110 foram escaneadas continuamente em intervalos de 1,5 segundos, deixando a matriz estacionária no local onde a água foi aplicada através de alguns tubos porosos, colocados na superfície. Reunir comum-deslocamento (COG) e dados de ponto médio comuns (CMP) foram reconstruídos a partir o cubo de dados lapso de tempo. Seria praticamente impossível coletar dados CMP na mesma taxa com sistemas GPR de bistatico convencional. É muito importante deixar a antena na gravação estacionária durante o experimento para obter dados de lapso de tempo significativos e reprodutíveis.

Embora dados GPR matriz tem sido usados para estimar EM onda velocidades¹⁵, existem apenas alguns estudos que analisaram dados GPR matriz lapso de tempo para estimar velocidades de onda EM para processos transitórios como infiltração de água. Neste estudo, a estrutura de velocidade de onda eletromagnética do (EM) foi estimada a partir dos dados CMP lapso de tempo. Em vez de realizar análise de semelhança, a curva de hipérbole para o tempo de viagem em dois sentidos foi montada heuristicamente os sinais refletidos em radargrams a CMP para estimar a velocidade média das ondas na zona de umectação por causa do baixo sinal-ruído (S / N) relação nos dados. Quando a relação S/N é baixa, análise de semelhança não pode ser usado para gerar um espectro de velocidade confiável. Um método de filtragem adequado precisa ser desenvolvido para usar o método de análise de semelhança. Um sensor de umidade do solo de haste-tipo foi instalado ao lado da antena GPR para medir as mudanças no teor de umidade do solo durante o experimento de infiltração; os sensores foram deslocados nas profundidades de 10, 20, 30, 40 e 60 cm e trabalhou de forma independente.

Usando a velocidade da onda EM estimado, a profundidade da frente de molhamento foi computadorizada em intervalos de 1 minuto do processo de infiltração. A evolução no tempo da molhadela estimado frente concorda bem com as observações de sensores de umidade do solo em profundidades abaixo de 20 cm. Em profundidades mais rasas, a estimativa GPR da profundidade da frente umectante mostra uma discrepância com as leituras dos sensores de umidade do solo.

Em geral, este estudo demonstra que a matriz sistema GPR é capaz de rastrear a evolução da profundidade da frente umectante durante a infiltração de água no solo, através da recolha de lapso de tempo comuns dados do ponto médio (CMP). Como este tipo de dados não foi facilmente coletado de GPR superfície convencional antes, os dados obtidos neste estudo são os primeiros que realmente mostram como a frente de molhamento evoluiu ao longo do tempo no subsolo. Trabalho futuro irá explorar a possibilidade de usar a inversão de dados para estimar os parâmetros hidráulicos do solo a partir dos dados obtidos durante o experimento.

Materials

List of materials used in this article

Name	Company	Catalog Number	Comments
Unidade de Radar	GeoScope Radar 3D Antena AS
DXG1820	Radar 3D Sonda
PR2/6	Delta-T