Sara Godoy del Olmo Calculation of contaminated soil volumes : Geostatistics applied to a hydrocarbons spill Lac Megantic Case

Similar documents
Radiation Physics Laboratory - Complementary Exercise Set MeBiom 2016/2017

Observer Bias and Reliability By Xunchi Pu

Higher order derivatives

MCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17)

Applied Statistics II - Categorical Data Analysis Data analysis using Genstat - Exercise 2 Logistic regression

Molecular Orbitals in Inorganic Chemistry

4.2 Design of Sections for Flexure

Dealing with quantitative data and problem solving life is a story problem! Attacking Quantitative Problems

Recursive Estimation of Dynamic Time-Varying Demand Models

INFLUENCE OF GROUND SUBSIDENCE IN THE DAMAGE TO MEXICO CITY S PRIMARY WATER SYSTEM DUE TO THE 1985 EARTHQUAKE

Two Products Manufacturer s Production Decisions with Carbon Constraint

Search sequence databases 3 10/25/2016

Estimation of apparent fraction defective: A mathematical approach

EXST Regression Techniques Page 1

Atomic and Laser Spectroscopy

Constants and Conversions:

Answer Homework 5 PHA5127 Fall 1999 Jeff Stark

Differentiation of Exponential Functions

Fourier Transforms and the Wave Equation. Key Mathematics: More Fourier transform theory, especially as applied to solving the wave equation.

ph People Grade Level: basic Duration: minutes Setting: classroom or field site

Extraction of Doping Density Distributions from C-V Curves

Collisions between electrons and ions

Hydrogen Atom and One Electron Ions

Linear Non-Gaussian Structural Equation Models

2. Laser physics - basics

MEASURING HEAT FLUX FROM A COMPONENT ON A PCB

Chapter 6: Polarization and Crystal Optics

15. Stress-Strain behavior of soils

That is, we start with a general matrix: And end with a simpler matrix:

Learning Spherical Convolution for Fast Features from 360 Imagery

Title: Vibrational structure of electronic transition

ARIMA Methods of Detecting Outliers in Time Series Periodic Processes

10. Limits involving infinity

First derivative analysis

Unfired pressure vessels- Part 3: Design

orbiting electron turns out to be wrong even though it Unfortunately, the classical visualization of the

COHORT MBA. Exponential function. MATH review (part2) by Lucian Mitroiu. The LOG and EXP functions. Properties: e e. lim.

AS 5850 Finite Element Analysis

Definition1: The ratio of the radiation intensity in a given direction from the antenna to the radiation intensity averaged over all directions.

PHA 5128 Answer CASE STUDY 3 Question #1: Model

Difference -Analytical Method of The One-Dimensional Convection-Diffusion Equation

Sec 2.3 Modeling with First Order Equations

Middle East Technical University Department of Mechanical Engineering ME 413 Introduction to Finite Element Analysis

Determination of Vibrational and Electronic Parameters From an Electronic Spectrum of I 2 and a Birge-Sponer Plot

Principles of Humidity Dalton s law

Simulated Analysis of Tooth Profile Error of Cycloid Steel Ball Planetary Transmission

Estimation of odds ratios in Logistic Regression models under different parameterizations and Design matrices

A=P=E M-A=N Alpha particle Beta Particle. Periodic table

Chapter 13 GMM for Linear Factor Models in Discount Factor form. GMM on the pricing errors gives a crosssectional

PHA 5127 Answers Homework 2 Fall 2001

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

Status of LAr TPC R&D (2) 2014/Dec./23 Neutrino frontier workshop 2014 Ryosuke Sasaki (Iwate U.)

Outline. Image processing includes. Edge detection. Advanced Multimedia Signal Processing #8:Image Processing 2 processing

1 Isoparametric Concept

Interactive spicpms data treatment using Nanocount

Supplementary Materials

Laboratory work # 8 (14) EXPERIMENTAL ESTIMATION OF CRITICAL STRESSES IN STRINGER UNDER COMPRESSION

Data Assimilation 1. Alan O Neill National Centre for Earth Observation UK

Chapter 8: Electron Configurations and Periodicity

4037 ADDITIONAL MATHEMATICS

