BRPI0803719B1 - Método para identificar uma posição de uma fonte de ruído em uma gravação sísmica marinha - Google Patents

Método para identificar uma posição de uma fonte de ruído em uma gravação sísmica marinha Download PDF

Info

Publication number
BRPI0803719B1
BRPI0803719B1 BRPI0803719-1A BRPI0803719A BRPI0803719B1 BR PI0803719 B1 BRPI0803719 B1 BR PI0803719B1 BR PI0803719 A BRPI0803719 A BR PI0803719A BR PI0803719 B1 BRPI0803719 B1 BR PI0803719B1
Authority
BR
Brazil
Prior art keywords
noise
source
signals
noise source
recording
Prior art date
Application number
BRPI0803719-1A
Other languages
English (en)
Inventor
Jack Dewayne Kinkead
Original Assignee
Pgs Geophysical As
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Pgs Geophysical As filed Critical Pgs Geophysical As
Publication of BRPI0803719A2 publication Critical patent/BRPI0803719A2/pt
Publication of BRPI0803719B1 publication Critical patent/BRPI0803719B1/pt

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Environmental & Geological Engineering (AREA)
  • Acoustics & Sound (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Oceanography (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

método para supressão de ruído em sinais sísmicos utilizando transformações espaciais. a presente invenção refere-se a um método para identificar uma posição de uma fonte de ruído em uma gravação sísmica marinha que inclui definir, para pelo menos uma gravação de explosão, um conjunto de posições possíveis de fonte de ruído. uma diferença entre o tempo de percurso do ruído a partir de cada posição possível de fonte de ruído até cada uma dentre várias posições de receptor para a pelo menos uma gravação de explosão é determinada. os sinais a partir de pelo menos um subconjunto de posições de receptores são alinhados no tempo com relação à diferença entre os tempos de percurso para cada posição possível de fonte de ruído. os sinais alinhados no tempo são então superpostos. a posição da fonte de ruído é determinada a partir dos sinais superpostos alinhados no tempo. este resultado pode então ser utilizado para construir um modelo de ruído na posição do receptor.

Description

PARA IDENTIFICAR UMA POSIÇÃO DE UMA FONTE DE RUÍDO EM
UMA GRAVAÇÃO SÍSMICA MARINHA.
Antecedentes da Invenção
Campo da Invenção
A presente invenção refere-se em geral ao campo de aquisição e processamento de dados sísmicos. Mais especificamente, a invenção se refere ao método para processar sinais sísmicos para atenuar os efeitos de certos tipos de ruído nos sinais sísmicos.
Técnica de Anterior
A agrimensura sísmica é conhecida na técnica para determinar estruturas de formações rochosas abaixo da superfície da terra. A agrimensura sísmica geralmente inclui dispor uma série de sensores sísmicos na superfície da terra em um padrão selecionado, e seletivamente atuar uma fonte de energia sísmica posicionada próxima dos sensores sísmicos. A fonte de energia pode ser um explosivo, um vibrador, ou no caso de agrimensura sísmica executada em ambientes marinhos, um ou mais canhões de ar ou canhões de água.
A energia sísmica que emana a partir da fonte percorre através das formações de subsuperfície da terra até que ela alcance um limite de impedância acústica nas formações. Os limites de impedância acústica tipicamente ocorrem onde a composição e/ou as propriedades mecânicas da formação terrestre se alteram. Tais limites tipicamente são referidos como interfaces. Em uma interface, alguma energia sísmica é refletida de volta para a superfície da terra. A energia refletida pode ser detectada pelos um ou mais sensores sísmicos dispostos na superfície. O processamento de sinal sísmico conhecido na técnica possui como um, dentre vários objetivos,
Petição 870190010080, de 30/01/2019, pág. 4/13
a determinação das profundidades e das localizações geográficas de interfaces abaixo da superfície da terra. A profundidade e a localização das interfaces são inferidas a partir do tempo de percursos da energia sísmica até as interfaces e de volta até os sensores na superfície.
A agrimensura sísmica é executada no oceano e em outros corpos de água (agrimensura sísmica marinha) para determinar a estrutura da formação terrestre abaixo do leito do mar. Os sistemas de agrimensura sísmica marinha conhecidos na técnica incluem uma embarcação de agrimensura sísmica que reboca uma ou mais fontes de energia sísmica, e a mesma 10 ou uma embarcação de agrimensura diferente que reboca uma ou mais fitas. As fitas são arranjos de receptores ou sensores sísmicos dispostos ao longo de um cabo que são rebocados pela embarcação. Tipicamente, uma embarcação sísmica irá rebocar várias fitas dispostas para ficarem separadas por uma distância lateral selecionada uma da outra, em um padrão sele15 cionado para permitir a determinação relativamente completa de estruturas geológicas da subsuperfície em três dimensões. Também é conhecido na técnica colocar cabos possuindo sensores sísmicos (cabos de leito de oceano) ao longo do leito marinho, e atuar uma fonte de energia sísmica na água. Tipicamente, a fonte de energia sísmica será rebocada por uma em20 barcação tal como na agrimensura do tipo fita.
Nas interfaces, como anteriormente explicado, alguma energia a partir da fonte é refletida e em última análise detectada pelos sensores sísmicos. Em adição à energia sísmica refletida, tanto ruído coerente como ruído incoerente podem estar presentes na energia sísmica detectada. A pre25 sença de ruído na energia detectada pelos sensores sísmicos reduz a relação de sinal para ruído (SNR) dos sinais sísmicos de interesse. Portanto, um objetivo dos especialistas em sismologia é buscar métodos para reduzir os efeitos do ruído sobre os sinais detectados pelos sensores sem reduzir de forma apreciável os componentes verdadeiros do sinal sísmico dos sinais 30 detectados.
Os métodos da técnica anterior que têm sido utilizados para reduzir os efeitos do ruído e adquirir uma representação de melhor qualidade
Vi de uma estrutura particular de subsuperfície incluem utilizar várias atuações da fonte sísmica (vários disparos ou explosões) para gravar várias medidas de sensor a partir substancialmente da mesma estrutura de subsuperfície, e então somar ou sobrepor tais medições para acentuar a potência do 5 sinal enquanto substancialmente reduzindo os efeitos de ruído aleatório ou incoerente.
A Patente U.S. 5.818.795 que é designada para um associado do cessionário da presente invenção e que proporciona um resumo detalhado dos métodos da técnica anterior e dos sistemas endereçando o problema 10 de supressão de ruído em sinais sísmicos descreve um método para reduzir o efeito de um ruído de rajada em gravações de sinal sísmico sem eliminar sinais sísmicos de interesse.
A Patente U.S. 5.761.152, que é designada para um associado do cessionário da presente invenção, descreve um método e sistema para 15 agrimensura sísmica marinha. O método descrito na patente '152 inclui aumentar a dobra (número de reflexões gravadas a partir do mesmo refletor), e como resultado a relação de sinal para ruído de sinais sísmicos coerentes, sem incorrer nos problemas de arrasto, entrelaçamento, manipulação complicada de convés associada com comprimento aumentado da fita, número 20 aumentado de fitas, e distância aumentada entre as fitas. Os deslocamentos da fonte e da fita e o tempo de disparo das fontes dianteiras e traseiras da embarcação em uma seqüência de retardo de tempo são otimizados para aumentar a dobra enquanto evitando influência substancial pelos sinais sísmicos resultantes da fonte de uma embarcação em relação aos sinais sísmi25 cos resultantes da fonte da outra embarcação.
A Patente U.S. 6.751.999 expedida para Fookes e outros, a qual é designada para um associado do cessionário da presente invenção, descreve um método para atenuar o ruído a partir de sinais sísmicos marinhos causados por um ruído na água. O método inclui determinar um tempo de 30 chegada de um evento de ruído em cada um dos vários sensores sísmicos, estimar uma posição da fonte de ruído a partir dos tempos de chegada, e atenuar o evento de ruído a partir dos sinais detectados pelos sensores sís
micos.
A descrição precedente não pretende ser uma explicação exaustiva dos tipos de ruídos e dos métodos para reduzir os efeitos dos mesmos em sinais sísmicos. Existem dois tipos de ruído, entretanto, que continuam a 5 apresentar uma necessidade de métodos aperfeiçoados de processamento. Um destes tipos de ruído é gerado por embarcações ou plataformas na vizinhança da embarcação sísmica (aqui, denominado ruído contínuo de navio). O outro tipo de ruído é ruído impulsivo refletido que se origina como resultado da atuação da fonte de energia sísmica e da subseqüente reflexão 10 da energia sísmica percorrendo lateralmente através da água até um refletor na água e de volta para os sensores no arranho (aqui, denominado ruído de retorno dispersado). Continua a existir uma necessidade por métodos aperfeiçoados para atenuar o ruído de retorno dispersado e o ruído contínuo de navio nos dados sísmicos marinhos.
Sumário da Invenção
Um método para identificar uma posição de uma fonte de ruído em uma gravação sísmica marinha inclui definir, para pelo menos uma gravação de explosão, um conjunto de posições possíveis de fonte de ruído. Uma diferença entre o tempo de percurso de ruído a partir de cada posição 20 possível de fonte de ruído até cada uma dentre várias posições receptoras para pelo menos uma gravação de explosão é determinada. Os sinais a partir de pelo menos um subconjunto das posições receptoras são alinhadas no tempo com respeito à diferença entre os tempos de percurso para cada posição possível de fonte de ruído. Os sinais alinhados no tempo são então 25 superpostos. A posição da fonte de ruído é determinada a partir dos sinais superpostos alinhados no tempo.
Um método para agrimensura sísmica marinha de acordo com outro aspecto da invenção inclui rebocar uma fonte de energia sísmica e várias fitas sísmicas em um corpo de água. Em tempos selecionados, a fonte é 30 atuada e os sinais detectados por cada um dos vários receptores em cada fita são gravados. Para cada atuação, uma posição geodésica da fonte e uma posição geodésica de cada receptor, são determinadas. Para cada atu
ação da fonte, um conjunto de posições possíveis de fonte de ruído com respeito à posição da fonte em cada atuação, é definido. Uma diferença entre o tempo de percursos do ruído a partir de cada posição possível de fonte de ruído até cada posição de receptor para cada gravação de explosão é 5 determinada. Os sinais a partir de cada um dos vários subconjuntos das posições de receptor com relação à diferença entre os tempos de percurso para cada posição possível de fonte de ruído para cada gravação de explosão são determinados. Os sinais alinhados no tempo para cada gravação de explosão são superpostos. A posição da fonte de ruído para cada gravação de 10 explosão é determinada a partir dos sinais superpostos.
Outros aspectos e vantagens da invenção serão aparentes a partir da descrição seguinte e das reivindicações anexas.
Breve Descrição dos Desenhos
A figura 1 apresenta um sistema ilustrativo de aquisição sísmica 15 como relacionado a um sistema de coordenadas de fonte de ruído de acordo com a invenção.
A figura 2 apresenta duas das fitas a partir do sistema da figura em maiores detalhes para ilustrar um cálculo de diferença de tempo de percurso de ruído de acordo com um exemplo de um método de acordo com 20 a invenção.
A figura 3 é um fluxograma de uma implementação ilustrativa de um método de acordo com a invenção.
As figuras 4A até 4C apresentam exemplos de dados sísmicos processados de acordo com uma concretização de um método de acordo 25 com a invenção.
As figuras 5A até 5C apresentam exemplos de dados sísmicos processados de acordo com uma concretização de um método de acordo com a invenção.
Descrição Detalhada
Um sistema de aquisição sísmica marinha ilustrativo que pode adquirir dados sísmicos para uso com os métodos de acordo com vários aspectos da invenção é apresentado esquematicamente na figura 1. Uma em-
barcação sísmica 10 é apresentada rebocando várias fitas lateralmente separadas S1, S2, S3, S4 através de um corpo de água 11 tal como um lago ou oceano. A embarcação sísmica 10 inclui equipamento na mesma apresentado geralmente por 12 e referido por conveniência como um sistema de gravação que inclui dispositivos (nenhum dos seguintes apresentados separadamente) para navegação, controle de fonte de energia sísmica, temporização e gravação de sinais sísmicos detectados por cada um dos vários receptores sísmicos longitudinalmente separados 16 ao longo de cada fita S1, S2, S3, S4. O sistema de gravação 12 nos tempos selecionados causa a atuação de uma fonte de energia sísmica 14, tal como um canhão de ar ou uma série de tais canhões de ar. A fonte de energia sísmica 14 é apresentada sendo rebocada pela embarcação sísmica 10, mas como será apreciado pelos versados na técnica, a embarcação que realmente reboca a fonte de energia sísmica pode ser uma embarcação diferente. Adicionalmente, o número de tais fontes ou de arranjos das mesmas e o número de fitas e seu espaçamento lateral não são limites em relação ao escopo desta invenção.
Os receptores sísmicos 16 detectam energia sísmica refletida pelos limites de impedância acústica na formação abaixo do leitor do corpo de água. Os sinais gerados pelos receptores 16 em resposta a detectar a energia sísmica refletida são sinais de interesse ao se interpretar a estrutura e a composição das formações abaixo do leito do mar.
Ao gravar tais sinais sísmicos, na hora de cada atuação da fonte de energia sísmica 14, o sistema de gravação inicia a geração de uma gravação indexada por tempo (tipicamente com respeito ao tempo de atuação da fonte) dos sinais produzidos por cada receptor 16. A posição geodésica da fonte 14 e de cada receptor 16 na hora de cada atuação de fonte, também é indexada com as gravações de sinal. A informação anterior gravada para cada atuação de fonte individual é conhecida como uma gravação de explosão.
Podem existir uma ou mais fontes de ruído na água 11, por exemplo, outras embarcações, instalações de produção de óleo e gás e outros dispositivos, que produzem energia acústica que também é detectada pelos receptores 16. Tal ruído, anteriormente identificado como ruído contínuo de navio e como ruído de retorno dispersado, pode interferir com a análise dos sinais sísmicos detectados. Em um método de acordo com a invenção, tal ruído pode ser identificado nos sinais sísmicos gravados e atenuado por um processo executado como dito a seguir.
Um conjunto de localizações possíveis ou a posição para a fonte de ruído, para cada gravação de explosão, apresentado como Ni até Nn na figura 1 pode ser predefinido. Em um exemplo, as localizações predefinidas Ni até Nn podem ser definidas nos círculos, cada círculo possuindo um raio selecionado, r, a partir de uma origem predefinida O. No presente exemplo, a origem predefinida O pode ser a posição geodésica da fonte de energia sísmica 14 na hora da atuação. A origem real utilizada em qualquer exemplo é uma questão de escolha para o projetista do sistema e não é pretendida para limitar o escopo da invenção.
Um ângulo Θ que define a segunda coordenada das posições predefinidas de fonte de ruído ao longo de um círculo definido pode ser definido com respeito a qualquer referência selecionada. No presente exemplo, a referência selecionada R pode ser uma referência geodésica, tal como o Norte geográfico, ou seja, quando Θ = 0, a posição da fonte de ruído está a uma distância r da fonte de energia sísmica 14 em uma direção Norte da posição da fonte de energia sísmica 14 (e da origem O). A referência R também pode, por exemplo, ser a direção geodésica de percurso da embarcação 10 ou qualquer outra referência que permita a definição azimutal das posições possíveis de fonte de ruído com respeito a uma direção. Para cada gravação de explosão, um conjunto de tais posições possíveis de fonte de ruído (apresentado ao longo de um único círculo de raio r como até Nn) pode ser definido ao longo dos círculos utilizando um incremento selecionado do ângulo de referência Θ, por exemplo, cinco ou dez graus, e um incremento selecionado do raio r.
Referindo-se à figura 2, duas das fitas S3 e S4 são apresentadas em maiores detalhes com uma das posições possíveis definidas de fonte de ruído Nx para ilustrar um cálculo utilizado em um método de acordo com a invenção. Para cada posição possível definida de fonte de ruído Nx, existe uma distância particular entre tal posição definida de fonte de ruído Nx e cada receptor sísmico 16. Tais distâncias são apresentadas na figura 2 como d! até dw- Uma quantidade de tempo para o ruído percorrer a partir de cada posição possível definida de fonte de ruído Nx até cada receptor sísmico 16 será relacionada com a distância (d-! até dw) entre cada receptor e cada posição possível de fonte de ruído Nx e a velocidade sísmica do ruído. Assim, para cada receptor 16, é possível definir um tempo de chegada esperado do ruído para cada posição possível de fonte de ruído Nx. Também é possível definir uma diferença entre os tempos esperados de chegada de ruído em cada receptor com respeito a um receptor de referência (ou qualquer outro receptor). No exemplo apresentado na figura 2, se o receptor de referência for apresentado em r5, então uma distância d5 entre o receptor de referência e a posição possível de fonte de ruído Nx pode ser calculada como a seguir:
d5 = [(x5 - XN )2 + (?5 - yN )2]'/2 onde x e y representam as coordenadas geodésicas para cada um dos receptores de referência r5 e a posição possível da fonte de ruído Nx, respectivamente, em um sistema Cartesiano de coordenadas, definido com relação à origem (O na figura 1). As coordenadas geodésicas de cada um dos receptores 16 e da fonte (14 na figurai) se utilizada como a origem O, podem ser determinadas por certo equipamento no sistema de gravação (12 na figura 1) para cada gravação de explosão. As coordenadas geodésicas de cada posição possível de ruído Nx podem ser definidas ao longo de círculos, como anteriormente explicado. As coordenadas de cada posição possível de ruído são definidas pelas coordenadas (rx, θχ) como anteriormente explicado. Um cálculo similar para a distância entre a posição possível da fonte de ruído Nx e cada um dos outros receptores 16 pode ser feito para cada gravação de explosão, e para cada posição possível definida de fonte de ruído. As coordenadas Cartesianas com relação à origem O de cada posição possível definida de ruído (definida acima utilizando coordenadas polares) podem ser calculadas pela expressão:
xN = rcosd yN = rsin# e as coordenadas Cartesianas anteriores para cada uma das posições definidas de fonte de ruído Nx podem então ser utilizadas para determinar as distâncias entre cada posição de receptor e cada posição possível definida de fonte de ruído para cada gravação de explosão.
Tendo assim determinada uma diferença entre as distâncias para cada posição possível de fonte de ruído para cada receptor com respeito a um receptor de referência, uma diferença entre os tempos de percurso do ruído a partir de cada uma das posições possíveis definidas de fonte de ruído Nx até cada receptor 16 pode ser calculada para cada receptor com respeito ao receptor de referência (o qual pode ser qualquer receptor selecionado no sistema de aquisição) utilizando a seguinte expressão:
Δζ = (ίζ -d,)/v0 onde di representa a distância entre a posição possível definida de fonte de ruído Nx e o receptor de referência, e dq representa a distância entre o q-ésimo receptor e a posição possível definida de fonte de ruído Nx. v0 representa a velocidade sísmica do ruído. A velocidade sísmica do ruído para agrimensuras sísmicas marinhas tipicamente será igual à velocidade sísmica da água, apesar da velocidade do ruído poder depender da profundidade da água, da distância total até a fonte de ruído a partir dos receptores e da velocidade sísmica das formações próximas do leito marinho. Em alguns exemplos, a velocidade do ruído v0 pode ser ajustada e os cálculos anteriores de diferença de tempo repetidos para a velocidade ajustada. O conjunto original e os cálculos de diferença de tempo repetidos, podem ser utilizados como explicado abaixo.
As diferenças de tempo de percurso de ruído anteriores podem ser utilizadas para alinhar no tempo os sinais gravados a partir de cada um dos receptores. O alinhamento no tempo pode ser executado pelo deslocamento no tempo de forma linear dos sinais gravados a partir de cada um dos receptores de modo que o ruído, se for originado a partir de uma posição predefinida de fonte de ruído Nx, estaria presente nos sinais de todos os receptores substancialmente ao mesmo tempo. O dito anteriormente é repetido para todas as posições possíveis definidas de fonte de ruído Nx.
Os sinais sísmicos alinhados no tempo para cada posição possível definida de fonte de ruído podem então ser somados ou superpostos. A superposição pode ser executada para todos os sinais alinhados no tempo a partir de todos os receptores no sistema de aquisição sísmica, ou pode em alguns exemplos ser executado para cada um dentre vários subconjuntos selecionados de todos os receptores, apesar de que o receptor de referência utilizado para calcular as diferenças de tempo será o mesmo. Em um exemplo, os sinais alinhados no tempo a partir de cada fita podem ser somados ou superpostos. O resultado é uma gravação de sinal superposta para cada fita relacionada com cada posição predefinida possível de fonte de ruído Nx.
Os sinais superpostos (rastos) para cada posição predefinida de fonte de ruído, para cada fita, então podem ser avaliados para determinar uma posição mais provável de fonte de ruído. Para a posição predefinida de fonte de ruído Nx que mais provavelmente representa a posição real da fonte de ruído, um evento de sinal de alta amplitude pode ser observado nos rastos superpostos correspondendo à posição mais provável de fonte de ruído. Em um exemplo, os sinais superpostos alinhados no tempo são avaliados pela varredura com respeito ao ângulo azimutal Θ para cada raio r de cada posição possível de fonte de ruído Nx. Um benefício possível de superpor sinais alinhados no tempo em subconjuntos, tal como todos os receptores nas fitas individuais como explicado acima é que em alguns caos, mais do que uma posição possível de fonte de ruído irá indicar um evento de alta amplitude nos sinais superpostos alinhados no tempo. Se os sinais superpostos alinhados no tempo forem processados em subconjuntos como explicado acima, é mais provável que somente uma posição possível de fonte de ruído (para cada fonte de ruído real) irá indicar o mesmo evento de alta amplitude em cada um dos vários subconjuntos de sinais superpostos alinhados no tempo. Por exemplo, os sinais superpostos alinhados no tempo a partir de cada fita, processados como explicado acima, podem então ser filtrados
pela média para identificar a única mais provável posição de fonte de ruído.
A figura 3 apresenta um fluxograma de um processo ilustrativo como explicado acima. Para cada gravação de explosão, em 30, um conjunto de posições possíveis de fonte de ruído é definido. Um tempo de percurso diferencial de ruído a partir de cada uma das posições de fonte de ruído para cada receptor, é calculado, para cada gravação de explosão, em 32. Os sinais são alinhados no tempo em 34. Os sinais alinhados no tempo são superpostos para cada posição possível de fonte de ruído em 36. A posição possível de fonte de ruído N possuindo a energia máxima é então selecionada para cada subconjunto de receptores em 38. Os procedimentos 32, 34, 36 e 38 acima são repetidos para cada subconjunto de receptores (por exemplo, para cada fita). Os resultados a partir de todas as fitas (subconjuntos de receptores) podem ser alinhados pelo azimute comum (ângulo Θ) das posições possíveis de fonte de ruído e então filtrados pela média para determinar uma ou mais posições de fonte de ruído em 44.
Para atenuar os efeitos do ruído nos sinais gravados, um tempo de chegada esperado para o ruído, com respeito à posição determinada de fonte de ruído e às posições de receptor para cada gravação de explosão para cada receptor é então determinado, em 46. O evento de ruído é então de forma adaptável subtraído do sinal gravado de cada receptor no tempo de chegada de ruído calculado, em 48. O resultado são dados sísmicos possuindo efeito reduzido de ruído coerente, por exemplo, a partir das fontes na água tal como outras embarcações sísmicas e plataformas de perfuração.
Exemplos de dados sísmicos processados como explicado acima são apresentados nas figuras 4A, 4B e 4C. A figura 4A apresenta sinais a partir de uma fita em um arranjo de fitas (tal como apresentado na figura 1) que são alinhados no tempo e superpostos como explicado acima. A exibição na figura 4A inclui vários rastos superpostos alinhados no tempo lado a lado para cada direção azimute definida da posição possível de fonte de ruído ao longo de um círculo (definido como explicado acima). Em 50, um evento de alta amplitude pode ser observado nos rastos superpostos alinhados no tempo para uma série de azimutes adjacentes definidos. O evento de
3,3 alta amplitude 50 indica alguma incerteza no cálculo de azimute da posição da fonte de ruído. Uma exibição similar para os sinais sísmicos obtidos a partir de uma fita diferente no mesmo sistema de aquisição é apresentada na figura 4B. Um evento de alta amplitude 52 nos sinais superpostos alinhados no tempo na figura 4B também apresenta alguma incerteza quanto ao azimute da posição da fonte de ruído. A figura 4C apresenta os sinais superpostos alinhados no tempo a partir das figuras 4A e 4B após serem filtrados pela média. Observe que um evento de alta amplitude 54 apresenta consideravelmente menos incerteza do azimute quando à posição da fonte de ruído.
Exibições similares de sinais superpostos alinhados no tempo são apresentadas na figura 5A, onde para uma das fitas em um arranjo, um primeiro evento de alta amplitude 56 é observável ao longo de um primeiro azimute 56 e em um segundo evento de alta amplitude 58 é observável ao longo de um segundo azimute. Os dois eventos de alta amplitude 56, 58 podem representar a mesma posição possível de fonte de ruído, mas em algumas circunstâncias, duas posições possíveis de fonte de ruído irão resultar em um evento de alta amplitude nos rastos superpostos alinhados no tempo. Uma exibição similar para traços a partir de outra fita no mesmo arranjo é apresentada na figura 5B, onde eventos de alta amplitude 60, 62 são observáveis ao longo de mais do que um azimute. Pela filtragem pela média dos rastos a partir da figura 5A e 5B, os rastos superpostos alinhados no tempo filtrados pela média irão produzir uma exibição que possui somente um evento de alta amplitude 64. Assim, a filtragem pela média dos rastos a partir de diferentes subconjuntos de receptores (por exemplos, diferentes fitas) pode levar a uma solução única para a posição da fonte de ruído.
Enquanto a invenção foi descrita com relação a um número limitado de concretizações, os versados na técnica, tendo o benefício desta revelação, irão apreciar que outras concretizações podem ser planejadas, as quais não se afastem do escopo da invenção como descrito neste documento. Por conseqüência, o escopo da invenção somente deve ser limitado pelas reivindicações anexas.

Claims (20)

  1. REIVINDICAÇÕES
    1. Método para identificar uma posição de uma fonte de ruído em uma gravação sísmica marinha, caracterizado pelo fato de que compreende:
    definir para pelo menos uma gravação de explosão um conjunto de posições possíveis de fonte de ruído (Nn), a referida pelo menos uma gravação de explosão compreendendo gravações indexadas por tempo de atuação de uma fonte de energia sísmica (14) de sinais detectados por cada um de uma pluralidade de receptores sísmicos (16) lançados em um corpo de água, os sinais detectados em resposta à atuação da fonte de energia sísmica (14);
    calcular uma diferença (32) entre o tempo de percurso (At) de ruído a partir de cada posição possível de fonte de ruído (Nn) para cada uma dentre uma pluralidade de posições de receptor (16) para a referida pelo menos uma gravação de explosão;
    alinhar no tempo sinais (34) a partir de pelo menos um subconjunto de posições de receptor (16) com relação à diferença entre os tempos de percurso para cada posição possível de fonte de ruído (Nn);
    superpor os sinais alinhados no tempo (36); e determinar a posição da fonte de ruído (44) a partir dos sinais superpostos.
  2. 2. Método de acordo com a reivindicação 1, caracterizado pelo fato de que a diferença entre os tempos de percurso é determinada pelo cálculo de uma distância (d) entre cada posição possível de fonte de ruído (Nn) e cada receptor no pelo menos um subconjunto.
  3. 3. Método de acordo com a reivindicação 2, caracterizado pelo fato de que adicionalmente compreende ajustar a velocidade do ruído (38) e repetir (40) o cálculo de diferença, o alinhamento no tempo, a superposição e a determinação da posição da fonte de ruído.
  4. 4. Método de acordo com a reivindicação 1, caracterizado pelo fato de que o conjunto de posições possíveis de fonte de ruído é definido ao longo de pelo menos um círculo possuindo um raio selecionado (r) a partir
    Petição 870190010080, de 30/01/2019, pág. 5/13 de uma origem (O).
  5. 5. Método de acordo com a reivindicação 1, caracterizado pelo fato de que a origem (O) é definida como uma posição de uma fonte de energia sísmica (14) para pelo menos uma gravação de explosão.
  6. 6. Método de acordo com a reivindicação 1, caracterizado pelo fato de que o pelo menos um subconjunto compreende receptores ao longo de uma primeira fita sísmica selecionada (S1) em um sistema de aquisição sísmica possuindo uma pluralidade fitas (S1-S4).
  7. 7. Método de acordo com a reivindicação 6, caracterizado pelo fato de que o cálculo da diferença, o alinhamento no tempo, a superposição e a determinação da posição da fonte de ruído são repetidos (40) para pelo menos um subconjunto adicional de receptores.
  8. 8. Método de acordo com a reivindicação 7, caracterizado pelo fato de que pelo menos um subconjunto adicional compreende uma segunda fita selecionada (S2) no sistema de aquisição sísmica.
  9. 9. Método de acordo com a reivindicação 7, caracterizado pelo fato de que adicionalmente compreende alinhar os sinais superpostos a partir da primeira fita (S1) e os sinais superpostos a partir da segunda fita (S2) com relação a um azimute comum para cada posição possível de fonte de ruído e filtrar pela média (42) os sinais superpostos azimutalmente alinhados para identificar a posição da fonte de ruído.
  10. 10. Método de acordo com a reivindicação 1, caracterizado pelo fato de que adicionalmente compreende determinar um tempo de chegada esperado de um evento de ruído (46) em sinais a partir de cada receptor a partir da posição determinada da fonte de ruído e de forma adaptável subtrair (48) ruído a partir dos sinais a partir de cada receptor sísmico baseado no tempo de percurso esperado do ruído até cada receptor.
  11. 11. Método de acordo com a reivindicação 1, caracterizado pelo fato de ainda compreende:
    rebocar uma fonte de energia sísmica (14) e uma pluralidade de fitas sísmicas (S1-S4) em um corpo de água;
    em tempos selecionados, atuar a fonte (14) e gravar sinais dePetição 870190010080, de 30/01/2019, pág. 6/13 tectados por cada um dentre uma pluralidade de receptores (16) em cada fita (S1-S4);
    determinar para cada atuação, uma posição geodésica da fonte (14) e uma posição geodésica de cada receptor (16) e executar;
    a definição para cada atuação de fonte, o conjunto de posições possíveis de fonte de ruído (Nn) com relação à posição da fonte (O) em cada atuação, o cálculo da diferença (32) entre o tempo de percurso do ruído a partir de cada posição possível de fonte de ruído (Nn) até cada posição de receptor (16) para cada gravação de explosão, o alinhamento no tempo (34) os sinais a partir de cada um da pluralidade de subconjuntos das posições de receptor com relação à diferença entre os tempos de percurso para cada posição possível de fonte de ruído para cada gravação de explosão, e a superporsição (36) os sinaisalinhados no tempo (36) para cada gravação de explosão; e a determinação da posição de fonte de ruído (44) para cada gravação de explosão a partir dos sinais superpostos.
  12. 12. Método de acordo com a reivindicação 11, caracterizado pelo fato de que a diferença entre os tempos de percurso é determinada pelo cálculo de uma distância (d) entre cada posição possível de fonte de ruído e cada receptor em cada subconjunto.
  13. 13. Método de acordo com a reivindicação 12, caracterizado pelo fato de que adicionalmente compreende ajustar a velocidade do ruído (38) e repetir (40) o cálculo de diferença, o alinhamento no tempo, a superposição e a determinação da posição da fonte de ruído para cada gravação de explosão.
  14. 14. Método de acordo com a reivindicação 11, caracterizado pelo fato de que o conjunto de posições possíveis de fonte de ruído é definido ao longo de pelo menos um círculo possuindo um raio selecionado (r) a partir de uma origem (O) para cada gravação de explosão.
  15. 15. Método de acordo com a reivindicação 11, caracterizado
    Petição 870190010080, de 30/01/2019, pág. 7/13
    4 pelo fato de que a origem (O) é definida como a posição da fonte de energia sísmica (14) para cada gravação de explosão.
  16. 16. Método de acordo com a reivindicação 11, caracterizado pelo fato de que cada subconjunto compreende os receptores (16) em uma fita selecionada das fitas (S1-S4).
  17. 17. Método de acordo com a reivindicação 16, caracterizado pelo fato de que o cálculo da diferença, o alinhamento no tempo, a superposição e a determinação da posição da fonte de ruído são repetidos (40) para pelo menos um subconjunto adicional de receptores.
  18. 18. Método de acordo com a reivindicação 17, caracterizado pelo fato de que o pelo menos um subconjunto adicional compreende uma segunda fita (S2) das fitas.
  19. 19. Método de acordo com a reivindicação 18, caracterizado pelo fato de que adicionalmente compreende alinhar os sinais superpostos a partir da primeira fita (S1) e os sinais superpostos a partir da segunda fita (S2) para cada gravação de explosão com relação a um azimute comum para cada posição possível de fonte de ruído e filtrar pela média (42) os sinais superpostos azimutalmente alinhados para identificar a posição da fonte de ruído para cada gravação de explosão.
  20. 20. Método de acordo com a reivindicação 11, caracterizado pelo fato de que adicionalmente compreende determinar um tempo de chegada esperado de um evento de ruído (46) em sinais a partir de cada receptor a partir da posição determinada da fonte de ruído para cada gravação de explosão e de forma adaptável subtrair (48) ruído a partir dos sinais a partir de cada receptor sísmico baseado no tempo de percurso esperado do ruído até cada receptor para cada gravação de explosão.
BRPI0803719-1A 2007-08-14 2008-08-14 Método para identificar uma posição de uma fonte de ruído em uma gravação sísmica marinha BRPI0803719B1 (pt)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US11/893,032 US7869303B2 (en) 2007-08-14 2007-08-14 Method for noise suppression in seismic signals using spatial transforms
US11/893,032 2007-08-14

Publications (2)

Publication Number Publication Date
BRPI0803719A2 BRPI0803719A2 (pt) 2009-05-19
BRPI0803719B1 true BRPI0803719B1 (pt) 2019-05-21

Family

ID=39790545

Family Applications (1)

Application Number Title Priority Date Filing Date
BRPI0803719-1A BRPI0803719B1 (pt) 2007-08-14 2008-08-14 Método para identificar uma posição de uma fonte de ruído em uma gravação sísmica marinha

Country Status (8)

Country Link
US (1) US7869303B2 (pt)
AU (1) AU2008203152B2 (pt)
BR (1) BRPI0803719B1 (pt)
EG (1) EG25319A (pt)
GB (1) GB2451945B (pt)
MY (1) MY148207A (pt)
NO (1) NO341944B1 (pt)
SG (1) SG150432A1 (pt)

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9310503B2 (en) * 2009-10-23 2016-04-12 Westerngeco L.L.C. Methods to process seismic data contaminated by coherent energy radiated from more than one source
US8639442B2 (en) * 2010-11-23 2014-01-28 Westerngeco L.L.C. Identifying invalid seismic data
US9188693B2 (en) * 2012-03-26 2015-11-17 Apache Corporation Method for acquiring marine seismic data
FR2990028B1 (fr) * 2012-04-25 2014-05-16 Kietta Acquisition de donnees sismiques
US9329293B2 (en) 2013-03-12 2016-05-03 Pgs Geophysical, As Systems and methods for removing acquisition related effects from seismic data
US10338250B2 (en) 2013-03-14 2019-07-02 Pgs Geophysical As Method of removing incoherent noise
US10359528B2 (en) * 2013-03-15 2019-07-23 Pgs Geophysical As Systems and methods for randomizing firing times of simultaneous sources in marine surveys
US11391857B2 (en) 2013-12-30 2022-07-19 Pgs Geophysical As Methods and systems for attenuating residual acoustic energy in seismic data
US10598807B2 (en) 2014-02-18 2020-03-24 Pgs Geophysical As Correction of sea surface state
US9874646B2 (en) 2014-04-14 2018-01-23 Pgs Geophysical As Seismic data processing
US10132946B2 (en) 2014-08-13 2018-11-20 Pgs Geophysical As Methods and systems that combine wavefields associated with generalized source activation times and near-continuously recorded seismic data
US10317553B2 (en) 2014-08-13 2019-06-11 Pgs Geophysical As Methods and systems of wavefield separation applied to near-continuously recorded wavefields
CN105717543B (zh) * 2016-01-25 2018-07-13 浪潮(北京)电子信息产业有限公司 一种噪声压制方法及系统
US10267936B2 (en) 2016-04-19 2019-04-23 Pgs Geophysical As Estimating an earth response
FR3054890B1 (fr) 2016-08-02 2019-07-05 Kietta Controle de la position horizontale d’un cable sismique

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5241514A (en) 1992-05-06 1993-08-31 Amoco Corporation Identifying, characterizing, and mapping of seismic noise
US5448531A (en) * 1994-05-05 1995-09-05 Western Atlas International Method for attenuating coherent noise in marine seismic data
US5572483A (en) * 1995-07-25 1996-11-05 Western Atlas International, Inc. Method of reducing noise in seismic signals by adaptive filtering of a noise reference
US5761152A (en) 1996-10-29 1998-06-02 Pgs Exploration (Us), Inc. Method and system for increasing fold to streamer length ratio
US5818795A (en) 1996-10-30 1998-10-06 Pgs Tensor, Inc. Method of reduction of noise from seismic data traces
GB9726928D0 (en) * 1997-12-19 1998-02-18 Geco Prakla Uk Ltd Method of stacking seismic signals
US6751559B2 (en) 2002-09-10 2004-06-15 Pgs Exploration (Uk) Limited Method for suppressing noise from seismic signals by source position determination
US7197399B2 (en) * 2005-01-13 2007-03-27 Bp Corporation North America, Inc. Method of multiple attenuation

Also Published As

Publication number Publication date
EG25319A (en) 2011-12-12
AU2008203152A2 (en) 2012-03-22
GB0814590D0 (en) 2008-09-17
MY148207A (en) 2013-03-29
AU2008203152A1 (en) 2009-03-05
US7869303B2 (en) 2011-01-11
AU2008203152B2 (en) 2013-03-14
SG150432A1 (en) 2009-03-30
GB2451945A (en) 2009-02-18
NO20083162L (no) 2009-02-16
US20090046536A1 (en) 2009-02-19
NO341944B1 (no) 2018-02-26
GB2451945B (en) 2011-11-23
BRPI0803719A2 (pt) 2009-05-19

Similar Documents

Publication Publication Date Title
US7869303B2 (en) Method for noise suppression in seismic signals using spatial transforms
AU2003207039B2 (en) Method and system for acquiring marine seismic data using multiple seismic sources
US6751559B2 (en) Method for suppressing noise from seismic signals by source position determination
EP3324215B1 (en) Method of first arrival picking of multi-channel seismic survey data
US9188693B2 (en) Method for acquiring marine seismic data
US20050027454A1 (en) Method for separating seismic signals from two or more distinct sources
DK2326972T3 (en) DIMENSION OF SEISMIC INTERFERENCE NOISE USING A DUAL SENSOR RECORDING SYSTEM
US20110096625A1 (en) Methods to Process Seismic Data Contaminated By Coherent Energy Radiated From More Than One Source
BRPI1001848B1 (pt) método para adquirir e processar sinais sísmicos marinhos para extrair campos de onda ascendentes e descendentes de uma fonte de energia sísmica
GB2589011A (en) Modified simultaneous long-offset acquistion with improved low frequency performance for full wavefield inversion
BR102014015913A2 (pt) técnicas de pesquisa utilizando cordões em profundidades diferentes
BR102013019400A2 (pt) dispositivo e método para aquisição marinha sincronizada com reduzido ruído de interferência
BR102015013658A2 (pt) criação de imagem sísmica utilizando reflexos de ordem superior
US10962666B2 (en) Using seabed sensors and sea-surface reflections for structural imaging of a subsurface location in a geological formation
GB2494974A (en) Method of determining the relative position of two detectors at the bottom of the sea
EP2271953B1 (en) Using towed seismic surveys that do not have coinciding streamer positions in the time lapse analysis of a producing field
CN1290351A (zh) 地震数据获得方法
Minshull et al. Analysis of gas-hydrate provinces by ocean-bottom seismic methods
BR112018005726B1 (pt) Método e sistema para aquisição de um levantamento sísmico marinho, e, meio legível por computador não transitório

Legal Events

Date Code Title Description
B03A Publication of a patent application or of a certificate of addition of invention [chapter 3.1 patent gazette]
B03H Publication of an application: rectification [chapter 3.8 patent gazette]

Free format text: REFERENTE A RPI 2002 DE 19/05/2009, QUANTO AO ITEM (54).

B06A Patent application procedure suspended [chapter 6.1 patent gazette]
B09A Decision: intention to grant [chapter 9.1 patent gazette]
B16A Patent or certificate of addition of invention granted [chapter 16.1 patent gazette]

Free format text: PRAZO DE VALIDADE: 10 (DEZ) ANOS CONTADOS A PARTIR DE 21/05/2019, OBSERVADAS AS CONDICOES LEGAIS. (CO) 10 (DEZ) ANOS CONTADOS A PARTIR DE 21/05/2019, OBSERVADAS AS CONDICOES LEGAIS

B21F Lapse acc. art. 78, item iv - on non-payment of the annual fees in time

Free format text: REFERENTE A 13A ANUIDADE.

B24J Lapse because of non-payment of annual fees (definitively: art 78 iv lpi, resolution 113/2013 art. 12)

Free format text: EM VIRTUDE DA EXTINCAO PUBLICADA NA RPI 2632 DE 15-06-2021 E CONSIDERANDO AUSENCIA DE MANIFESTACAO DENTRO DOS PRAZOS LEGAIS, INFORMO QUE CABE SER MANTIDA A EXTINCAO DA PATENTE E SEUS CERTIFICADOS, CONFORME O DISPOSTO NO ARTIGO 12, DA RESOLUCAO 113/2013.