Exploração sísmica em terrenos complexos

Prévia do material em texto

Exploração sísmica em terrenos complexos: uma experiência de processamento
nos Apeninos do Sul
Abstract 
	Discutimos uma seqüência de processamento de dados adotada para reprocessar uma linha sísmica que atravessa o sul do sul da Itália Apeninos do mar Tirreno até a margem adriática e investigar as áreas de overthrust e foreland. Primeiro determinamos as principais causas do muito baixo a proporção S / N nos dados de campo e, em seguida, propor um processamento sequência destinada a explorar o conteúdo do sinal, também fazendo uso do conhecimento geológico a priori desta área. Nosso trabalho indica uma combinação de causas para a baixa qualidade dos dados sísmicos. Estes incluem o comprimento de a propagação (cerca de 20 km) que é desfavorável por causa de A rápida variação na geologia da superfície próxima, tectônica complexidade, aquisição de linha torta e topografia em bruto associado a rochas de afloramento caracterizadas por velocidades altamente variáveis.
	Com base no resultado desta análise de dados, apresentamos uma sequência de processamento conduzida pelo conhecimento da região configuração tectônica e pelo conhecimento do raso geologia subterrânea. O principal esforço é remover o Componentes de ruído de grande porte, perto da superfície. O menor A relação S / N torna impossível ou quase impossível aplique técnicas altamente sofisticadas como como migração em profundidade ou datação da equação de onda. Portanto, usamos técnicas robustas especificamente projetadas para resolver cada problema que degradou a qualidade dos dados. O mais relevante dessas técnicas foram a remoção de maus traços onde a qualidade inaceitavelmente baixa foi detectada pela energia e critérios de degradação de frequência; estimativa e correção para mudanças de tempo estático atribuíveis a condições de superfície próxima; Otimização da classificação intermediária comum (CMP) para atenuar os efeitos deletérios da linha torta aquisição; Aplicação de uma técnica de empilhamento ponderada para maximizar o poder da pilha e a aplicação de prestack f -x deconvolução para atenuar o ruído não correlacionado. O resultado dessa seqüência de processamento é comparado com o resultado de uma seqüência mais padrão que foi anteriormente aplicado aos mesmos dados e também é discutido em termos do possível modelo geológico que possa provar. A realização de uma seção sísmica mostrando bastante contínua e eventos estruturados até os oito, dependendo na interpretação, podem estar relacionados à descontinuidade de Moho ou para camadas sedimentares muito profundas a eficácia da abordagem de processamento que propomos.