Chapter 13 Aggregate Supply

Lecture 37 (Schrödinger Equation) Physics Spring 2018 Douglas Fields

Brief Introduction to Statistical Mechanics

SUMMER 17 EXAMINATION

Transitional Probability Model for a Serial Phases in Production

Studies of Turbulence and Transport in Alcator C-Mod Ohmic Plasmas with Phase Contrast Imaging and Comparisons with GYRO*

Cramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter

What are those βs anyway? Understanding Design Matrix & Odds ratios

5.80 Small-Molecule Spectroscopy and Dynamics

Design Guidelines for Quartz Crystal Oscillators. R 1 Motional Resistance L 1 Motional Inductance C 1 Motional Capacitance C 0 Shunt Capacitance

de/dx Effectively all charged particles except electrons

The van der Waals interaction 1 D. E. Soper 2 University of Oregon 20 April 2012

MA 262, Spring 2018, Final exam Version 01 (Green)

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002

Solution: APPM 1360 Final (150 pts) Spring (60 pts total) The following parts are not related, justify your answers:

Unit 6: Solving Exponential Equations and More

Eigenvalue Distributions of Quark Matrix at Finite Isospin Chemical Potential

u x v x dx u x v x v x u x dx d u x v x u x v x dx u x v x dx Integration by Parts Formula

Application of Vague Soft Sets in students evaluation

1 Minimum Cut Problem

CO-ORDINATION OF FAST NUMERICAL RELAYS AND CURRENT TRANSFORMERS OVERDIMENSIONING FACTORS AND INFLUENCING PARAMETERS

Overview. Flood Risk in NSW FLOODING, FLOOD MODELLING AND FLOOD MAPPING COMMUNICATION CHALLENGES. Floodplain Risk Management Process

Dynamic Modelling of Hoisting Steel Wire Rope. Da-zhi CAO, Wen-zheng DU, Bao-zhu MA *

Problem Set #2 Due: Friday April 20, 2018 at 5 PM.


GEOMETRICAL PHENOMENA IN THE PHYSICS OF SUBATOMIC PARTICLES. Eduard N. Klenov* Rostov-on-Don, Russia

CHAPTER 16 HW: CONJUGATED SYSTEMS

Solution of Assignment #2

Brief Notes on the Fermi-Dirac and Bose-Einstein Distributions, Bose-Einstein Condensates and Degenerate Fermi Gases Last Update: 28 th December 2008

Note If the candidate believes that e x = 0 solves to x = 0 or gives an extra solution of x = 0, then withhold the final accuracy mark.

Calculus II (MAC )

The pn junction: 2 Current vs Voltage (IV) characteristics

Alpha and beta decay equation practice

Engineering 323 Beautiful HW #13 Page 1 of 6 Brown Problem 5-12

Instantaneous Cutting Force Model in High-Speed Milling Process with Gyroscopic Effect

Exam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator.

1 1 1 p q p q. 2ln x x. in simplest form. in simplest form in terms of x and h.

3-2-1 ANN Architecture

Estimation of the two-photon QED background in Belle II

Outline. Why speech processing? Speech signal processing. Advanced Multimedia Signal Processing #5:Speech Signal Processing 2 -Processing-

Transcription:

wwwnvisol-canadaca Sara Godoy dl Olmo Calculation of contaminatd soil volums : Gostatistics applid to a hydrocarbons spill Lac Mgantic Cas

Gostatistics: study of a PH contamination CONTEXT OF THE STUDY This study was conductd for pdagogical purposs In Jun 24 a gostatistical training was prformd for th Environmntal Ministry of Qubc and for othr public and privat companis in ordr to bttr undrstand th Gostatistical mthod and its applications to th contaminatd sits Th data st obtaind in th Lac Mgantic first charactrizations was usd for prforming th study cas during this training 2

Gostatistics: study of a PH contamination Ptrolum Hydrocarbons 3

Gostatistics: study of a PH contamination OBJECTIVES OF THE GEOSTATISTICAL STUDY Th gostatistical tratmnt has bn prformd in ordr to: Accuratly undrstand th contaminants bhaviour in soil (corrlations btwn contamination and dpth/topography/gology ) Onc th bhaviour of th contamination is stablishd, w prform a contamination mapping in ordr to visualiz th pollution migration Dlinat, if possibl, th pollution xtnt and th volums of contaminatd soil affctd by th PH spill 4

