Extraindo cabeçalho de arquivo formato SEGY pelo RED HAT sem software de interpretação

Publicado por Rodrigo Paiva em 28/04/2019

[ Hits: 2.786 ]

22 0

Denuncie Favoritos Indicar Impressora

Extraindo cabeçalho de arquivo formato SEGY pelo RED HAT sem software de interpretação

O formato SEGY é arquivo que possui o conteúdo de um sísmica. Para abrir esta sísmica é necessário ter algum software da indústria de Petróleo como, por exemplo, o Petrel ou GeoFrame da Schlumberger, ou outros diversos softwares que necessitam de licenciamento.

Como este arquivo é binário e precisamos da informação do header para identificação do dado, saber de onde veio o dado, informações básicas como País, Estado, Bacia sedimentar, levantamento sísmico, tipo de processamento e empresa que gerou o processamento.

Estas informações são importantes para estudo da área estuda em questão. Utilizamos um comando pelo Linux, sem a necessidade de ter a licença do software para extrair esta informação para um arquivo de texto.

Neste procedimento, será capaz de extrair o cabeçalho do SEGY sem algum programa de interpretação.

Procedimento: Workstation Red Hat Linux

Passos

1. Abra um Terminal Linux.

2. Digite:

bash

3. Digite o comando:

dd se = nome.segy bs = 3200 count = 1 conv = ascii >> /temp/test_header.txt

4. Abra o editor de texto (vi, nedit ou outros)

Foi gerado um arquivo "test_header.txt" com todo conteúdo do cabeçalho do SEGY.

Exemplo:

C 1 CLIENT PETROLEO COMPANY CREW NO
C 2 LINE 0001-0322 AREA ESPIRITO SANTO
C 3 REEL NO DAY-START OF REEL YEAR OBSERVER
C 4 INSTRUMENT MFG MODEL SERIAL NO
C 5 DATA TRACES/RECORD 24 AUXILIARY TRACES/RECORD 0 CDP FOLD 24
C 6 SAMPLE INTERVAL 4000 SAMPLES/TRACE 4500 BITS/IN BYTES/SAMPLE 4
C 7 RECORDING FORMAT FORMAT THIS REEL BRV04 MEASUREMENT SYSTEM
C 8 SAMPLE CODE FLOATING PT
C 9 GAIN TYPE
C10 FILTERS
C11 SOURCE TYPE NUMBER/POINT POINT INTERVAL
C12 PATTERN LENGTH WIDTH
C13 SWEEP START HZ END HZ LENGTH MS CHANNEL NO TYPE
C14 TAPER START LENGTH MS END LENGTH MS TYPE
C15 SPREAD OFFSET MAX DISTANCE GROUP INTERVAL
C16 GEOPHONES PER GROUP SPACING FREQUENCY MFG MODEL
C17 TYPE LENGTH WIDTH
C18 TRACES SORTED BY RECORD PROJECT LINE ID
C19 AMPLITUDE RECOVERY
C20 MAP PROJECTION ZONE ID COORDINATE UNITS
C21 PROCESSING: E&P-SEDE CLAU S070 - CONVERSAO PARA SEGY FEITA EM 1998
C22 PROCESSING: PROJETO MASTER - ORGANIZACAO DO MASTER
C23 I3001V05V 12- 8-94 9.08.05 SYS94224.T090610.RA000.S062B352.TEMP01 D0
C24 S575AC02V 12- 8-94 9.09.52 S0TD.M2S013.E0030494.E239.L0352.MIG01 D0
C25 + INPUT 303 12- 8-94 20.37.33 D0
C26 + RNA 202 FPS 12- 8-94 20.37.33 D0
C27 + TVFILTER 301 FPS 12- 8-94 20.37.33 D0
C28 + GAIN 305 FPS 12- 8-94 20.37.34 D0
C29 + OUTFIX 2A2 12- 8-94 20.37.35 D0
C30 I3001V05V 12- 8-94 20.37.35 S0TD.M2S013.E0030494.E239.L0352.FLT01 D0
C31 + INPUT 303 16- 8-94 8.58.34 D0
C32 + SHALLAVG 201 16- 8-94 8.58.34 D0
C33 + OUTFIX 2A2 16- 8-94 8.58.34 D0
C34 I3001V05V 16- 8-94 8.58.34 S0TD.M2S013.E0030494.E239.L0352.FLT02 D0
C35 P280D14 16- 8-94 12.04.07 S0TD.M2S013.E0030494.E239.L0352.MIGCORD D0
C36 + INPUT 303 23- 9-98 16.35.50 D0
C37 + SEGYSTKW 102 FPS 23- 9-98 16.35.50 D0
C38 + OUTPUT 3C3 23- 9-98 16.35.50 D0
C39 I3001V05V 23- 9-98 16.35.50 SYS98266.T163546.RA000.S0705352.TEMP D0
C40 END EBCDIC 47940039

Desta maneira, foi alcançado a identificação das informações do dado, sem a utilização de um software que requer licenciamento pago da indústria de Petróleo.

Grato,
Rodrigo Paiva
Schlumberger SIS - Brasil

Outras dicas deste autor

Como descobrir qual distribuição Linux que está instalada

Realizar cópia mantendo mesma permissão da pasta

Configuração de uma uma pequena rede RedHat Linux

Leitura recomendada

Ouça o que diz a Vaca

Instalando o LimeWire no Conectiva 10

Instalando o aMSN 0.97-1 no Famelix 2.0

Google Earth 7.1 no Ubuntu 14.04 64 bits (BETA 1)

Rodando o Hamachi como serviço no Linux

Comentários

[1] Comentário enviado por eldermarco em 29/04/2019 - 14:22h

Show! Na época eu fiz isso em Python, mas da maneira como tu fez é bem mais simples. Apenas complementando, o arquivo SEG-Y não possui as quebras de linha no cabeçalho EBCDIC, mas é uma boa ideia inserir elas a cada 80 caracteres, para ficar mais legível. Pode ser feito com o comando:

$ dd if=arquivo.sgy bs=3200 count=1 conv=ascii 2>/dev/null | sed -r 's/(.{80})/\1\n/g'

Valeu pela dica.

3 0