Introdução
	Exploração sísmica em terrenos complexos, como cintos de pressão e regimes tectônicos complexos são dificultados por vários fatores que muitas vezes levam a dados de baixa qualidade. Diferentes abordagens de processamento oferece formas possíveis de melhorar a qualidade do final seções em áreas de sobretensão. Wu et al. (1998) exploram as vantagens de migração de prestack de topografia no Husky conjunto de dados estruturais e propor um procedimento iterativo onde A forte interação entre geofísicos e intérpretes ajuda otimize o campo de velocidade para a migração de prestack da topografia. Zhu et al. (1998) afirmam que a migração de profundidade de prestack de um dado horizontal após o tomo-datuming, isto é, após o pré-cálculo datação de equação de onda com velocidade próxima da superfície estimada da tomografia de tomografia computadorizada, produz imagens melhores do que direto migração da topografia. Eles testam o procedimento em dados reais de várias áreas complexas, incluindo o Husky Conjunto de dados estruturais e conclui que a sua suposição é particularmente verdade quando topografia em bruto e forte superfície próxima A variação de velocidade existe simultaneamente. Infelizmente, ambas abordagens tendem a ser inaplicáveis ​​quando a relação S / N do os dados do campo são muito baixos; assim, técnicas que são menos sensíveis ao ruído do que a migração de prestack são necessários. Este é o caso de uma linha sísmica regional que atravessa os apeninos do sul (Figura 1) perto de importantes descobertas de hidrocarbonetos. Conforme mostrado na Figura 2, os domínios estruturais investigados por Esta linha inclui as principais áreas de sobreavaliação e foreland. Extenso exploração de hidrocarbonetos na bacia do foredeep e no o cinto de impulso levou a um conhecimento de primeiro plano muito maior da arquitetura interna dos apeninos do sul (ver Balduzzi et al., 1982a, b; Mostardini e Merlini, 1986; Casnedi, 1988a, b; Sella et al., 1988; Casero et al., 1991; D'Andrea et al., 1993; Roure e Sassi, 1995; La Bella et al., 1996; Mattavelli et al. 1993). A espinha dorsal da cadeia montanhosa tem uma profundidade sistema dúplex compressional - um alvo para exploração de petróleo, sobreposto por uma espessa pilha de nappes sem raízes. O enterrado O sistema duplex é uma pilha de cavalos de carbonato mesozoico-terciário
Desmontado de um reino de plataforma de águas rasas (Plataforma de Apúlia) ao longo de evaporites triássicos. As nappas sem raízes consistem em Sementes sedimentares mesozoicas-terciárias da bacia e plataforma overlain inconformably por Miocene-Pleistocene pushtop-sheettop
depósitos. Processos de violação de duplex e fora de sequência A propagação de impulso contribuiu para produzir um complexo bastante complexo. arquitetura estrutural (por exemplo, grandes pilhas informais dentro as unidades do telhado). Linhas sísmicas comerciais e registros de poços são um ferramenta fundamental para decifrar a estrutura dos Apeninos. Mesmo assim, devido ao comprimento de registro limitado (geralmente 5 s TWT) das seções sísmicas, as estruturas profundas do cinto de impulso são geometrias estruturais de primeira ordem críticas e desconhecidas ainda não foram resolvidos. Por exemplo, não sabemos a localização da borda traseira do sistema duplex na subsuperfície ou a profundidade do único impulso em correspondência com o Margem tirrenia da cadeia de montanhas. Além disso, é Desconhece se o porão está envolvido no encurtamento crustal ou até que ponto a crosta da Púlia não deformada se estende por baixo o cinto de impulso. Então, criamos e adquirimos uma nova região linha de reflexão sísmica que atravessa todo o Apenninic cadeia para responder a essas perguntas. A linha foi adquirida no âmbito do projeto de RCOP (o projeto italiano de exploração de crosta profunda) com os parâmetros ilustrado na Tabela 1. A mesma linha também foi adquirida com uma fonte de energia explosiva, usando a mesma gravação espalhada como a aquisição da Vibroseis e, assim, permitir uma quantitativa comparação das duas fontes. Nós fomos designados para reprocessar Os primeiros 100 km de dados de Vibroseis a partir do Tirreno costa para o site de Fiumara Venosa (da CMP 5-2600, Figura1). Em vez de se concentrar em cada etapa do processamento, nosso artigo enfatiza os principais passos que visivelmente melhoraram os dados qualidade. Assim, descrevemos os problemas mais importantes afetando os dados de campo e depois apresenta o processamento adotado soluções.
Análise da Qualidade de Dados
Para efetivamente definir uma estratégia de processamento apropriada, primeiro realizamos uma análise detalhada dos dados para identificar os diferentes tipos de problemas e o ruído a serem abordados e atenuados pelas fases de processamento subsequentes. Os dados disponíveis para esta análise foram
FIG. 1. Mapa de localização que se refere à área indicada no quadro superior; alguns locais de CMP ao longo da linha são indicados.
1) o conjunto completo dos registros de campo Vibroseis e dinamite (empilhados e correlacionados), relatório do observador, dados do uphole, mapa de elevação, resultados de testes de origem e gráfico de empilhamento;
2) um mapa geológico (escala 1:25 000) detalhando o litológico unidades; e
3) uma seção transversal geológica coincidente com o caminho do linha sísmica.
Condições próximas à superfície
	Conforme mostradona Figura 3, a rugosidade da topografia e a sucessão abrupta, dentro de um único comprimento, de Litologias duras e superficiais caracterizam o curso da linha sísmica . Essas condições exigem etapas de processamento muito cuidadosas como computação estática e análise de velocidade; além disso, sob tais circunstâncias, esquemas de muting convencionais de primeias pausas não são aplicáveis. Abaixo da seção transversal geológica da Figura 3, também planejamos as velocidades estimadas do refractor calculado para correções estáticas de dados. O subtraculo A velocidade é comparada com a geologia na superfície: a final modelo de velocidade mostra boa correspondência entre litologia e velocidade, resultando em correlação ótima quando as falhas são presente. Dentro de um tamanho de propagação, pode haver várias centenas diferença de metros na elevação e mais de 1000 m / s de variação da velocidade da superfície próxima. Estas condições severamente teste muitos dos pressupostos simplificadores feitos pelo convencional algoritmos e métodos de processamento.
	Outra consequência da topografia em bruto e as dificuldades A acessibilidade para os caminhões Vibroseis é a linha crooked aquisição dos dados sísmicos. É mostrado um exemplo típico
