31:T570,{ var data = unika_hourly_mw.map((month, month_no)=>{ return { y:month, x:[...Array(24).keys()], mode: 'lines', line: {shape:'spline',color: d3.quantize(d3.interpolateHcl("#d66000", "#a9a9b4"), 12)[month_no],width:2}, fill:'', type: 'scatter', hoverinfo: 'name+y', showlegend:true, name: monthNames[month_no] } }) var layout = { height: 270, width: width, font: font, hovermode:'closest', showlegend: true, legend: {orientation: 'h',xanchor: 'left', x:0, y: 1.15}, margin: {l: 90,r: 90,b: 30,t: 40}, xaxis: { showgrid: true, zeroline: false, ticks:'outside', dtick:1, }, yaxis: { title: 'MW', range:[0,210], showgrid: true, zeroline: false, tickformat: '.0f', ticks:'outside' } }; const div=DOM.element('div') plt.newPlot(div, data, layout,{displayModeBar: false}); div.on('plotly_hover', d => { var curveNumber = d.points[0].curveNumber var vals = div.data.map((_, i) => i === curveNumber ? 1 : 0.2) plt.restyle(div, 'opacity', vals) }) div.on('plotly_unhover', d => { var curveNumber = d.points[0].curveNumber var vals = div.data.map((_, i) => i === curveNumber ? 1 : 1) plt.restyle(div, 'opacity', vals) }) return html`

${div}

Unika wind P50 daily power profiles by calendar month.

` }32:T7b5,{ var data = [{ y:unika_diurnal_summary.map(v => v.mean_off), x:[...Array(12).keys()].map(v=>v+1), mode: 'lines', line: {shape:'hvh',color: 'black',width:0.5, dash:'dot'}, fill:'', type: 'scatter', hoverinfo: 'name+y', showlegend:false, name:'' },{ y:unika_diurnal_summary.map(v => v.mean_peak), x:[...Array(12).keys()].map(v=>v+1), mode: 'lines', line: {shape:'hvh',color: 'black',width:0.5, dash:'dot'}, fill:'', type: 'scatter', hoverinfo: 'name+y', showlegend:false, name:'' },{ y:unika_diurnal_summary.map(v => v.mean), x:[...Array(12).keys()].map(v=>v+1), mode: 'lines', line: {shape:'hvh',color: 'black',width:0.5}, fill:'', type: 'scatter', hoverinfo: 'name+y', showlegend:false, name:'' }, { y:unika_diurnal_summary.map(v => v.min), x: [...Array(12).keys()].map(v=>v+1), marker: {color: '#2E8B57', opacity:0.0 }, type: 'bar', showlegend:false, hoverinfo: 'y', }, { y:unika_diurnal_summary.map(v => v.max-v.min), x: [...Array(12).keys()].map(v=>v+1), marker: {color: '#2E8B57', opacity:0.7 }, type: 'bar', showlegend:false, hoverinfo: 'y', }] var layout = { height: 300, width: width*0.7, font: font, hovermode:'closest', barmode: 'stack', showlegend: true, legend: {orientation: 'h',xanchor: 'left', x:0, y: 1.1}, margin: {l: 90,r: 90,b: 30,t: 40}, xaxis: { tickvals: [ 1,2,3,4,5,6,7,8,9,10,11,12 ], ticktext:['Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec'], showgrid: true, zeroline: false, ticks:'outside', dtick:1, }, yaxis: { title: 'MW', showgrid: true, zeroline: false, tickformat: '.0f', ticks:'outside' } }; const div=DOM.element('div') plt.newPlot(div, data, layout,{displayModeBar: false}); return html`

${div}

Unika maximum, minimum and mean average diurnal power by calendar month.

` }33:T45d,{ var data = [ { y:unika_daily_diurnal_requirement_mwh, x: [...Array(12).keys()].map(v=>v+1), text:unika_daily_diurnal_requirement_mwh.map(v => d3.format('.0f')(v)), mode: 'text', textposition: 'inside', marker: {color: '#2E8B57', opacity:0.7 }, type: 'bar', showlegend:false } ] var layout = { height: 250, width: width*0.7, font: font, hovermode:'closest', showlegend: true, legend: {orientation: 'h',xanchor: 'left', x:0, y: 1.1}, margin: {l: 90,r: 90,b: 30,t: 40}, xaxis: { tickvals: [ 1,2,3,4,5,6,7,8,9,10,11,12 ], ticktext:['Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec'], showgrid: true, zeroline: false, ticks:'outside', dtick:1, }, yaxis: { title: 'MWh', showgrid: true, zeroline: false, tickformat: '.0f', ticks:'outside' } }; const div=DOM.element('div') plt.newPlot(div, data, layout,{displayModeBar: false}); return html`

${div}

Unika daily storage requirement to produce flat daily output based on P50 .

` }34:T9be,{ function trace(name, set, col, dash) { let max = set[d3.maxIndex(set, v => (v.energyTransfer))] console.log(max) return [{ x: [max.exceedance], y: [max.energyTransfer], text: [d3.format('.2r')(max.energyTransfer) + 'MWh'], type: 'scatter', mode: 'markers+text', textposition: 'top-right',showlegend: false, marker: {size:15, color:'grey'}, hoverinfo: 'none', },{ name: name, x: set.map(v => v.exceedance), y: movingAverage(set.map(v => (v.energyTransfer)),5), type: 'scatter', mode: 'lines', hoverinfo: 'x+y', showlegend:true, line: {color: col, shape: '', dash: dash, width:1.5} }] } let mon=6 var data = [] data=data.concat(trace(scenarios[0].name, modelledFdc[0], colr(1, 1), '')) data=data.concat(trace(scenarios[1].name, modelledFdc[1], colr(3, 1), '')) data=data.concat(trace(scenarios[2].name, modelledFdc[2], colr(5, 1), '')) data=data.concat(trace(scenarios[3].name, modelledFdc[3], colr(7, 1), '')) data=data.concat([{ x: [0.27,0,0,0.27], y: [1000,1000,0,0], fill: 'toself', mode:'none', fillcolor: 'rgba(220,120,120,0.3)', showlegend:false, type: 'scatter' }]) data.push({ y:unika_daily_diurnal_requirement_mwh, x: hydrology_calmonths.map(m => flowExceedance(m.Flow_median)), text: [...Array(12).keys()].map(v => monthNames[v]), width: 0.05, mode: 'text', textposition: 'inside', marker: {color: '#2E8B57', opacity:0.7 }, type: 'bar', showlegend:false }) var layout = { width: width*0.9, height:360, font: font, margin: {l: 65,r: 10,b: 45,t: 40}, annotations: [ { x: 0.15, y: 200, text: 'No storage', showarrow: false, }], showlegend: true, legend: {orientation: 'h', x: 0,xanchor: 'left',y: 1.1}, xaxis: { showgrid: false, zeroline: false, ticks:'outside', range: [1,0], tickformat: '%.0f', title: 'Flow exceedance' }, yaxis: { title: 'MWh / day', range: [0, 1300], zeroline: true, showline: false, separatethousands:true, tickformat: '.0f' } }; const div=DOM.element('div') plt.newPlot(div, data, layout,{displayModeBar: false}); return html`

${div}

Total daily energy storage and transfer achievable depending on river flow for four storage scenarios (lines) and median daily storage requirement to smooth diurnal wind output (bars).

` }35:T4b9,{ let pcts = modelledMonths[2].map((m,i) => (m.energyTransfer / unika_daily_diurnal_requirement_mwh[i]) + 0.1) var data = [{ y: pcts, text: pcts.map(m => d3.format('.0%')(m)), x: modelledMonths[2].map((v,i) => monthNames[i]), textfont: {color:'white'}, mode: 'text', textposition: 'auto', marker: {color: pcts.map(v => v>1?'#2E8B57':''), opacity:0.7 }, type: 'bar', showlegend:false },{ y: modelledMonths[2].map((m,i) => 1), x: modelledMonths[2].map((v,i) => monthNames[i]), type: 'scatter', mode: 'lines', line: {color: 'black', width: 1}, showlegend: false }] var layout = { width: width*0.9, height:360, font: font, margin: {l: 65,r: 10,b: 45,t: 40}, showlegend: true, legend: {orientation: 'h', x: 0,xanchor: 'left',y: 1.1}, xaxis: { showgrid: false, zeroline: false, ticks:'outside', title: 'Month' }, yaxis: { title: '', zeroline: true, showline: false, separatethousands:true, tickformat: '%.0f', range: [0,1.1] } }; const div=DOM.element('div') plt.newPlot(div, data, layout,{displayModeBar: false}); return html`

${div}

` }36:Td91,// Prepare storage and energy profiles for a day modelDay = function(scenario, flow) { let model = {} let ngo = fdc.filter(v=>v.Flow<=flow)[0] // Copy over model.flow = ngo.Flow model.flowEWR = ngo.EWRTotal model.flowSpill = ngo.FlowSpill model.levelTailwater = ngo.LevelTailwater model.levelHeadpond = ngo.LevelHeadpond model.nethead = ngo.HeadTurbine1 model.energy = ngo.Energy model.power = ngo.Power model.exceedance = ngo.Exceedance // Total daily generation volume Mm3 model.generationVolume = ngo.FlowCanal * (24 * 60 * 60) / 1e6 // Rough calculation of the maximum flow during peak time - flow that produced maximum output at current head and high efficiency model.maxGenerationFlow = params.maxPlantPower / ngo.HeadTurbine1 / 9.81 / params.effectiveEfficiency * 1000 // Limit maximum flow to canal capacity model.maxGenerationFlow = d3.min([model.maxGenerationFlow, params.canalCapacity]) // Maximum flow that could be used during the peak period from storage m3/s model.maxFlowFromStore = model.maxGenerationFlow - ngo.FlowCanal // Max volume that could be used by the peaks Mm3 if (scenario.ref==='nostore') { model.maxVolFromStore = 0 } else { model.maxVolFromStore = model.maxFlowFromStore * (params.storeHours * 60 * 60) / 1e6 } if (scenario.ref==='nostore') { model.maxVolToStore = 0 } else { // Max volume that could be taken from the night period for storage Mm3 model.maxVolToStore = d3.max([((ngo.FlowAvailableForGeneration - scenario.flowMin) * ((24-params.storeHours) * 60 * 60) / 1e6),0]) } // Volume per night that will be stored (min of available or needed) model.storeVol = d3.min([model.maxVolFromStore, model.maxVolToStore]) // Not enough head for generation so no point storing if (ngo.LowHeadShutoff === 1) { model.storeVol = 0 } //Flows model.extraFlow_Off = model.storeVol / (60 * 60 * params.storeHours / 1e6) model.extraFlow_Peak = model.storeVol/ (60 * 60 * (24-params.storeHours) / 1e6) model.flowCanal_Off = ngo.FlowCanal - model.extraFlow_Off model.flowCanal_Peak = ngo.FlowCanal + model.extraFlow_Peak //Total flow model.flowTotal_Off = model.flowCanal_Off + ngo.FlowSpill + ngo.EWRTotal model.flowTotal_Peak = model.flowCanal_Peak + ngo.FlowSpill + ngo.EWRTotal // Total volume model.volume_Off = model.flowCanal_Off * params.storeHours * 60 * 60 / 1e6 model.volume_Peak = model.flowCanal_Peak * (24-params.storeHours) * 60 * 60 / 1e6 model.volumeTotal = model.volume_Off + model.volume_Peak // Energy if (model.generationVolume>0) { model.energy_Off = ngo.Energy / model.generationVolume * model.volume_Off model.energy_Peak = ngo.Energy / model.generationVolume * model.volume_Peak model.energyTransfer = ngo.Energy / model.generationVolume * model.storeVol } else { model.energy_Off = 0 model.energy_Peak = 0 model.energyTransfer = 0 } model.energyTotal = model.energy_Off + model.energy_Peak // Power model.power_Off = model.energy_Off / params.storeHours model.power_Peak = model.energy_Peak / (24-params.storeHours) //Headpond levels model.headpondLevelChange = model.storeVol / params.storageCurve //Tailwater levels model.tailwaterLevel_Off = tailwaterLevel(model.flowTotal_Off) model.tailwaterLevel_Peak = tailwaterLevel(model.flowTotal_Peak) model.tailwaterChange = model.tailwaterLevel_Peak - model.tailwaterLevel_Off return model }37:[{"id":"c454c304237cb8f85550eb2ceaad3284679841e19331e79cf71323e43ee22d583e252e9d2c8b5c467714a865a81396cac89a93c5617396b12bf180a95e951011","href":"https://static.observableusercontent.com/files/c454c304237cb8f85550eb2ceaad3284679841e19331e79cf71323e43ee22d583e252e9d2c8b5c467714a865a81396cac89a93c5617396b12bf180a95e951011","name":"unika_p50-1.csv","type":"text/csv","lastModified":1635661122772.268,"size":1267}] 2f:{"title":"Diurnal smooting of other renewables","cells":[{"id":0,"value":"# Diurnal smooting of other renewables","pinned":false,"hidden":false,"mode":"md"},{"id":80,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":96,"value":"$32","pinned":false,"hidden":false,"mode":"ojs"},{"id":102,"value":"$33","pinned":false,"hidden":false,"mode":"ojs"},{"id":57,"value":"$34","pinned":false,"hidden":false,"mode":"ojs"},{"id":145,"value":"$35","pinned":false,"hidden":false,"mode":"ojs"},{"id":36,"value":"html`\n

