Report on iteration 1

started at Thu 11 May 2017 08:14:50 BST in folder /Users/nickc/xval/sample_sprint_scripts/contoured_vint hostname NJCMBP15 process ID: 10319 command line sprint

01_Horizon_1_TVDSS_all_xyz.csv corresponds to 01_Horizon_1_Time.dat.
02_Horizon_2_TVDSS_all_xyz.csv corresponds to 02_Horizon_2_Time.dat.
03_Horizon_3_TVDSS_all_xyz.csv corresponds to 03_Horizon_3_Time.dat.
04_Horizon_4_TVDSS_all_xyz.csv corresponds to 04_Horizon_4_Time.dat.
05_Horizon_5_TVDSS_all_xyz.csv corresponds to 05_Horizon_5_Time.dat.
06_Top_Reservoir_TVDSS_all_xyz.csv corresponds to 06_Top_Reservoir_Time.dat.

Using 5 of 5 horizons
SLD: Layer array: 1 1 1 1 1

Data loading File 1: 01_Horizon_1_TVDSS_all_xyz.csv contains 5 wells
File 2: 02_Horizon_2_TVDSS_all_xyz.csv contains 5 wells
File 3: 03_Horizon_3_TVDSS_all_xyz.csv contains 5 wells
File 4: 04_Horizon_4_TVDSS_all_xyz.csv contains 6 wells
File 5: 05_Horizon_5_TVDSS_all_xyz.csv contains 7 wells

Subset wells not being treated separately.

Minimum isochron for horizon 2 is -27.513. There are 1 negative points.

Horizon 1 "01_Horizon_1_Time.dat"

Time map "01 Horizon 1 Time.dat"

V₀,k to minimise max depth error (top horizon -1 = Datum ):

Best point found is v0 = 1618.60716158007 k = 0.456900948527121 with error 5.8796220695416:
Controlling wells (note +/- sign of error not reported here)
Well_C with 5.8796220695416;
Well_A with 5.87962206521001;
Well_G with 5.87962204691325;
Well_F with 4.0024032553489

Depth (not forced to tie to wells):

Average velocity:

Range from 1730.44 to 1889.93.

Horizon 2 "02_Horizon_2_Time.dat"

Time map "02 Horizon 2 Time.dat"

Isochron from "01 Horizon 1 Time.dat" to "02 Horizon 2 Time.dat"

V₀,k with uplift (p) to minimise max depth error (top horizon 0 = 01 Horizon 1 Time.dat ):

Best point found is v0 = 2205.27179269258 k = 0.0167449143438167 p = 0.998718148159038 with error 8.71517510703347:
Controlling wells (note +/- sign of error not reported here)
Well_A with 8.71517510703347;
Well_C with 8.71517510480908;
Well_G with 8.71517509334899;
Well_F with 2.25033895998195

Depth (not forced to tie to wells):

Isochore

Interval velocity:

Range from 2223.78 to 2249.38.

Average velocity:

Range from 1772.10 to 1914.40.

Horizon 3 "03_Horizon_3_Time.dat"

Time map "03 Horizon 3 Time.dat"

Isochron from "02 Horizon 2 Time.dat" to "03 Horizon 3 Time.dat"

Iteration 1: horizon 2 method 0: Mapped (to tie) V₀ at k = 0.5 using V = Vo + k Z referenced to seafloor:

Well name	Vo	Isochron	Well TVDSS	Top Depth
Well_A	1996.0	258.0	1382.0	1111.5
Well_C	2169.4	178.3	1247.2	1043.7
Well_E	2011.3	89.5	1020.0	926.8
Well_F	2079.0	80.2	1094.0	1007.7
Well_G	2110.8	69.0	1371.0	1294.3

Variogram range is 5000

Variogram model:

$VAR1 = { 'variogram_range' => 5000 };

Depth (not forced to tie to wells):

Isochore

The isochore grid appears to have zero defined values!

Interval velocity:

Range from 2077.65 to 2281.88.

Average velocity:

Range from 1799.44 to 1933.71.

Horizon 4 "04_Horizon_4_Time.dat"

Time map "04 Horizon 4 Time.dat"

Isochron from "03 Horizon 3 Time.dat" to "04 Horizon 4 Time.dat"

Iteration 1: horizon 3 method 0: Mapped (to tie) Vint to tie with isochron drift:

Well name	Vo	Drift	Isochron	Well TVDSS	Top Depth
Well_A	2083.3	0.1	149.4	1537.7	1382.1
Well_B	1418.0	0.1	145.0	1421.0	1318.2
Well_C	2196.4	0.1	161.2	1424.0	1247.0
Well_E	1682.7	0.1	124.8	1125.0	1020.0
Well_F	1695.5	0.0	83.8	1165.0	1094.0
Well_G	2076.9	0.0	62.6	1436.0	1371.0

to tie Vint kriged with external drift.

Rank correlation coefficient is 0.445676274944568

R-squared correlation coefficient is 0.00115250349100461

Vint kriged with external drift.

Well name	Vint	Drift
Well_A	2083.25095867034	0.0747015827862207
Well_B	1418.03476939699	0.0724988948804029
Well_C	2196.43459262282	0.0805827638210448
Well_E	1682.74138061586	0.0623805369420649
Well_F	1695.4618174791	0.0418803354364369
Well_G	2076.89508960349	0.0313008019071277

Variogram range is 5000

Variogram model:

$VAR1 = { 'variogram_range' => 5000 };

Depth (not forced to tie to wells):

Isochore

Interval velocity:

Range from 1236.53 to 2636.28.

Average velocity:

Range from 1791.11 to 2020.37.

Horizon 5 "05_Horizon_5_Time.dat"

Time map "05 Horizon 5 Time.dat"