na Figura 4 na área entre Albanella e Muro Lucano localizado no lado ocidental da cadeia dos Apeninos. Tão alto tortuosidade em muitas partes da linha provoca problemas visíveis na classificação média comum (CMP) e degrada o imagem sísmica, particularmente em períodos pouco profundos. Assim, o slalom A linha deve ser projetada cuidadosamente para evitar a degradação.
Contaminação do ruído incoerente
	Muitas coleções de Vibroseis estão dramaticamente contaminadas por ruído incoerente que aparece de repente em grande deslocamento(4000 m). Esse ruído geralmente parece ruído não correlacionado (branco na banda de Vibroseis), afetando todo o registro de rastreamento comprimento. Em alguns casos, o espectro de ruído é colorido, apresentandouma frequência dominante entre 20 e 30 Hz. Em qualquer caso, Tanto o poder de ruído quanto o espectro de ruído são, basicamente, o tempo invariante. Esse ruído pode ser de configurações de ganho inadequadas
dos amplificadores de ponto flutuante do gravador de grupo sísmico (SGR) e / ou a outros motivos ainda por determinar. Em alguns casos, A ocorrência na propagação deste ruído peculiar coincide
com a transição para rochas de superfície dura. Assim, o receptor pobre O acoplamento pode ser um fator concorrente na degradação do sinal, embora um padrão convincente consistente em superfície foi encontrado apenas ao longo de alguns segmentos de toda a linha. Também descobrimos que o 
FIG. 2. Seção esquemática geológica em todo o sul dos Apeninos ao longo do traço da linha sísmica. Resultados de perfuração para hidrocarbonetos A exploração forneceu dados importantes para restringir as isóbuas base-de-Plioceno / Pleistoceno na área de Foreland, a arquitetura interna das unidades do telhado da correia de impulso e a parte superior do sistema duplex duplex na metade direita da figura. As estruturas profundas da cadeia montanhosa sob Monte Soprano e Monte Alburno, pelo contrário, nunca foram diretamente explorado uma vez que a profundidade máxima dos poços na área excede ligeiramente 4000 m. Portanto, a seção geológica nesta parte de
O perfil é baseado na interpretação da linha sísmica.
conteúdo de ruído no Vibroseis e os dados explosivos aumentam significativamente em áreas de estresse tectônico intenso, o que pode ser a causa da dispersão de campos de onda difusa. Isto é particularmente evidente em alguns segmentos da linha, como o paralelo uma falha normal regional que separa o carbonato de Alburno maciço da planície costeira de Sele (Figura 5) e em áreas de superação. Os dados de melhor qualidade foram registrados em áreas de menor complexidade estrutural (Figura 6). A má qualidade das quantidades significativas de dados nos forçou. encontrar uma técnica capaz de detectar automaticamente vestígios onde o A qualidade está abaixo de um limite predefinido. Para este fim, usamos dois indicadores de qualidade: variação de amplitude média com o tempo e variação de freqüência com o tempo. Um bom traço deve exibir ambos decadência de amplitude e decadência de frequência com o tempo. Pouco traço é um que mostra um status constante de amplitude e freqüência conteúdo com o tempo. As amplitudes de rms e os espectros de freqüência foram computados em intervalos de tempo de cerca de 0,5 s, deslizando ao longo da eixo temporal dos dados. Após alguns testes, limiares adequados para As amplitudes de rms e as variações de freqüência de pico foram selecionadas, e todos os traços de campo foram avaliados por esses atributos. Dependendo sobre a boa ou má natureza resultante de cada traço, um código apropriado foi inserido em um local do seu cabeçalho. Traços de natureza incerta, ou seja, com indicações contrastantes da dois indicadores, ou com características muito próximas dos limites foram marcados com um código distinto que habilitou a classificação separada e inspeção visual para definir seu status final. A Figura 7 mostra a porcentagem de traços ruins em relação a o número total de vestígios em cada fonte comum se reúnem. O As litologias de superfície na localização da fonte também são mostradas. O poucas áreas com 0% de traços ruins foram inacessíveis para a Vroserose caminhões; assim, nenhum dado foi gravado. A litologia na fonte A localização não mostra uma relação convincente com a ocorrência de maus traços. Podem estabelecer-se correlações mais fortes com a complexidade estrutural da superfície próxima. Quanto maior a ocorrência de traços ruins entre as fontes 500 e 750 está associada com a posição de uma falha regional (veja também a Figura 3). A porcentagem de traços ruins permanece constantemente alta (40% como uma média aproximada) ao longo de toda a cadeia apenina
e começa a diminuir em correspondência com a transição para áreas menos perturbadas da bacia do Foreland (área do Sub- Apenninic Hills, Figura 2). A Figura 8 mostra a forte dependência
de qualidade de dados na fonte para o deslocamento do receptor. Em compensações maior do que ± 4000 m, mais de 30% dos dados são classificados tão mau.
Tabela 1. Parâmetros de aquisição da linha CROP 04.
	Os resultados desta análise também foram armazenados em um banco de dados que pode ser usado para estudos adicionais sobre possíveis causas e efeitos relações entre qualidade de dados e parâmetros sísmicos (por exemplo, tortuosidade da linha, tipo de fonte de energia, fonte-receptor offset, etc.) ou entre a qualidade dos dados e as características geológicas (por exemplo, litologia superficial na fonte e / ou no receptor, configuração tectônica local, etc.). O resultado deste específico estudo, não detalhado neste artigo, é útil para o planejamento aquisição sísmica em áreas próximas ou em zonas com características.
