- Seisunix - mailman.seismic-unix.org

Re: Problems with sucdpbin
by Dominique ROUSSET 07 Apr '23

07 Apr '23

Hi, good... but not yet very good ! I suspect a problem of numerical accuracy with such high values of X and Y coordinates. You can test this hypothesis by subtracting the mean of X and Y coordinates to have 0 centered data... If it's true, it should be better to correct the algorithm, I know... D. Le 07/04/2023 à 08:41, JUAN LORENZO DIAZ NAVEAS a écrit : > Dear Glenn and Dominique, > > Thanks a lot for your help. Indeed the keyword scalco=-10000 was > guilty for scaling sx,sy,gx and gy coordinates. So now I am able to > produce an su file with cdp headers. > > But now unfortunately I found another problem. The su file with cdp > headers has about 66% of cdp=0. The cdp=0 traces are at the > "beginning" of the su file. > > And the intriguing aspect is that I have a crooked line with 10 > points. See the attached map. There the black dots are the source > coordinates, the red dots are the receiver coordinates and the green > circles are the points of the crooked line. > > As you can guess I smoothed both the source and receiver coordinates > using the GMT trend1d command. > > The command line for my binning looks like this: > > *sucdpbin > xline=723193.2,744776.011111,766358.822222,787941.633333,809524.444444,831107.255556,852690.066667,874272.877778,895855.688889,920934 > \ > yline=7242472.4,7242099.3532,7241693.29184,7241254.1793,7240782.11164,7240277.08676,7239739.07882,7239168.07407,7238564.07513,7237820 > \ > < SO104-17_geom.su verbose=2 dcdp=25 cdpmin=1 distmax=100 > > SO104-17_geom-cmp.su > * > > See below the ranges of the file SO104-17_geom-cmp.su and the > attached log file with the verbosity of the sucdpbin command. > > I tried also with distmax less than 100m, the smallest value being > 15m, but the result was similar. > > Any idea on what is going on? > > Best wishes > > Juan Díaz Naveas > ------------------------------------------------------------------- > Juan Diaz Naveas (Dr. rer. nat.) > Escuela de Ciencias del Mar www.cienciasdelmar.pucv.cl > <http://www.cienciasdelmar.pucv.cl> > Pontificia Universidad Catolica de Valparaiso Tel.: +56-32-2274269 > Avenida Universidad 330, Curauma Cel.: +56-9-91630936 > Valparaiso, CHILE E-Mail: juan.diaz(a)pucv.cl > ------------------------------------------------------------------- > *466560 traces: > tracl 1 513204 (1 - 513204) > tracr 1 513204 (1 - 513204) > fldr 1 3888 (1 - 3888) > tracf 1 120 (1 - 120) > ep 1 3888 (1 - 3888) > cdp 0 2572 (0 - 14) > trid 1 > nhs 1 > offset 413 3388 (3388 - 413) > scalel -10000 > scalco 1 > sx 723302 917547 (917547 - 723302) > sy 7237924 7242471 (7237924 - 7242471) > gx 723715 920934 (920934 - 723715) > gy 7237823 7242464 (7237823 - 7242464) > counit 1 > ns 3585 > dt 4000 > gain 9 > afilf 102 > afils 209 > lcf 3 > hcf 102 > lcs 6 > hcs 209 > year 1995 > > Shot coordinate limits: > North(723302,7.24247e+06) South(917547,7.23792e+06) > East(917547,7.23792e+06) West(723302,7.24247e+06) > > Receiver coordinate limits: > North(723740,7.24246e+06) South(920934,7.23782e+06) > East(920934,7.23782e+06) West(723715,7.24246e+06) > > Midpoint coordinate limits: > North(723521,7.24247e+06) South(919240,7.23787e+06) > East(919240,7.23787e+06) West(723508,7.24247e+06)* > > El jue, 6 abr 2023 a las 3:10, Dominique ROUSSET > (<dominique.rousset(a)univ-pau.fr>) escribió: > > Hi Juan, > > after a quick look at sucdpbin.c, i note that this tools honors > the scalco keyword. So, since your data bear a -10000 scalco, you > have to divide xline and yline parameters by 10000 to work with > actual units, not with sx and sy units. > > Let me know it it works > Good luck ! > > D. > > Le 06/04/2023 à 02:47, JUAN LORENZO DIAZ NAVEAS a écrit : >> Dear SU users, >> >> I am trying to make the CMP binning of a seismic line (see the >> headers ranges below). >> >> The command line is: >> >> *sucdpbin < so104-17_geom.su <http://so104-17_geom.su> >> xline=723715,920883 yline=7242464,7237825 verbose=2 dcdp=25 >> distmax=500 > so104-17_cdp.su <http://so104-17_cdp.su>* >> >> As you see I have used a distmax of 500m and still I am not able >> to get the binning. >> >> Just in case I also tried again the same command, but sorting the >> date so that the first traces in the su file correspond to the >> minimum x coordinate and the last traces to the maximum >> coordinate (in the field the line was shot from east to west), >> but it failed as well. >> >> Any help would be very much appreciated. >> >> Best wishes >> >> Juan >> >> *466560 traces: >> tracl 1 513204 (1 - 513204) >> tracr 1 513204 (1 - 513204) >> fldr 1 3888 (1 - 3888) >> tracf 1 120 (1 - 120) >> ep 1 3888 (1 - 3888) >> trid 1 >> nhs 1 >> offset 413 3388 (3388 - 413) >> scalel -10000 >> scalco -10000 >> sx 723302 917547 (917547 - 723302) >> sy 7237924 7242471 (7237924 - 7242471) >> gx 723715 920934 (920934 - 723715) >> gy 7237823 7242464 (7237823 - 7242464) >> counit 1 >> ns 3585 >> dt 4000 >> gain 9 >> afilf 102 >> afils 209 >> lcf 3 >> hcf 102 >> lcs 6 >> hcs 209 >> year 1995 >> >> Shot coordinate limits: >> North(-91.7547,-723.792) South(-72.3302,-724.247) >> East(-72.3302,-724.247) West(-91.7547,-723.792) >> >> Receiver coordinate limits: >> North(-92.0934,-723.782) South(-72.374,-724.246) >> East(-72.3715,-724.246) West(-92.0934,-723.782) >> >> Midpoint coordinate limits: >> North(-91.924,-723.787) South(-72.3521,-724.247) >> East(-72.3509,-724.247) West(-91.924,-723.787)* >> ------------------------------------------------------------------- >> Juan Diaz Naveas (Dr. rer. nat.) >> Escuela de Ciencias del Mar www.cienciasdelmar.pucv.cl >> <http://www.cienciasdelmar.pucv.cl> >> Pontificia Universidad Catolica de Valparaiso Tel.: +56-32-2274269 >> Avenida Universidad 330, Curauma Cel.: +56-9-91630936 >> Valparaiso, CHILE E-Mail: juan.diaz(a)pucv.cl >> ------------------------------------------------------------------- >> >> _______________________________________________ >> Seisunix mailing list --seisunix(a)mailman.seismic-unix.org >> To unsubscribe send an email toseisunix-leave(a)mailman.seismic-unix.org > > -- > |Dominique Rousset LFCR UMR5150 / E2S UPPA / Dpt des géosciences| > |Univ. Pau et des Pays de l'Adour BP 1155 F-64013 Pau Cedex - France| > |mailto:dominique.rousset@univ-pau.fr <mailto:dominique.rousset@univ-pau.fr> http://lfc.univ-pau.fr| > | Off: +33 559 407 423 Cell: +33 660 861 391 | > > _______________________________________________ > Seisunix mailing list -- seisunix(a)mailman.seismic-unix.org > To unsubscribe send an email to > seisunix-leave(a)mailman.seismic-unix.org > > > _______________________________________________ > Seisunix mailing list --seisunix(a)mailman.seismic-unix.org > To unsubscribe send an email toseisunix-leave(a)mailman.seismic-unix.org -- |Dominique Rousset LFCR UMR5150 / E2S UPPA / Dpt des géosciences| |Univ. Pau et des Pays de l'Adour BP 1155 F-64013 Pau Cedex - France| |mailto:dominique.rousset@univ-pau.fr http://lfc.univ-pau.fr| | Off: +33 559 407 423 Cell: +33 660 861 391 |