Gostatistics: study of a PH contamination Rmmbr: A gostatistical study is conductd as follows: Data fild intgration 2 Exploratory Data Analysis 3 Variogram fitting 4 Intrpolation mthods: Krigging / Simulations 5 3D mapping / Estimat of contaminatd soil volums 5

Gostatistics: study of a PH contamination Intgrating th fild data Laboratory data - In total, 45 PH laboratory analyss - Prfrntial sampling prformd in surfac fills - Maximum sampling dpth = 8 m Fild Data: gology - Fills layr (-2m) and natural soil (silt and clay) until 8 m dpth Rmdiation thrshold ppm for th PH Fild data must b go-rfrncd in a coordinat X Y Z systm 6

Gostatistics: study of a PH contamination 2 Exploratory Data Analysis Objctivs: Gological study, Contamination bas map, Quick statistics, histograms Corrlations btwn concntration and dpth/gology, Exprimntal variogram construction

Gostatistics: study of a PH contamination 2 Exploratory Data Analysis Contamination bas map 9 Nb Sampls: Minimum: Maximum: Man: Std Dv: 8 s i c n u q r F 45 9 2556222 864959 6 5 4 3 2 8 2 3 4 5 6 8 Contamin C-C5

Gostatistics: study of a PH contamination 2 Exploratory Data Analysis Corrlation btwn concntration and dpth Th highst concntrations ar masurd btwn 5 and m dpth Concntration moynn n C-C5 (ppm) 6 529,5 5 4 3 9 935,5 53,9 2 9,5 6,8 9 533,33 3,25,38 33,33

Gostatistics: study of a PH contamination 2 Exploratory Data Analysis Sampling fills Sand, gravl, Silt-F Natural soil Silt-NS, clay, organic mattr

Gostatistics: study of a PH contamination 2 Exploratory Data Analysis Corrlation btwn concntration and gology Man Minimum Lithology Cod Frquncy Effctiv 2 Gravl 2,8 53 584,53 428,9 3 Sand 3 5,2 225 2539,8 99,9 4 Silt_F 4,6 34 94,4 35,4 5 Silt_NS 5 28, 26 8,95 254, Clay, 3 6 Org Mattr 6,6 3,4 28, (PH ppm) concntration Maximum concntration

Gostatistics: study of a PH contamination 2 Exploratory Data Analysis Exprimntal Variogram Construction B/ Isotrop Dirctional Variogram NO/D9 5 C C n i m a t n o C : m a r g o i r a V 2 3424 382 36 8 2 466 3929 2349 29235 399 N 6 5 9 4 D-9 3 2 54 3 Contamin C-C5 2 3 Distanc (m) 5 C C n i m a t n o C : m a r g o i r a V 6 5 4 3 2-6 85 6 3 D-9 Contamin C-C5 2 3 Distanc (m) 4 5

Géostatistiqu: étud d un contamination aux HP t au Ni 3 Variogram Fitting B/ Isotrop Dirctional Variogram NO/D9 5 C C n i m a t n o C : m a r g o i r a V 3 3424 382 36 8 2 4663929 2349 29235 399 6 3 5 9 4 3 3 2 9 2 Distanc (m) N D-9 3 5 C C n i m a t n o C : m a r g o i r a V 6 5 4 3 2 2 Distanc (m) 3

Gostatistics: study of a PH contamination 3 Variogram Fitting Cross Validation Error 3D Variogram 4 Standard rror Man Varianc -44558 2988668826-355 96465 Prcntag of robust data Corrlation Cofficint Z/Z* 986 68

Gostatistics: study of a PH contamination 4 Intrpolation mthod: ordinary krigging Ptrolum Hydrocarbons -2 m 5

Gostatistics: study of a PH contamination 4 Intrpolation mthod: ordinary krigging Ptrolum Hydrocarbons 2-25 m 6

Gostatistics: study of a PH contamination 4 4 Intrpolation mthod: ordinary krigging Ptrolum Hydrocarbons 25-5 m