SEQUÊNCIA DE PROCESSAMENTO
A sequência de processamento aplicada aos dados está descrita em Figura 9. Em vez de descrever cada passo em ordem seqüencial, discutimos grupos de operações de processamento, não necessariamente seguindo uma ordem de fluxo de dados. As etapas de processamento em cada grupo têm o objetivo comum de resolver um problema específico detectado na análise de dados anterior. Focamos nossa discussão sobre o seguintes grupos de etapas de processamento:
	1) operações de atenuação do ruído aplicadas em toda a parte
 	sequência de processamento, mas tentamos remover o
	maiores componentes de ruído o mais cedo possível;
	2) remoção de turnos de tempo estático devido a condições de superfície próxima,
	o que implica a aplicação de duas formas diferentes
	métodos de cálculo de estática de refração e uma técnica
	de estática residual consistente em superfície;
	3) otimização da classificação e empilhamento CMP para atenuar
	os efeitos deltérios da aquisição da linha torta; e
	4) análise de velocidade orientada por modelo para explorar a priori geologicalconhecimento da área.
	Para as outras operações de processamento, não é tão essencial para alcançar
	o resultado final como mencionado acima, mencionamos apenas o
	Parâmetros de configuração pertinentes (Figura 9).
FIG. 4. Um exemplo da aquisição da linha torta. Note o dispersão relevante dos locais reais do CMP.
FIG. 3. Seção geológica das formações rasas ao longo do perfil sísmico. O gráfico abaixo da seção mostra o refrator velocidade estimada para o cálculo das correções estáticas. Observe a correlação entre alterações de litologia e velocidade
FIG. 5. Um coletor de fonte comum (fonte 530) localizado em uma área de estresse tectônico 
intenso e paralelo a um falha regional.
FIG. 6. Recolhimentos de fonte de informações localizados na área de foreland perto de CMP2600 (fonte 1579). A qualidade dos dados é melhor do que o exemplo da Figura 5.
Atenuação do ruído
	Com base na classificação de rastreamento realizada na análise de dados fase, os traços onde a qualidade, conforme definido pela amplitude e os indicadores de frequência, caiu abaixo de um determinado limite foram removidos imediatamente do conjunto de dados. Isso levou a rejeitar cerca de 30% dos traços de campo. Em algumas áreas, removemos até 70% dos dados (Figura 7). Esta foi nas zonas de perturbação tectônica alta e complexa geologia próxima à superfície - particularmente em deslocamentos de fonte-receptor elevados. Recolhimentos de compensação comum, que imitam o refletor (deslocado no tempo) A geometria do coletor de dobra única de incidência normal, foram considerados o domínio mais apropriado para se candidatar f-x deconvolução. De fato, esse tipo de agrupamento de dados permite avaliação fácil da eficácia da operação e umaControle de qualidade importante verifica o risco de criar falsos eventos aumentando o ruído correlacionado. A desvolução f-x foi realizada em coletores de compensação comum depois de aplicar estática correções para o datum final (plano) e depois de preencher as lacunas de vestígios faltantes com traços nulos. Uma equalização de rastreamento de amplitude também foi aplicado antes da deconvolução f-x para evitar artefatos que resultam de dados não balanceados. Caixas de compensação de 160 m foram selecionados de modo que, em condições nominais, dois vestígios de cada um registro de origem pertencia ao mesmo agrupamento de compensação comum. Comprimento total dividido no comprimento do cabo ativo (20 km) e a complexidade da área sugeriu a classificação das compensações positivas e os negativos; assim, obtivemos um total de 120 comutadores Seções. A Figura 10 mostra o agrupamento em um deslocamento constante de 3000m antes e depois de aplicar o filtro f-x. Os parâmetros do filtro (20 traços de correlação em intervalos de tempo de 0,4 s, filtro de seis rastros comprimento) foram cuidadosamente selecionados para evitar a introdução excessiva coerência nos dados. Os mesmos parâmetros foram adotados para 