Model

\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":147,"value":"modelledMonths = {\n let ret = []\n ret=ret.concat(scenarios.map(s => {\n return hydrology_calmonths.map((m, i) => {\n let d=modelDay(s, m.Flow_median)\n d.Flow_median = m.Flow_median\n return d\n })\n }))\n return ret\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":50,"value":"modelledFdc = scenarios.map(s => {\n return fdc.map(v => {\n return modelDay(s, v.Flow)\n })\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":29,"value":"$36","pinned":false,"hidden":false,"mode":"ojs"},{"id":47,"value":"tailwaterLevel = function(flow) {\n let logFlow = Math.log10(flow)\n if (logFlow >= 3.8) {\n return (-21.429 * logFlow**2) + (182.03 * logFlow) + 597.15\n }\n else {\n if (logFlow >= 2.3) {\n return (6.5299 * logFlow**2) - (30.155 * logFlow) + 999.52\n } else {\n return (-2.6828 * logFlow**2) + (11.265 * logFlow) + 952.68\n }\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":42,"value":"scenarios = [ \n {\n ref: 'unitrated',\n name: 'Unit Rated Flow',\n flowMin: params.unitRatedFlow\n },\n {\n ref: 'halfrated',\n name: '50% Unit Rated',\n flowMin: params.unitRatedFlow * 0.5\n },\n {\n ref: 'fiftycumec',\n name: formatFlow(50),\n flowMin: 50\n },\n {\n ref: 'zero',\n name: 'Zero Off-peak',\n flowMin: 0\n }\n]","pinned":false,"hidden":false,"mode":"ojs"},{"id":40,"value":"params = { return {\n canalCapacity: 1100,\n unitRatedFlow: 250,\n maxPlantPower: 185,\n storageCurve: 10, // Mm3/m headpond rise - From Feasibility Study Rev D Appendix H - Mott MacDonald\n effectiveEfficiency: 0.9212, // P80 efficiency. Used to determine maximum storage flows\n storeHours: 13\n}}","pinned":false,"hidden":false,"mode":"ojs"},{"id":83,"value":"html`\n

Unika

\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":104,"value":"unika_daily_diurnal_requirement_mwh = unika_diurnal_difference.map((m, i) => {\n return d3.sum(m.filter(v => v>0).map(v => v))\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":106,"value":"unika_diurnal_difference = [...Array(12).keys()].map(month_no => {\n return unika_hourly_mw[month_no].map(v => v - unika_monthly_mw[month_no])\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":108,"value":"unika_monthly_mw=unika_monthly_gwh.map((v,i)=>v*1000/(monthDays(i+1)*24))","pinned":false,"hidden":false,"mode":"ojs"},{"id":110,"value":"unika_monthly_gwh=[...Array(12).keys()].map(month_no=>{\n return d3.sum(unika_p50.map(v=>v[month_no]))\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":98,"value":"unika_diurnal_summary = [...Array(12).keys()].map(month_no => {\n return {\n month_no: month_no +1,\n month: monthNames[month_no],\n min: d3.min(unika_hourly_mw[month_no]),\n max: d3.max(unika_hourly_mw[month_no]),\n range: d3.max(unika_hourly_mw[month_no]) - d3.min(unika_hourly_mw[month_no]),\n mean: d3.mean(unika_hourly_mw[month_no]),\n hours_off: unika_hourly_mw[month_no].filter(v => v < d3.mean(unika_hourly_mw[month_no])).length,\n hours_peak: unika_hourly_mw[month_no].filter(v => v >= d3.mean(unika_hourly_mw[month_no])).length,\n mean_off: d3.mean(unika_hourly_mw[month_no].filter(v => v < d3.mean(unika_hourly_mw[month_no]))),\n mean_peak: d3.mean(unika_hourly_mw[month_no].filter(v => v >= d3.mean(unika_hourly_mw[month_no]))), \n }\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":87,"value":"unika_hourly_mw = [...Array(12).keys()].map(month_no=>{\n return [...Array(24).keys()].map(hour_no=> unika_p50[hour_no][month_no]/monthDays(month_no+1)*1000)\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":92,"value":"unika_p50=d3.csvParse(await FileAttachment(\"unika_p50-1.csv\").text()).map(v => {\n return [...Array(12).keys()].map(k => v[k+1])\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":34,"value":"html`\n

Ngonye

\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":68,"value":"flowExceedance = function(flow) {\n let filt = fdc.filter(v=>v.Flow<=flow)\n if (filt.length > 0) {\n return filt[0].Exceedance\n } else {\n return 0\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":71,"value":"hydrology_calmonths = (await d3.csv(\"https://gis.westernpower.org/data/output_data/ngonye_flow_calmonthly.csv\",d3.autoType)).sort((a, b) => a.Month-b.Month)","pinned":false,"hidden":false,"mode":"ojs"},{"id":31,"value":"fdc=(await d3.csv(\"https://gis.westernpower.org/data/output_data/base_pe_fdc.csv?\" + Math.floor(Math.random() * 1000), d3.autoType))","pinned":false,"hidden":false,"mode":"ojs"},{"id":27,"value":"html`\n

Common

\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":25,"value":"function formatFlow(flow){\n return d3.format(',.0f')(flow) + ' m\\u00b3/s'\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":22,"value":"makeTrans=function(col,opacity) {\n let c=d3.color(col)\n c.opacity=opacity\n return c.formatRgb()\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":90,"value":"function monthDays(i) {\n if(i === 4 || i === 6 || i === 9 || i === 11) {\n\t\treturn 30\n\t} else if(i === 2) {\n\t\treturn 28\n\t} else {\n\t\treturn 31\n\t}\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":21,"value":"monthNames={\n return {\n 0:'Jan',\n 1:'Feb',\n 2:'Mar',\n 3:'Apr',\n 4:'May',\n 5:'Jun',\n 6:'Jul',\n 7:'Aug',\n 8:'Sep',\n 9:'Oct',\n 10:'Nov',\n 11:'Dec',\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":19,"value":"monthDict={\n return {\n 'Jan':{index: 0, group: 1},\n 'Feb':{index: 1, group: 1},\n 'Mar':{index: 2, group: 6},\n 'Apr':{index: 3, group: 6},\n 'May':{index: 4, group: 2},\n 'Jun':{index: 5, group: 2},\n 'Jul':{index: 6, group: 3},\n 'Aug':{index: 7, group: 3},\n 'Sep':{index: 8, group: 4},\n 'Oct':{index: 9, group: 4},\n 'Nov':{index: 10, group: 4},\n 'Dec':{index: 11, group: 1},\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":17,"value":"function colr(i = 1, a = 1) {\n return d3hsv.hsv(26, 0.85, ((6 - i + 1.5) * 15)/100, a).toString();\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":15,"value":"function movingAverage(values, N) {\n let i = 0;\n let sum = 0;\n const means = new Float64Array(values.length).fill(NaN);\n for (let n = Math.min(N - 1, values.length); i < n; ++i) {\n sum += values[i];\n }\n for (let n = values.length; i < n; ++i) {\n sum += values[i];\n means[i - Math.floor(N / 2)] = sum / N;\n sum -= values[i - N + 1];\n }\n return means;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":13,"value":"font = ({\n family: 'Courier New, monospace',\n size: 13,\n color: '#6f6f6f'\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":9,"value":"luxon = require(\"luxon@1/build/amd/luxon.js\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":11,"value":"DateTime = luxon.DateTime","pinned":false,"hidden":false,"mode":"ojs"},{"id":7,"value":"d3hsv = require(\"d3-hsv@0.1\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":5,"value":"plt = require(\"https://cdn.plot.ly/plotly-latest.min.js\")","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W37","stdlib":"1","titling":"auto"}