Gostatistics: study of a PH contamination 4 Intrpolation mthod: ordinary krigging Ptrolum Hydrocarbons 5-8 m 8

Gostatistics: study of a PH contamination 4 Intrpolation mthod: ordinary krigging Ptrolum Hydrocarbons 3D viw 9

Gostatistics: study of a PH contamination 4 Intrpolation mthod: ordinary krigging Ptrolum Hydrocarbons Concntrations > ppm - 3D viw Is th whol sit contaminatd? Why? 2

Gostatistics: study of a PH contamination 4 Intrpolation mthod: ordinary krigging Krigging rproducs th xpctd valu of th distribution: stimats ar strongly influncd by th man valu of th distribution Original data histogramm - PH Nb Sampls: Minimum: Maximum: Man: Std Dv: 5 s i c n u q r F Nb Sampls: Minimum: Maximum: Man: Std Dv: 5 s i c n u q r F 5 2 2 6 36 9 9868 5862 5 Krigging histogramm- PH 564-3564 62923 685365 2428222 4 3 2 2 3 4 5 6 8 Contamin C-C5 2 3 4 5 6 Contam HP stimatd Sara 3D omni

Gostatistics: study of a PH contamination 5 Intrpolation Mthods: Conditional simulations Krigging Conditional Simulations Krigging allows rproducing th xpctd Each simulation rspcts th valus valu of th distribution: stimats ar obsrvd in th sampl points Man and strongly influncd by th man valu varianc ar prsrvd in ach ralization Svral random fonctions whit th sam man and varianc can hav th sam variogram Allows only on contamination Quantifying possibl 22 th Simulations can rproduc th statistics of ordr (histogram) and 2 (variogram) mapping of th Prforming a larg numbr of simulations allows to bttr rproduc th ral divrsity of possibl cass, spcially if w ar considring a random variabl uncrtainty is not Th uncrtainty is quantifid so it is possibl to asss th risk of th projct

Gostatistics: study of a PH contamination 5 Intrpolation Mthod: Conditional Simulations Conditional simulations allow us to prform svral charts (with qual probability) to bttr rflct th divrsity of ral possibl cass Simulation : grid 5*5*5 m 23

Gostatistics: study of a PH contamination 5 Intrpolation Mthod: Conditional Simulations Simulation 3 : grid 5*5*5 m 24

Gostatistics: study of a PH contamination 5 Intrpolation Mthod: Conditional Simulations Simulation 25: grid 5*5*5 m tc 25

Gostatistics: study of a PH contamination 5 Intrpolation Mthod: Conditional Simulations Prforming a post-tratmnt of conditional simulations allows to ask th following qustions: Whr ar th contaminatd aras? This qustion is usd to locally stimat th risk of xcding th rmdiation thrshold What is th volum of contaminatd soil? This qustion is usd to stimat th ovrall lvl of contamination in th study ara How to answr ths qustions??? 26

Gostatistics: study of a PH contamination 5 Intrpolation Mthod: Conditional Simulations Whr ar th contaminatd aras? This qustion is usd to locally stimat th risk of xcding th rmdiation thrshold Stting a rmdiation thrshold probability map for th occurrnc of PH in soil (for ach bloc) in concntrations highr than 5 ppm 2

Gostatistics: study of a PH contamination 5 Intrpolation Mthod: Conditional Simulations Blocks showing a probability of xcding th rmdiation thrshold gratr than 6% 28

Gostatistics: study of a PH contamination 5 Intrpolation Mthod: Conditional Simulations Which risk lvl ar w willing to accpt? How to manag blocks showing an uncrtain stimat? Blocks having a probability of xcding th rmdiation thrshold (5 ppm) btwn 3 and 6% 29