FIG. 7. Classificação da qualidade do traço de Vibroseis ao longo do perfil com base em amplitude e em indicadores de freqüência versus tempo. O maior número de traços ruins, entre as estações 500 e 750, 
outra deconvolução f-x aplicada nos dados da ponderada seção de pilha no final da seqüência proposta.
 Outra operação efetiva na redução do ruído é a pilha ponderada, onde os pesos são as estimativas locais da relação S / N. O cálculo do local, ou seja, amostra por amostra e trace-by-trace, a relação S / N dos dados é realizada por meio de de decomposição de valor singular (SVD) ou, alternativamente, por correlacionando vestígios adjacentes (Grion e Mazzotti, 1998). Estes os pesos foram calculados em coletores de CMP corrigidos de NMO após análise de velocidade final e aplicação de estabilidade residual estimando a coerência dos dados dentro de um pequeno espaço de tempo deslizante janela. As dimensões da janela são três ou cinco vestígios por um Comprimento de tempo escolhido para conter pelo menos uma onda. Figura 11 mostra o CMP original (à esquerda) e os pesos pertinentes (à direita) calculado com um método baseado na soma dos valores de desfasamento zero da correlação cruzada entre três vestígios adjacentes e com um tempo de janela de 112 ms. Pesos maiores são atribuídos para amostras de dados que mostram um espaço espacial mais pronunciado continuidade, enquanto as zonas ruidosas são penalizadas. Estes coeficientes pode ser usado diretamente para pesar as amostras de dados, conforme feito neste caso, ou pode ser processado posteriormente (por exemplo, alisado ou obtido) para produzir os pesos de empilhamento final. Pesos semelhantes foram também obtido com um método mais exigente da CPU baseado em a técnica SVD. Como os resultados foram muito semelhantes, o método de crosscorrelation, duas vezes mais rápido do que o método SVD, foi preferido.
Remoção de turnos de tempo estáticos
O cálculo das correções estáticas exigiu muita força de trabalho. A escolha da primeira quebra foi realizada manualmente em ambos o Vibroseis e os dados de dinamite. Escolha inicial de os dados de Vibroseis são conhecidos por ser problemáticos, principalmente por causa de os eventos precursores que são paralelos ao primeiro intervalo verdadeiro. Nisso caso a redução foi dificultada pelo ruído que muitas vezes obscureceu as primeiras pausas. Para tornar esta operação mais confiável, incluímos os dados explosivos (baixos) que mostram muito mais limpo primeiro quebra e guiou a escolha no lado adjacente Vibroseis se reúne. Tentativas de operação com dois refractores ficaram frustrados com a insuficiente cobertura de primeiros intervalos, especialmente no deslocamento distante, onde a escolha precisa foi raramente possível devido à contaminação do ruído. Assim, trabalhamos com um único refrator. Em média, cortes confiáveis foram retirados a 3000 m.
FIG. 8. Classificação da qualidade de rastreio de Vibroseis versus offset de fonte-receptor. A porcentagem de traços ruins aumenta com o aumento offset de fonte-receptor. Para compensações superiores a ± 4000 m, a porcentagem de traços ruins excede 30%.
FIG. 9. Fluxograma da sequência de processamento. Os passos marcados por uma caixa sombreada são discutidos em detalhes no texto. O Parâmetros de configuração relevantes das outras etapas são indicados.
A velocidade da primeira camada do modelo foi reconstruída interpolando as informações derivadas de tempo de burracos acima vezes. Como conseqüência da complexidade superficial e a aquisição da linha torta, a profundidade da velocidade da superfície
 