1 0

Problems with sucdpbin
by JUAN LORENZO DIAZ NAVEAS 06 Apr '23

06 Apr '23

Dear SU users, I am trying to make the CMP binning of a seismic line (see the headers ranges below). The command line is: *sucdpbin < so104-17_geom.su <http://so104-17_geom.su> xline=723715,920883 yline=7242464,7237825 verbose=2 dcdp=25 distmax=500 > so104-17_cdp.su <http://so104-17_cdp.su>* As you see I have used a distmax of 500m and still I am not able to get the binning. Just in case I also tried again the same command, but sorting the date so that the first traces in the su file correspond to the minimum x coordinate and the last traces to the maximum coordinate (in the field the line was shot from east to west), but it failed as well. Any help would be very much appreciated. Best wishes Juan *466560 traces:tracl 1 513204 (1 - 513204)tracr 1 513204 (1 - 513204)fldr 1 3888 (1 - 3888)tracf 1 120 (1 - 120)ep 1 3888 (1 - 3888)trid 1nhs 1offset 413 3388 (3388 - 413)scalel -10000scalco -10000sx 723302 917547 (917547 - 723302)sy 7237924 7242471 (7237924 - 7242471)gx 723715 920934 (920934 - 723715)gy 7237823 7242464 (7237823 - 7242464)counit 1ns 3585dt 4000gain 9afilf 102afils 209lcf 3hcf 102lcs 6hcs 209year 1995Shot coordinate limits: North(-91.7547,-723.792) South(-72.3302,-724.247) East(-72.3302,-724.247) West(-91.7547,-723.792)Receiver coordinate limits: North(-92.0934,-723.782) South(-72.374,-724.246) East(-72.3715,-724.246) West(-92.0934,-723.782)Midpoint coordinate limits: North(-91.924,-723.787) South(-72.3521,-724.247) East(-72.3509,-724.247) West(-91.924,-723.787)* ------------------------------------------------------------------- Juan Diaz Naveas (Dr. rer. nat.) Escuela de Ciencias del Mar www.cienciasdelmar.pucv.cl Pontificia Universidad Catolica de Valparaiso Tel.: +56-32-2274269 Avenida Universidad 330, Curauma Cel.: +56-9-91630936 Valparaiso, CHILE E-Mail: juan.diaz(a)pucv.cl -------------------------------------------------------------------