Gostatistics: study of a PH contamination 5 Intrpolation Mthod: Conditional Simulations What is th volum of contaminatd soil? This qustion is usd to stimat th ovrall lvl of contamination in th study ara 9 8 P95 ( 395635) P9 ( 44235) P85 ( 432985) Optimistic volum s i c n u q r F 6 5 P5 ( 486935) Most probabl volum 4 3 2 P5 ( 539 P ( 554 P 5 ( 54 35 4 45 5 55 6 65 Volum (m3) Risk curv 3 Pssimistic Volum

Gostatistics: study of a PH contamination 5 Intrpolation Mthod: Conditional Simulations Prforming conditional simulations shows that: Th stimatd volum of contaminatd soil is btwn 54 m3 (pssimistic volum) and 395 63 m3 (optimistic volum) Th most probabl volum is 486 93 m3, which corrsponds to 25% probability of xcding th rmdiation thrshold (low risk) Th uncrtain volum (probability btwn 3% and 6%) rprsnts a total volum of 23 8 m3 This volum rprsnts 48% of th most probabl contaminatd volum announcd by th risk curv 3

Gostatistics: study of a PH contamination 5 Intrpolation Mthod: Conditional Simulations Th calculations show a grat uncrtainty of th contaminatd soil stimat 2 ) m ( Z 3 4 5 - Sampling barly prformd from 5 m dpth 6 8-32 9 2 X (m) 3 4 5

Gostatistics: study of a PH contamination 5 Intrpolation Mthod: Conditional Simulations If th stimat is stoppd at 5 m dpth: 9 8 P95 ( 2855) P9 ( 24) P85 ( 29) s i c n u q r F Th diffrnc btwn th pssimistic and optimistic volums is lss important, Th most probabl volum is 248,25 m3, which is half of th uncrtain volum calculatd abov, and corrsponds to approximatly 25% probability of xcding th rmdiation thrshold Th volum of contaminatd soil showing uncrtain stimation (3% to 6% probability of xcding th thrshold) is 66 35 m3 This volum rprsnts 25% of th most probabl contaminatd volum (P5) announcd by th risk curv 6 5 P5 ( 24825) 4 3 2 P5 ( 285) P ( 28345 P 5 ( 2984 33 2 25 Volum (m3) 3 35

Gostatistics: study of a PH contamination 5 Intrpolation Mthod: Conditional Simulations Mapping : Soil having 25% probability or mor of bing abov th rmdiation thrshold and corrsponding to th most probabl contaminatd volum 34

Gostatistics: study of a PH contamination 5 Intrpolation Mthod: Conditional Simulations Changing th grid: blocks siz = 25*25* m Blocks having a probability gratr than 6% of xcding 5 ppm 35

Gostatistics: study of a PH contamination 5 Intrpolation Mthod: Conditional Simulations Changing th grid: blocks siz = 25*25* m It is possibl to produc maps suitabl for dcontamination which will b basd on th risk w ar willing to tak in ordr to dfin a ral stratgy for th rmdiation projct Blocks with 3 to 4% probability of xcding th critrion Blocks with 4 to 6% probability of xcding th critrion Blocks with mor than 6% probability of xcding th critrion 36

Gostatistics: study of a PH contamination Comparison btwn diffrnt mthods: Polygons vs Gostatistics Polygons mthod prformd according to th obsrvd gology (thicknss of th fills layr) Total volum 3 m3 to 26 3 m3 3

Gostatistics: study of a PH contamination Comparison btwn diffrnt mthods: Polygons vs Gostatistics Risk curv Risk curv prformd on th fills for th north ara, along th railways P95 ( 6525) P9 ( 685) P85 ( 25) 9 8 s i c n u q r F 6 5 P5 ( 95) 4 3 2 P5 ( 8595) P ( 8865) P 5 ( 925 38 6 8 9 Volum (m3)

Gostatistics: study of a PH contamination Conclusions Prforming th gostatistical study highligths th following information : Gostatistics provid a cohrnt analysis of data collctd at ach stp of a contaminatd sit charactrization by prforming statistical procssing of spatial data Gostatistics ultimat goal is to dtrmin th uncrtaintis associatd to intrpolation mthods in ordr to rduc th financial risks of rmdiation sits projcts Th choic of th gostatistical mthod usd for ach projct has a larg influnc whn making stimats, spcially if th original data dos not show a clar spatial structur To mak rlvant calculations, it is ncssary to conduct a comprhnsiv analysis of th input data bfor intgrating it into th gostatistical tratmnt Each choic must b justifid by th gostatistical thoris but also judgd by th xprts working on contaminatd sits 39

4 Thank you for your attntion