FIG. 10. Recolhimento de compensação comum (3000 m) em relação à área ocidental (CMP 5-565) com correções estáticas ao ponto de referência plano (a) antes e (b) após a deconvolução f-x.
O modelo é estimado através de códigos desenvolvidos internamente baseados nos métodos WIM (Zanzi, 1996) e GRM (Palmer, 1980). O GRM foi aplicado em tiros explosivos integrados com um Seleção de registros da Vibroseis onde fonte explosiva Os dados eram muito escassos. A correspondência de fase do explosivo e Os dados da Vibroseis foram realizados calibrando a escolha procedimento em locais onde ambas as fontes foram usadas. Onde as primeiras chegadas da dinamite reversa e reversa Os tiros não se sobrepuseram suficientemente, o 
FIG. 11. Recolha CMP corrigida de NMO com estatística residual aplicada (esquerda) e os pesos de empilhamento correspondentes (à direita) calculados por meio de crosscorrelation. O valor dos pesos é dada pela excursão das wavelets. Os pesos mais altos são localizado ao longo de reflexões coerentes. Baixas zonas de qualidade de dados são destacado pelos pesos correspondentes de baixa amplitude.
método WIM foi aplicado nos dados da Vibroseis e os resultados foram mesclados. Estática foram calculados para um ponto de referência flutuante projetado por suavização a topografia. Foi utilizada uma velocidade de substituição variável para explicar as rápidas variações laterais da superfície próxima litologia. Uma lição que aprendemos é que o projetoda flutuação O ponto de referência deve ser calibrado com o comprimento da extensão espalhada. Neste caso específico, dada a propagação bastante longa, uma O alisamento da topografia foi necessário para definir um adequado datum flutuante. O valor final (plano) foi estabelecido em uma elevação de 400 m acima do nível do mar, que é a elevação média do perfil, resultando em um intervalo estático de -200 a +100 ms. (Para um exemplo da tendência estática, veja a elevação e a estática perfis no topo da Figura 16.) 
A correspondência entre a variação da velocidade do refractor e as mudanças na litologia superficial são visíveis (Figura 3). Observe também que a repetição de unidades litológicas ao longo de a linha combina perfeitamente com os valores de velocidade do refractor (por exemplo, o Carbonatos de Alburno-Cervati). Essa característica ajudou a determinar as correções estáticas ótimas para o Monte Soprano área. Nesta zona, onde a taxa de repetição de tiro é muito baixa porque de problemas de acessibilidade, as velocidades do refractor diretamente estimado a partir dos primeiros intervalos (linha pontilhada, Figura 3) não são confiáveis e não são comparáveis ​​aos obtidos para o mesmo litologia em outras partes da linha (4000-4500 m / s). Assim, para Esta seção específica, correções estáticas foram determinadas pelo teste e erro e observando visualmente os resultados. Cada computação estática O teste foi realizado utilizando uma velocidade de refractor diferente, incluindo aqueles pertencentes à mesma litologia ao longo da linha. 
A boa correlação entre o modelo derivado do refractário e os elementos geológicos indicam boa confiabilidade em nosso estimativas da velocidade do refractor. Além disso, a eficácia de as correcções estáticas de dados que aplicamos são evidentes na comparação as pilhas de velocidade constante antes e depois da aplicação (Figura 12). 
Os turnos de tempo residual foram ainda corrigidos por meio de estabilidade residual consistente em superfície em NMOcorrectedCMPgathers. Estes foram calculados ao longo dos portões do tempo centrado nos principais eventos sísmicos, com uma largura temporal > 1 s. Em geral, impomos pelo menos três portas horárias para cada CMPlocation e deixe o algoritmo procurar a solução que maximizou simultaneamente o poder da pilha de todo o tempo janelas consideradas. O turno máximo de tempo permitido em cada a estação era de 30 ms.
Otimização da classificação CMP
A alta tortuosidade da linha degrada significativamente a pilha qualidade, especialmente para as reflexões rasas. Para limitar este negativo efeito tanto quanto possível, testámos três critérios diferentes para a classificação CMP, como ilustrado na Figura 13. Curvatura mínima, o método utilizado na sequência de processamento anterior, favorece a distribuição de CMP que minimiza a curvatura do linha. Os resultados deste procedimento de classificação são mostrados na Figura 13a, onde apresentamos um supergather que é a soma de cinco adjacentes CMPs na área entre as Montanhas Albanella e Alburno (CMP 617). Com a densidade máxima, a linha segue o zonas do perfil onde os CMPs estão mais densamente espaçados. A Figura 13b ilustra os resultados desta reunião de CMP com o mesmo supergather. Finalmente, em curto-circuito, a linha resultante segue um caminho através das zonas mais densamente povoadas por CMPs compostas por traços com deslocamento de fonte-receptor até um valor dado. A Figura 13c mostra os resultados desta técnica até um desvio máximo de ± 2,5 km (cerca de um quarto do spread). A Figura 13 revela que o método de curto-deslocamento produz o resultados S / N ótimos. Na verdade, com este método, as formas de As reflexões são mais distintas e se aproximam de um quasihyperbolic deslocamento (por exemplo, as reflexões de 0,6-0,8 s). O As mesmas observações foram alcançadas ao examinar os resultados em outros segmentos da linha. Somente em segmentos limitados da linha, entre os CMPs 1271-1390 e 1470-1490, fez o máximo O critério de densidade produz os melhores resultados. Assim, com o Exceção desses dois segmentos onde usamos o máximo método de densidade, reunimos todos os CMPs da linha seguindo o critério de compensação curta.
Análise de velocidade orientada por modelo
Além disso, a determinação do campo ótimo de velocidade de empilhamento Esta fase foi consideravelmente melhorada em relação ao original dados, ainda havia áreas em que a escolha da velocidade era difícil. A variação da velocidade da superfície lateral e o baixo S / N
O índice nos levou a prosseguir com uma análise preliminar global do campo de velocidade. Pilhas de velocidade constante (CVS) para toda a linha com intervalos de velocidade de 250 a 500 m / s foramexigiu muito mão de obra. Embora a qualidade dos dados em produziu e estudou cuidadosamente para determinar a velocidade inicial campo. Isso foi refinado nos espectros de velocidade de CMP para uma definição mais detalhada. A pilha resultante foi comparada com o CVS, e o procedimento foi iterado até o resultado foi considerado satisfatório. Mesmo a menor evidência de lateral coerência que apareceu ou desapareceu durante essas iterações foi discutido com intérpretes e comparado com o bdados geológicos disponíveis da área. Em alguma parte da linha, os dados geológicos (afloramentos, dados de furos publicados) tornaram possível definir modelos plausíveis da subsuperfície que foram
usado como restrições suaves para orientar a análise de velocidade. Além disso, A revisão do campo de velocidade usando a análise de espectros de velocidade foi também realizada após a aplicação de estática residual.
FIG. 12. Resultados de correção estática: (a) CVS (5500 m / s) na área entre as Montanhas Albella e Alburno antes da aplicação as correções estáticas; (b) o mesmo CVS que em (a), mas com a estática ao ponto de referência final (400 m acima do nível do mar) aplicada.
FIG. 13. Resultados da otimização da reunião do CMP. Três diferentes critérios de coleta são comparados em um super CMP (soma de cinco CMPs adjacentes): (a) curvatura mínima
(aplicado na sequência de processamento anterior); (b) máximo densidade; e (c) deslocamentos curtos. A flecha branca aponta para as reflexões que se tornam particularmente evidentes quando a compensação curta o critério é usado.
Discursão dos Resultados
Nossos resultados finais em termos de seção de tempo empilhado são comparados com o resultado de uma sequência de processamento que teve já foi aplicado. A Tabela 2 mostra esta seqüência padrão. Difere da que propomos em muitos aspectos: 
não houve remoção baseada em medição de baixa qualidade traços;
2) a estabilidade da refração foi calculada diretamente no plano datum e com uma velocidade de substituição constante, enquanto a escolha das primeiras quebras foi realizada em vibroseis apenas dados;
3) a classificação CMP foi realizada com o menos efetivo método mínimo de curvatura; e
4) tenta aumentar a coerência dos dados envolvidos Estabilidade residual CMP-consistente (trimtics trim) e pós-recuperação técnicas de aprimoramento de mergulho.
FIG. 14. Comparação dos resultados finais na área Controne: (a) seção final da pilha da sequência de processamento anterior; (b) seção de pilha final com a seqüência usada neste trabalho.
FIG. 15. Comparação dos resultados finais na área de San Fele: (a) seção de pilha final da sequência de processamento anterior; (b) final seção de pilha produzida com a seqüência usada neste trabalho.	
Assim, além do filtro de simulação de matriz, não havia tentativa precoce de remover a grande quantidade de ruído que afeta os dados e a estática e a classificação CMP foram realizados de forma diferente. As tentativas de recuperar alguma coerência nos dados foram feito in extremis por meio do aprimoramento da coerência pós-cotação técnicas. Apresentamos a comparação entre os resultados dos dois Processando seqüências em duas áreas de todo o perfil. Um é localizado no início da linha, entre CMPs 500 e 790, perto da área de Controne. A nova seção depilha deste área tem eventos contínuos e estruturados às vezes até 8 s (Figura 14b); tais eventos não eram claramente discerníveis no versão anterior (Figura 14a). Dependendo da interpretação, esses refletores podem estar relacionados à descontinuidade do Moho ou a camadas sedimentares muito profundas. A outra área, perto de a área de San Fele (CMPs 1720-1960), atravessa um espesso antiformal Pilha de imbricatos de Lagonegro dentro das unidades de telhado de o sistema duplex (Sanformé antiforme). Este é o mais difícil segmento do perfil e, como mostrado na Figura 15a, originalmente teve uma resposta de reflexão sísmica muito baixa. Este tem claramente melhorado, e a pilha final agora apresenta muitos eventos (pelo menos na parte superior da seção) até 4 s (Figura 15b). A nova sequência de processamento produziu melhorias semelhantes nos outros segmentos da linha. Foi tomado muito cuidado verificando a confiabilidade dos resultados. Na verdade, como a qualidade dos dados do campo bruto eram tão pobres, estávamos constantemente preocupados sobre a possibilidade de criar eventos falsos, introduzindo processando artefatos ou aumentando um tipo particular de barulho. Para evitar este resultado mais indesejável, verificamos o resultado de cada etapa de processamento para verificar que nada extraordinário aconteceu e que os resultados intermediários estavam de acordo com os modelos geológicos plausíveis de a área. Além disso, sempre que possível, operações consistentes em superfície foram aplicadas, e nenhuma tentativa foi feita artificialmente ou excessivamente forçar a coerência dos dados. Finalmente, as operações multicanal que são projetados para atenuar o ruído aleatório e que melhorar implicitamente a continuidade das reflexões, por exemplo, f-x deconvolução, foram realizadas após uma cuidadosa escolha do parâmetros. No que diz respeito à qualidade, a seção sísmica reprocessada
é comparável às melhores linhas comerciais da região.
FIG. 16. Pilha final do segmento ocidental do perfil (das CMP 5 a 790). Os gráficos para elevação, estática para o datum final, e a dobra da pilha também são mostradas. O corpo livre de reflexão entre 2 e 4 s na extremidade esquerda da seção e entre 0 e 2 s em o lado direito corresponde a uma sequência espessa (> 5000 m) de carbonatos mesozoicos-terciários de águas rasas (unidade Alburno-Cervati) sobreposição tectônica (ver duas sinergias marcadas) pela nappa da Sicilida e pelos depósitos de topo da folha de impulso que cobrem de maneira inconformável o Nappe da Ligúria. A embalagem de refletores fortes que se aprofunda de 2 s (extremidade direita) para 4 s (extremidade esquerda) foi interpretada em Figura 2 como uma pilha de depósitos basais intercalados entre a nuca Alburno-Cervati e o sistema duplex duplex das Apúlia carbonatos. Observe o pacote de refletores fracos de escalada do nordeste no lado direito da seção entre 6 e 8 s; o geológico A interpretação desses refletores (estruturas tectônicas envolvidas no transporte orogênico de acordo com a nossa interpretação) tem forte influenciou a reconstrução geométrica geral da correia de pressão esboçada na Figura 2.
No entanto, registros adequados para interpretação geralmente não excede 5-6 s em linhas comerciais, enquanto contínua e bem estruturada os eventos são reconhecíveis na nova seção de pilha até 8 s (veja a Figura 16). Esta diferença é substancial para a reparação a profundidade do único impulso e para construir um equilíbrio seção geológica. A Figura 2 é uma seção transversal simplificada derivada a partir de uma interpretação preliminar da seção sísmica.
O ponto é representado pela presença de um pacote bem organizado de refletores sob as Montanhas Alburno a uma profundidade de 6-8 s (veja a extremidade direita da seção na Figura 16). A geometria desses reflectores sugere uma profundidade (20-25 km) único impulso da cadeia Apenninic na área. Trabalho contínuo, incluindo migração em profundidade do horizonte e balanceamento de seção, é sendo abordado para entender melhor a geometria real e o significado geológico desses eventos.
Conclusões
A organização deste projeto em duas fases distintas de análise de dados e processamento de dados, embora exijam uma quantidade de mão-de-obra e tempo de CPU, tem sido a base para a melhoria bem-sucedida da qualidade dos dados. Outra chave O fator tem sido a interação constante entre geofísicos, responsável pelo reprocessamento dos dados, e geólogos com um
conhecimento geológico preciso da área. Esta cooperação tem foi particularmente frutífero na fase de análise de dados, na verificação a consistência das velocidades do refractor com a geologia da superfície, e na análise / interpretação da velocidade. Nenhuma fase única dessa seqüência de processamento pode ser considerada A chave para o sucesso deste projeto. Em vez disso, o todo processo produziu uma melhoria significativa nos resultados finais.
O projeto levou cerca de um ano e foi realizado por um equipe que incluía geofísicos seniores e geólogos e dois bolsistas de pesquisa em tempo integral. As fases mais exigentes pelo tempo O projeto foi a análise preliminar da qualidade dos dados, a computação estática e a análise da velocidade. Este trabalho não só melhorou a qualidade e a interpretabilidade deste conjunto de dados específicos, mas seus resultados também podem servir como base para a futura aquisição e processamento de sísmicas dados em áreas próximas ou em regiões com características semelhantes.