2 1

Announcing Incremental release 44R28 of Seismic Unix
by John Stockwell 23 Mar '23

23 Mar '23

This release does contain a fix for the “types.h” (tirpc) installation problem seen by users on Ubuntu 22.04. I have also started adding platform specific installation instructions on https://wiki.seismic-unix.org <https://wiki.seismic-unix.org/> under Documentation. Right now this includes a list of all the packages that you need to install under Ubuntu 22.04 before you can install SU.

1 0

Compilation probles wis dist CWP-44.26
by Ramon Carbonell 23 Mar '23

23 Mar '23

Dear All, Has anybody struggled with compiling version 44.26 and having trouble when compiling with the xdr I am getting messages like src/su/lib/fgettr.c:148: undefined reference to `xdr_float' src/su/lib/fputtr.c:123: undefined reference to `xdr_vector' I am running Ubuntu 22.04 and somehow, although xdr.h is installed at rpc/xdr.h I am not sure why the compiler does not locate it. Anu suggestions. -- Ramon -------------------------------------------------------------------------------------------- Ramon Carbonell CSIC-GEO3BCN GeoSciences Barcelona (http://Geo3BCN.csic.es/) https://orcid.org/0000-0003-2019-1214

8 16

Fwd: Re: Compilation probles wis dist CWP-44.26
by Dominique Rousset 14 Mar '23

14 Mar '23

-------- Courriel d’origine -------- De : Paul Fertser <fercerpav(a)gmail.com> Envoyé : 14 mars 2023 21:09:01 GMT+01:00 À : Ramon Carbonell <ramon.carbonell(a)csic.es> Cc : seisunix(a)seismic-unix.org Objet : Re: [Seisunix] Compilation probles wis dist CWP-44.26 Hello, On Tue, Mar 14, 2023 at 06:16:14PM +0100, Ramon Carbonell wrote: > compiling with the xdr > I am getting messages like > src/su/lib/fgettr.c:148: undefined reference to `xdr_float' > src/su/lib/fputtr.c:123: undefined reference to `xdr_vector' > I am running Ubuntu 22.04 and somehow, although xdr.h is installed at rpc/xdr.h > I am not sure why the compiler does not locate it. That's because it's a linking error, the compiler found your include but the linker doesn't know it needs to use an additional library. This change makes it buildable on current Ubuntu LTS: diff --git a/src/Makefile.config b/src/Makefile.config index a1b8699c11e2..f7cdd4feab08 100644 --- a/src/Makefile.config +++ b/src/Makefile.config @@ -68,7 +68,7 @@ RANFLAGS = ICHMODLINE = chmod 644 $@ MCHMODLINE = chmod 755 $@ -POSTLFLAGS = +POSTLFLAGS = $(shell pkg-config --libs libtirpc) #----------------------------------------------------------------------- # use both X11 path conventions @@ -92,8 +92,8 @@ CPP = /usr/bin/cpp CC = /usr/bin/gcc #OPTC = -g -c90 -m64 -Wall -ansi -Wno-long-long #OPTC = -g -std=c90 -m64 -Wall -ansi -Wno-long-long -OPTC = -g -std=c99 -m64 -Wall -ansi -Wno-long-long -CFLAGS = -I$I $(OPTC) $(CWP_FLAGS) +OPTC = -g -std=c99 -m64 -Wall -ansi -Wno-long-long -D_ISOC99_SOURCE +CFLAGS = -I$I $(OPTC) $(CWP_FLAGS) $(shell pkg-config --cflags libtirpc) # if gfortran is made using Homebrew FC = gfortran -- Be free, use free (http://www.gnu.org/philosophy/free-sw.html) software! mailto:fercerpav@gmail.com -- Envoyé de mon appareil Android avec Courriel K-9 Mail. Veuillez excuser ma brièveté.

1 0

Help with fixing Trace number
by eloraafrin5＠gmail.com 08 Dec '22

08 Dec '22

Hello greetings! I am a beginner. I have a big su file with huge shot points, but I want to see a particular shot point in reduced velocity mode. I made a single file for those particular shot by using key=fldr. Now this file should have total 636 traces. I wrote a command suwind < indut su file key=tracl max=636 ....... However, in image I got values from2.84*10^5 to 2.85*10^5. Even when I want to see the image in reduced velocity mode, I need to set rv= >1000 whereas it should be rv=1 to 2. It seems like there are something wrong with units in the headers. Would anyone please help me with how can I fix these??? I am a new user; I will appreciate your help and time. Here I attached the header info. Thanks FYI 636 traces: tracl 1 636 (1-636) tracr 284929 285564 (284929-285564) fldr 2245 tracf 1 636 (1-636) cdp 2245 offset 112464 120401 (112464-120401) sx 376962 sy 5176643 gx 376962 gy 5176643 ns 2001 dt 4000

6 6

Inverse Scattering Series for inter-bed multiple attenuation
by Ameur Hamdane 05 Dec '22

05 Dec '22

Hello SU users; Is there a written Inverse Scattering Series interbed multiple attenuation in Seismic Unix, I haven't used the software in a decade but I don't remember seeing one. I need to try something using this method,please advise? Thanks alot

1 0

How to handle trace header information?
by tanahashi.m＠gmail.com 11 Nov '22

11 Nov '22

Dear all, I have performed a time resample of some 500ns sampled 3D seismic data to 1ms intervals using the suresamp command as shown below, but was unable to transcribe the 189-192 bytes of the Inline and the 192-196 bytes of the crossline numbers. segyread tape=Source3D.sgy|suresamp rf=0.5|segyhdrs|segywrite tape=Output3D_1ms.sgy I solved this problem by using the ASCII input/output function for trace header on SeiSee, a free Windows software, to write out the header information of the original file and then write it into the header of the output file. Is there a way to transcribe inline and crossline numbers using the SU function? Manabu Tanahashi Tsukuba, Japan

3 5

Help with Adding Offset Header in a SU File
by eloraafrin5＠gmail.com 27 Oct '22

27 Oct '22

Hello, I am a beginner. I need to apply sureduce command to a particular shot which is in su file format. I tried with different rv values but could not see any differences in the image. I checked with the headers of that file and found I have ns and dt headers but not offset. Can anyone please help me how can I add offset header to my su file? And please let me know if there are any other things to fix to get result. Thanks

2 1

test message. Discard. Sorry
by Dominique ROUSSET 19 Oct '22

19 Oct '22

testing archiving of new messages -- |Dominique Rousset LFCR UMR5150 / E2S UPPA / Dpt des géosciences| |Univ. Pau et des Pays de l'Adour BP 1155 F-64013 Pau Cedex - France| |mailto:dominique.rousset@univ-pau.fr http://lfc.univ-pau.fr| | Off: +33 559 407 423 Cell: +33 660 861 391 |

3 3

2025

2024

2023

2022

2021

2020

2019

Seisunix