Elastic wave propagation in fractured media with spring-type and frictional contact deformation laws

O solo sob nossos pés raramente é um bloco sólido e ininterrupto. Em muitos lugares, desde profundos reservatórios de petróleo até a própria crosta terrestre, a rocha é atravessada por fraturas — fissuras planas e finas que podem se estender por quilômetros. Quando a energia viaja através desse cenário fraturado, seja uma onda sísmica de um terremoto ou um pulso sonoro usado para inspecionar um material, as rachaduras mudam a história. Elas dispersam a energia, enfraquecem o sinal e, às vezes, até fazem com que as faces das rochas se esfreguem umas contra as outras ou deslizem para longe. Compreender exatamente como essas ondas interagem com as rachaduras é vital para geólogos que tentam mapear o que jaz no subsolo e para engenheiros que garantem a segurança de estruturas. No entanto, a física dessa interação é notoriamente difícil de definir. Uma simples rachadura não é apenas um vão; é um limite complexo onde as faces da rocha podem se pressionar, deslizar com fricção ou abrir e fechar bruscamente, tudo isso enquanto a onda passa em uma fração de segundo.

Por décadas, cientistas confiaram em modelos simplificados para prever esses comportamentos, muitas vezes tratando as rachaduras como se fossem feitas de molas invisíveis que se esticam e comprimem. Embora úteis, esses modelos falham quando as ondas são fortes o suficiente para forçar as faces das rochas a se tocarem, deslizarem ou se separarem completamente. Eles não conseguem capturar o movimento de "stick-slip" (aderência e deslizamento), onde a fricção mantém a rocha no lugar até que a força se torne grande demais, nem podem descrever com precisão o endurecimento não linear que ocorre quando uma rachadura é espremida até o fechamento. Um novo estudo realizado por pesquisadores da Universidade de Bergen, do Centro de Pesquisa da Noruega e da Universidade Técnica de Munique aborda essas lacunas ao construir um arcabouço computacional abrangente que simula ondas elásticas movendo-se através de rocha fraturada com um realismo sem precedentes. Em vez de depender de suposições simplificadas, a equipe desenvolveu um método que pode lidar com quatro níveis distintos de complexidade, variando desde o comportamento básico de mola até a mecânica de contato total, onde a fricção e a abertura e fechamento físico das rachaduras são calculados em tempo real.

Os pesquisadores construíram um ambiente digital unificado onde a rocha e as rachaduras são tratadas como um sistema único e interconectado. Em suas simulações, a rocha é representada como um volume tridimensional, enquanto as fraturas são modeladas como superfícies bidimensionais inseridas nela. Isso permite que o computador rastreie como uma onda viajando pela rocha atinge uma rachadura, como a rachadura se deforma e como essa deformação envia ondulações de volta para a rocha. A equipe testou quatro formas específicas de como as rachaduras poderiam se comportar. Os dois primeiros modelos usaram leis de mola, onde a força através da rachadura é proporcional ao quanto a rocha se move, com uma versão permitindo uma relação linear simples e outra usando uma regra não linear mais complexa que leva em conta o fato de que as rachaduras ficam mais difíceis de espremer conforme se fecham. Os outros dois modelos introduziram a realidade da fricção e do contato. Nesses cenários, as faces das rachaduras podem se tocar e deslizar umas contra as outras, governadas por leis de fricção semelhantes às que impedem os pneus de um carro de derrapar, enquanto também impõem uma regra de que as faces das rochas não podem passar uma pela outra.

