CVSPH0 - DYNAMIC TIME SHIFT AND GATHER 00010000 C*********************************************************************** C COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991 * C*********************************************************************** C 00020000 C DESIGNER J E GAISER 00030000 C AUTHOR J E GAISER 00040000 C LANGUAGE S/370 FORTRAN H 00050000 C DATE 06-06-82 00060000 C 00070000 C 00080000 C CALL VSPI0(OH,ICC,AUTO3,IABORT,RA) 00090000 C CALL VSPI1(OH,OTR,VEL,PASS,IABORT,RA,SA) 00100000 C CALL VSPI2(OH,OTR,VEL,PASS,IABORT,RA,SA) 00110000 C CALL VSPI3(OH,OTR,VEL,PASS,IABORT,RA,SA) 00120000 C 00130000 C 00140000 C THIS SPARC DEVELOPMENT PROGRAM COMPUTES AND APPLIES NMO TIME SHIFTS00150000 C TO CORRECT VSPS TO TWO-WAY VERTICAL TIME. SHIFTING IS BASED ON SOURCE00160000 C RECEIVER GEOMETRY AND STRATIGRAPHIC DIP. INPUT DATA MUST HAVE 00170000 C VERTICAL TIMES AND DEPTHS IN THE HEADERS AND THE TRACE INPUT ORDER 00180000 C MUST BE FROM TD TO SHALLOW. 00190000 C 00200000 C 00210000 C VSPI 00220000 C ===================================================================== 00230000 CU PROCESS VSPI -- VSP OFFSET STACKING 00240000 CU DATA CARD (1) -- DEFINES S OR D RANGE AND PLOT TYPE 00250000 CU 00260000 CU NO. OF CARDS: REQUIRED = 1 ALLOWED = 1 00270000 CU 00280000 CU REQ OR OPT 00290000 CU DF COLS DESCRIPTION OR DEFAULT 00300000 CU -- ----- ----------- -----------00310000 CU 1 1- 4 'VSPI' | REQ |00320000 CU 2 - 5 PROCESS NUMBER | 0 |00330000 CU 3 - 6 NOT USED | |00340000 CU 4 - 7 PROCESSING MODE |LINE CARD|00350000 CU 'S' = SHOT POINT MODE | |00360000 CU 'D' = DEPTH POINT MODE | |00370000 CU 5 8-10 FOR PROGRAMMING USE | |00380000 CU 6 11-15 STARTING SHOTPOINT, DEPTH POINT | REQ |00390000 CU 7 16-20 ENDING SHOTPOINT, DEPTH POINT | DF6 |00400000 CU 8 21-25 NOT USED | |00410000 CU 9 26-35 TEN COLUMNS, FIELDS 9 AND 10 | NEXT NO.|00420000 CU 10 DATA SET IDENTIFICATION | |00430000 CU NORMALLY BLANK--SYSTEM ASSIGNS NEXT NUMBER | |00440000 CU 11 36-40 DATA SET TYPE | T |00450000 CU 'T' = TAPE WITH IBM STANDARD LABEL | |00460000 CU 'N' = TAPE WITHOUT IBM STANDARD LABEL | |00470000 CU 12 41-45 TAPE RECORDING DENSITY (NOT USED FOR DISK) | H FOR T |00480000 CU 'L' = 800 BPI, 'M' = 1600 BPI, 'H' = 6250 BPI | L FOR N |00490000 CU 13 46-50 NOT USED | |00500000 CU 14 51-55 ONLINE PLOTTER CODE. SEE OPERATIONS FOR CODE. | |00510000 CU 15 56-60 NOT USED | |00520000 CU 16 61-65 NOT USED | |00530000 CU 17 66-70 NOT USED | |00540000 CU 18 71-75 NOT USED | |00550000 CU 19 76-80 'ESP' = NO DISPLAY | ESP |00560000 CU 'OFESP' = OFFLINE VARIAN DISPLAY | |00570000 CU 'ONESP' = ONLINE DISPLAY | |00580000 CU 'OFTAD' = OFFLINE TIADD DISPLAY | |00590000 CU -----------00600000 CU DF NOTES 00610000 CU -- ----- 00620000 CU 00630000 CU 1 ONLY ONE PROCESS OCCURRENCE IS ALLOWED PER JOB. 00640000 CU THIS PROCESS CORRESPONDS TO RESEARCH PROGRAM 'GEOA'. 00650000 CU 00660000 CU 19 OUTPUT IS IN VECTOR FORMAT. 'OFESP' IMPLIES THAT THE OUTPUT 00670000 CU WILL BE PLOTTED OFFLINE USING THE ELECTROSTATIC PAPER PLOTTER. 00680000 CU 'ONESP' IMPLIES THAT THE OUTPUT WILL GO TO A DISK FILE FOR 00690000 CU PLOTTING ON AN ONLINE VARIAN PLOTTER. 00700000 CU 'OFTAD' IMPLIES THAT THE OUTPUT WILL BE PLOTTED OFFLINE USING 00710000 CU THE TIADD PLOTTER. 00720000 CU 00730000 CU EJECT 00740000 C PROCESS VSPI -- VSP OFFSET STACKING 00750000 CU DATA CARD (2) -- PROCESSING PARAMETERS 00760000 CU 00770000 CU NO. OF CARDS: REQUIRED = 1 ALLOWED = 1 00780000 CU 00790000 CU REQ OR OPT 00800000 CU DF COLS DEFINITION OR DEFAULT 00810000 CU -- ----- ---------- -----------00820000 CU 1 1- 4 'VSPI' | REQ |00830000 CU 2 - 5 PROCESS NUMBER | 0 |00840000 CU 3 - 6 NOT USED | |00850000 CU 4 - 7 NOT USED | |00860000 CU 5 8-10 'ANA' | REQ |00870000 CU 6 11-15 NOT USED | |00880000 CU 7 16-20 NOT USED | |00890000 CU 8 21-25 CDP INTERVAL (FEET) | 110 |00900000 CU 9 26-30 CDP WIDTH (FEET) | DF8 |00910000 CU 10 31-35 CDP OVERLAP (PERCENT OF CDP INTERVAL) | 0 |00920000 CU 11 36-40 TAPER LENGTH (PERCENT OF CDP INTERVAL) | 0 |00930000 CU 41-80 NOT USED | |00940000 CU -----------00950000 CU 00960000 CU DF NOTES 00970000 CU -- ----- 00980000 CU 00990000 CU 8 DEFINES THE DEPTH POINT INTERVAL IN FEET. THE FIRST CDP WILL 01000000 CU BE HALF OF DF8 AWAY FROM THE B.O.L. POINT (ON CARD 3) IN THE 01010000 CU DIRECTION TOWARDS THE E.O.L. POINT. THE NEXT DEPTH POINT WILL 01020000 CU THEN BE DF8 FEET FROM THE FIRST AND SO ON FOR THE THIRD CDP. 01030000 CU 01040000 CU 9 THE DISTANCE PERPENDICULAR TO THE LINE TO INCLUDE IN THE CDP. 01050000 CU THIS ALSO DEFINES THE LINE SPACING WHEN STACKING MULTIPLE LINES. 01060000 CU THE DEFAULT IS THE CDP INTERVAL. 01070000 CU 01080000 CU 10 THIS AMOUNT WILL BE ADDED TO EACH SIDE OF THE CDP INTERVAL AND 01090000 CU CDP WIDTH. THE TOTAL LENGTH OF THE CDP INTERVAL OR CDP WIDTH 01100000 CU CAN NOT EXCEED 1000 FEET. 01110000 CU 01120000 CU 11 COSINE TAPER APPLIED TO THE SAMPLES TO SMOOTH CDP EDGE EFFECTS. 01130000 CU (THIS WILL REDUCE THE EFFECTIVE CDP INTERVAL BY ABOUT HALF 01140000 CU THE TAPER LENGTH). 01150000 C 01160000 C 01170000 CU EJECT 01180000 C ======================================================================01190000 C PROCESS VSPI -- VSP OFFSET STACKING 01200000 CU DATA CARD (3) -- DEFINES BEGINNING OF LINE (B.O.L.) 01210000 CU AND END OF LINE (E.O.L.) COORDINATES 01220000 CU 01230000 CU NO. OF CARDS: REQUIRED = 0 ALLOWED = 1 01240000 CU 01250000 CU REQ OR OPT 01260000 CU DF COLS DEFINITION OR DEFAULT 01270000 CU -- ----- ---------- -----------01280000 CU 1 1- 4 'VSPI' | REQ |01290000 CU 2 - 5 PROCESS NUMBER | 0 |01300000 CU 3 - 6 NOT USED | |01310000 CU 4 - 7 NOT USED | |01320000 CU 5 8-10 'LIN' | REQ |01330000 CU 6 11-15 NOT USED | |01340000 CU 7 16-20 NUMBER OF LINES TO STACK | 1 |01350000 CU 8 21-30 B.O.L. X COORDINATE OF LINE 1(NORTH IS POSITIVE)| TD COOR |01360000 CU |NOTE DF8|01370000 CU 10 31-40 B.O.L. Y COORDINATE OF LINE 1(EAST IS POSITIVE)| TD COOR |01380000 CU 12 41-50 E.O.L. X COORDINATE OF LINE 1(NORTH IS POSITIVE)| 0 |01390000 CU 14 51-60 E.O.L. Y COORDINATE OF LINE 1(EAST IS POSITIVE)| 0 |01400000 CU 16 61-80 NOT USED | |01410000 CU -----------01420000 CU 01430000 CU DF NOTES 01440000 CU -- ----- 01450000 CU 01460000 CU 7 THIS DEFINES THE NUMBER OF 3-D LINES TO STACK. EACH LINE IS 01470000 CU PARALLEL TO LINE 1 DEFINES IN DF8-DF14 AND FORMS A RECTANGULAR 01480000 CU GRID. THE ADDITIONAL LINES HAVE CDP COORDINATES THAT ARE 01490000 CU COUNTER CLOCKWISE FROM LINE 1 WHEN B.O.L. X AND B.O.L. Y ARE 01500000 CU VIEWED AS THE ORIGIN. 01510000 CU 01520000 CU 8- DEFINES START AND END OF THE FIRST LINE. CAN BE USED TO REVERSE 01530000 CU 14 THE DIRECTION OF THE LINE NUMBERS. IF CDP 1 IS TO BE NEAR T.D., 01540000 CU CODE B.O.L. COORDINATES NEAR T.D. X-Y COORDINATES; IF CDP 1 IS 01550000 CU TO BE NEAR THE WELLHEAD, CODE B.O.L. NEAR 0. THE DEFAULTS ARE 01560000 CU THE COORDINATES OF THE FIRST TRACE (USUALLY TD) FOR B.O.L. AND 01570000 CU ZERO (WELLHEAD) FOR E.O.L. 01580000 CU EJECT 01590000 C ======================================================================01600000 C PROCESS VSPI -- VSP OFFSET STACKING 01610000 CU DATA CARD (4) -- CDP PLOT PARAMETERS FOR DISPLAYING THE CDP 01620000 CU COVERAGE IN THE X-Y PLANE RELATIVE TO THE 01630000 CU WELLHEAD AND THE BOREHOLE 01640000 CU 01650000 CU NO. OF CARDS: REQUIRED = 0 ALLOWED = 1 01660000 CU 01670000 CU REQ OR OPT 01680000 CU DF COLS DEFINITION OR DEFAULT 01690000 CU -- ----- ---------- -----------01700000 CU 1 1- 4 'VSPI' | REQ |01710000 CU 2 - 5 PROCESS NUMBER | 0 |01720000 CU 3 - 6 NOT USED | |01730000 CU 4 - 7 NOT USED | |01740000 CU 5 8-10 'PLT' | REQ |01750000 CU 6 11-15 ' XY' | REQ |01760000 CU 7 16-20 HORIZONTAL SCALE (FT/INCH) FLOATING POINT | 100.0 |01770000 CU 8 21-25 GROUP INCREMENT FOR PLOTTING REFLECTION POINTS | 5 |01780000 CU 9 26-30 FIRST DEPTH TO ANNOTATE ON PLOT | TD |01790000 CU 10 31-35 NEXT DEPTH OF 2 | OPT |01800000 CU 11 36-40 NEXT DEPTH OF 3 | OPT |01810000 CU 12 41-45 NEXT DEPTH OF 4 | OPT |01820000 CU 13 46-50 NEXT DEPTH OF 5 | OPT |01830000 CU 51-80 NOT USED | |01840000 CU -----------01850000 CU 01860000 CU DF NOTES 01870000 CU -- ----- 01880000 CU 01890000 CU 8 THIS DEFINES THE INCREMENT FOR PLOTTING RECEIVER LOCATIONS AND 01900000 CU THE REFLECTION COVERAGE FOR EACH RECEIVER. E.G., A VALUE OF 2 01910000 CU CAUSES EVERY OTHER RECEIVER TO BE PLOTTED. 01920000 CU 01930000 CU 9 REFLECTION COVERAGE FOR 5 DIFFERENT HORIZONS CAN BE PLOTTED. 01940000 CU THE DEFAULT IS FOR THE DEEPEST RECEIVER (FIRST TRACE INPUT). 01950000 CU THE EXACT TVD DEPTH OF THE RECEIVER MUST BE CODED. 01960000 CU EJECT 01970000 C ======================================================================01980000 C 01990000 C PROCESS VSPI -- VSP OFFSET STACKING 02000000 CU DATA CARD (5) -- PLOT LABEL 02010000 CU 02020000 CU NO. OF CARDS: REQUIRED = 0 ALLOWED = 1 02030000 CU 02040000 CU REQ OR OPT 02050000 CU DF COLS DEFINITION OR DEFAULT 02060000 CU -- ----- ---------- -----------02070000 CU 1 1- 4 'VSPI' | REQ |02080000 CU 2 - 5 PROCESS NUMBER | 0 |02090000 CU 3 - 6 NOT USED | |02100000 CU 4 - 7 NOT USED | |02110000 CU 5 8-10 'PLT' | REQ |02120000 CU 6 11-15 NOT USED | |02130000 CU 7 16-20 NOT USED | |02140000 CU 8 21-68 LABEL FOR PLOT | |02150000 CU 69-80 NOT USED | |02160000 CU -----------02170000 CUEND 02180000 C ======================================================================02190000 C 02200000 C 02210000 C 02220000 C 02230000 C ======================================================================02240000 C EJECT 02250000 C 02260000 C 02270000 C VSPI 02280000 C MEMORY AND SPACE ALLOCATIONS FOR SPARC DEVELOPMENT 02290000 C PROCESS VSPI -- VSPI 02300000 C DATA CARD (6) -- DEFINES INPUT PARAMETERS MUST BE LAST CARD 02310000 C 02320000 C NO. OF CARDS: REQUIRED = 1 02330000 C 02340000 C 02350000 C DF COLS DEFINITION REQ OR OPT02360000 C -- ----- ----------- -----------02370000 C 1 1-4 'VSPI' | REQ |02380000 C 2 -5 PROCESS NUMBER | 0 |02390000 C 3 6-7 NOT USED | |02400000 C 4 8-10 'SPC' | REQ |02410000 C 5 11-15 PROGRAM MEMORY SIZE (K-BYTES) | 500 |02420000 C 6 16-20 BLANK COMMON RESERVED SPACE (WORDS) | 0 |02430000 C 7 21-25 BLANK COMMON SCRATCH SPACE (WORDS) | 0 |02440000 C 8 26-30 NUMBER OF GROUPS (WORKFILE 1) | REQ |02450000 C 9 31-35 NUMBER OF TRACES/GROUP (WORKFILE 1) | 1 |02460000 C 10 36-40 RECORD LENGTH (WORDS) (WORKFILE 1) |HDR + TR |02470000 C 11 41-45 NUMBER OF GROUPS (WORKFILE 2) | REQ |02480000 C MUST CODE A '1' IN THIS DATA FIELD | |02490000 C 12 46-50 NUMBER OF TRACES/GROUP (WORKFILE 2) | 1 |02500000 C 13 51-55 RECORD LENGTH (WORDS) (WORKFILE 2) |HDR + TR |02510000 C 14 56-60 NUMBER OF GROUPS (WORKFILE 3) | REQ |02520000 C MUST CODE A '1' IN THIS DATA FIELD | |02530000 C 15 61-65 NUMBER OF TRACES/GROUP (WORKFILE 3) | 1 |02540000 C 16 66-70 RECORD LENGTH (WORDS) (WORKFILE 3) |HDR + TR |02550000 C -----------02560000 C 02570000 C 02580000 C=======================================================================02590000 C 02600000 C VARIABLES IN RESERVED AREA RA 02610000 C 02620000 C WIENER INTERPOLATION FILTERS COEF(12,99) (R4) 02630000 C RA ( 1) . . . . . COEF ( 1, 1) 02640000 C RA ( 12) . . . . . COEF (12, 1) 02650000 C RA ( 13) . . . . . COEF ( 1, 2) 02660000 C RA (1188) . . . . . COEF (12,99) 02670000 C 02680000 C SHOTPOINTS/DEPTH POINTS FOR EDITING DELET(100) (I4) 02690000 C RA (1189) . . . . . DELET ( 1) 02700000 C RA (1288) . . . . . DELET (100) 02710000 C 02720000 C NOT USED 02730000 C RA (1289) 02740000 C RA (1300) 02750000 C 02760000 C PRINTOUT ARRAYS 02770000 C 02780000 C SHOTPOINT/DEPTH POINT NUMBER LSHT(12) (I4) 02790000 C RA (1301) . . . . . LSHT ( 1) 02800000 C RA (1312) . . . . . LSHT (12) 02810000 C 02820000 C TRACE NUMBER LTRC(12) (I4) 02830000 C RA (1313) . . . . . LTRC ( 1) 02840000 C RA (1324) . . . . . LTRC (12) 02850000 C 02860000 C STATIC SHIFT LSHF(12) (R4) 02870000 C RA (1325) . . . . . LSHF ( 1) 02880000 C RA (1336) . . . . . LSHF (12) 02890000 C 02900000 C CRITICAL ANGLE FLAG LGTN(12) (A4) 02910000 C RA (1337) . . . . . LGTN ( 1) 02920000 C RA (1348) . . . . . LGTN (12) 02930000 C 02940000 C VERTICAL TWO-WAY TRAVELTIME LVRT(12) (R4) 02950000 C RA (1349) . . . . . LVRT ( 1) 02960000 C RA (1360) . . . . . LVRT (12) 02970000 C 02980000 C FREQUENCY SCALAR LBP1(12) (R4) 02990000 C RA (1361) . . . . . LBP1 ( 1) 03000000 C RA (1372) . . . . . LBP1 (12) 03010000 C 03020000 C REFLECTION POINT OFFSET LRFO(12) (I4) 03030000 C RA (1373) . . . . . LRFO ( 1) 03040000 C RA (1384) . . . . . LRFO (12) 03050000 C 03060000 C VERTICAL DEPTH OF RECEIVER LDPT(12) (I4) 03070000 C RA (1385) . . . . . LDPT ( 1) 03080000 C RA (1396) . . . . . LDPT (12) 03090000 C 03100000 C NOT USED 03110000 C RA (1397) 03120000 C RA (1499) 03130000 C 03140000 C INTERVAL VELOCITIES INVEL(500) (R4) 03150000 C RA (1501) . . . . . INVEL ( 1) 03160000 C RA (2000) . . . . . INVEL (500) 03170000 C 03180000 C INTERVAL THICKNESSES INDPTH(500) (R4) 03190000 C RA (2001) . . . . . INDPTH ( 1) 03200000 C RA (2500) . . . . . INDPTH (500) 03210000 C 03220000 C VERTICAL TRAVELTIMES AND DEPTHS 03230000 C RA (2501) . . . . . VERT TT OF SP/DP 1, TRACE 1 03240000 C RA (2502) . . . . . VERT DP OF SP/DP 1, TRACE 1 03250000 C RA (2503) . . . . . VERT TT OF SP/DP 1, TRACE 2 03260000 C RA (2504) . . . . . VERT DP OF SP/DP 1, TRACE 2 03270000 C . 03280000 C . 03290000 C RA ( ) . . . . . VERT TT OF SP/DP LAST, TRACE LAST 03300000 C RA ( ) . . . . . VERT DP OF SP/DP LAST, TRACE LAST 03310000 C 03320000 SUBROUTINE VSPI0(OH,ICC,AUTO3,IABORT,RA) 03330000 C 03340000 RETURN 03350000 END 03360000