Isochron from "04 Horizon 4 Time.dat" to "05 Horizon 5 Time.dat"

multilinear apparent isochore regression.

Isochore = 37619.4432470802

+ Isochron (msTWT) * 2.21728734588195

+ Time (msTWT) * -1.39126693289038

+ X coordinate * 0.00797133554233838

+ Y coordinate * -0.0859881164712175

+ Depth of previous horizon * -0.260126957619227

with RMS error 0.217750116327888

Depth (not forced to tie to wells):

Isochore

Interval velocity:

Range from 1305.25 to 4630.48.

Average velocity:

Range from 1722.82 to 3020.89.

Depth errors gridded with a IDQ range:

Variogram range is 800

Final average velocity map:

Final depth map:

Cross validation RMS 42.07 for 8 wells 28 values with radius 500

Depth errors:

Horizon	01 Horizon 1 Time.dat	02 Horizon 2 Time.dat	03 Horizon 3 Time.dat	04 Horizon 4 Time.dat	05 Horizon 5 Time.dat
Mean depth error	-0.11	2.44	0.02	-0.04	0.00
Stdev of depth error	5.01	6.40	0.10	0.51	0.22
Number of points	5	5	5	6	7
Minimum	-5.88	-8.72	-0.10	-1.02	-0.33
Maximum	5.88	8.72	0.19	0.66	0.31
Horizon	01 Horizon 1 Time.dat	02 Horizon 2 Time.dat	03 Horizon 3 Time.dat	04 Horizon 4 Time.dat	05 Horizon 5 Time.dat
Well_A	-5.88	-8.72	-0.10	0.22	0.13
Well_B				-1.02
Well_C	5.88	8.72	0.19	0.66	0.18
Well_E	-2.43	1.24	-0.03	-0.19	0.31
Well_F	-4.00	2.25	0.01	0.00	-0.28
Well_G	5.88	8.72	0.01	0.10	0.03
Well_H					-0.03
Well_I					-0.33

Cross-validation residuals:

Horizon	01 Horizon 1 Time.dat	02 Horizon 2 Time.dat	03 Horizon 3 Time.dat	04 Horizon 4 Time.dat	05 Horizon 5 Time.dat		All
RMS	9.65	15.04	10.94	44.15	71.39		42.07
Mean	0.28	3.23	4.29	12.60	16.63		8.25
Stdev	9.65	14.69	10.06	42.31	69.42		41.25
Number of points	5	5	5	6	7		28
Minimum	-16.05	-19.15	-7.15	-53.15	-31.13		-53.15
Maximum	11.13	24.65	19.84	63.49	181.73		181.73
Horizon	01 Horizon 1 Time.dat	02 Horizon 2 Time.dat	03 Horizon 3 Time.dat	04 Horizon 4 Time.dat	05 Horizon 5 Time.dat	n	RMS
Well_A	-16.05	-19.15	-7.15	36.58	-29.39	5	23.99
Well_B	-	-	-	-53.15	-	1	53.15
Well_C	11.13	12.39	11.76	63.49	19.19	5	31.03
Well_E	0.34	-1.73	-3.79	-26.17	-31.13	5	18.28
Well_F	-2.76	-0.00	0.78	1.85	-0.92	5	1.58
Well_G	8.76	24.65	19.84	53.01	181.73	5	85.92
Well_H	-	-	-	-	-20.94	1	20.94
Well_I	-	-	-	-	-2.13	1	2.13

Report on iteration 1

Horizon 1 "01_Horizon_1_Time.dat"

Time map "01 Horizon 1 Time.dat"

V0,k to minimise max depth error (top horizon -1 = Datum ):

Depth (not forced to tie to wells):

Average velocity:

Horizon 2 "02_Horizon_2_Time.dat"

Time map "02 Horizon 2 Time.dat"

Isochron from "01 Horizon 1 Time.dat" to "02 Horizon 2 Time.dat"

V0,k with uplift (p) to minimise max depth error (top horizon 0 = 01 Horizon 1 Time.dat ):

Depth (not forced to tie to wells):

Isochore

Interval velocity:

Average velocity:

Horizon 3 "03_Horizon_3_Time.dat"

Time map "03 Horizon 3 Time.dat"

Isochron from "02 Horizon 2 Time.dat" to "03 Horizon 3 Time.dat"

Iteration 1: horizon 2 method 0: Mapped (to tie) V0 at k = 0.5 using V = Vo + k Z referenced to seafloor:

Variogram model:

Depth (not forced to tie to wells):

Isochore

Interval velocity:

Average velocity:

Horizon 4 "04_Horizon_4_Time.dat"

Time map "04 Horizon 4 Time.dat"

Isochron from "03 Horizon 3 Time.dat" to "04 Horizon 4 Time.dat"

Iteration 1: horizon 3 method 0: Mapped (to tie) Vint to tie with isochron drift:

to tie Vint kriged with external drift.

Vint kriged with external drift.

Variogram model:

Depth (not forced to tie to wells):

Isochore

Interval velocity:

Average velocity:

Horizon 5 "05_Horizon_5_Time.dat"

Time map "05 Horizon 5 Time.dat"

Isochron from "04 Horizon 4 Time.dat" to "05 Horizon 5 Time.dat"

multilinear apparent isochore regression.

Depth (not forced to tie to wells):

Isochore

Interval velocity:

Average velocity:

Depth errors gridded with a IDQ range:

Final average velocity map:

Final depth map:

Depth errors:

Cross-validation residuals:

V₀,k to minimise max depth error (top horizon -1 = Datum ):

V₀,k with uplift (p) to minimise max depth error (top horizon 0 = 01 Horizon 1 Time.dat ):

Iteration 1: horizon 2 method 0: Mapped (to tie) V₀ at k = 0.5 using V = Vo + k Z referenced to seafloor: