gisis-uff / rtm_python-fortran Goto Github PK

Esperar para ver se ele realmente melhora a nossa imagem final

Relembrar os modelos usados, condição de imagem utilizada, etc

• Teste de performace - Quanto tempo demora para rodar 2, 4,8, 16 tiros, Comparar com uma modelagem em serial.

• Checar funcionamento do registro das matrizes de tempo de trânsito

• Checar o funcionamento dos snap shots

Ler as condições de tempo de excitação e correlação cruzada

• Modelo real repitido

• Sismograma modelo homogeneo

• Sismograma sem onda direta

• Matriz de tempo de trânsito

• Imagem Final

• Modelo suave

@fetim To colocando aqui os resultados para 12 tiros dos três modelos que testamos: Marmousi, Plano Paralelo e Sigsbee. O que você acha das imagens?

Marmousi 383x141: Fx0 = 10, Fz0 = 2 e SpaFonte = 30

Plano Paralelo(3 camadas) 383x141: Fx0 = 10, Fz0 = 2 e SpaFonte = 30

Sigsbee 500x500: Fx0 = 10, Fz0 = 2 e SpaFonte = 40

Observações:

1. O modelo Marmousi está sendo bem imageado. Acho que com mais tiros ele vai ficar bom.

2. O modelo plano paralelo está com uns artefatos que eu não sei da onde podem estar vindo. O sismograma dele não está mostrando somente as camadas esperadas. Segue abaixo alguns exemplos:

tiro 6:

tiro 12:

3. O modelo Sigsbee está demorando muito em cada etapa (com relação aos outros dois). Acho que pros próximos testes será melhor modificar o tamanho dele. E como você pode ver, não tem imagem. Olhando os sismogramas que não deveriam ter onda direta:

tiro 6:

tiro 12:

Eu não faço ideia do que pode estar causando esse problema no sigsbee.

Eu preciso resolver esse problema do sigsbee no maximo ate terça que vem pq quero deixar o relatorio pronto o mais rapido possivel. E eu preciso de pelo menos 2 modelos complexos. :(

O que você acha que podemos fazer??

Para os mesmos modelos que a condição de tempo de excitação

O cetale pediu para testar se fortransubroutines.so funciona sem ser compilado na máquina

Atualizar os scripts do Python na pasta compartilhada

descobrir pq a injeção não está sendo no campo inteiro!

Até agora temos 8 tiros rodando em paralelo. Pensar em como aumentar isso

Fazer script que grava posição das fontes em *.dat. Baseado no número de fontes e espaçamento

Armazenar esse valores nos vetore Fz e Fx para serem usados no algoritmo de modelagem

Praticamente pronto. Assim que tivermos as imagens finais para todos os modelos podemos fechar a issue.

Adicionar as bordas para fora!

@Caicau

• Realizar testes com o script atual

• Definir melhor estratégia para distribuição dos tiros

Resultados obtidos para o Marmousi filtro passa banda utilizando as frequencias de 10Hz a 50Hz, com slope de 2Hz

criar outro programa para RTM

• Criar uma função de plot para o modelo de velocidade

• Criar uma função de plot para o Sismograma

• Finalizar a subrotina de Snapshot

Pensar em implementar isso para o artigo

@fetim Tive problemas ao tentar compilar a F2Py. Ela não indicou nenhum problema específico, mas acredito que esteja relacionado com a declaração das variáveis strings. Você pode ver o que eu modifiquei no último commit que eu fiz no meu repositório. Qualquer dúvida só falar que eu tento responder. Se você não conseguir resolver, não tem problema. Quinta eu vou pro laboratório e tento verificar

obc - oceam bottom cable.

Implementar posicionamento dos receptores

Sugestão: Plano paralelo 2 e 4 camadas, Marmousi
O Segbee não apresentou uma boa imagem final mas podemos testar mais um pouco

Pensar em como será feita.

Use no ínicio do script em python

# -*- coding: utf-8 -*-

Fiz um teste no python 2 e funcionou a acentuação.

• Concluido ?

@fetim Como te falei, o modelo homogêneo não está só com a camada de água como esperado. Acredito que seja isto que esteja causando um sismograma com mais reflexões além da onda direta.

Tentei colocar as imagens abaixo, não sei se vai funcionar

• Instruções para execução do programa
• Inserir instruções para execução rápida da modelagem e migração
• Explicar como configurar o modelo de velocidade para o correto funcionando

Verificar quais são os códigos que estão sendo usados e separá-los

Oi @Caicau. Mais tarefas 🤙

• Criar subrotina em fortran para realizar a modelagem e registrar os sismogramas.

• Registrar os snapshots (opcional - o usuário decide)

• Registrar matriz de tempo de trânsito

• Começar a estudar a teoria de como é feita a condição de Imagem por correlação cruzada.

Implementar as bordas de Cerjan

@Caicau, Essa é a issue da migração. Se precisar de alguma ajuda coloque aqui as dúvidas

trecho para ser alterado: linhas 100 a 110
if (regTTM == 0) then

 if (**caminho_modelo  == "../modelos_utilizados/marmousi_vp_383x141.bin"**) then
    CALL Seismogram(Nt,Nx,shot,"Marmousi","../sismograma/",Seism)
 end if

 if (**caminho_modelo == '../modelos_utilizados/velocitymodel_Hmgns_wtrly.bin'**) then   
    CALL Seismogram(Nt,Nx,shot,"Homogeneo","../sismograma_modelo_camada_de_agua/",Seism)
 end if

end if

@Caicau próxima tarefa!

para remover a onda direta e a reflexão do fundo do mar

• implementação em script diferente

• Modelagem com o modelo que somente tem a camada de água e registrar os sismogramas

• Ler os sismogramas observados (gerados na primeira modelagam)

• Realizar diferençca entre sismograma observado e sismograma com somente a camada de agua

• Registrar sismograma sem onda direta

@Caicau
definir quais condições iremos usar

@fetim Felipe, o script rodou mas, eu to achando a matriz de tempo de trânsito com um artefato estranho na parte superior. Você pode dar uma olhada pra mim, por favor? Seguem abaixo as imagens da matriz e do modelo suavizado utilizado.