Este curso exemplificará o processamento sísmico de um dado sintético, modelado por traçado de raios. Esta linha abriga todos os fluxos construídos ao longo do curso para modelar e processar este dado sintético.
Mesmo que um dado sintético não seja suficiente para representar ou exemplificar toda a complexidade do processamento de uma dado real, ele já será suficiente para ilustrar os passos principais do processamento sísmico 2D em tempo.
Nota: Todos os arquivos criados no processamento desta linha são armazenados em ~/IVSimBGf. Alguns arquivos auxiliares ou temporários são armazenados em ~/IVSimBGf/temp. Caso você não tenha estes diretórios você deve cria-los antes de executar os fluxos desta linha.
Modelamento
Simulação de uma aquição marinha 2D
Por Ricardo Biloti biloti@gebrproject.com
Fluxo composto pelo(s) programa(s):
Este fluxo produz um dado sintético através do método de traçamento de raios. Ele utiliza o programa S88 Modeling, que por sua vez faz uso do conceituado programa de traçamento de raios Seis88.
O S88 Modeling simula uma aquisição ao longo de uma linha em 2D. Nesse fluxo, simulamos uma seqüência de seções de tiro comum, semelhante a uma aquisição de marinha, no qual os receptores estão distribuídos apenas para um dos lados da fonte (para direita, neste exemplo).
O modelo criado é 3D, mas todas as suas propriedades não variam ao longo da direção transversal da linha (situação conhecida como 2.5D). Portanto, apenas uma fatia do modelo, paralela à linha de aquisição, precisa ser descrita. Além disso, as interfaces dentro de cada fatia tornam-se apenas por curvas, representadas por splines cúbicas interpolantes.
O modelo definido neste fluxo é composta por quatro camadas. 251 tiros foram simulados. O deslocamento entre tiros consecutivos é de 50m. A resposta de cada tiro é registrada em 69 receptores, distribuídos regularmente, a partir de 100m à frente a posição de tiro. A distância entre receptores adjacentes também é de 50m. Cada receptor registra 6 segundos, a uma taxa de amostragem de 4ms.
I/O table
Input | (none) |
Output | /home/biloti/IVSimBGf/data.bin |
Error | (none) |
Parameters for S88 Modeling program
Parameter | Value |
Interfaces |
Instance 1 |
X coordinates (km) |
0,16 |
Z coordinates (km) |
0,0 |
Point types (see Seis88 doc for details) |
auto |
Instance 2 |
X coordinates (km) |
0,1.035,3.375,7.2,10.665,14.85,16 |
Z coordinates (km) |
1.1,1.055,0.785,1.51,1.055,1.685,1.715 |
Point types (see Seis88 doc for details) |
auto |
Instance 3 |
X coordinates (km) |
0,1.95,3.105,5.040,6.974,9,11.475,13.180,15.250,16 |
Z coordinates (km) |
3.400,3.535,3.355,2.590,3.245,2.275,2.335,2.140,2.635,2.725 |
Point types (see Seis88 doc for details) |
auto |
Instance 4 |
X coordinates (km) |
0.000, 3.880,8.250,13.600,16.000 |
Z coordinates (km) |
4.975,4.840, 3.800, 5.015, 5.275 |
Point types (see Seis88 doc for details) |
auto |
Instance 5 |
X coordinates (km) |
0,16 |
Z coordinates (km) |
5.5,5.5 |
Point types (see Seis88 doc for details) |
auto |
Velocities |
Velocities at layers' top (km/s) |
1.48,2,3.4,5.5 |
Velocities at layers' bottom (km/s) |
1.48,2,3.4,5.5 |
Densities |
Experiment |
Number of shots |
251 |
First shot x coordinate (km) |
0 |
Shot displacement in x direction (km) |
0.05 |
Source depth (km) |
0.0 |
Number of receivers |
69 |
First receiver x coordinate relative to shot position (km) |
0.1 |
Receiver displacement in x direction (km) |
0.05 |
Initial time |
0.0 |
Source type |
3 |
Modeled wave specification |
Primary P-waves |
1 |
Primary S-waves |
0 |
Seismograms specification |
Minimum modeled time |
0.0 |
Sampling rate |
0.004 |
Maximum modeled time |
6 |
Control parameters |
Seis88 binary command |
seis |
Syntpl binary command |
syntpl |
Ray shooting method |
0 |
Radius of the vicinity of a receveir |
0.0005 |
Tolerance for boundary rays |
0.00002 |
Time step in integration |
0.1 |
Initial shooting angle for primary reflections |
0.0 |
Angle step for primary reflections |
0.0010 |
Final shooting angle for primary reflections |
3.1415 |
Initial shooting angle for direct wave |
0.0 |
Angle step for direct wave |
0.0010 |
Final shooting angle for direct wave |
3.1415 |
Accuracy for ray tracing integration |
0.00001 |
Source magnitude |
1.0 |
Source frequence |
25 |
Source parameter 1 |
3.5 |
Source time shift |
0.0 |
Implosive source |
off |
Working directory |
/home/biloti/IVSimBGf/temp |
Preserve ray information |
on |
Verbose |
on |
Show debug information |
off |
Dry run (do not perform anything, just test) |
off |
Exibição de raios
Exibe os raios traçados em um experimento de tiro comum
Por Ricardo Biloti biloti@gebrproject.com
Fluxo composto pelo(s) programa(s):
O programa S88 Modeling tem uma opção para preservar a informação dos raios traçados. Quando esta opção estiver ativa, o programa Rays, em conjunto com o programa de visualização científica Grace, podem ser utilizados para exibir o modelo e os raios traçados para um experimento de tiro comum específico.
I/O table
Input | /home/biloti/IVSimBGf/temp/lu1-0070.dat |
Output | (none) |
Error | (none) |
Parameters for Rays program
Parameter | Value |
Maximum x coordinate (right corner) |
16 |
Turn off filling in of layers |
off |
Land color palette |
off |
Display all rays in black |
off |
Suppress rays |
off |
Parameters for Grace program
Parameter | Value |
Free page layout |
on |
Hide all tool bars |
on |
Add timestamp to graph |
off |
Input data in format (X Y1 Y2 Y3 ...) |
off |
Export device |
PNG |
Geometria
Acrescenta informação de geometria ao dado binário
Por Ricardo Biloti biloti@gebrproject.com
Fluxo composto pelo(s) programa(s):
- SU Add Header
- SU Set Header
O fluxo de modelagem, com base no programa S88 Modeling, produz um conjunto de dados binário, sem qualquer informação adicional. A tarefa deste fluxo é adicionar a geometria de aquisição ao conjunto de dados. Isso é feito em duas etapas. Primeiro, o programa SU Add Header cria o espaço necessário para o armazenamento de um cabeçalho para cada traço do dado, ou seja, entre cada conjunto de amostras que representa um traço será incluído espaço para o armazenamento de 240 bytes. Em seguida, o programa SU Set Header preenche vários campos do cabeçalho com informações que você especifique.
Nesse fluxo, apenas as cinco palavras-chave mais importante estão sendo definidas. Elas são a taxa de amostragem (em μs), as coordenadas X da fonte e do grupo de receptores, e o índice CDP.
O programa SU Set Header é capaz de definir a geometria de aquisição apenas para dados regularmente adquiridos, que é justamente o caso dos dados gerados pelo S88 Modeling. Para utiliza-lo é necessário prescrever os valores de cada palavra-chave do cabeçalho para o primeiro traço do dado e os incrementos para estes valores entre traços de um mesmo grupo e entre grupos (geralmente um tiro comum reunir).
Por exemplo, normalmente o conjunto de dados é adquirido com a mesma taxa de amostragem, 4ms no nosso caso. Portanto, o valor da taxa de amostragem para o primeiro traço é definida como 4000 (lembre-se que é medida em μs), e os incrementos, tanto entre traços de um mesmo tiro como entre traços de tiros distintos são zero, uma vez que tal propriedade é constante para todos os traços.
A posição do tiro se comporta de forma diferente. Definindo a origem do sistema de coordenadas coincidente com a posição da primeira fonte, o valor da propriedade X source coordinate para o primeiro traço é zero. O incremento nesta propriedade de um traço para outro, dentro do mesmo grupo, ou seja, dentro de um tiro comum, também é zero. No entanto, entre os grupos de traços, ou seja, entre famílias de tiro comum, o incremento na posição de tiro é de 50m, no caso deste exemplo.
I/O table
Input | /home/biloti/IVSimBGf/data.bin |
Output | /home/biloti/IVSimBGf/data.su |
Error | (none) |
Parameters for SU Add Header program
Parameter | Value |
Binary file format |
0 |
Trace sorting code |
3 |
Number of data traces per record |
1 |
Parameters for SU Set Header program
Parameter | Value |
List of header key words |
dt,sx,gx,offset,cdp,d2 |
Value(s) on first trace |
4000,0,100,100,1,25 |
Increment(s) within group |
0,0,50,50,1,0 |
Group increment(s) |
0,50,50,0,2,0 |
Number of elements in group |
69,69,69,69,69,69 |
Fatia
Extrai seções do dado
Por Ricardo Biloti biloti@gebrproject.com
Fluxo composto pelo(s) programa(s):
Este fluxo extrai subconjuntos do dado sísmico. Para tanto utiliza o program SU Window, que é capaz de selecionar os traços para os quais o valores de uma determinada palavras-chave do cabeçalho do traço está dentro de um intervalo prescrito.
Para extrair uma seção de afastamento comum, basta definir o parâmetro Header keyword to windon on como Offset (signed distance from source to receiver group), e os parâmetros Minimum value of the key header word to pass e Maximum value of the key header word to pass iguais ao valor do afastamento desejado.
I/O table
Input | /home/biloti/IVSimBGf/data.su |
Output | /home/biloti/IVSimBGf/co-smallest.su |
Error | (none) |
Parameters for SU Window program
Parameter | Value |
Header keyword to windown on |
offset |
Time window control |
Special selection options |
Take absolute value of the key header word |
off |
Echoes information |
off |
Análise espectral
Exibe o conteúdo de freqüência do dado
Por Ricardo Biloti biloti@gebrproject.com
Fluxo composto pelo(s) programa(s):
- SU Spec FX
- SU Simple X Image
O dado sísmico é gravado digitalmente, o que significa que as freqüências que podem ser adequadamente representados estão limitadas pela freqüência de Nyquist, relacionada com a taxa de amostragem. Durante o processamento sísmico, o conteúdo de freqüência dos dados pode eventualmente ser afetados.
Este fluxo inspeciona o conteúdo de freqüência dos dados. O programa SU Spec FX substitui cada traço pela amplitude da sua transformada de Fourier.
I/O table
Input | /home/biloti/IVSimBGf/data.su |
Output | (none) |
Error | (none) |
Parameters for SU Spec FX program
Parameter | Value |
Parameters for SU Simple X Image program
Parameter | Value |
Title of the plot |
Clean data |
Colormap |
hsv2 |
Axis 1 properties (fast dimension) |
Axis label |
Frequency [Hz] |
Axis 2 properties (slow dimension) |
Colorscale bar properties |
Display the color scale |
on |
Clipping properties |
Percentile used to determine clip value |
100.0 |
Portion of the lower range to blank out |
0.0 |
Curves |
Appearance properties |
Swap axes |
off |
Block interpolation |
off |
Quiet mode |
on |
File to save mouse picks |
/dev/stdout |
Acrescenta ruído
Acrescenta ruído branco ao dado
Por Ricardo Biloti biloti@gebrproject.com
Fluxo composto pelo(s) programa(s):
Uma nítida diferença entre os dados reais e sintéticos, é a presença de ruído. Todo dado real está inexoravelmente contaminado com o ruído coerente e incoerente. Portanto, para usar um dado sintético como exemplo no processamento sísmico é preciso acrescentar um pouco de ruído para torná-lo um pouco mais realista.
O programa SU Add Noise acrescenta ruído aleatório a um dado. Este fluxo de emprega-o para acrescentar ruído ao dado criado no fluxo de modelagem.
Adicionar ruído aleatório aos dados tem o efeito colateral de alterar o conteúdo espectral dos dados. O ruído aleatório possui uma distribuição de freqüências uniforme, enquanto os dados sísmicos tem banda limitada.
Esta discrepância pode ser inspecionada com a ajuda do fluxo de Análise espectral. Se este fluxo for executado sobre o dado sintético original o todo o conteúdo espectral ficará confinado entre as freqüências de 5 Hz e 50 Hz. No entanto, um conjunto de dados com ruído aleatório tem um espectro plano.
Para corrigir o conteúdo espectral dos dados é preciso usar um filtro passa-banda. Isto pode ser realizado através da definição do parâmetro Filter frequencies. Definindo este parâmetro como 5,10,40,50 o usuário está de fato definindo um filtro trapezoidal. Isso significa que todas as freqüências menores que 5 Hz e maior do que 50 Hz são completamente amortecidas, e as freqüências de 10Hz a 40Hz são completamente preservada.
I/O table
Input | /home/biloti/IVSimBGf/data.su |
Output | /home/biloti/IVSimBGf/data-sn20.su |
Error | (none) |
Parameters for SU Add Noise program
Parameter | Value |
Signal to noise ratio |
20 |
Flat noise probability distribution |
off |
Random number seed |
2010 |
Filter frequencies |
5,10,40,50 |
Echos information |
off |
Autocorrelação
Calcula a autocorrelação do pulso
Por Ricardo Biloti biloti@gebrproject.com
Fluxo composto pelo(s) programa(s):
- SU Auto Correlation
- SU Simple X Wiggle
Este fluxo calcula a autocorrelação dos traços sísmicos, para estimar o comprimento do pulso, que é um parâmetro necessário para realizar a deconvolução do pulso.
Como entrada, ao invés do todo conjunto de dados, apenas uma pequena parte foi usada, já que apenas uma estimativa grosseira é necessária. No nosso caso, a pequena porção de dados que alimenta o fluxo é a seção afastamento comum de menor afastamento presente no dado.
I/O table
Input | /home/biloti/AGU2010/co-smallest-sn20.su |
Output | (none) |
Error | (none) |
Parameters for SU Auto Correlation program
Parameter | Value |
Do not normalize maximum absolute output to 1 |
off |
Turn off symmetric output |
off |
Parameters for SU Simple X Wiggle program
Parameter | Value |
Axis 1 properties (fast dimension) |
Axis 2 properties (slow dimension) |
Value at which axis ends |
400 |
Wiggle properties |
Do not draw wiggle trace |
off |
Variable area fill |
1 |
Wiggle excursion in traces corresponding to clip |
1.0 |
Clip large wiggles at the boundary of the plot box |
off |
Clipping properties |
Percentile used to determine clip value |
95.0 |
Curves |
Appearance properties |
Swap axes |
off |
Use 8 point sinc interpolation |
off |
Echoes information |
off |
File to save mouse picks |
/dev/stdout |
Deconvolução do pulso
Deconvolução para aumentar a resolução vertical
Por Ricardo Biloti biloti@gebrproject.com
Fluxo composto pelo(s) programa(s):
O dado de reflexão, registrados durante o experimento sísmico, é a resposta impulsiva da Terra, conhecida como refletividade, convoluída com o pulso sísmico. O objetivo da deconvolução do pulso é remover a influência do pulso. Quando isto é feito com sucesso, o conteúdo de freqüência do dado é aumentado e a resolução vertical é melhorada.
O filtro de Wiener-Levinson, calculado pelo programa SU PEF, para deconvoluir o pulso pressupõe duas hipóteses. Primeiro, o pulso é zero-lag, ou seja, tem fase mínima, o que significa que a maior parte de sua energia está concentrada no início do pulso. Em segundo lugar, o pulso é causal. Dado que raramente isto é o caso, a deconvolução não é ótima. No entanto, a experiência mostra que os resultados deste processamento são muito melhores com ela.
Para calcular o filtro de Wiener-Levinson é necessário estimar a autocorrelação do pulso, geralmente desconhecido. Com base nas hipóteses acima, a autocorrelação do pulso pode ser razoavelmente estimada através da autocorrelação do traço sísmico registrado. Isso é feito pelo fluxo de autocorrelação.
I/O table
Input | /home/biloti/IVSimBGf/data-sn20.su |
Output | /home/biloti/IVSimBGf/data-sn20-decon.su |
Error | (none) |
Parameters for SU PEF program
Parameter | Value |
Maximum lag of prediction filter |
0.07 |
Show Wiener filter on each trace |
off |
CDP interpolation method |
linear |
Ordenação por CDP
Reordena o dado em famílias CDP
Por Ricardo Biloti biloti@gebrproject.com
Fluxo composto pelo(s) programa(s):
O conjunto de dados original é composto por famílias de tiro comum. Mas, para a maior parte do processamento sísmico, os dados devem estar organizados por famílias CDP e, em segundo lugar, por asfatamento fonte–receptor.
O programa SU Sort é capaz de ordenar os dados pelo valor de qualquer palavra-chave do cabeçalho, seja em ordem ascendente ou descendente. Ele é usado aqui para produzir uma cópia do conjunto de dados ordenados pela CDP e Offset.
Uma limitação deste programa é que ou a entrada ou a saída deve ser um arquivo em disco. Isso significa que ele deve ser colocado no início ou no fim de um fluxo de processamento.
I/O table
Input | /home/biloti/IVSimBGf/data-sn20-decon.su |
Output | /home/biloti/IVSimBGf/data-sn20-decon-CDP.su |
Error | (none) |
Parameters for SU Sort program
Parameter | Value |
Sort keys |
Instance 1 |
Sort ascending by |
cdp |
Sort descending by |
cdp |
Instance 2 |
Sort ascending by |
offset |
Sort descending by |
offset |
CDP fold
Conta o número de traços em cada família CDP
Por Ricardo Biloti biloti@gebrproject.com
Fluxo composto pelo(s) programa(s):
- SU Key Count
- Parse SU Key Count
- Grace
CDP fold é a quantidade de traços em cada família CDP. Via de regra, as famílias CDP no ínicio e fim da linha sísmica têm menos traços.
Este fluxo le o dado sísmico, ordenado por CDP, e produz um gráfico com a quantidade de traços em cada família CDP. Para tanto, três programa são encadeados: SU Key Count, Parse SU Key Count, e Grace.
O primeiro deles, SU Key Count, conta o número de traços consecutivos no dado, com um mesmo valor para uma palavra-chave específica do cabeçalho do dado. Como estamos supondo que o dado sísmico já está ordenado por CDP, basta contar o número de repetições consecutivas da palavra-chave CDP ensemble number, selecionada pelo parâmetro Header keyword to count changings do programa.
A saída produzida pelo programa SU Key Count é adequada para ser lida, mas não para ser utilizada por um programa de gráfico. O programa Parse SU Key Count processa esta saída, convertendo-a em uma tabela de valores adequadamente formatada para manipulações computacionais
O programa Grace é uma ferramenta utilizada para a exibição de gráficos 2D geral. Ele é empregado neste fluxo para gerar o gráfico final.
I/O table
Input | /home/biloti/IVSimBGf/data-sn20-decon-CDP.su |
Output | (none) |
Error | (none) |
Parameters for SU Key Count program
Parameter | Value |
Header keyword to count changings |
cdp |
Echoes information |
on |
Parameters for Parse SU Key Count program
This program has no parameters. |
Parameters for Grace program
Parameter | Value |
Free page layout |
on |
Hide all tool bars |
off |
Add timestamp to graph |
off |
Input data in format (X Y1 Y2 Y3 ...) |
off |
Export device |
PNG |
Análise de velocidade
Estima perfis de velocidades para a correção NMO
Por Ricardo Biloti biloti@gebrproject.com
Fluxo composto pelo(s) programa(s):
Análise de velocidade é um processo para estimar um modelo de velocidade a partir do conjunto de dados sísmicos. Este processo é aplicado em cada seção CMP dos dados e produz um modelo de velocidade RMS no domínio do tempo.
Este é um fluxo interativo que utiliza o programa Velan CMP, desenvolvido na Universidade Federal da Bahia, com base em programas do pacote Seismic Un*x. Ele permite ao usuário escolher os tempos e velocidades que provavelmente representam eventos de reflexão em cada família CMP.
Para cada família CMP uma seção de coerência é produzida. Os valores mais elevados nesta seção indicam a presença de potenciais eventos de reflexão. É de responsabilidade do intérprete selecionar os pontos associados à reflexões primárias.
Para definir a lista de famílias CMP a serem analisadas o usuário deve fornecer o valor do primeiro índice CDP as ser considerado, o incremento neste índice, e a quantidade de CDPs a serem processados. Quanto maior a quantidade de CDPs analisados, mais refinado do modelo de velocidade final será.
Como o programa produz uma gráfico de coerência para cada família CMP, o arquivo de entrada devem estar previamente ordenado por famílias CMP.
A seleção dos pares tempo-velocidade deve ser feita na janela de coerência, posicionando o ponteiro do mouse sobre os melhores pontos e pressionando tecla "S" no teclado. O usuário deve selecionar os pontos com tempos crescentes. Para concluir o picking de eventos em uma seção CMP o usuário deve fechar a janela de coerência, pressionando a tecla "Q" no teclado.
Após o picking em uma seção CMP, o programa mostra esta mesma seção CMP corrigida pelo efeito de NMO e o resultado do empilhamento de seus traços. Neste ponto, o usuário pode aceitar a escolha e avançar para a próxima seção CMP ou optar por descartar esta escolha e refazer o picking para a mesma seção CMP.
O resultado final deste fluxo é um arquivo de parâmetros contendo os pares tempo—velocidade, no formato esperado pelo programa SU NMO.
I/O table
Input | (none) |
Output | (none) |
Error | (none) |
Parameters for Velan CMP program
Parameter | Value |
Input data file |
/home/biloti/IVSimBGf/data-sn20-decon.su |
Output file for picked values |
/home/biloti/IVSimBGf/picked-velocities.dat |
First CDP used in the velocity analisys |
65 |
Interval between sucessives CDP |
50 |
Quantity of CDP used |
10 |
First velocity |
1400 |
Last velocity |
3000 |
Quantity of CVS panels |
10 |
Size of each CVS panel (odd numbers only) |
11 |
Number of velocities |
75 |
Semblance on power |
1 |
Strech mute applied to NMO and semblance plots |
1.5 |
Number of traces in each CDP |
0 |
Offset increment |
1 |
Threshold for significance values |
0.4 |
Plot type |
suvelan |
Percentile to clip |
99 |
Temporary directory |
/home/biloti/IVSimBGf/temp |
NMO + Empilhamento
Aplica a correção NMO e empilha o dado para produzir uma seção de afastamento nulo
Por Ricardo Biloti biloti@gebrproject.com
Fluxo composto pelo(s) programa(s):
A correção NMO é uma transformação que mapeia eventos hiperbólica de linhas horizontais. Esta correcção, quando aplicada corretamente, desconta o sobretempo de reflexão devido ao afastamento entre a fonte e o receptor, remapeando o evento para o tempo de reflexão associado a configuração de afastamento nulo. Com o evento horizontalizado, os traços podem ser somados para produzir uma traço de afastamento nulo simulado.
Este fluxo utiliza o programa SU NMO para aplicar esta transformação em cada CMP do dados de entrada. Mas para tanto é necessário um modelo de velocidade RMS, como o estimado pelo fluxo de análise velocidade.
Tendo o modelo de velocidade, corrigir os dados pelo efeito NMO e empilhar os traços para obter a seção de afastamento nulo é feito por fluxo com os programas SU NMO e SU Stack, em seqüência.
O único parâmetro necessário para o SU NMO é o arquivo de parâmetros. O usuário deve fornecer o arquivo criado pelo fluxo de análise de velocidade, contendo os eventos selecionados. O dado de entrada para esse fluxo é a todo conjunto de dados ordenados por CDP.
Após a aplicação da correção NMO, o programa SU Stack empilha os traços de cada família CMP e, portanto, produz a seção de afastamento nulo simulada.
I/O table
Input | /home/biloti/IVSimBGf/data-sn20-decon-CDP.su |
Output | /home/biloti/IVSimBGf/data-sn20-stacked.su |
Error | (none) |
Parameters for SU NMO program
Parameter | Value |
Times and velocities |
Traveltime curve |
Do not divide output samples by NMO stretch factor |
off |
Perform (approximate) inverse NMO |
off |
Read cross-line offset from trace header |
off |
Scan upward to find first sample to kill |
off |
Parameters file |
/home/biloti/IVSimBGf/picked-velocities-10.dat |
Parameters for SU Stack program
Parameter | Value |
Header word to stack on |
cdp |
Repeat the stack trace |
off |
Number of repetitions of the stack trace |
10 |
Echoes information |
off |
Migração em tempo
Migra o dado para resolver complicações moderadas
Por Ricardo Biloti biloti@gebrproject.com
Fluxo composto pelo(s) programa(s):
Apenas para estruturas geológicas muito simples, a seção empilhada pode ser diretamente interpretada. A técnica de migração é planejada para resolver situações mais complexas, facilitando a interpretação dos dados.
A migração move o evento de reflexão do local onde foi registrado para o local onde a reflexão realmente aconteceu.
A migração pode ser feita no domínio do tempo ou no domínio de profundidade. O migração em tempo é destina-se a situações onde a variação lateral de velocidade é de pequena a moderada. Para situações mais complexas, com fortes variações laterais de velocidade, a migração em profundidade é recomendada. No entanto, a migração em profundidade necessita um modelo de velocidade em profundidade, que é difícil de ser estimado. Além disso, o processo de migração em profundidade também é muito mais sensível a erros no modelo de velocidade.
Este fluxo migra os dados em tempo com a ajuda do programa SU Stolt. Apesar de sua simplicidade, devido à robustez da migração em tempo, mesmo com uma estimativa pobre do modelo de velocidade, a seção migrada é mais facilmente interpretada. Nós exploramos esse fato, migrando os dados com um modelo de velocidade constante.
I/O table
Input | /home/biloti/IVSimBGf/data-sn20-stacked.su |
Output | /home/biloti/IVSimBGf/data-sn20-time-migrated.su |
Error | (none) |
Parameters for SU Stolt program
Parameter | Value |
CDP properties |
Minimum CDP in dataset |
1 |
Maximum CDP in dataset |
569 |
Distance between adjacent CDP bins (m) |
25 |
Number of offsets to mix (for unstacked data only) |
1 |
RMS velocities |
Instance 1 |
Times corresponding to RMS velocities below (s) |
0 |
RMS velocities corresponding to times above (m/s) |
1480 |
Instance 2 |
Times corresponding to RMS velocities below (s) |
3.8 |
RMS velocities corresponding to times above (m/s) |
1800 |
Stretch factor |
1.0 |
Scale factor to velocities |
1.0 |
Length of side tapers (# traces) |
0 |
Length of bottom tapers (# traces) |
0 |
Echoes extra information |
off |
Parameters for SU X Image program
Parameter | Value |
Title of the plot |
Time migrated section |
Colormap |
hsv6 |
Axis 1 properties (fast dimension) |
Grid line style |
none |
Axis 2 properties (slow dimension) |
Axis label |
CDP |
First sample |
1 |
Sampling interval |
1 |
Grid line style |
none |
Colorscale bar properties |
Display the color scale |
off |
Colorscale bar width in pixels |
16 |
Colorscale x-position in pixels |
3 |
Clipping properties |
Percentile used to determine clip value |
100.0 |
Balance clip equally |
off |
Portion of the lower range to blank out |
0.0 |
Curves |
Appearance properties |
Upper left corner x coordinate in pixels |
50 |
Upper left corner y coordinate in pixels |
50 |
Window width in pixels |
550 |
Window height in pixels |
700 |
Swap axes |
off |
Block interpolation |
off |
Plot title color |
red |
Axes labels color |
blue |
Font name for colorscale bar labels |
times_roman10 |
Grid color |
blue |
Quiet mode |
on |
File to save mouse picks |
/dev/stdout |