2b:null 2c:["$","title",null,{"children":"Overly Simplified Stratigraphic Modeling | Justin Gosses | Observable"}] 31:T580,html `

Simple Stratigrahpic Concepts We Might Be Able To Put into an Explorable Explanation, in order

Superposition: oldest things on bottom: show age goes up as depth goes up as measured from top.
- How to show: Plot of depth vs. Age. Assume continuous deposition & no post-deposition tectonic movement.
Walters law 1: neighboring facies on a map are neighboring facies in a well log
- How to show: Demonstrate that facies next to one another on the map are the same neighboring facies in the well.
Walkers law 2: same as above, with demonstrated exception of when there is a lot of erosion.
- How to show: Demonstrate non-neighbors can be next to one another in well if there is erosion that takes out the facies that would have been between. Show erosion where missing facies would have been.
Relative vs. Absolute time are not 1to1: As most time not preserved.
- How to show:Plot preserved time vs. all time per facies part of GDE. Need ability to turn erosion, rate deposition, non-deposition on or off during some facies changes.
Unconformity, disconformity, confirmity, paraconformity
- How to show: add lines to the previous plots that label changes in preservation state.

`32:T526,html `

How to put in code?

A 2D grid that represents a map surface of depositional environments. Each grid block represents a facies (or type of depositional environment) within an overall GDE (gross depositional environment) map. Initially, we'll just use a 1D array within the 2D grid for the GDE and move up and down in dip direction only.
Series of time points in an array. Holds the movement of the GDE map relative to the location of the well at each time point. We'll assume the GDE stays constant and just moves up or down dip.
Small arrays of 10 values, one array for each typical log shape pattern Log shapes include upwards coarsening, finding, serrated sand, block sand, shale, etc. Values go from 0 to 10. There are 10 values in each array.
Object that relates each facies type to the probability each log shape pattern will occur during deposition of that facies.There can be one or more log shapes per facies.
Object that relates each facies type & value of recent movement to the probability of deposition, non-deposition, or erosion.This will get recorded at each time step like movement amount, position, facies, and log shape. <= It is currently not implimented!

`33:Te88,different_length_movements = [{ "shortSortaRandom":[1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,-2,9,7,0,0,0,1,-8,-5,-1,1,-1,10,10,-4,5,-1,1,1,0,0,5,0,0,0,0,10,-2], "mediumSortaRandom":[1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,-2,9,7,0,0,0,1,-8,-5,-1,1,-1,10,10,-4,5,-1,1,1,0,0,5,0,0,0,0,10,-2,-2,0,0,4,5,6,-8,1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,10,1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,-2,9,7,0,0,0,1,-8,-5,-1,1,-1,10,10,-4,5,-1,1,1,0,0,5,0,0,0,0,10,-2,-2,0,0,4,5,6,-8,1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,10,1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,-2,9,7,0,0,0,1,-8,-5,-1,1,-1,10,10,-4,5], "longSortaRandom":[1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,-2,9,7,0,0,0,1,-8,-5,-1,1,-1,10,10,-4,5,-1,1,1,0,0,5,0,0,0,0,10,-2,-2,0,0,4,5,6,-8,1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,10,1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,-2,9,7,0,0,0,1,-8,-5,-1,1,-1,10,10,-4,5,-1,1,1,0,0,5,0,0,0,0,10,-2,-2,0,0,4,5,6,-8,1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,10,1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,-2,9,7,0,0,0,1,-8,-5,-1,1,-1,10,10,-4,5,1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,-2,9,7,0,0,0,1,-8,-5,-1,1,-1,10,10,-4,5,-1,1,1,0,0,5,0,0,0,0,10,-2,-2,0,0,4,5,6,-8,1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,10,1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,-2,9,7,0,0,0,1,-8,-5,-1,1,-1,10,10,-4,5,-1,1,1,0,0,5,0,0,0,0,10,-2,-2,0,0,4,5,6,-8,1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,10,1,5,4,10,10,-2,-7,-6,2,2,15,-6,-3,-2,5,0,-7,-9,-9,-6,-2,9,7,0,0,0,1,-8,-5,-1,1,-1,10,10,-4,5], "rhythmicSimpleShort":[0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1], "rhythmicSimpleMedium":[0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1,-1], "rhythmicshort":[0,1,2,3,4,5,6,7,6,5,4,3,2,1,0,1,-1,-2,-3,-4,-5,-6,-7,-6,-5,-4,-3,-2,-1,0,0,1,2,3,4,5,6,7,6,5,4,3,2,1,0,1,-1,-2,-3,-4,-5,-6,-7,-6,-5,-4,-3,-2,-1,0,0,1,2,3,4,5,6,7,6,5,4,3,2,1,0,1,-1,-2,-3,-4,-5,-6,-7,-6,-5,-4,-3,-2,-1,0,0,1,2,3,4,5,6,7,6,5,4,3,2,1,0,1,-1,-2,-3,-4,-5,-6,-7,-6,-5,-4,-3,-2,-1,0], "rhythmicmedium":[0,1,2,3,4,5,6,7,6,5,4,3,2,1,0,1,-1,-2,-3,-4,-5,-6,-7,-6,-5,-4,-3,-2,-1,0,0,1,2,3,4,5,6,7,6,5,4,3,2,1,0,1,-1,-2,-3,-4,-5,-6,-7,-6,-5,-4,-3,-2,-1,0,0,1,2,3,4,5,6,7,6,5,4,3,2,1,0,1,-1,-2,-3,-4,-5,-6,-7,-6,-5,-4,-3,-2,-1,0,0,1,2,3,4,5,6,7,6,5,4,3,2,1,0,1,-1,-2,-3,-4,-5,-6,-7,-6,-5,-4,-3,-2,-1,0,0,1,2,3,4,5,6,7,6,5,4,3,2,1,0,1,-1,-2,-3,-4,-5,-6,-7,-6,-5,-4,-3,-2,-1,0,0,1,2,3,4,5,6,7,6,5,4,3,2,1,0,1,-1,-2,-3,-4,-5,-6,-7,-6,-5,-4,-3,-2,-1,0,0,1,2,3,4,5,6,7,6,5,4,3,2,1,0,1,-1,-2,-3,-4,-5,-6,-7,-6,-5,-4,-3,-2,-1,0,0,1,2,3,4,5,6,7,6,5,4,3,2,1,0,1,-1,-2,-3,-4,-5,-6,-7,-6,-5,-4,-3,-2,-1,0] }]34:T757,function createLog(probability_dep_or_eros_per_faces,facies_array,fraction_log_shapes_within_depositional_env,names_depositional_env,starting_position_in_2D_grid, x_movement_of_GDE_over_time){ //let position_x let position_x = starting_position_in_2D_grid[1] let last_position_x = starting_position_in_2D_grid[1] let log = [] let positions = [] let facies_at_time = [] //let position_x = starting_position_in_2D_grid[1] for (let i = 0; i < x_movement_of_GDE_over_time.length; i++){ //// Get position given last or starting position, /// apply in turn each movement in x_movement_of_GDE_over_time and get new position let facies_at_position = 0 /// position_x = position_x+x_movement_of_GDE_over_time[i] let max_facies_array_position = facies_array[1].length if (position_x >= max_facies_array_position){ position_x = max_facies_array_position-1 } //// facies_at_position = facies_array[starting_position_in_2D_grid[0]][position_x] //// let log_shape_key = fraction_log_shapes_within_depositional_env[0][String(facies_at_position)] let random01 = Math.random() let logShape let log_shape_key_keysArray = Object.keys(log_shape_key) for (let i = 0; i < log_shape_key_keysArray.length; i++){ if(random01 < log_shape_key[log_shape_key_keysArray[i]]){ logShape = log_shape_key_keysArray[i] } } let log_shape_array = log_shapes[0][logShape] log.push(log_shape_array); positions.push(position_x) facies_at_time.push(names_depositional_env[0][facies_at_position]) } let final_log = log.flat().reverse() let final_positions = positions.reverse() let final_facies_at_time = facies_at_time.reverse() let final_x_movement_of_GDE_over_time =x_movement_of_GDE_over_time.reverse() return [final_log,positions,facies_at_time,x_movement_of_GDE_over_time] }35:T743,function addDepthToLogArray (log,positions,facies_at_time,x_movement_of_GDE_over_time){ //// x_movement_of_GDE_over_time = x_movement_of_GDE_over_time.reverse() positions = positions.reverse() //// let new_log = [] let facies_at_time_counter = 0 let facies_position = 0 let movement_position = 0 let movement = 0 let last_facies_position = facies_at_time[0] let last_movement_position = positions[0] let last_movement = x_movement_of_GDE_over_time[0] ///// reflection coefficient calculation from bottom of notebook let log_rfe = find_reflectionCoEfficient_forEachItemInArray(log) ///// for (let m = 0; m < log.length; m++){ let new_log_obj = {"depth": "0", "curveGR": "0", "curveName":"fakeSand","Facies":0,"Location":0,"Movement":0, "UWI": "fake"} new_log_obj["depth"] = m new_log_obj["curveGR"] = log[m] new_log_obj["reflectionCoEf"] = log_rfe[m] if(facies_at_time_counter == 10){ //// facies_position=facies_position+1 movement_position=movement_position+1 movement=movement+1 //// new_log_obj["Facies"]= facies_at_time[m/10] new_log_obj["Location"]= positions[m/10] new_log_obj["Movement"]= x_movement_of_GDE_over_time[m/10] ///// last_facies_position = facies_position last_movement_position = movement_position last_movement = movement /// facies_at_time_counter = 0 } else if(m<=10){ new_log_obj["Facies"]= facies_at_time[0] new_log_obj["Location"]= positions[0] new_log_obj["Movement"]= 0 } else{ new_log_obj["Facies"]= facies_at_time[last_facies_position] new_log_obj["Location"]= positions[last_movement_position] new_log_obj["Movement"]= 0 } new_log.push(new_log_obj) facies_at_time_counter += 1 } return new_log }36:T487,function find_reflectionCoEfficient_forEachItemInArray(gr_array){ let sand = [8,9,10] let silt = [4,5,6,7] let shale = [0,1,2,3] let gr_to_facies = {1:"shale",2:"shale",3:"shale",4:"shale",5:"siltySand",6:"siltySand",7:"sand",8:"sand",9:"sand",10:"sand"} // let {"shale":"","siltySand":"p_siltySan","sand":""} let reflectionCoEfficient_array = [] let p_array = [] let v_array = [] for (let m = 0; m < gr_array.length; m++){ if(gr_array[m]<=shale[shale.length-1]){ p_array.push(p_shale_sat); v_array.push(V_shale_sat)} else if(silt[0]< gr_array[m]<=silt[silt.length-1]){ p_array.push(p_silt_sat); v_array.push(V_silt_sat); } else if(sand[0]< gr_array[m]<=sand[sand.length-1]){ p_array.push(p_sand_sat); v_array.push(V_sandstone_sat); } else{p_array.push(0); v_array.push(0);} } let rce for (let n = 0; n < gr_array.length; n++){ if(n==0){ reflectionCoEfficient_array.push(0) } else{ rce = reflectionCoEfficient(p_array[n-1],p_array[n],v_array[n-1],v_array[n]) reflectionCoEfficient_array.push(rce) } } return reflectionCoEfficient_array }37:T1996,// function CurveBox(well_curve_config_template,curve_stroke_d3){ // //// Getting the variables out of the input json template // well_curve_config_template = well_curve_config_template[0] // let multipleLines = well_curve_config_template["multipleLines"] // let curveNames = well_curve_config_template["curveNames"] // let curveColors = well_curve_config_template["curveColors"] // let curveName = curveNames[0] // let curveColor = curveColors[0] // //// // let data = well_curve_config_template["data"] // let width = well_curve_config_template["width"] // let height = well_curve_config_template["height"] // let margin = well_curve_config_template["margin"] // let depth_curve_name = well_curve_config_template["depth_curve_name"] // //// Calculate depth min and max // let depth_min // if(!depth_min){depth_min = d3.min(data, function(d) { return +d[depth_curve_name];});} // let depth_max // if(!depth_max){depth_max = d3.max(data, function(d) { return +d[depth_curve_name];});} // // Calculate x domain extent for one or more than one curve // let mins = [] // let maxes = [] // for (let i = 0; i < curveNames.length; i++) { // let min_this = d3.min(data, function(d) { return +d[curveNames[i]]}) // let max_this = d3.max(data, function(d) { return +d[curveNames[i]]}) // mins.push(min_this) // maxes.push(max_this) // } // let min_all_curves = d3.min(mins) // let max_all_curves = d3.max(maxes) // //// Calculate Axis & Scales // let x = d3.scaleLinear().domain([min_all_curves,max_all_curves]).nice().range([margin.left, width - margin.right]) // let y = d3.scaleLinear().domain([depth_max, depth_min]).nice().range([height - margin.bottom, margin.top]) // let xAxis = g => g.attr("transform", `translate(0,${height - margin.bottom})`).call(d3.axisBottom(x).ticks(width / 80).tickSizeOuter(0)) // let yAxis = g => g.attr("transform", `translate(${margin.left},0)`).call(d3.axisLeft(y)).call(g => g.select(".domain").remove()) // /////// All variables created above now they are applied below // /////// , some in for loops for multiple curves and areas // /////// attempt at function for color used in ..... // function colors(data_array){ // let color_array = [] // for (let i =0;i 0){ // color_array.push =("blue") // } // else { // color_array.push("green") // } // } // return color_array // } // //// START to make the firt SVG // const svg = d3.create("svg") // svg.attr("width",width) // .attr("height",height); // svg.append("g") // .call(xAxis); // svg.append("g") // .call(yAxis); // /////// throw away code for single curve to plot that will be deleted soon // ///// was here: // //// Code that assumes multiple curves are plotted in same curvebox // if(multipleLines == "yes"){ // for (let k = 0; k < curveNames.length; k++) { // ///// code that creates a line for each Curve in order provided and applies // ///// the color in the color array in order provided // let another_line = d3.line().x(d => x(d[curveNames[k]])).y(d => y(d[depth_curve_name])); // svg.append("path") // .datum(data) // .attr("fill", "none") // .attr("stroke", curveColors[k]) // .attr("stroke-width", curve_stroke_d3) // .attr("stroke-linejoin", "round") // .attr("stroke-linecap", "round") // .attr("d", another_line); // } // // define the area filled under the curve // let two_curve_fill_flag = "no" // for (let i = 0; i < well_curve_config_template["fill"].length; i++) { // //// // if (well_curve_config_template["fill"][i]["fill"] == "yes"){ // let number_colors = well_curve_config_template["fill"][i]["fillColors"].length // let curveName1 = well_curve_config_template["fill"][i]["curveName"] // let threshold = -99999999 // let fillColor = "gray" // console.log("got inside fill,",curveName1,'number_colors',number_colors) // ////// // for (let j = 0; j < number_colors; j++) { // console.log("got to start of J loop",j) // let area1 = d3.area() // if (number_colors != 0){ // threshold = well_curve_config_template["fill"][i]["cutoffs"][j] // fillColor = well_curve_config_template["fill"][i]["fillColors"][j] // } // if(well_curve_config_template["fill"][i]["fillDirection"] == "left"){ // let startFromLeft = well_curve_config_template["margin"]["left"] // area1 // .x1(d => x(d[curveName1])) // .x0(d => startFromLeft) // .defined(d => ((d[curveName1])>threshold)) // .y(d => y(d[depth_curve_name])); // } // if(well_curve_config_template["fill"][i]["fillDirection"] == "right"){ // let startFromRight = well_curve_config_template["margin"]["right"] // area1 // .x0(d => x(d[curveName1])) // .x1(d => startFromRight) // .defined(d => ((d[curveName1]) y(d[depth_curve_name])); // } // if(well_curve_config_template["fill"][i]["fillDirection"] == "between"){ // let between2Curve = well_curve_config_template["fill"][i]["curve2"] // area1 // .x1(d => x(d[curveName1])) // .x0(d => x(d[between2Curve])) // .defined(d => ((d[curveName1])>threshold)) // .y(d => y(d[depth_curve_name])); // } // svg.append("path") // .datum(data) // .attr("class", "area") // .attr("d", area1) // .attr("stroke", fillColor) // .attr("fill",fillColor) // .attr("fill-opacity",0.8); // ///// // console.log("got to end of J loop",j) // } // } // } // } // return svg.node(); // }38:T91b,// from https://observablehq.com/@jashkenas/inputs#sliderDemo function input(config) { let { form, type = "text", attributes = {}, action, getValue, title, description, format, display, submit, options } = config; const wrapper = html`

`; if (!form) form = html``; Object.keys(attributes).forEach(key => { const val = attributes[key]; if (val != null) form.input.setAttribute(key, val); }); if (submit) form.append( html`` ); form.append( html`` ); if (title) form.prepend( html`

${title}

` ); if (description) form.append( html`

${description}

` ); if (format) format = typeof format === "function" ? format : d3format.format(format); if (action) { action(form); } else { const verb = submit ? "onsubmit" : type == "button" ? "onclick" : type == "checkbox" || type == "radio" ? "onchange" : "oninput"; form[verb] = e => { e && e.preventDefault(); const value = getValue ? getValue(form.input) : form.input.value; if (form.output) { const out = display ? display(value) : format ? format(value) : value; if (out instanceof window.Element) { while (form.output.hasChildNodes()) { form.output.removeChild(form.output.lastChild); } form.output.append(out); } else { form.output.value = out; } } form.value = value; if (verb !== "oninput") form.dispatchEvent(new CustomEvent("input", { bubbles: true })); }; if (verb !== "oninput") wrapper.oninput = e => e && e.stopPropagation() && e.preventDefault(); if (verb !== "onsubmit") form.onsubmit = e => e && e.preventDefault(); form[verb](); } while (form.childNodes.length) { wrapper.appendChild(form.childNodes[0]); } form.append(wrapper); return form; }39:T4b6,function select(config = {}) { let { value: formValue, title, description, disabled, submit, multiple, size, options } = Array.isArray(config) ? { options: config } : config; options = options.map(o => typeof o === "object" ? o : { value: o, label: o } ); const form = input({ type: "select", title, description, submit, attributes: { disabled }, getValue: input => { const selected = Array.prototype.filter .call(input.options, i => i.selected) .map(i => i.value); return multiple ? selected : selected[0]; }, form: html` ` }); form.output.remove(); return form; }3a:[] 2f:{"id":"4dc9bf67a311fcb7","slug":"overly-simplified-stratigraphic-modeling","title":"Overly Simplified Stratigraphic Modeling","thumbnail":"b445f5e48b0569931e71a87a7924c990725525dce97128421445d7d4dd7ec0a5","access_level":"public","create_time":"$D2019-10-19T19:26:58.316Z","update_time":"$D2022-02-20T01:05:24.010Z","view_time":null,"first_public_time":"$D2019-10-19T20:29:57.143Z","fork":true,"stars":2,"starred":false,"owner":{"type":"team","id":"3168f49a68f07e98","login":"justingosses","name":"Justin Gosses","avatar_url":"https://avatars-next.observableusercontent.com/eba66f9b9de3134c0c22c453a1d562b450a8cf7f2274c4ed1ff16b7455312520","tier":"starter_2024"},"version":1908,"publish_version":1908,"role":null,"authored":false,"authors":[{"type":"individual","id":"44bc5e36abcab180","login":"justingosses","name":"Justin Gosses","avatar_url":"https://avatars-next.observableusercontent.com/eba66f9b9de3134c0c22c453a1d562b450a8cf7f2274c4ed1ff16b7455312520","tier":null}],"members":[{"type":"individual","id":"44bc5e36abcab180","login":"justingosses","name":"Justin Gosses","avatar_url":"https://avatars-next.observableusercontent.com/eba66f9b9de3134c0c22c453a1d562b450a8cf7f2274c4ed1ff16b7455312520","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"subscribed":null,"forks":0,"collection":{"id":"60edc6e9d1104035","slug":"geology","title":"Geology","owner":{"id":"3168f49a68f07e98","login":"justingosses","name":"Justin Gosses"}},"type":"notebook","body":{"title":"Overly Simplified Stratigraphic Modeling","cells":[{"id":0,"value":"md`# Overly Simplified Stratigraphic Modeling`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1865,"value":"md `Simplistic playground for how cyclic changes in sea-level move different depositional environments across a location where deposition is occuring creating the patterns often seen in a geologic well log or cliff face.`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1645,"value":"html `

`","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":808,"value":"html `draft 2, work in progress, vertical curve with D3.js, horiztonal with vega`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1219,"value":"plot1 = CurveBox(well_curve_config_template,curve_stroke_d3)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1305,"value":"viewof short_or_long_time_plotted = select({\n title: \"Select how the 1D GDE (gross depositional environment) moves relative to the well location over time. This also impacts the amount of time modeled.\",\n description: \"\",\n options: [\"shortSortaRandom\",\"mediumSortaRandom\",\"longSortaRandom\",\"rhythmicSimpleShort\",\"rhythmicSimpleMedium\",\"rhythmicshort\",\"rhythmicmedium\"],\n value: \"short\"\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1449,"value":"html `Change the length of different facies in the GDE`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1342,"value":"rgbForm = html``","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1882,"value":"html ` The blue vertical line is where the well is. The colored dots represent different environments that will move left or right as sea level goes up or down. `","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":934,"value":"vegalite({\n \"width\": screen.width-screen.width*.450,\n \"height\": 120,\n title:\"new_log2 - Facies Distribution on Gross Depositional Map. Well location relative to map starting at position defined above\",\n data: {values: array_facies_on_gde_original},\n \"layer\": [\n {\n mark: \"point\",\n encoding: {\n x: {field: \"PositionOnGDEarray\", type: \"quantitative\"},\n y: {field: \"Facies\",\"type\": \"ordinal\", \"sort\": [\"coastal_plain\", \"shoreface\", \"shelf\",\"deepwater\"]},\n color: {field: \"Facies\", \"type\": \"nominal\"},\n }, \n },\n { \n \"mark\": \"rule\",\n \"encoding\": {\n // \"x\":{\"value\":starting_x_position*13.8},\n \"x\": {\n \"field\": \"starting_x_position\",\n \"type\": \"quantitative\",\n \"aggregate\": \"mean\"\n },\n \"size\": {\"value\": 3},\n \"color\": {\"field\": \"symbol\", \"type\": \"nominal\"}\n }\n }\n// },\n]\n \n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1460,"value":"viewof starting_x_position = slider({\n min: 10, \n max: facies_array[1].length-10, \n step: 1, \n value: 24, \n title: \"Define starting position of well relative to GDE before GDE moves over time \", \n description: \"Minimum of 10 km. from furthest up-hill point in depositional enviroments above.\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1630,"value":"html `

Now lets define some log shape probabilities for each facies belt in the GDE

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1539,"value":"viewof prob_serratedSand_vs_UpCSand_in_CoastalPlain = slider({\n min: 0.1, \n max: 1, \n step: 0.1, \n value: 0.3, \n title: \"Probabilty serrated vs. upwards coarsening log shape in COASTAL PLAIN facies\", \n description: \"greater number results in more change of second facies\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1543,"value":"viewof prob_upF_vs_blocky_in_shoreface = slider({\n min: 0.1, \n max: 1, \n step: 0.1, \n value: 0.3, \n title: \"Probabilty upwards fining vs. blocky log shape in SHOREFACE facies\", \n description: \"greater number results in more change of second facies\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1548,"value":"viewof prob_upC_vs_serratedShale_in_shelf = slider({\n min: 0.1, \n max: 1, \n step: 0.1, \n value: 0.3, \n title: \"Probabilty upwards coarsening vs. blocky log shape in SHELF facies\", \n description: \"greater number results in more change of second facies\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1557,"value":"viewof prob_serratedSand_vs_shale_in_dw = slider({\n min: 0.1, \n max: 1, \n step: 0.1, \n value: 0.3, \n title: \"Probabilty upwards fining vs. almost flat shale log shape in DEEPWATER facies\", \n description: \"greater number results in more change of second facies\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":439,"value":"vegalite({\n \"width\": screen.width-screen.width*.450,\n \"height\": 120,\n data: {values: new_log2},\n title:\"new_log2 - How GDE map moved over time in dip direction relative to fixed well locatin. \",\n mark: \"line\",\n encoding: {\n y: {field: \"Location\", type: \"quantitative\", axis: {}},\n x: {field: \"depth\", type: \"ordinal\"\n },\"color\": {\"value\": \"brown\"}\n },\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":483,"value":"vegalite({\n \"width\": screen.width-screen.width*.450,\n \"height\": 120,\n data: {values: new_log2},\n title:\"new_log2 - The magnitude of movement within each time point of the GDE map relative to fixed core location\",\n mark: \"line\",\n encoding: {\n y: {field: \"Movement\", type: \"quantitative\", axis: {}},\n x: {field: \"depth\", type: \"ordinal\"\n },\"color\": {\"value\": \"brown\"}\n },\n})\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":591,"value":"vegalite({\n \"width\": screen.width-screen.width*.450,\n \"height\": 120,\n data: {values: new_log2},\n title:\"new_log2 - Which GDE Elements Deposited Which Part Of Logs\",\n mark: \"point\",\n encoding: {\n y: {field: \"Facies\", type: \"nominal\",\"sort\": [\"coastal_plain\", \"shoreface\", \"shelf\",\"deepwater\"]},\n x: {field: \"depth\", type: \"ordinal\"},\n color: {field: \"Facies\", \"type\": \"nominal\"},\n },\n },\n)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":173,"value":"vegalite({\n \"width\": screen.width-screen.width*.450,\n \"height\": 120,\n data: {values: new_log2},\n title:\"new_log2 - pseudo GR (gamma ray)\",\n mark: \"line\",\n encoding: {\n y: {field: \"curveGR\", type: \"quantitative\", axis: {\n \t}\n },\n x: {field: \"depth\", type: \"ordinal\", axis: {\n \tmaxExtent:110,\n \tminExtent:0,\n \t}\n }\n },\n})\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1174,"value":"// vegalite({\n// \"width\": screen.width-screen.width*.450,\n// \"height\": 120,\n// data: {values: new_log2},\n// title:\"new_log2 - Reflection Coefficient on point by point basis...needs more work\",\n// mark: \"line\",\n// encoding: {\n// y: {field: \"reflectionCoEf\", type: \"quantitative\", axis: {}},\n// x: {field: \"depth\", type: \"ordinal\"\n// },\"color\": {\"value\": \"orange\"}\n// }\n// })\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":775,"value":"// vegalite({\n// \"width\": screen.width-screen.width*.450,\n// \"height\": 120,\n// data: {values: array_facies_on_gde_original},\n// title:\"new_log2 - Original Facies Distribution on 1D Gross Depositional Map\",\n// mark: \"point\",\n// encoding: {\n// x: {field: \"PositionOnGDEarray\", type: \"quantitative\"},\n// y: {field: \"Facies\",\"type\": \"ordinal\", \"sort\": [\"coastal_plain\", \"shoreface\", \"shelf\",\"deepwater\"]},\n// color: {field: \"Facies\", \"type\": \"nominal\"},\n// },\n// },\n// )","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":692,"value":"html `

Premise:

Stratigraphy 101 is filled with basic laws and concepts, why not put them into a very simple explorable, changable, fixable model? Although most geologists learn Python as its great for working with data, JavaScript in Observable notebooks are great way to build interactive explanations!

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":313,"value":"$31","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":32,"value":"$32","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":728,"value":"html `

The basic idea...

The basic idea is that, there is a distribution of facies within a GDE map, over time the position of the map relative to the well changes and different facies can be deposited. Each facies can deposited one or more log shape pattern. Each time deposition occurs, a log shape pattern array is added to the well.

If there are 48 time steps, and each time step only deposits and deposits one log shape pattern array of 10 values, then the final well will be 480 units long.

\n`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":15,"value":"html `

Future Things to Add in

Erosion!
Erosion based on rate of GDE movement & over which facies.
More than just two log curve shapes per facies.
Interaction of GDE movement and a topographical surface that affects deposition or erosion given accommodation space?
Get in some real logs and try to create fake versions by twiddling with the knobs of our log creation machine?

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":46,"value":"html `

Potentially additional Stratigraphic Concepts:

Group, formation, member, bed
Cyclostratigraphy
Sequence stratigraphy / sea-level rise & fall
Unconformity, Angular unconformity,Disconformity,Paraconformity,Nonconformity
Interpretation of depositional environments from rock record.
Asymmetry of Transgressive and\nRegressive Cycles
Lithologic vs. chronostratigraphic correlation surfaces

Most concepts in this slide presentation saved as PDF

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":308,"value":"html `

https://people.ucsc.edu/~jzachos/migrated/eart120/lectures_08/LECT12.pdf

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":892,"value":"html `

INPUTS

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":989,"value":"html `

Changable Defininitions

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":206,"value":"html `

The Gross Depositional Environment Map is defined as an array of array where each facies element at that part of the GDE has a number. You can check what number = what GDE element/facies below:

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":885,"value":"names_depositional_env = [{\"0\":\"deepwater\",\"1\":\"shelf\",\"2\":\"shoreface\",\"3\":\"coastal_plain\"}]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1698,"value":"html `The variable facies_array_original is the GDE that we start with. It is used with the sliders above and results in the sliders for facies lengths to create the facies_array GDE variable, which is then used by the rest of the code`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":59,"value":"facies_array_original = [[3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],[3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0],[3,3,3,3,3,3,3,3,3,3,3,3,3,3,2,2,2,2,2,2,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0]]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1343,"value":"coastal_plain_length = Generators.input(rgbForm.a)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1345,"value":"shoreface_length = Generators.input(rgbForm.b)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1347,"value":"shelf_length = Generators.input(rgbForm.c)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1349,"value":"deepwater_length = Generators.input(rgbForm.d)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1375,"value":"facies_lenth_obj = [{\"coastal_plain_length\":coastal_plain_length,\"shoreface_length\":shoreface_length,\"shelf_length\":shelf_length,\"deepwater_length\":deepwater_length}]","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1400,"value":"slider_controlled_facies_single_row = sliderChangeGDEarray(facies_lenth_obj,names_depositional_env)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1415,"value":"html `Now let's replace the second row of our facies array with the new one controlled by the slider`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1393,"value":"function swap(json){\n var ret = {};\n for(var key in json){\n ret[json[key]] = key;\n }\n return ret;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1422,"value":"function merge_facies_array_with_slider_controlled_facies(slider_controlled_facies_single_row,facies_array_original){\n facies_array_original[1] = slider_controlled_facies_single_row\n return facies_array_original\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1368,"value":"///////// WORK IN PROGRESS\nfunction sliderChangeGDEarray(facies_lenth_obj,names_depositional_env){\n let facies_array_strip = []\n let names_depositional_env_swapped = swap(names_depositional_env[0])\n let facies_keys = Object.keys(facies_lenth_obj[0])\n let number_of_facies = facies_keys.length\n for (let i = 0; i < number_of_facies; i++){\n for (let j = 0; j < facies_lenth_obj[0][facies_keys[i]]; j++){\n let facies_trimmedOf__length = facies_keys[i].replace(\"_length\",\"\")\n facies_array_strip.push(parseInt(names_depositional_env_swapped[facies_trimmedOf__length]))\n }\n }\n return facies_array_strip\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1430,"value":"facies_array = merge_facies_array_with_slider_controlled_facies(slider_controlled_facies_single_row,facies_array_original)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1703,"value":"html `\"starting_position_in_2D_grid\" is the starting position of the well in the GDE`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":92,"value":"/// The position of the well in the GDE array of arrays. Over time steps, the map moves relative to the well.\nstarting_position_in_2D_grid = [1,starting_x_position]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1718,"value":"html `This defines the probability of each log shape for a given facies area of the GDE. The probabilities are in part defined by the slider up top.`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":345,"value":"fraction_log_shapes_within_depositional_env = [{\"0\":{\"shale\":1,\"serrated_shale\":prob_serratedSand_vs_shale_in_dw},\"1\":{\"up_coarsen\":1,\"serrated_shale\":prob_upC_vs_serratedShale_in_shelf},\"2\":{\"blocky\":1,\"up_fine\":prob_upF_vs_blocky_in_shoreface},\"3\":{\"serrated_sand\":1,\"up_fine\":prob_serratedSand_vs_UpCSand_in_CoastalPlain}}]","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":981,"value":"html `

Values in second item must be less than first! For each facies number deposited at a given time step, a random number between 0 and 1 is drawn. If number is less than last log shape for that facies number, that log shape is used for that time step. If not, it progresses to the second to last log shape, etc. until a log shape is found for that time step of deposition based on the facies element in the GDE map over the well at that time step.

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":209,"value":"html `

Define How Map Moves Over Time Compared to Well Position

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1727,"value":"html `The \"different_length_movements\" is a JSON that defines a few choices for how the GDE moves in dip direction over time. Basically, we're taking a simple approach where the GDE relative to the well position can only move up or down dip. You can pick between these choices up above.`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1309,"value":"$33","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1732,"value":"html `The \"switchOutDifferentLengthsOfTimeAndMovement\" function takes in the choice for movement array picked via the dropdown up top and places it into the \"x_movement_of_GDE_over_time\" variable which is then used further.`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1314,"value":"function switchOutDifferentLengthsOfTimeAndMovement(short_or_long_time_plotted){\n let keys = Object.keys(different_length_movements[0])\n for (let i = 0; iNOTE: Log shapes currently most be an array of length 10 as the code assumes that when calculating time. This should be changed in future!`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":228,"value":"log_shapes = [{\"up_coarsen\":[2,2,3,3,4,5,6,7,8,8],\"up_fine\":[7,8,7,6,5,4,3,2,1,1],\"blocky\":[5,7,9,10,10,10,10,10,10,10],\"serrated_sand\":[1,6,2,8,4,9,9,1,2,4],\"serrated_shale\":[1,4,1,3,1,1,3,1,5,1],\"shale\":[0,1,1,1,1,1,1,1,1,0]}]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":890,"value":"html `

Fixed Defininitions

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":62,"value":"probability_dep_or_eros_per_faces = [{\"0\":{\"erosion\":0.0,\"deposition\":0.9,\"nondeposition\":0.1},\"1\":{\"erosion\":0.1,\"deposition\":0.8,\"nondeposition\":0.1},\"2\":{\"erosion\":0.3,\"deposition\":0.5,\"nondeposition\":0.2},\"3\":{\"erosion\":0.5,\"deposition\":0.2,\"nondeposition\":0.3}}]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":211,"value":"html `

Functions that take in the inputs above & make intermediate data objects before visualizing

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":342,"value":"function randomHit(modifier) {\n return Math.random() < modifier;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":83,"value":"$34","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":105,"value":"log1_and_position = createLog(probability_dep_or_eros_per_faces,facies_array,fraction_log_shapes_within_depositional_env,names_depositional_env,starting_position_in_2D_grid, x_movement_of_GDE_over_time)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":278,"value":"log1 = log1_and_position[0]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":279,"value":"positions = log1_and_position[1]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":399,"value":"facies_at_time = log1_and_position[2]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":177,"value":"$35","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":182,"value":"new_log2 = addDepthToLogArray(log1,positions,facies_at_time,x_movement_of_GDE_over_time)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1055,"value":"new_log2[1]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1064,"value":"new_log2[476]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":813,"value":"function replaceNumbersWithFacies(facies_array,names_depositional_env){\n let new_facies_array_with_words = []\n for (let i = 0; i < facies_array.length; i++){\n new_facies_array_with_words[i]= []\n for (let j = 0; j < facies_array[i].length; j++){\n new_facies_array_with_words[i].push(names_depositional_env[0][facies_array[i][j]])\n }\n }\n return new_facies_array_with_words\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":778,"value":"function turnFaciesArraysToObj(facies_array,GDE_y_value){\n let array_facies_on_gde_original = []\n \n for (let m = 0; m < facies_array[GDE_y_value].length; m++){\n let template = {\"PositionOnGDEarray\":0,\"Facies\":0}\n let facies = facies_array[GDE_y_value][m]\n template[\"PositionOnGDEarray\"]=m\n template[\"Facies\"]=names_depositional_env[0][facies]\n template[\"starting_x_position\"] = starting_x_position\n array_facies_on_gde_original.push(template)\n }\n\n return array_facies_on_gde_original\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":785,"value":"array_facies_on_gde_original = turnFaciesArraysToObj(facies_array,1)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":823,"value":"facies_array_1 = replaceNumbersWithFacies(facies_array,names_depositional_env)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1070,"value":"html `

Reflection Coefficient Calculations

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1072,"value":"viewof V_sandstone_sat = slider({\n min: 2000,\n max: 3500,\n step: 100,\n value: 2200,\n format: \",\",\n description:\n \"P wave velocity(m/s) - 2000 to 3500, in steps of one hundred,\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1096,"value":"viewof V_silt_sat = slider({\n min: 1500,\n max: 2200,\n step: 100,\n value: 1800,\n format: \",\",\n description:\n \"P wave velocity(m/s) - 1100 to 2200, in steps of one hundred,\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1089,"value":"viewof V_shale_sat = slider({\n min: 1100,\n max: 2500,\n step: 100,\n value: 1500,\n format: \",\",\n description:\n \"P wave velocity(m/s) - 1100 to 2500, in steps of one hundred,\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1106,"value":"viewof p_shale_sat = slider({\n min: 2.0,\n max: 3.4,\n step: 0.1,\n value: 2.5,\n format: \",\",\n description:\n \"p density (g/cm^3) - 2.0 to 3.4, in steps of 0.1,\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1110,"value":"p_sand_sat","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1108,"value":"viewof p_silt_sat = slider({\n min: 2.0,\n max: 3.4,\n step: 0.05,\n value: 2.3,\n format: \",\",\n description:\n \"p density (g/cm^3) - 2.0 to 2.4, in steps of 0.1,\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1115,"value":"viewof p_sand_sat = slider({\n min: 2.0,\n max: 3.4,\n step: 0.05,\n value: 2.2,\n format: \",\",\n description:\n \"p density (g/cm^3) - 2.0 to 3.4, in steps of 0.1,\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1131,"value":"function reflectionCoEfficient(p1,p2,v1,v2){\n let ref_coef = ((p2*v2)-(p1*v1))/((p2*v2)+(p1*v1))\n return ref_coef\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1134,"value":"$36","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1152,"value":"test_gr = [1,1,1,6,6,6,6,7,8,9,1,9,9,1,4,5,6,1,5,5,5]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1202,"value":"md `not sure this looks right?`","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1196,"value":"find_reflectionCoEfficient_forEachItemInArray(test_gr)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1154,"value":"find_reflectionCoEfficient_forEachItemInArray(log1)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1169,"value":"new_log2","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1208,"value":"md `## D3.js Vertical Log Visualization attempt`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1234,"value":"html `Borrowing some things trialed from here: https://observablehq.com/@justingosses/well-log-in-d3-js-v5`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1286,"value":"curve_stroke_d3 =1.6","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1849,"value":"CurveBox = wellioviz.curveBox","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1214,"value":"$37","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1252,"value":"html `This makes the vertical log at the very top of this notebook.`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1238,"value":"md `-----------------------------------`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1073,"value":"html `

Borrowed functions from other notebooks

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1075,"value":"// from https://observablehq.com/@jashkenas/inputs#sliderDemo\nfunction slider(config = {}) {\n let {\n min = 0, max = 1, value = (max + min) / 2, step = \"any\", precision = 2,\n title, description, getValue, format, display, submit,\n } = typeof config === \"number\" ? {value: config} : config;\n precision = Math.pow(10, precision);\n if (!getValue) getValue = input => Math.round(input.valueAsNumber * precision) / precision;\n return input({\n type: \"range\", title, description, submit, format, display,\n attributes: {min, max, step, value},\n getValue\n });\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1078,"value":"$38","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":220,"value":"html `

Imported Libraries

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1081,"value":"d3format = require(\"d3-format@1\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1217,"value":"// d3 = require(\"https://d3js.org/d3.v5.min.js\")\nd3 = require(\"d3@5\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":171,"value":"vegalite = require(\"@observablehq/vega-lite@0.1\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1307,"value":"// import {select} from \"@jashkenas/inputs\"","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1898,"value":"$39","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1638,"value":"import {multipleLogPlot} from \"@justingosses/well-log-in-d3-js-v5\"","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1771,"value":"// wellioviz = require('https://bundle.run/wellioviz@0.0.11')","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1904,"value":"wellioviz = require('https://bundle.run/wellioviz@0.0.32')","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1845,"value":"wellioviz.define_wellioviz()","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1842,"value":"wellioviz.help()","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1661,"value":"html `

IN PROGRESS AREA....

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1779,"value":"example_template = wellioviz.curveBoxTemplateExamples(\"example\")","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1817,"value":"example_template[0][\"components\"][0][\"rectangles\"][0] = []","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1871,"value":"example_template[0][\"curve_box\"][\"height_multiplier_components\"] = 1","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1828,"value":"example_template","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1211,"value":"well_curve_config_template_old = [\n {\"multipleLines\":\"yes\",\"curveNames\":[\"curveGR\"],\n \"curveColors\":[\"black\"],\n \"fill\":[\n {\"curveName\":\"curveGR\",\"fill\":\"yes\",\"fillDirection\":\"left\",\"cutoffs\":[0,4,6],\"fillColors\":[\"gray\",\"orange\",\"yellow\"],\"curve2\":\"\"},\n {\"curveName\":\"RESD\",\"fill\":\"no\",\"fillDirection\":\"left\",\"cutoffs\":[],\"fillColors\":[],\"curve2\":\"\"}\n ],\n \"data\":new_log2,\n \"width\":200,\n \"height\":600,\n \"margin\":({top: 20, right: 3, bottom: 30, left: 30}),\n \"depth_curve_name\":\"depth\"}]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1650,"value":"well_curve_config_template_location_old = [\n {\"multipleLines\":\"yes\",\"curveNames\":[\"Location\"],\n \"curveColors\":[\"black\"],\n \"fill\":[\n {\"curveName\":\"curveGR\",\"fill\":\"yes\",\"fillDirection\":\"left\",\"cutoffs\":[0,4,6],\"fillColors\":[\"gray\",\"orange\",\"yellow\"],\"curve2\":\"\"},\n {\"curveName\":\"RESD\",\"fill\":\"no\",\"fillDirection\":\"left\",\"cutoffs\":[],\"fillColors\":[],\"curve2\":\"\"}\n ],\n \"data\":new_log2,\n \"width\":200,\n \"height\":600,\n \"margin\":({top: 20, right: 3, bottom: 30, left: 30}),\n \"depth_curve_name\":\"depth\"}]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1665,"value":"well_curve_config_template_Movement_old = [\n {\"multipleLines\":\"yes\",\"curveNames\":[\"Movement\"],\n \"curveColors\":[\"black\"],\n \"fill\":[\n {\"curveName\":\"curveGR\",\"fill\":\"no\",\"fillDirection\":\"left\",\"cutoffs\":[0,4,6],\"fillColors\":[\"gray\",\"orange\",\"yellow\"],\"curve2\":\"\"},\n {\"curveName\":\"RESD\",\"fill\":\"no\",\"fillDirection\":\"left\",\"cutoffs\":[],\"fillColors\":[],\"curve2\":\"\"}\n ],\n \"data\":new_log2,\n \"width\":200,\n \"height\":600,\n \"margin\":({top: 20, right: 3, bottom: 30, left: 30}),\n \"depth_curve_name\":\"depth\"}]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1781,"value":"well_curve_config_template = wellioviz.minimumDataIntoTemplateFunc(example_template,new_log2,[\"made up well\"],[\"curveGR\"],[\"black\"],[\"\"],[[\n {\"curve_name\":\"curveGR\",\"fill\":\"yes\",\"fill_direction\":\"left\",\"cutoffs\":[0,4,6],\"fill_colors\":[\"gray\",\"orange\",\"yellow\"],\"curve2\":\"\"}]],\"well_holder\",200,600,\"depth\")","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1790,"value":"well_curve_config_template_location = wellioviz.minimumDataIntoTemplateFunc(example_template,new_log2,[\"made up well\"],[\"Location\"],[\"black\"],[\"\"],[[\n {\"curve_name\":\"Location\",\"fill\":\"no\",\"fill_direction\":\"left\",\"cutoffs\":[0,4,6],\"fill_colors\":[\"gray\",\"orange\",\"yellow\"],\"curve2\":\"\"}]],\"well_holder\",200,600,\"depth\")","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1797,"value":"well_curve_config_template_Movement = wellioviz.minimumDataIntoTemplateFunc(example_template,new_log2,[\"made up well\"],[\"Movement\"],[\"black\"],[\"\"],[[\n {\"curve_name\":\"Location\",\"fill\":\"no\",\"fill_direction\":\"left\",\"cutoffs\":[],\"fill_colors\":[],\"curve2\":\"\"}]],\"well_holder\",200,600,\"depth\")","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1775,"value":"result = wellioviz.multipleLogPlot(\"well_holder\",[well_curve_config_template,well_curve_config_template_location,well_curve_config_template_Movement])","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1643,"value":"// result_old = multipleLogPlot(\"well_holder\",[well_curve_config_template,well_curve_config_template_location,well_curve_config_template_Movement])","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1244,"value":"md `-------------------------------`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":1676,"value":"html `

THINGS TO DO

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":214,"value":"html `

The final log created by moving the depositional environment and its facies over time relative to the fixed position of the well.

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":216,"value":"html `

Currently not taking into account deposition/erosion/etc.

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"},{"id":868,"value":"html `

Currently not plotting preserved time

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined"}],"files":"$W3a","stdlib":"1","titling":"auto"}} 30:"$2f:body"