Seismický Unix - Seismic Unix
Seismický Unix je otevřené seismické nástroje balíček, který byl podpořen Centrem pro vlnové jevy (CWP) při Colorado School of Mines (CSM). V současné době ji podporuje John Stockwell.
 Analýza rychlosti pomocí SU | |
| Vývojáři | CWP, John Stockwell |
|---|---|
| Stabilní uvolnění | SU 43R1 / 3. ledna 2012 |
| Operační systém | Unix, Linux, Solaris, Mac OS X, Cygwin |
| Typ | Technické výpočty |
| Licence | Nová licence BSD[1] |
| webová stránka | Seismický unixový domov |
Dějiny
Einar Kjartansson začal psát to, co se nyní nazývá SU (balíček SY), koncem sedmdesátých let, ještě jako postgraduální student Jon Claerbout je Stanfordský průzkumný projekt (SEP). Pokračoval v rozšiřování balíčku, když byl na počátku osmdesátých let profesorem na univerzitě v Utahu. V roce 1984, během delší návštěvy SEP, představil Einar SY Shuki Ronen, tehdejší postgraduální student na Stanfordu. Ronen dále rozvíjel SY v letech 1984 až 1986. Ostatní studenti SEP ho začali používat a přispěli kódem a nápady. SY byl inspirován mnoha dalšími softwary vyvinutými na SEP a těžil ze základů, které položil Claerbout a mnoho jeho studentů; Rob Clayton, Stew Levin, Dave Hale, JeffThorson, Chuck Sword a další, kteří propagovali seismické zpracování na Unixu v sedmdesátých a na počátku osmdesátých let.
V roce 1986 přinesl Shuki Ronen tuto práci na CWP v Colorado School of Mines během svého jednoročního postdoktorského jmenování tam Ronen pomohl Cohenovi přeměnit SU na nepodporovatelný a exportovatelný produkt.
Chris Liner (domovská stránka ), zatímco student ve středisku, přispěl k mnoha grafickým kódům používaným v předpracovní stanici (tj. grafický terminál) věku SU[když? ]. Liner nadále podporuje využívání SU při výzkumu svých studentů na University of Houston.
Craig Artley, nyní s divizí Landmark společnosti Halliburton, významně přispěl ke grafickým kódům ještě jako student CWP a nadále významně přispívá k obecnému balíčku[když? ].
Dave Hale napsal několik kódů zpracování těžkého zvedání a také většinu klíčových vědeckých a grafických knihoven[když? ].
Spolupráce Johna Stockwella se SU začala v roce 1989. Byl do značné míry zodpovědný za Makefile v balíčku. Byl hlavním kontaktem projektu od prvního veřejného vydání SU v září 1992 (vydání 17). Po smrti Jacka Cohena v roce 1996 převzal Stockwell roli hlavního řešitele projektu SU a od té doby v ní zůstal. Počet řádků kódu se za 11 let více než ztrojnásobil.
V uplynulých dvou desetiletích bylo do SU mnoho přispěvatelů.
Syntax
Rutiny Seismic Unix běží pod terminálem Unix a mohou dosáhnout maximální efektivity při použití skriptovacích technik Bourne Shell (sh) nebo Bourne-again Shell (bash).
Jednoduché rutiny
Mnoho programů se spouští jednoduše příkazem na terminálu, například pro vizualizaci seismogramu, jako stopy po kroucení
$ suxwigb nebo jako obrázek
$ suximage Propracovanější rutiny
Je také možné použít funkce bash k vypracování složitějších struktur zpracování:
$ pro ((i=1;i <=100;i ++)); dělat> sufdmod2 output.mov nx=200 nz=300 tmax=5 xs=$ i zs=0 hsfile=seismogram.$ i.su> Hotovo Ve výše uvedeném příkladu Seismic Unix vytvoří 100 seismogramů na 100 různých pozicích zdroje
SU Data
Zde bude vysvětlení toho, jak jsou data SU, jejich záhlaví a jak jsou uspořádána ve velkém souboru SU s více než jedním shromážděním:
- hlavička — data — hlavička — data --...
Programy SU
Seismic Unix má mnoho procesů potřebných pro geofyzikální zpracování. Je možné jej použít k manipulaci a vytváření vlastních seismogramů a také k jejich převodu mezi standardním souborem SU a průmyslovým standardem, SEG Y.
Zde naleznete seznam programů, které balíček SU obsahuje, se stručným popisem a odkazem na jeho stránku nápovědy.[2]
Komprese dat
- Diskrétní kosinová transformace
- dctcomp Komprese diskrétní kosinovou transformací
- dctuncomp Diskrétní kosinová transformace dekomprese
- Balení
- supack1 Zabalte sledovaná data do znaků
- suunpack1 Rozbalte nepříjemná sledovací data od znaků po plováky
- supack2 Sbalte trasovací data do 2 bajtových šortek
- suunpack2 Rozbalte nepříjemná trasovací data od šortek po plováky
- Waveletové transformace
- wpc1comp2 Komprimujte 2D seizmický řez trasování po trasách pomocí Wavelet Packets
- wpc1uncomp2 Dekomprimujte 2D seismický řez, který byl komprimován pomocí Wavelet Packets
- wpccompress Komprimujte 2D sekci pomocí Wavelet paketů
- wpcuncompress Dekomprimujte 2D sekci
- wptcomp Komprese pomocí komprese paketů Wavelet
- wptuncomp Dekomprimujte komprimovaná data WPT
- wtcomp Komprese pomocí vlnkové transformace
- wtuncomp Dekomprese komprimovaných dat WT
Úpravy, třídění a manipulace
- Upravit + nástroje
- suabshw Nahraďte klíčové slovo záhlaví jeho absolutní hodnotou
- suazimuth Vypočítat trasování AZIMUTH vzhledem k polím záhlaví sx, sy, gx, gy a nastavit uživatelem definované pole záhlaví na tuto hodnotu
- podmnožina Vyberte SUBSET vzorků ze trojrozměrného souboru
- sucho Změňte slovo záhlaví pomocí jednoho nebo dvou polí slova záhlaví
- sucountkey Spočítejte počet jedinečných hodnot pro dané klíčové slovo
- suedit Prozkoumejte segy diskety a upravte záhlaví
- sugethw Sugethw zapisuje hodnoty vybraných klíčových slov
- sukill Nulové stopy
- Sunan odebrat NaN & Inf ze vstupního proudu
- suquantile zobrazit některé kvantily nebo řady datové sady
- surange získejte maximální a minimální hodnoty pro nenulové položky záhlaví
- utišit Nastavte jedno nebo více slov v záhlaví pomocí trasovacího čísla, modu a celočíselného dělení pro výpočet hodnot slov v záhlaví nebo zadání hodnot v záhlaví ze souboru
- sutab tisknout nenulové hodnoty záhlaví a data pro negrafické terminály
- suwind stopy okna podle klíčového slova
- suxedit prozkoumat segy diskety a upravit záhlaví
- Třídit
- susort seřadit podle všech klíčových slov segy záhlaví
- susorty vytvořit malou 2-D společnou výstřelovou koncovou datovou sadu, ve které data zobrazují hodnoty geometrie, aby pomohla vizualizovat třídění dat
- Manipulovat
- fcat rychlá kočka s 1 čtením na soubor
- maxdiff najít absolutní maximální rozdíl ve dvou segy datových sadách
- segyhdrmod nahraďte záhlaví textu v souboru SEGY
- suaddnoise přidejte stopám hluk
- jímka CoMPare dvě seismické datové sady, vrátí 0 do shellu, pokud je stejný a 1, pokud se liší
- sudiff, susum, suprod, suquo rozdíl, součet, součin, kvocient dvou souborů dat SU prostřednictvím suop2
- suflip převrátit datovou sadu různými způsoby
- suhtmath proveďte unární aritmetické operace na segy trasách s hodnotami záhlaví
- suinterp interpolovat stopy pomocí automatického výběru událostí
- sumixgathers smíchejte dva shromáždění
- sunull vytvořit nulové (všechny nuly) stopy
- suop dělat unární aritmetické operace na segys
- suop2 proveďte binární operaci na dvou souborech dat
- supermute permutovat nebo transponovat 3D datovou trubku
- suramp Začátek a / nebo konec stop lineárně zužujte na nulu
- surecip součet protichůdných kompenzací v připravených datech
- recip součet protichůdných (vzájemných) offsetů v cdp seřazených datech
- suresamp Převzorkování v čase
- resamp Převzorkování 1. dimenze dvourozměrné funkce f (x1, x2)
- suswapbytes Zaměňte bajty v datech SU a převádějte data z velkého endianového na malé endianové pořadí bajtů a naopak
- šicí stroj Zúžte okrajové stopy datového panelu na nulu
- sutxtaper Zúžte (X, T) okraje datového panelu na nulu
- suvcat připojit jednu sadu dat k jiné, s překrývající se oblastí nebo bez ní. Data v překrytí lze určit jednou z několika metod
- suzero nulová data v časovém okně
- swapbytes Vyměňte bajty různých datových typů
- transp Transponujte matici prvků n1 x n2
Filtrování, transformace a atributy
- Jednorozměrné filtrování
- suband Lichoběžníkovitý Sin hranatý zúžený Bandpass filtr přes SUFILTER
- subfiltrovat použít pásmový filtr Butterworth
- suconv konvoluce s filtrem dodaným uživatelem
- sueipofi Polarizační filtr založený na vlastních obrázcích (SVD) pro třísložková data
- sufilter aplikuje nulovou fázi, sine-square zúžený filtr
- sufrac vezměte obecnou (zlomkovou) časovou derivaci nebo integrál dat plus fázový posun. Vstupem jsou data TIME DOMAIN
- supef Wienerovo prediktivní filtrování chyb
- supofilt POlarizační FILTER pro třísložková data
- sushape Wienerův tvarovací filtr
- sutvband časově variabilní pásmový filtr (sinusový kvadrát)
- suxcor korelace s uživatelem dodaným filtrem
- Dvourozměrné filtrování
- sudipfilt DIP - nebo lepší - SLOPE Filter v doméně f-k
- sufxdecon náhodný útlum šumu pomocí FX-DECONvolution
- filtr suk1k2 symetrický krabicový filtr K-domény definovaný kartézským součinem dvou sin ^ 2-kónických polygonálních filtrů definovaných v k1 a k2
- sumedian MEDIÁNSKÝ filtr o uživatelem definované polygonální křivce se vzdáleností podél křivky určenou klíčovým slovem v záhlaví
- sukfilter radiálně symetrický K-doména, sin ^ 2-zúžený, polygonální filtr
- sukfrac aplikovat FRACionální síly i | k | k datům, s fázovým posunem
- Transformace a atributy
- entropie vypočítat ENTROPY signálu
- mrafxzwt Analýza s více rozlišeními funkce F (X, Z) pomocí vlnkové transformace
- suamp výstupní zesilovač, fáze, skutečné nebo obrazové trasování z dat domény
- doplňkové atributy stopové atributy okamžitá amplituda, fáze nebo frekvence
- suenv Okamžitá amplituda, frekvence a fáze prostřednictvím: atributů
- suhilb Hilbertova transformace
- suhrot Horizontální rotace třísložkových dat
- dost stopové stopy v reálném čase až po složité frekvenční stopy
- suifft fft komplexní frekvenční stopy ke stopám v reálném čase
- sugabor Výstupem je časově-frekvenční reprezentace seismických dat pomocí techniky Gaborovy transformační multifilterové analýzy
- suharlan separace signálu od šumu metodou invertibilní lineární transformace Harlana
- sulog časová osa log-úsek seismických stop
- suilog časová osa inverzní logaritmus seizmických stop
- supolární POLarizační analýza třísložkových dat
- suradon vypočítat dopředu nebo dozadu radonovou transformaci nebo odstranit násobky pomocí parabolické radonové transformace k odhadu násobků a odečtení
- sutaup dopředné a inverzní globální nakloněné zásobníky T-X a F-K
- sutsq časová osa časově čtvercový úsek seismických stop
- přece převést stopy na zobrazení ve zkráceném čase
- suspecfk F-K Fourier SPECTRum souboru dat
- suspecfx Fourierovo SPECTRUM (T -> F) stop
- suspeck1k2 2D (K1, K2) Fourierovo SPECTRUM souboru dat (x1, x2)
Zisk, NMO, zásobník a standardní procesy
- Standardní procesy
- suagc provádět agc na datech SU
- sudipdivcor Oprava divergence (šíření) závislá na dipu
- sudivcor Korekce divergence (šíření)
- sugain použít různé typy zisku k zobrazení stop
- grm Zobecněná vzájemná lomová analýza pro jednu vrstvu
- sumix vypočítat vážený klouzavý průměr (trace MIX) na panelu seismických dat
- součet ztlumit nad (nebo pod) uživatelem definovanou polygonální křivku se vzdáleností podél křivky určenou klíčovým slovem v záhlaví
- sunmo NMO pro libovolnou rychlostní funkci času a CDP
- supgc Programmed Gain Control - apply agc like function but the same function to all traces zachování relativní amplitudy prostorově
- prosby Fázový zásobník nebo fázově vážený zásobník (PWS) sousedních stop, které mají stejné slovo záhlaví klíče
- suresstat Výpočet statiky zdroje a přijímače konzistentní s povrchem
- sustack skládat sousední stopy mající stejné slovo záhlaví klíče
- sustatický Statické korekce výšky, použít korekce ze záhlaví nebo ze souboru statiky zdroje a přijímače
- sustaticrrs Korekce STATIC nadmořské výšky, použití korekcí ze záhlaví nebo ze souboru statiky zdroje a přijímače, zahrnuje použití statiky Residual Refraction Statics
- rozčarování Clip outliers in data
- Smíšený
- suacor Automatická korelace
- suttoz Převzorkování z času do hloubky
- suvibro Generuje rozmítání Vibroseis (lineární, lineární segment, dB na oktávu, dB na hertz, výkon T)
- délka Upravte stopy proměnné délky na běžnou délku
SU komunita
Seismic Unix má velmi velkou komunitu, kterou využívá mnoho laboratoří a výzkumných pracovníků po celém světě. Existuje skupina seznamových serverů, kde můžete získat pomoc a zveřejnit své dotazy týkající se SU. Chcete-li se přihlásit k jeho odběru, klikněte na tady.
Můžete také zobrazit staré příspěvky, o kterých uživatelé již hovořili, kliknout tady
Nový web Seismic Unix od Dominique Rousset v roce 2019, klikněte sem tady
Ocenění
2002 - Společnost geofyziků průzkumu Zvláštní uznání [3]
1994 - Cena University to Industry od kapitoly Colorado společnosti Transfer technologií [4]
Reference
Viz také
- Průvodce Instalační průvodce seismickým Unixem (hostovaný na brazilské stránce)
- Seismická Unix wiki
- Nová seismická Unix wiki 2019