BRPI1103858A2 - removedor de fantasma de campo de onda de dados sísmicos registrados usando múltiplas fontes sísmicas em diferentes profundidades de água - Google Patents
removedor de fantasma de campo de onda de dados sísmicos registrados usando múltiplas fontes sísmicas em diferentes profundidades de água Download PDFInfo
- Publication number
- BRPI1103858A2 BRPI1103858A2 BRPI1103858A BRPI1103858A BRPI1103858A2 BR PI1103858 A2 BRPI1103858 A2 BR PI1103858A2 BR PI1103858 A BRPI1103858 A BR PI1103858A BR PI1103858 A BRPI1103858 A BR PI1103858A BR PI1103858 A2 BRPI1103858 A2 BR PI1103858A2
- Authority
- BR
- Brazil
- Prior art keywords
- sources
- seismic
- source
- additional
- wave field
- Prior art date
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 50
- 230000001960 triggered effect Effects 0.000 claims abstract description 20
- 230000000977 initiatory effect Effects 0.000 claims abstract description 13
- 230000004913 activation Effects 0.000 claims abstract description 3
- 238000000034 method Methods 0.000 claims description 65
- 230000000694 effects Effects 0.000 claims description 10
- 238000001914 filtration Methods 0.000 claims description 7
- 238000005314 correlation function Methods 0.000 claims description 5
- 238000001228 spectrum Methods 0.000 description 20
- 230000015572 biosynthetic process Effects 0.000 description 10
- 238000005755 formation reaction Methods 0.000 description 10
- 230000014509 gene expression Effects 0.000 description 10
- 239000011435 rock Substances 0.000 description 6
- 238000000354 decomposition reaction Methods 0.000 description 5
- 239000011159 matrix material Substances 0.000 description 5
- 238000012545 processing Methods 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 4
- 230000035508 accumulation Effects 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000001934 delay Effects 0.000 description 3
- 239000000203 mixture Substances 0.000 description 3
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000003213 activating effect Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 2
- 230000001066 destructive effect Effects 0.000 description 2
- 238000005553 drilling Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000010304 firing Methods 0.000 description 2
- 230000001976 improved effect Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000003199 nucleic acid amplification method Methods 0.000 description 2
- 230000001902 propagating effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 230000001131 transforming effect Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000002706 hydrostatic effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004181 pedogenesis Methods 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
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
- G01V1/364—Seismic filtering
-
- 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
- G01V1/3808—Seismic data acquisition, e.g. survey design
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/003—Seismic data acquisition in general, e.g. survey design
- G01V1/005—Seismic data acquisition in general, e.g. survey design with exploration systems emitting special signals, e.g. frequency swept signals, pulse sequences or slip sweep arrangements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/20—Trace signal pre-filtering to select, remove or transform specific events or signal components, i.e. trace-in/trace-out
- G01V2210/21—Frequency-domain filtering, e.g. band pass
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/56—De-ghosting; Reverberation compensation
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Oceanography (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
removedor de fantasma de campo de onda de dados sísmicos registrados usando múltiplas fontes sísmicas em diferentes profundidades de água dados sísmicos são adquiridos através de acionar uma primeira fonte em um primeiro tempo e uma ou mais fontes adicionais sísmicas cada uma com seus tempos característicos com relação a um tempo de registro de sinal, as fontes substancialmente colocadas em diferentes profundidades. um primeiro campo de onda é determinado para ocorrer se a primeira fonte for acionada em um tempo selecionado com relação tempo de iniciação dos registros e sendo o tempo ajustado para a profundidade de água. um ou mais campos de onda adicionais são determinados para ocorrer se uma ou mais fontes adicionais fossem cada uma acionadas no dito tempo selecionado com relação ao dito tempo de iniciação, e sendo o tempo ajustado para profundidades de água de uma ou mais fontes adicionais. o primeiro campo de onda e o ou os campos adicionais de onda são unidos para determinar uma fonte campo de onda sem fantasma, que corresponde ao acionamento de uma única fonte de energia sísmica.
Description
Relatório Descritivo da Patente de Invenção para "REMOVE-DOR DE FANTASMA DE CAMPO DE ONDA DE DADOS SÍSMICOS REGISTRADOS USANDO MÚLTIPLAS FONTES SÍSMICAS EM DIFERENTES PROFUNDIDADES DE ÁGUA".
Referência remissiva a aplicações relacionadas Não aplicável.
Declaração Relativa à pesquisa ou ao desenvolvimento com patrocínio Federal Não aplicável.
Antecedentes da invenção Campo Da Invenção A invenção diz respeito geralmente ao campo da aquisição e processamento de dados sísmicos marinhos. Mais especificamente, a invenção diz respeito a métodos para a redução do efeito de reflexões na superfície da água na energia emitida pelas fontes de energia sísmica marinha.
Antecedentes da técnica Em exploração sísmica, os dados sísmicos são adquiridos pela transmissão de energia acústica dentro da Terra próximo da sua superfície, e a detecção da energia acústica que é refletida dos limites entre diferentes camadas de formações rochosas de subsuperfície. A energia acústica é refletida quando há uma diferença na impedância acústica entre camadas como uma adjacente a um limite. Sinais representando a energia acústica detectada são interpretados para supor estruturas e a composição das estruturas de formação rochosa de subsuperfície.
Em exploração sísmica marinha, uma fonte de energia sísmica, como um canhão de ar, ou uma matriz de canhão de ar, é tipicamente usada para transmitir a energia acústica nas formações abaixo do fundo da água. O canhão de ar ou matriz é ativado em uma profundidade selecionada na água, tipicamente enquanto o canhão de ar ou matriz é rebocado por uma embarcação. A mesma ou uma embarcação diferente reboca um ou mais cabos de sensor sísmico, chamados “cabos flutuadores”, na água. Geralmente o cabo flutuador se estende atrás da embarcação junto com a direção em que o cabo flutuador é rebocado. Tipicamente, um cabo flutuador inclui uma pluralidade de hidrófonos dispostos no cabo afastados espaçadamente, em posições conhecidas junto com o cabo. Hidrófonos, conforme conhecidos na técnica, são sensores que geram um sinal ótico ou elétrico correspondendo à pressão da água ou o tempo gradual (dp/dt) de pressão na água. A embarcação que reboca um ou mais cabos flutuadores tipicamente inclui equipamento de registro para produzir um registro, desenvolvido com respeito ao tempo dos sinais gerados pelos hidrófonos em resposta a energia acústica detectada. O registro de sinais é processado, conforme previamente explicado, para supor estruturas de e composições das formações do solo abaixo dos lugares em que a análise sísmica é executada.
Dados sísmicos marinhos incluem um efeito que limita a exatidão da dedução sobre a estrutura e a composição das formações rochosas de subsuperfície. Esse efeito, conhecido como fonte de fantasma, aparece porque a água tem uma densidade substancialmente diferente e uma velocidade de propagação das ondas de pressão do que o ar acima da superfície da água. A fonte de fantasma pode ser entendida conforme se segue. Quando o canhão de ar ou a matriz de canhão de ar é ativado, a energia acústica irradia geralmente voltada para fora do canhão de ar ou matriz. Metade da energia viaja para baixo onde ela passa através do fundo da água e para o interior das formações rochosas de subsuperfície. A outra metade da energia acústica viaja para cima de um canhão ou matriz e maior parte dessa energia se reflete da superfície da água e, portanto ela viaja para baixo. A energia acústica refletida irá ser atrasada no tempo e também ser trocada em fase por cerca de 180 graus da energia acústica se propagando diretamente para baixo. A superfície-refletida, energia acústica viajando para baixo é comumente conhecida como um sinal “fantasma”. O sinal fantasma interfere com o campo de onda se propagando diretamente para baixo causando uma interferência construtiva em algumas partes da faixa de frequências e interferência destrutiva em outras partes da faixa de frequências. A interferência destrutiva causa a sequência de descontinuidades no espectro, igualmente espaçadas em frequência incluindo uma descontinuidade na fre- quência zero (0 Hz). As frequências dessas descontinuidades no sinal acústico detectado são relacionadas com a profundidade em que o canhão de ar ou matriz de canhão é disposto, conforme é bem conhecido na técnica. O efeito da fonte de fantasma é tipicamente referido como uma “fonte fantasma. A energia sísmica emitida pela fonte é atenuada com a distância de propagação por causa do espalhamento geométrico, perda de transmissão, e absorção. A absorção de energia de frequência superior em uma taxa maior que a energia de frequência inferior é bem conhecida na técnica. Então, para uma penetração profunda, é desejável maximizar a energia emitida por uma fonte em frequências menores. Desde que a fonte fantasma causa uma descontinuidade no sinal registrado a 0 Hz, o fantasma está limitando a energia na extremidade de baixa frequência do sinal. A quantia de energia nas baixas frequências pode ser melhorada pelo rebocamento das fontes em profundidades maiores. Contudo, isso causa com que as descontinuidades fantasmas no espectro ocorram em frequências menores, e, portanto limita as partes de alta frequência do espectro necessárias para uma imagem de alta resolução de alvos menos profundos. Também, quando usando um ou mais canhões de ar como uma fonte de energia sísmica, a frequência fundamental do(s) canhão(ões) aumenta com uma profundidade crescente. Portanto, o aumentar da energia na extremidade de baixa frequência quando rebocando o(s) canhão(ões) de ar mais profundamente devido à fonte fantasma, é contrariado pelo aumentar na frequência fundamental do(s) ca-nhão(ões) de ar.
Uma maneira conhecida de aumentar o do nível do sinal emitido pela fonte através a largura de banda quando usando canhão(ões) de ar é a de aumentar o volume total de ar liberado pelo(s) canhão(ões) de ar e/ou para aumentar a pressão operante. Contudo, o volume máximo de ar que pode ser liberado para todo disparo e a pressão máxima de ar é limitada pelo equipamento de fonte e sistema de fornecimento de ar disponível. Para trocar os volumes máximos e a pressão pode ser algo muito dispendioso e moroso. Também, aumentando a força da fonte pode ter um impacto na vida marinha. Então, maximizando o uso do sinal emitido pela fonte pode ser de grande valor e reduzir a necessidade de aumentar o nível de energia emitido pela fonte. Pela extração de uma ascendente (com fantasma) e os campos de ondas que se propagam diretamente para baixo vindos da fonte, os efeitos da fonte fantasma são eliminados e o sinal em torno de todas as descon-tinuidades fantasma é aumentado, incluindo o sinal nas frequências menores próximas a descontinuidade em 0 Hz. Esses campos de ondas separados podem também ser trocados em tempo para a superfície do mar ou uma referência comum de profundidade usando a já conhecida profundidade(s) da fonte. Portanto, pela aplicação de uma 180 mudança de fase para o sinal com fantasma, os campos de ondas separados podem ser sintetizados juntos construtivamente. Desta forma quase toda energia emitida pela fonte é utilizada. A técnica conhecida na técnica para a extração da fonte fantasma é descrita na obra M. Egan et al., Full deghosting of OBC data with o-ver/under source acquisition, 2007 Annual Meeting, San Antonio, TX, Society of Exploration Geophysicists. A técnica descrita na publicação de Eram ET al. inclui o rebocamento de uma primeira fonte de energia sísmica em uma primeira profundidade na água, e rebocamento uma segunda fonte de energia sísmica em uma segunda profundidade da água. As fontes são canhões de ar ou matrizes disso. A segunda fonte é também rebocada a uma distância selecionada atrás da primeira fonte. A primeira fonte é ativada e os sinais sísmicos são registrados correspondendo a atuações da primeira fonte. Depois do rebocamento a embarcação se moveu de forma que a segunda fonte esteja disposta consideravelmente na mesma posição geodésica como a primeira fonte foi no tempo de sua atuação, A segunda fonte é ativada e sinais sísmicos são novamente registrados. Um conjunto de dados sísmicos “sem fantasma” é obtido usando a técnica descrita mais completamente na publicação Egan et al.
Uma das questões principais sobre/sob a fonte técnica descrita na publicação Egan et al. Se referindo ao tratado acima é que o número de posições de disparo é metade quando comparada com as técnicas conven- cionais de atuação da fonte, causando uma cobertura de dobra a ser a metade. Outra questão com essa técnica é se os receptores sísmicos são rebocados atrás de uma embarcação e, portanto se movendo de disparo a disparo, é que os receptores se moveram a uma distância considerável entre quando as fontes em diferentes profundidades são acionadas. Para manter o número de posições de disparo e cobertura de dobra como na aquisição sísmica marinha convencional, e para minimizar a diferença nas posições do receptor quando as fontes em diferentes profundidades são acionadas, é desejável ter um método para uma extração da fonte fantasma que permitas com que as fontes rebocadas em diferentes profundidades serem acionadas durante a gravação de cada registro de disparo.
Isto é também conhecido na técnica, o uso de dois ou mais cabos flutuadores rebocados em diferentes profundidades para obter um sinal gravado que é consideravelmente livre de fantasmas causados pelas reflexões na superfície da água próximas aos cabos flutuadores. Veja, por exemplo, Posthumus, B. [1993] Deghosting of twin streamer configuration, Geophysical Prospecting, vol. 41, pp. 267-286. Adaptação direta do método Pos- thiimiiQ nara rpmnrân Ha fantacmac “cnhrA/Qnh” ac fnntAc annHo fnntAC Am II IUI 111*155 |·#dl d I v?l 11Ibr\^ÍAKJ UC Idl 114312511IGI25 55vJUI v?/ d25 I ví11lví%5 j dUl Ivlv3 IUI IΙν?ϊ5 ví111 diferentes profundidades são acionadas com atrasos de tempo dentro dos mesmos registros sísmicos provou ser impreciso.
Portanto, ainda continua a necessidade por técnicas melhoradas para a redução dos efeitos das fontes fantasmas nos dados sísmicos marinhos.
Sumário da Invenção Em uma modalidade, a invenção é um método para determinar um campo de onda de uma fonte de energia sísmica marinha sem fantasma dos registros de sinal dos dados sísmicos adquiridos pela ativação de uma primeira fonte sísmica em um primeiro tempo e uma ou mais fontes sísmicas adicionais em seus próprios tempos característicos com respeito a um índice de tempo dos registros de sinal, em que as fontes são consideravelmente estabelecidas e dispostas em diferentes profundidades de um corpo de água. Um primeiro campo de onda é determinado a partir dos registros de sinal que ocorreríam se uma primeira fonte fosse acionada em um tempo selecionado com respeito a um tempo de iniciação dos registros de sinal. O primeiro campo de onda é ajustado pelo tempo com respeito à profundidade na água da primeira fonte. Um ou mais campos de ondas adicionais são determinados a partir dos registros de sinal que ocorreríam se uma ou mais fontes adicionais estivessem, cada uma, acionadas em um determinado tempo selecionado com respeito ao dito tempo de iniciação dos registros de sinal. Um ou mais campos de onda adicionais são ajustados pelo tempo com respeito às profundidades na água de uma ou mais fontes adicionais. O primeiro campo de onda é combinado com um ou mais campos de ondas adicionais para determinar um campo de onda sem fantasma correspondendo à atuação de uma única fonte sísmica.
Outros aspectos e vantagens da invenção irão se tornar aparentes a partir da seguinte descrição e das reivindicações em apêndice.
Breve descrição dos desenhos A invenção e suas vantagens podem ser mais facilmente entendidas com referência a seguinte descrição detalhada e os desenhos anexados, em que: Figura 1 mostra a aquisição de dados sísmicos em um corte transversal para mostrar um exemplo de um arranjo das fontes de energia sísmica;
Figura 2 mostra uma vista em planta para adquirir dados sísmicos para mostrar o exemplo de um arranjo de cabos flutuadores de receptor sísmico;
As Figuras 3A, 3B, e 3C mostram, respectivamente, exemplos de assinaturas de duas fontes operadas em diferentes momentos em diferentes profundidades, a decomposição por uma técnica anterior (o método descrito em Posthumus [1993]), e a decomposição por um método de acordo com a invenção;
As Figuras 4A e 4B mostram, respectivamente, espectros de campos de ondas sem fantasma in Figura 3 usando o método Posthumus e o método de uma invenção;
As Figuras 5A e 5B mostram, respectivamente, erro espectral em campos de ondas sem fantasma quando usando o método de Posthumus e o método de uma invenção;
As Figuras 6A e 6B mostram, respectivamente, o espectro de um denominador não filtrado e exato em uma equação de acordo com o método de uma invenção, e o espectro de um filtro de passagem de banda aplicado aos dois últimos termos de tal denominador;
As Figuras 7A e 7B mostram, respectivamente, o espectro do denominador nas equações usadas no método de Posthumus, e o espectro de uma versão do filtro de passagem de banda do denominador em um método de uma invenção;
As Figuras 8A e 8B mostram, respectivamente, o espectro dos últimos dois termos do denominador na equação do método de uma invenção sem filtro de passagem de banda e com um filtro de passagem de banda anlirarlrv a apiiuauu, c Figura 9 shows a flowchart ilustrando uma modalidade do método de uma invenção para determinar a campo de onda de uma fonte de energia sísmica marinha sem fantasma.
Enquanto a invenção irá ser descrita em conexão com suas modalidades preferidas, isso irá ser entendido que uma invenção não está limitada a estes. Por outro lado, a invenção tem a intenção para cobrir todas as alternativas, modificações, e equivalentes que podem ser incluídos dentro do escopo de uma invenção, conforme definido pelas reivindicações em apêndice.
Descrição detalhada Figura 1 mostra em uma vista em corte transversal um exemplo de um arranjo da adquirir dados sísmicos de acordo com a invenção. A análise sísmica da embarcação 10 se move junto com a superfície 11A de um corpo de água 11 como um lago ou o oceano. A embarcaçãolO tipicamente inclui equipamento mostrado geralmente em 12 e é referido como conveniência como um “sistema de gravação”. Um sistema de gravação 12 pode incluir dispositivos (nenhum dos quais é mostrado separadamente) para a- cionar seletivamente as fontes de energia sísmica 14, 16 (explicado abaixo), para gravar os sinais gerados pelos sensores ou receptores 20 (explicado abaixo) em resposta a energia sísmica conferida na água 11 (pelas fontes 14, 16) e através disso nas formações rochosas 19, 21 abaixo do fundo da água 13, e para determinar a posição geodésica de uma embarcaçãolO, as fontes de energia sísmicas 14, 16 e cada um de uma pluralidade de sensores sísmicos ou receptores 20 em qualquer tempo. A embarcaçãolO é mostrada rebocando duas fontes de energia sísmicas 14, 16, para propósitos de clareza ilustrativa apenas. A invenção tem a intenção de ser empregada com dois ou mais fonte sísmicas e não é restrita a apenas duas fontes. As fontes de energia sísmicas 14, 16 podem ser de qualquer tipo de fonte de energia marinha incluindo, mas não limitado a, canhões de ar e canhões de água, ou matrizes de tais fontes de energia. No exemplo mostrado na Figura 1, as fontes 14, 16 são rebocadas a consideravelmente a mesma distância atrás da embarcaçãolO e em diferentes profundidades na água 11. Tal arranjo das fontes é referido como “so-bre/sob.” Em outros exemplos, as fontes 14, 16 podem ser rebocadas por uma embarcação diferente (não mostrado), ou pode ser em uma posição fixa (desde que as profundidades são diferentes conforme mostrado na Figura 1). Então, tendo uma embarcação de pesquisa 10 rebocando as fontes 14, 16 não é um limite no escopo da presente invenção. A embarcaçãolO é também mostrada rebocando um cabo flutu-ador sísmico 18 tendo uma pluralidade de receptores sísmicos espaçados longitudinalmente 20 nisso. Contudo, essa invenção é geralmente relacionada com fontes de energia sísmicas, e então pode ser usada em conjunto com qualquer tipo de arranjo de receptor sísmico. Em modalidades diferentes, uma invenção pode ser usada como um reboque de cabos flutuadores sísmicos em qualquer configuração, cabos de fundo de oceano, sensores posicionados em poços, etc., e com qualquer tipo de sensor de recepção, incluindo, mas não limitado a, sensores de pressão, sensores graduais de tempo de pressão, sensores de velocidade, medidores de aceleração, etc., ou qualquer combinação disso.
Em um tempo selecionado durante a operação de um sistema de aquisição mostrado na Figura 1, o sistema de aquisição 12 aciona um primeiro de uma das fontes de energia sísmicas, por exemplo, a fonte 14. O tempo selecionado pode ser desenvolvido com respeito ao in início da gravação respostas do receptor ou qualquer outro tempo referência relacionado com o tempo de gravação do sinal, e tipicamente como um primeiro atraso selecionado do início da gravação. Energia da primeira fonte 14 viaja voltada para fora a partir desse ponto conforme mostrado em 22 e 24. Alguma energia viaja para baixo, 22, aonde ela é refletida nos limites da impedância acústica, por exemplo, o fundo da água 13 e nos limites 15, 17 entre diferentes formações rochosas 19, 21. Apenas o fundo das reflexões da água é mostrado na Figura 1 para a clareza da ilustração. Porções indo para cima da energia da primeira fonte 14, mostrada em 24, são refletidas da superfície da água 11A conforme mostrado na Figura 1. Um sistema de gravação 12 é configurado para acionar uma segunda fonte de energia sísmica, por exemplo, a fonte 16, na extremidade de um segundo tempo relativo ao início das gravações dos dados sísmicos, ou, alternativamente, depois um tempo selecionado antes ou depois da atuação da primeira fonte 14. A Energia viajando voltada para fora de uma segunda fonte 16 se move junto caminhos similares conforme a energia da primeira fonte 14. Na presente invenção, os atrasos de tempo acima descritos, são selecionados de modo que a energia de ambas as fontes seja detectada pelos receptores 20 em cada registro de dados (chamado de um registro de disparo). É esperado que o tempo de acionar as fontes em diferentes profundidades possa variar entre os registros de disparos em uma maneira aleatória, semi-aleatória, ou sistemática. A Figura 2 mostra o arranjo da Figura 1 em vista em planta para ilustrar o rebocamento de uma pluralidade cabos flutuadores lateralmente espaçados18. Os cabos flutuadores 18 podem ser mantidos em suas posições relativamente laterais e longitudinais com respeito a uma embarca-ção10 usando o equipamento de rebocamento 23 de tipos bem conhecidos na técnica. Novamente, uma embarcaçãolO é mostrada rebocando duas fontes de energia sísmicas 14, 16, para propósitos de clareza ilustrativa a- penas. A invenção tem a intenção de ser empregada com dois ou mais fonte sísmicas e não é restrita to a apenas duas fontes. O que é também mostrado na Figura 2 é que a primeira fonte 14 e uma segunda fonte 16 podem ser lateralmente deslocadas (e/ou longitudinalmente deslocadas em outras modalidades) para evitar, no caso das fontes 14, 16 serem canhões de ar ou matrizes disso, tendo dispersado ar na água 11 da primeira fonte 14 afetando a energia sísmica viajando para cima de uma segunda fonte 16, e vice versa. O deslocamento lateral e/ou longitudinal é contemplado como sendo apenas a poucos metros, de modo que as fontes 14, 16 provenham energia equivalente a sendo elas que ocorreríam se as fontes 14, 16 estivessem no mesmo plano vertical e na mesma distância longitudinal atrás da embarcação, ou expressada diferentemente, na essencialmente à mesma posição geodésica. Evitando a dispersão de ar acima das fontes quando acionadas, os efeitos da superfície da água (11A in Figura 1) na energia emitida de cada fonte 14, 16 irá ser, ajustada pela profundidade da água, consideravelmente a mesma que o efeito disso na outra fonte 14, 16. A atuação da fonte e a gravação de sinal explicado acima podem ser repetidas uma pluralidade de vezes enquanto uma embarcaçãolO, as fontes 14, 16 e os cabos flutuadores 18 se movem através da água 11. Conforme explicado acima, cada registro de disparo irá incluir, para cada receptor 20, sinais correspondendo à energia sísmica produzida por ambas as fontes: a primeira fonte 14 e uma segunda fonte 16. A Figura 1 e a 2 ilustram duas fontes cada uma operando em uma profundidade diferente. Isso é apenas para propósitos ilustrativos. O método não está limitado a duas fontes operando em diferentes profundidades, mas pode incluir qualquer multiplicidade de fontes operando em diferentes profundidades.
Os dados sísmicos obtidos pela execução de uma análise sísmica, representativa da subsuperfície terrestre, são processados para gerar informação relacionada com a estrutura geológica e as propriedades das formações de subsuperfície do solo na área sendo pesquisada. Os dados sísmicos processados são processados para exibição e análise de conteúdo em potencial de hidrocarbonetos destas formações subterrâneas. O objetivo do processamento dos dados sísmicos é para extrair dos dados sísmicos tanta informação quanto possível no que diz respeito às formações subterrâneas de forma a visualizar adequadamente a subsuperfície geológica. De forma a identificar lugares na subsuperfície terrestre aonde há uma probabilidade pra achar acúmulos de petróleo, grandes quantidades de dinheiros estão sendo gastos em reunir, processar, e interpretar os dados sísmicos. O processo de construir as superfícies de refletor definindo a camadas subterrâneas da terra de interesse dos dados sísmicos registrados que geram uma imagem da terra em profundidade ou tempo. A imagem da estrutura de subsuperfície terrestre é produzida de maneira a permitir que um interpretador selecione lugares com a maior probabilidade de ter acúmulos de petróleo. Para verificar a presença de petróleo, um poço precisa se perfurado. A perfuração de poços para determinar onde os depósitos de petróleo estão presentes ou não, é uma tarefa extremamente dispendiosa e que consome muito tempo. Por essa razão, existe uma necessidade contínua de melhorar o processamento e a exibição de dados sísmicos, de forma a produzir uma imagem da estrutura de subsuperfície terrestre que irá melhorar a habilidade de um interpretador, quando a interpretação for feita por um computador ou um humano, para acessar a probabilidade que um acúmulo de petróleo exista em um local em especial na subsuperfície terrestre. O processamento e a exibição de dados sísmicos adquiridos facilitam mais as decisões com precisão no caso de quando e onde perfurar, e através disso reduzindo o risco de se perfurarem buracos secos.
Tendo explicado uma aquisição técnica que pode ser usada com a presente invenção, um método de uma invenção para decomposição de campo de onda irá agora ser explicado. As seguintes modalidades exemplares empregam duas fontes em duas profundidades para simplicidade ilustrativa apenas. A invenção abrange o uso de qualquer número de fontes em diferentes profundidades.
Primeiro, no domínio de tempo, os operadores fantasma gn (fun- ções que determinam a assinatura fantasma de suas respectivas fontes sísmicas) podem ser definidos para cada fonte conforme se segue. Para a primeira fonte sísmica, e para uma segunda fonte sísmica, Nas expressões mencionadas anteriormente, di é a profundidade de operação da primeira fonte; d2 é a profundidade de operação da segunda fonte; a é o ângulo de emissão da propagação de energia da fonte, relativa à linha vertical; v é velocidade acústica na água, que é cerca de 1500 m/s; e <5é o operador Dirac delta.
Deixar w(t) representar o campo de onda que seria registrado em um receptor selecionado a partir da primeira fonte se = 0 segundo (no início do registro), se não houve reflexão fantasma de fonte, e se a assinatura de fonte for um reforço perfeito. O sinal realmente registrado r(t) que resulta do acionamento da primeira fonte sendo acionado em um tempo real após o registro ser iniciado t = U e a segunda fonte acionada no tempo t = t2 pode ser expresso como: onde n(t) representa ruído, Si(t) e s2(t) representam a assinatura sem fantasma das primeira e segunda fontes, respectivamente, e (*) representa convolação do domínio de tempo.
Mudar r(t) mais cedo no tempo por ti segundos na Equação (1) e, portanto, corrigir este para o tempo real de acionamento da primeira fonte com relação ao início do registro (tempo t = 0) no registro de projeção fornece a expressão para um primeiro sinal de tempo corrigido Oi(t): Transformar a Equação (2a) no domínio de frequência fornece a seguinte expressão: (2b) em que Νι(ω) é a transformação de n(t+ti) e Na Equação (2b) e expressões adicionais, capitalização de cada variável representa sua transformação a partir do domínio de tempo ao domínio de frequência.
Mudar r(t) mais cedo no tempo por t2 - {(d2 - di)cos(a)/v} segundos na Equação (1) e, portanto corrigindo este para o tempo de acionamento da segunda fonte e compensando por ser estar mais perto dos refletores do que os resultados da primeira fonte na expressão seguinte para um segundo sinal de tempo corrigido o2(t): Transformar a Equação (3a) ao domínio de frequência fornece a expressão: (3b) em que Ν2(ω) é a transformação de e Se duas cópias (como ilustrado aqui, embora mais cópias fossem feitas) dos campos de onda registrados forem feitas, e Οι(ω) é definido como os dados registrados corrigidos para o tempo de acionamento da primeira fonte, e 02(ω) é definido como os dados registrados corrigidos para o tempo de acionamento da segunda fonte (explicado acima), e ambos os re- gistros são corrigidos para a pequena diferença entre os tempos de chegada devido às profundidades de operação diferentes de cada uma das duas fontes, o resultado pode ser expresso como explicado acima: onde AT2é o tempo de acionamento da segunda fonte menos o tempo de acionamento da primeira fonte, mais a diferença de tempo causada pelas profundidades de operação diferentes para as fontes, e ΔTV é, de modo similar, o tempo de acionamento da primeira fonte menos o tempo de acionamento da segunda fonte, mais a diferença de tempo causada pelas diferentes fontes de profundidade de operação. Essencialmente, a Equação (4a) representa a primeira fonte de dados do tempo corrigido e pode ser descrita como o campo de onda da primeira fonte mais o campo de onda da segunda fonte com um atraso de tempo que depende na diferença de tempos de acionamento entre as duas fontes e a diferença nas fontes de profundidade, mais o ruído. A Equação (4b) representa segunda fonte de dados registrados de tempo corrigido de uma maneira similar.
Através de multiplicar as Equações (4a) e (4b) com o conjugado complexo das funções fantasma pelas respectivas fontes de profundidade, os resultados são as seguintes expressões: em que * representa o conjugado complexo dos respectivos operadores fantasma em um domínio de frequência.
Adicionar as Equações (5a) e (5b) fornece a expressão seguinte: A partir da Equação (6), a expressão seguinte pode ser derivada para calcular a assinatura de fonte corrigida, remoção de fantasma da fonte de campo de onda: (7) onde O método proposto por Posthumus (referência feita na seção de Antecedentes no presente documento) para remover os fantasmas usando receptores superiores e inferiores, em vez de fontes superiores e inferiores acionadas nos mesmos registros sísmicos em tempos diferentes, usa a expressão seguinte para calcular o campo de onda sem fantasma de fita. (8) Ao comparar a Equação (7) com a Equação (8), é evidente que as duas Equações são muito similares. A diferença entre as duas Equações está no denominador, em que calcular a assinatura de fonte corrigido da fonte sem fantasma do campo de onda inclui a assinatura de fontes sem fantasma, e dois termos adicionais. Os dois termos adicionais, Gi(cü)G2(u))* e G2(u))Gi(cú)* são equivalentes a correlação cruzada das duas funções fantasma em um domínio de tempo, e expf-iujATrf e βχρ(-ίωΔΤ2) são operadores de atraso de tempo para corrigir a diferença entre os tempos de acionamento e as fontes de profundidade de operação (7\Γί e ΔΤ2, como explicado acima). Foi determinado que o uso de termos extras, como explicado acima pode reduzir substancialmente erros na decomposição de sinal como contrastado com o uso do método Posthumus diretamente para fontes superiores e inferiores acionadas para serem registradas no mesmo registro de projeção. Adicionalmente, a menos que as duas (ou mais) fontes operadas em diferentes profundidades sejam especificamente projetadas para ser quase idênticas, a assinatura de fontes precisa ser incluída na expressão para computar as diferenças nas assinaturas como descrito na Equação (7).
Porque estas funções fantasma correlacionadas de forma cruzada são mudadas para frente e para trás no tempo, o denominador terá picos e descontinuidades no espectro relacionado às diferenças em tempo de acionamento entre duas fontes. Os picos e descontinuidades podem potencialmente causar uma amplificação de ruído nas partes do espectro, especialmente em relação às extremidades de frequência alta e baixa. Portanto, é desejável amortecer ou filtrar estes dois termos adicionais quando calculando o campo de onda sem fantasma. As Equações (7) e (8) são similares exceto por estes dois termos adicionais no denominador da Equação 7. O amortecimento de dois termos adicionais na Equação (7) faz com que seja possível abordar o resultado de aplicar a Equação (8) a dados de cima/ debaixo, onde duas fontes são ativadas nos mesmos registro, com diferentes tempos de ativação. Entretanto, mesmo se alguma filtragem de passagem de banda dos dois termos adicionais na Equação (7) for necessária para evitar a amplificação do ruído em direção as extremidades superior e inferior do espectro, os erros introduzidos são muito pequenos, sabendo-se que a filtra- gem é feita nas partes do espectro onde a correlação cruzada das duas funções fantasma têm energia limitada e/ou nas partes do espectro com energia sísmica limitada. Tal versão filtrada da Equação (7) pode ser expressa como: (9) onde e F(ω) é um filtro de passagem de banda. As Equações (7) e (9) podem ser prontamente expandidas para os casos de ativar três ou mais fontes em três ou mais profundidades nos mesmos registros com atrasos de tempo diferentes.
Porque é normalmente possível separa perfeitamente o ruído do sinal, apenas a primeira parte do lado à direita das Equações (7) ou (9) seria normalmente usada para calcular a assinatura de fonte corrigida, remoção de fantasma de campo de onda.
Observe que as modalidades da invenção ilustradas acima incluem determinar as funções de correlação cruzada de operadores fantasma da primeira e uma ou mais fontes adicionais. Em geral, o método da invenção pode ser estendido para quaisquer números de fontes, cada uma operando em diferentes profundidades de água. Além disso, as assinaturas das fontes individuais são computadas. A Figura 3A mostra exemplos de assinaturas de duas fontes operadas em diferentes times e em diferentes profundidades na curva 40, onde as assinaturas sem fantasma das duas fontes são idênticas e, portanto não computadas. A Curva 42 na Figura 3B mostra uma assinatura de fonte sem fantasma calculada através do uso de acordo com a técnica Posthumus. Os erros na assinatura de fonte são claramente visíveis na curva 42 tanto em um tempo mais cedo quanto mais tarde do que o sinal sem fantasma. A diferença de tempo é relacionada à diferença nos tempos de acionamento das fontes em diferentes profundidades, e a diferença nos tempos de chegada devido às diferenças nas fontes de profundidade. A decomposição por um método de acordo com a invenção é mostrada na curva 44 na Figura 3C. Muito pouco erro na assinatura de fonte sem fantasma é visível usando o método de acordo com a invenção.
As Figuras 4A e 4B mostram, respectivamente, o espectro dos campos de onda sem fantasma mostrados na Figura 3 usando o método Posthumus, e o método da invenção.
As Figuras 5A e 5B mostram, respectivamente, o erro de espectro nos campos de onda sem fantasma quando é usado o método Posthumus e o método da invenção.
As Figuras 6A e 6B mostram, respectivamente, o espectro do denominador exato e não filtrado em uma equação (por exemplo, Equação (7)) de acordo com o método da invenção, e o espectro de um filtro de passagem de banda aplicado aos dois últimos termos em tal denominador como explicado acima.
As Figuras 7A e 7B mostram, respectivamente, o espectro do denominador nas Equações usadas no método Posthumus, e do espectro da versão filtrada da passagem de banda do denominador no método da invenção.
As Figuras 8A e 8B mostram, respectivamente, o espectro dos dois últimos termos do denominador na equação do método da invenção sem a filtragem de passagem de banda e com um filtro de passagem de banda aplicado. A Figura 9 mostra um fluxograma que ilustra uma modalidade do método da invenção para determinar campo de onda sem fantasma de fonte de energia sísmica marinha.
No bloco 90, os dados sísmicos são adquiridos a partir dos re- gistros de sinais através do acionamento de uma primeira fonte sísmica em um primeiro tempo e uma ou mais fontes adicionais sísmicas em seus próprios tempos característicos com relação a um tempo índice dos registros de sinais, as fontes substancialmente colocadas e dispostas em diferentes profundidades em um corpo de água.
No bloco 91, um primeiro campo de onda é determinado a partir de registros de sinais no bloco 90 que ocorrería se a primeira fonte fosse acionada em um tempo selecionado com relação tempo de iniciação dos registros de sinais, o primeiro campo de onda sendo o tempo ajustado com relação à profundidade na água da primeira fonte. Em uma modalidade da invenção, o tempo selecionado pode ser o tempo de iniciação.
No bloco 92, um ou mais campos de onda adicionais são determinados a partir de registros de sinais no bloco 90 que ocorrería se uma ou mais fontes adicionais fossem cada uma acionadas no dito tempo selecionado com relação ao dito tempo de iniciação dos registros de sinais, um ou mais campos de onda adicionais sendo tempo ajustado com relação à profundidade na água de uma ou mais fontes adicionais.
No bloco 93, o primeiro campo de onda a partir do bloco 91 é combinado com um ou mais campos de onda adicionais a partir do bloco 92 para determinar um campo de onda sem fantasma que corresponde ao acionamento de uma única fonte sísmica.
Os métodos de acordo com a invenção podem fornecer um uma amplitude de sinal em dados sísmicos marinhos sema necessidade de aumentar o tamanho das fontes de energia sísmica usadas. A invenção fornece uma técnica para decompor sinais sísmicos registrados com esta e que são registrados usando apenas uma fonte, enquanto faz o uso de duas ou mais fontes operadas em diferentes profundidades e são acionadas para serem registradas nos mesmo registros de projeção. O campo de onda fontes marinhas de energia sísmica é substancialmente sensível a pressão hidrostática, que por sua vez é uma função de profundidade de fonte. Adicionalmente, pode ser desejável usar diferentes configurações de fontes para as profundidades de operação diferentes.
Portanto, uma correção para a assinatura de respostas de fonte pode necessitar ser aplicada como mostrado nas Equações (7) e (9). Observe que tal correlação seria necessária se as respostas de fonte individuais, excluindo os fantasmas, fossem especificamente projetadas para serem quase idênticas. Há uma variedade de técnicas conhecidas para projetar, medir ou calcular os campos de onda das fontes sísmicas, que têm diferentes níveis de precisão. O campo de onda ou as posições selecionadas no campo de onda podem ser medidos diretamente (por exemplo, medição de campo afastado) ou o campo de onda pode ser calculado com base nos modelos físicos da fonte. Há também vários métodos de monitoramento de fonte, que determinam o campo de onda do conjunto de fonte de projeção a projeção, usando diversos sensores dispostos no conjunto de fonte sísmica. Isto inclui o chamado Notional source method, descrito em Ziolkowski, A., Parkes, G., Hatton, L. e Haugland, T., Γ19821. The sianature ofan airaun arrav: Comouta-tion from near-field measurements including interactions, Geophysics, vol. 47, p. 1413-1421, e, por exemplo, Method of Sismic Source Monitoríng Using Modeled Assinatura de fontes with Calibration Function, Patente Ns U.S 7.218.572 concedida à Parkes e relacionado com a presente invenção.
Enquanto a invenção foi descrita com relação a um número limitado de modalidades, aqueles versados na técnica, tendo o benefício desta descrição, vão apreciar que outras modalidades podem ser desenvolvidas sem se distanciar do escopo invenção, como descrito no presente documento. Consequentemente, o escopo da invenção deve ser limitado apenas pelas reivindicações em anexo.
REIVINDICAÇÕES
Claims (20)
1. Método para determinar o campo de onda sem fantasma de fonte de energia sísmica marinha a partir de registros de sinais de dados sísmicos adquiridos através de acionar uma primeira fonte sísmica em um primeiro tempo e uma ou mais fontes adicionais sísmicas em seus próprios tempos característicos com relação a um tempo índice dos registros de sinais, as fontes substancialmente colocadas e dispostas em diferentes profundidades em um corpo de água, que compreende: determinar a partir dos registros de sinais um primeiro campo de onda que iria ocorrer se a primeira fonte foi acionada em um tempo selecionado com relação tempo de iniciação dos registros de sinais, o primeiro campo de onda sendo o tempo ajustado com relação à profundidade na água da primeira fonte; determinar a partir dos registros de sinais um ou mais campos de onda adicionais que iriam ocorrer se uma ou mais fontes adicionais fossem cada uma acionadas no dito tempo selecionado com relação ao dito tempo de iniciação dos registros de sinais, um ou mais campos de onda adicionais sendo o tempo ajustado com relação à profundidade na água de uma ou mais fontes adicionais; e unir o primeiro campo de onda com um ou mais campos de onda adicionais para determinar o campo de onda sem fantasma que corresponde ao acionamento de uma fonte sísmica única.
2. Método, de acordo com a reivindicação 1, em que a primeira fonte e uma ou mais fontes adicionais compreendem, cada uma. pelo menos, um canhão de ar.
3. Método, de acordo com a reivindicação 1, que compreende adicionalmente a correção para as assinaturas das fontes rebocadas em diferentes profundidades.
4. Método, de acordo com a reivindicação 1, que compreende adicionalmente determinar funções de correlação cruzada de operadores fantasma da primeira fonte e uma ou mais fontes adicionais.
5. Método, de acordo com a reivindicação 4, que compreende adicionalmente filtragem de passagem de banda das funções de correlação cruzada.
6. Método, de acordo com a reivindicação 1, em que a primeira fonte e uma ou mais fontes adicionais são dispostas em diferentes planos verticais para evitar efeitos de ar na determinação dos campos de onda individuais.
7. Método, de acordo com a reivindicação 1, que compreende adicionalmente projetar fontes rebocadas em diferentes profundidades, de modo que, tenham assinaturas sem fantasmas substancialmente idênticas.
8. Método, de acordo com a reivindicação 1, em que o tempo selecionado é o tempo de iniciação.
9. Método, de acordo com a reivindicação 3, em que a união do primeiro campo de onda com um ou mais campos de onda adicionais para determinar o campo de onda sem fantasma compreende aplicar a seguinte equação: onde e onde \Ν(ω) é o campo de onda sem fantasma; Si(uj) e S2(gj) são assinaturas sem fantasma para as primeira e segunda fontes, respectivamente; Gi(üj) e G2(uj) são operadores fantasma para as primeira e segunda fontes sísmicas, respectivamente; Νι(ω) e Ν2(ω) são termos de ruídos para as primeira e segunda fontes sísmicas, respectivamente; e onde ti e t2 são tempos de acionamento das primeira e segunda fontes sísmicas, respectivamente; di e d2 são profundidades de operação das primeira e segunda fontes sísmicas, respectivamente; aéo ângulo de emissão de propagação de energia a partir de fontes, relativas à linha vertical; e v é uma velocidade acústica na água.
10. Método, de acordo com a reivindicação 5, em que unir o primeiro campo de onda com um ou mais campos de onda adicionais para determinar o campo de onda sem fantasma, usando filtragem de passagem de banda, compreende aplicar a equação seguinte: onde e onde W(w) é o campo de onda sem fantasma; Si(u)) e S2(u)) são as assinaturas sem fantasma para as primeira e segunda fontes, respectivamente; Gi(ω) e G2(cú) são operadores fantasma para as primeira e segunda fontes sísmicas, respectivamente; Νι(ω) e Ν2(ω) são termos de ruídos para as primeira e segunda fontes sísmicas, respectivamente; F(ω) é o filtro de passagem de banda; e onde ti e t2 são tempos de acionamento das primeira e segunda fontes sísmicas, respectivamente; di e d2 são profundidades de operação das primeira e segunda fontes sísmicas, respectivamente; a é o ângulo de emissão de propagação de energia das fontes, relativo à linha vertical; evé uma velocidade acústica na água.
11. Método para determinar uma fonte de campo de onda sem fantasma de energia sísmica marinha em um corpo de água, que compreende; acionar uma primeira fonte de energia sísmica em uma primeira profundidade na água em um primeiro tempo selecionado com relação à iniciação do registro de sinais sísmicos; acionar uma ou mais fontes adicionais de energia sísmica colocadas substancialmente com a primeira fonte, cada fonte adicional com suas próprias profundidades de água características e cada fonte adicional com seus próprios tempos característicos com relação à iniciação do registro de sinais sísmicos, os tempos de acionamento são selecionados, de modo que, a energia das fontes adicionais está presente nos registros de sinais sísmicos; registrar sinais sísmicos que resultam do acionamento da primeira fonte e uma ou mais fontes adicionais; determinar a partir dos registros de sinais um primeiro campo de onda que iria ocorrer se a primeira fonte foi acionada em um tempo selecionado com relação tempo de iniciação de registro de sinais sísmicos, o primeiro campo de onda sendo o tempo ajustado com relação à profundidade na água da primeira fonte; determinar a partir dos registros de sinais um ou mais campos de onda adicionais que iria ocorrer se uma ou mais fontes adicionais fossem cada uma acionadas no dito tempo selecionado com relação tempo de iniciação de registro de sinais sísmicos, e cada campo de onda adicional sendo o tempo ajustado com relação à profundidade na água de sua fonte associada; e unir o primeiro campo de onda com um ou mais campos de onda adicionais para determinar o campo de onda sem fantasma correspondente ao acionamento de uma única fonte sísmica.
12. Método, de acordo com a reivindicação 11, em que a primeira fonte e uma ou mais fontes adicionais compreendem, cada uma, pelo menos, um canhão de ar.
13. Método, de acordo com a reivindicação 11, em que a primeira fonte e uma ou mais fontes adicionais são dispostas em planos verticais diferente para evitar efeitos de ar, na determinação dos campos de onda individuais.
14. Método, de acordo com a reivindicação 11, que compreende adicionalmente corrigir as assinaturas das fontes rebocadas em diferentes profundidades.
15. Método, de acordo com a reivindicação 11, que compreende adicionalmente determinar funções de correlação cruzada de operadores fantasma da primeira fonte e uma ou mais fontes adicionais.
16. Método, de acordo com a reivindicação 11, que compreende adicionalmente filtragem de passagem de banda de funções de correlação cruzada.
17. Método, de acordo com a reivindicação 11, que compreende adicionalmente projetar fontes rebocadas em diferentes profundidades, de modo que tenham substancialmente assinaturas sem fantasma substancialmente idênticas.
18. Método, de acordo com a reivindicação 11, em que o tempo selecionado é a tempo de iniciação.
19. Método, de acordo com a reivindicação 14, em que a união do primeiro campo de onda com um ou mais campos de onda adicionais para determinar o campo de onda sem fantasma, compreende aplicar a seguinte equação: onde e onde W(uS) é o campo de onda sem fantasma; Si(oü) e S2(cv) são assinaturas sem fantasma para as primeira e segunda fontes, respectivamente; Gi(üj) e Θ2(ω) são operadores fantasma para as primeira e segunda fontes sísmicas, respectivamente; Ν^ω) e Ν2(ω) são termos de ruídos para as primeira e segunda fontes sísmicas, respectivamente; e onde ti e t2 são tempos de acionamento da primeira e segunda fontes sísmicas, respectivamente; di e d2 são profundidades de operação da primeira e segunda fontes sísmicas, respectivamente; σ é o ângulo de emissão de propagação de energia das fontes, com relação à linha vertical; e v é a velocidade acústica na água.
20. Método, de acordo com a reivindicação 16, em que a união do primeiro campo de onda com um ou mais campos de onda adicionais para determinar os campos de onda sem fantasma, usando filtragem de passagem de banda, compreende aplicar a sequinte equação: onde e onde \Ν(ω) é o campo de onda sem fantasma; S^üj) e S2(ü)) são assinaturas sem fantasma para as primeira e segunda fontes, respectivamente; Θ-ι(ω) e G2(co) são operadores fantasma para as primeira e segunda fontes sísmicas, respectivamente; Νι(ω) e Ν2(ω) são termos de ruídos para as primeira e segunda fontes sísmicas, respectivamente; F(üj) é o filtro de passagem de banda; e onde ti e t2 são tempos de acionamento das primeira e segunda fontes sísmicas, respectivamente; c/? e d2 são profundidades de operação das primeira e segunda fontes sísmicas, respectivamente; cr é o ângulo de emissão de propagação de energia de fontes, com relação à linha vertical; e vé a velocidade acústica na água.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| US12/806,098 | 2010-08-05 | ||
| US12/806,098 US10838095B2 (en) | 2010-08-05 | 2010-08-05 | Wavefield deghosting of seismic data recorded using multiple seismic sources at different water depths |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| BRPI1103858A2 true BRPI1103858A2 (pt) | 2017-02-14 |
| BRPI1103858B1 BRPI1103858B1 (pt) | 2020-02-18 |
Family
ID=44532652
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| BRPI1103858-6A BRPI1103858B1 (pt) | 2010-08-05 | 2011-08-05 | Removedor de fantasma de campo de onda de dados sísmicos registrados usando múltiplas fontes sísmicas em diferentes profundidades de água |
Country Status (10)
| Country | Link |
|---|---|
| US (1) | US10838095B2 (pt) |
| EP (1) | EP2416179A3 (pt) |
| CN (1) | CN102375155B (pt) |
| AU (1) | AU2011203527B2 (pt) |
| BR (1) | BRPI1103858B1 (pt) |
| CA (1) | CA2747028A1 (pt) |
| EA (1) | EA021738B1 (pt) |
| EG (1) | EG26195A (pt) |
| MX (1) | MX2011008219A (pt) |
| MY (1) | MY180475A (pt) |
Families Citing this family (27)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2087313A4 (en) * | 2006-09-22 | 2015-04-01 | Sercel Inc | SEISMIC NETWORK HAVING VARIABLE PRESSURE SPACED SOURCES |
| FR2981758B1 (fr) * | 2011-10-19 | 2013-12-06 | Cggveritas Services Sa | . |
| AU2013208683A1 (en) * | 2012-01-12 | 2014-07-31 | Geco Technology B.V. | Simultaneous marine vibrators |
| US9411062B2 (en) * | 2012-02-07 | 2016-08-09 | Pgs Geophysical As | Method and system for determining source signatures after source ghost removal |
| US9360576B2 (en) | 2012-04-27 | 2016-06-07 | Pgs Geophysical As | Methods and apparatus for generating deghosted seismic data |
| FR3000224B1 (fr) * | 2012-12-26 | 2015-01-09 | Cggveritas Services Sa | Method and device to deghost seismic data |
| US10338250B2 (en) | 2013-03-14 | 2019-07-02 | Pgs Geophysical As | Method of removing incoherent noise |
| US9322944B2 (en) | 2013-03-15 | 2016-04-26 | Pgs Geophysical As | Wavefield regularization by 3-D wavefield decomposition for geophysical data |
| US9618636B2 (en) | 2013-03-15 | 2017-04-11 | Westerngeco L.L.C. | Seismic acquisition using phase-shifted sweeps |
| FR3006451B1 (fr) * | 2013-05-30 | 2017-07-28 | Total Sa | Procedes et dispositifs d'acquisition de donnees sismiques relatives a une zone du sous-sol situee sous la mer. |
| US10436923B2 (en) * | 2013-06-11 | 2019-10-08 | Cgg Services Sas | Method and apparatus for receiver-side deghosting of seismic data |
| WO2014209728A1 (en) * | 2013-06-25 | 2014-12-31 | Westerngeco Llc | Seismic wavefield deghosting and noise attenuation |
| US11092710B2 (en) * | 2013-06-27 | 2021-08-17 | Pgs Geophysical As | Inversion techniques using streamers at different depths |
| US10459100B2 (en) * | 2013-06-27 | 2019-10-29 | Pgs Geophysical As | Survey techniques using streamers at different depths |
| US10598807B2 (en) * | 2014-02-18 | 2020-03-24 | Pgs Geophysical As | Correction of sea surface state |
| MX388342B (es) * | 2014-03-20 | 2025-03-19 | Reflection Marine Norge As | Generacion de un campo de ondas mediante el uso de un arreglo de vibradores sismicos. |
| US9903966B2 (en) | 2014-04-14 | 2018-02-27 | Pgs Geophysical As | Seismic data acquisition |
| US20170204705A1 (en) * | 2014-08-01 | 2017-07-20 | Nexen Data Solutions, Inc. | Method and system for measuring non-drilling times and their application to improve drilling unit efficiency |
| US10317553B2 (en) | 2014-08-13 | 2019-06-11 | Pgs Geophysical As | Methods and systems of wavefield separation applied to near-continuously recorded wavefields |
| 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 |
| US10073183B2 (en) | 2014-10-20 | 2018-09-11 | Pgs Geophysical As | Methods and systems that attenuate noise in seismic data |
| US9910063B2 (en) * | 2015-04-21 | 2018-03-06 | Pgs Geophysical As | Magnetometer as an orientation sensor |
| US10267936B2 (en) | 2016-04-19 | 2019-04-23 | Pgs Geophysical As | Estimating an earth response |
| US10838096B2 (en) | 2016-12-13 | 2020-11-17 | Pgs Geophysical As | Method and apparatus for dynamic roll compensation of streamer for marine geophysical surveying |
| US11119230B2 (en) * | 2017-08-16 | 2021-09-14 | Pgs Geophysical As | Geophysical survey techniques using selective-length processing |
| EP3788409B1 (en) * | 2018-05-02 | 2024-03-20 | TGS-NOPEC Geophysical Company | Seismic source operation at low frequencies |
| CN109633737B (zh) * | 2018-11-15 | 2019-10-18 | 国家海洋局第二海洋研究所 | 基于虚拟和真实深度组合的海上宽频带气枪震源 |
Family Cites Families (28)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| DE3378233D1 (en) | 1982-08-18 | 1988-11-17 | Horizon Exploration Ltd | Underwater seismic testing |
| US4658384A (en) * | 1985-01-07 | 1987-04-14 | Western Geophysical Co. Of America | Method for determining the far-field signature of an air gun array |
| FR2577686B1 (fr) | 1985-02-14 | 1988-05-20 | Inst Francais Du Petrole | Methode de prospection sismique utilisant des vehicules se deplacant en sens inverses |
| GB9226729D0 (en) | 1992-12-22 | 1993-02-17 | Wellcome Found | Therapeutic combination |
| US5924049A (en) | 1995-04-18 | 1999-07-13 | Western Atlas International, Inc. | Methods for acquiring and processing seismic data |
| EP0766836B1 (en) | 1995-04-18 | 2003-01-29 | Western Atlas International, Inc. | Uniform subsurface coverage at steep dips |
| GB0019054D0 (en) | 2000-04-03 | 2000-09-27 | Schlumberger Technology Corp | A seismic source,a marine seismic surveying arrangement,a method of operating a marine seismic source,and a method of de-ghosting seismic data |
| US6480440B2 (en) | 2001-03-07 | 2002-11-12 | Westerngeco, L.L.C. | Seismic receiver motion compensation |
| US6545944B2 (en) | 2001-05-30 | 2003-04-08 | Westerngeco L.L.C. | Method for acquiring and processing of data from two or more simultaneously fired sources |
| GB2389183B (en) * | 2002-05-28 | 2006-07-26 | Westerngeco Ltd | Processing seismic data |
| US6906981B2 (en) | 2002-07-17 | 2005-06-14 | Pgs Americas, Inc. | Method and system for acquiring marine seismic data using multiple seismic sources |
| CN100401778C (zh) | 2002-09-17 | 2008-07-09 | 弗拉迪米尔·切佩尔科维奇 | 具有高压缩率和最小必需资源的快速codec |
| MXPA05010458A (es) * | 2003-04-01 | 2006-03-21 | Exxonmobil Upstream Res Co | Fuente vibratoria de alta frecuencia conformada. |
| US7123543B2 (en) * | 2003-07-16 | 2006-10-17 | Pgs Americas, Inc. | Method for seismic exploration utilizing motion sensor and pressure sensor data |
| US6882938B2 (en) | 2003-07-30 | 2005-04-19 | Pgs Americas, Inc. | Method for separating seismic signals from two or more distinct sources |
| US7042801B1 (en) * | 2004-02-04 | 2006-05-09 | Seismoelectric Soundings, Inc. | System for geophysical prospecting using induce electrokinetic effect |
| US7359283B2 (en) | 2004-03-03 | 2008-04-15 | Pgs Americas, Inc. | System for combining signals of pressure sensors and particle motion sensors in marine seismic streamers |
| US7218572B2 (en) | 2004-05-20 | 2007-05-15 | Pgs Exploration (Uk) Limited | Method of seismic source monitoring using modeled source signatures with calibration functions |
| US7489590B2 (en) * | 2005-04-08 | 2009-02-10 | Westerngeco L.L.C. | Method and apparatus for source and receiver side wave field separation |
| US7196969B1 (en) * | 2006-02-09 | 2007-03-27 | Pgs Geophysical As | Three-dimensional two-way acoustic wave equation pre-stack imaging systems and methods |
| US7379386B2 (en) * | 2006-07-12 | 2008-05-27 | Westerngeco L.L.C. | Workflow for processing streamer seismic data |
| US7492665B2 (en) * | 2006-07-24 | 2009-02-17 | Westerngeco L.L.C. | Seismic data acquisition and source-side derivatives generation and application |
| US8116166B2 (en) * | 2007-09-10 | 2012-02-14 | Westerngeco L.L.C. | 3D deghosting of multicomponent or over / under streamer recordings using cross-line wavenumber spectra of hydrophone data |
| US8077544B2 (en) * | 2008-03-28 | 2011-12-13 | Westerngeco L.L.C. | Dual-wavefield multiple attenuation |
| US8345510B2 (en) | 2008-06-02 | 2013-01-01 | Pgs Geophysical As | Method for aquiring and processing marine seismic data to extract and constructively use the up-going and down-going wave-fields emitted by the source(s) |
| US8218393B2 (en) | 2008-06-30 | 2012-07-10 | Westerngeco L.L.C. | Technique and system to increase the length of a seismic shot record |
| US8811115B2 (en) * | 2008-08-14 | 2014-08-19 | Pgs Geophysical As | Attenuating seismic interference noise using a dual sensor recording system |
| WO2010019957A1 (en) | 2008-08-15 | 2010-02-18 | Bp Corporation North America Inc. | Method for separating independent simultaneous sources |
-
2010
- 2010-08-05 US US12/806,098 patent/US10838095B2/en not_active Expired - Fee Related
-
2011
- 2011-07-12 AU AU2011203527A patent/AU2011203527B2/en not_active Ceased
- 2011-07-21 CA CA2747028A patent/CA2747028A1/en not_active Abandoned
- 2011-07-25 MY MYPI2011003446A patent/MY180475A/en unknown
- 2011-07-31 EG EG2011071289A patent/EG26195A/en active
- 2011-08-01 EA EA201170912A patent/EA021738B1/ru not_active IP Right Cessation
- 2011-08-04 EP EP11176547.5A patent/EP2416179A3/en not_active Withdrawn
- 2011-08-04 MX MX2011008219A patent/MX2011008219A/es active IP Right Grant
- 2011-08-05 BR BRPI1103858-6A patent/BRPI1103858B1/pt not_active IP Right Cessation
- 2011-08-05 CN CN201110223583.0A patent/CN102375155B/zh not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| CN102375155A (zh) | 2012-03-14 |
| EP2416179A2 (en) | 2012-02-08 |
| BRPI1103858B1 (pt) | 2020-02-18 |
| EA021738B1 (ru) | 2015-08-31 |
| AU2011203527B2 (en) | 2016-04-21 |
| MY180475A (en) | 2020-11-30 |
| US10838095B2 (en) | 2020-11-17 |
| AU2011203527A1 (en) | 2012-02-23 |
| US20120033526A1 (en) | 2012-02-09 |
| EG26195A (en) | 2013-04-10 |
| CA2747028A1 (en) | 2012-02-05 |
| MX2011008219A (es) | 2012-02-20 |
| CN102375155B (zh) | 2017-04-12 |
| EA201170912A1 (ru) | 2012-03-30 |
| EP2416179A3 (en) | 2013-07-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| BRPI1103858A2 (pt) | removedor de fantasma de campo de onda de dados sísmicos registrados usando múltiplas fontes sísmicas em diferentes profundidades de água | |
| US8345510B2 (en) | Method for aquiring and processing marine seismic data to extract and constructively use the up-going and down-going wave-fields emitted by the source(s) | |
| US9684086B2 (en) | Method for acquiring and processing marine seismic data to extract and constructively use the up-going and down-going wave-fields emitted by the source | |
| BR102013002915A2 (pt) | Método e sistema para determinar assinaturas da fonte após remoção de fantasma da fonte | |
| EP2180346A1 (en) | Method for Combining Pressure and Motion Seismic Signals from Streamers where Sensors are not at a Common Depth | |
| NO343375B1 (no) | Seismiske marine langtidsundersøkelser som benytter interpolerte multikomponents streamer-trykkdata | |
| US7336560B2 (en) | Method for determination of sufficient acquisition coverage for a marine seismic streamer survey | |
| NO333602B1 (no) | Dannelse av vertikale seismiske profiler i et boreverktoy | |
| BR102014023204B1 (pt) | Processo tecnológico para identificar uma estrutura subsuperficial potencial, aparelho para realizar formação de imagem sísmica de uma estrutura subsuperficial potencial, memória legível em computador e processo tecnológico para gerar um produto de dados geofísicos | |
| BRPI1104046B1 (pt) | método para criação de imagem da subsuperfície da terra | |
| BRPI0904049A2 (pt) | método para determinar a formação de fator de qualidade de sensor duplo de sinais sìsmicos marinhos | |
| BRPI1006934B1 (pt) | método para separação de campos de onda para dados de sensores de função dupla usando tratamento de dados e migração do tipo kirchhoff | |
| BRPI1001377A2 (pt) | mÉtodo para cÁlculo de atributos sÍsmicos a partir de sinais sÍsmicos | |
| NO341944B1 (no) | Fremgangsmåte for støyundertrykking av støy i seismiske signaler ved bruk av romlige transformasjoner | |
| BRPI0902827A2 (pt) | método para somar sinais de cabo rebocado de sensor duplo usando análise de efeito fantasma de cruzamento | |
| BRPI0902474A2 (pt) | método para somar sinais de cabo sìsmico rebocável com duplo sensor utilizando velocidades de reflexão sìsmica | |
| BR102012007972B1 (pt) | Método e sistema para mapeamento da geologia terrestre e meio legível em computador não transitório | |
| BR102013009382A2 (pt) | Métodos e aparelho para gerar dados sísmicos com fantasma retirado | |
| CN103907032A (zh) | 采用基于相关性的层叠的地震成像系统和方法 | |
| NO20140741A1 (no) | Inversjonsteknikker ved bruk av streamere på ulike dybder | |
| US10393902B2 (en) | Method and system for acquisition of seismic data | |
| Tikhonov et al. | Borehole vicinity acoustic properties study using offset VSP data |
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] | ||
| B06F | Objections, documents and/or translations needed after an examination request according [chapter 6.6 patent gazette] | ||
| B06U | Preliminary requirement: requests with searches performed by other patent offices: procedure suspended [chapter 6.21 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: 20 (VINTE) ANOS CONTADOS A PARTIR DE 05/08/2011, 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 10A 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 2631 DE 08-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. |