Para garantir que seu novo arcabouço fosse preciso, a equipe realizou uma série de testes rigorosos. Eles simularam ondas atingindo uma única rachadura e compararam seus resultados computacionais com soluções matemáticas conhecidas para casos simples. Descobriram que seu método podia prever o comportamento da onda com alta precisão, correspondendo às expectativas teóricas tanto para a velocidade da onda quanto para a quantidade de energia que passava ou ricocheteava. Eles então levaram o sistema além, simulando ondas que eram fortes o suficiente para fazer as rachaduras abrirem, aderirem e deslizarem. Nesses cenários mais complexos, onde as faces das rachaduras estavam se esfregando com um coeficiente de fricção de 1,0, as simulações mostraram que o método permanecia estável e produzia resultados consistentes. Os pesquisadores verificaram que sua abordagem poderia lidar com as mudanças súbitas na física que ocorrem quando uma rachadura muda de aderência para deslizamento, uma transição que frequentemente atrapalha modelos mais simples.

Quando a equipe comparou os quatro modelos diferentes lado a lado, as diferenças nos padrões de onda resultantes foram marcantes. As simulações revelaram que a escolha do modelo importa significativamente. Quando uma rachadura era permitida a fechar e endurecer de acordo com a lei não linear, ela transmitia muito mais energia de onda do que uma rachadura modelada com uma mola linear simples. Isso ocorre porque o modelo não linear considera o fato de que, à medida que a rachadura fecha, as faces das rochas se pressionam mais fortemente, tornando a rachadura efetivamente mais rígida e mais transparente à onda. Além disso, a inclusão da fricção alterou a forma das frentes de onda. Nos modelos friccionais, a onda que emergia do outro lado da rachadura era mais larga e distorcida do que nas versões sem fricção. Os pesquisadores observaram que a fricção fazia com que as faces das rachaduras ficassem presas por períodos mais longos, atrasando a liberação de energia e alterando a maneira como a onda se dispersava.

O estudo não parou em rachaduras individuais. Os pesquisadores aplicaram seu arcabouço a um domínio tridimensional contendo seis fraturas que se interceptam, criando uma pequena rede de rachaduras com diferentes propriedades. Nesse cenário complexo, eles atribuíram diferentes limites de fechamento máximo a diferentes conjuntos de fraturas, o que significa que algumas rachaduras podiam fechar mais apertadamente do que outras. A simulação mostrou que o campo de onda tornou-se altamente heterogêneo. As fraturas que eram mais rígidas e podiam fechar mais apertadamente refletiam mais energia e transmitiam menos, criando sombras distintas no padrão da onda. Os pesquisadores foram capazes de rastrear exatamente o quanto cada rachadura se fechou e o quanto ela deslizou, revelando que mesmo pequenas diferenças nas propriedades físicas das rachaduras levaram a respostas de deformação visivelmente diferentes. Por exemplo, uma rachadura que era permitida a fechar mais apertadamente exibiu mais comportamento de "aderência", mantendo suas faces unidas por mais tempo antes de deslizar, o que, por sua vez, afetou a onda que passava por ela.

As implicações dessas descobertas estendem-se além dos resultados imediatos. Ao demonstrar que um único arcabouço computacional pode lidar com tudo, desde o comportamento simples de mola até o complexo contato friccional, os pesquisadores forneceram uma ferramenta que pode ser usada para interpretar dados sísmicos do mundo real com maior fidelidade. A capacidade de simular como as ondas se dispersam e atenuam na presença de rachaduras friccionais, de abertura e de fechamento significa que os geocientistas podem distinguir melhor entre diferentes tipos de formações rochosas e fraturas preenchidas por fluidos. O estudo confirma que ignorar os aspectos não lineares e friccionais da mecânica de fraturas pode levar a erros significativos na previsão de como a energia se move pelo subsolo. Embora o trabalho atual se concentre em interações puramente mecânicas, o arcabouço é construído sobre uma base que é compatível com modelos de fluxo de fluidos, sugerindo um caminho para simulações ainda mais realistas onde a água ou o petróleo movendo-se através das rachaduras interage com a passagem das ondas. O trabalho serve como uma demonstração de que, ao abraçar a complexidade total de como as rachaduras das rochas se comportam, podemos construir uma imagem mais clara do mundo oculto sob nossos pés.