kalman filter in the real time urbis model
DESCRIPTION
Kalman Filter in the Real Time URBIS model. Richard Kranenburg Master scriptie June 11, 2010. Kalman Filter in the Real Time URBIS model. Introduction Real Time URBIS model Problems and Goals Method Kalman filter equations Results Extensions on the Kalman Filter Conclusions. - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/1.jpg)
Kalman Filter in the Real Time URBIS modelRichard KranenburgMaster scriptie June 11, 2010
![Page 2: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/2.jpg)
Kalman Filter in the Real Time URBIS model Introduction Real Time URBIS model Problems and Goals Method Kalman filter equations Results Extensions on the Kalman Filter Conclusions
![Page 3: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/3.jpg)
Introduction
Company: TNO Business Unit: ‘BenO’ Department: ‘ULO’ Accompanists:
Michiel Roemer (TNO) Jan Duyzer (TNO) Arjo Segers (TNO) Kees Vuik (TUDelft)
![Page 4: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/4.jpg)
Real Time URBIS model
![Page 5: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/5.jpg)
Real Time URBIS Model
URBIS Model, standard concentration fields 11 sources, 4 wind directions, 2 wind speeds
![Page 6: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/6.jpg)
Real Time URBIS model Every hour interpolation between standard
concentration fields Correction for background and traffic fields
μ is the weight function dependent of wind direction (φ), wind speed (v), temperature (T), hour (h), day (d), month (m)
: standard concentration fields
88
1,
iiki
mk mc
im
![Page 7: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/7.jpg)
Real Time URBIS model
![Page 8: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/8.jpg)
Real Time URBIS model
Linear correction used by DCMR Average concentration of three stations
Schiedam Hoogvliet Maassluis
mmslkmslk
mhgvkhgvk
msdmksdmk
mk
DCMRk cycycycc ,,,,,,3
1
![Page 9: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/9.jpg)
Real Time URBIS model
![Page 10: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/10.jpg)
Problems and Goals
The model simulation can become negative No information about the uncertainty of the
simulation
Goal: Find an uncertainty interval for the concentration NO , which does not contain negative concentrations
x
![Page 11: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/11.jpg)
Method
Kalman filter connects the model simulations with a series of measurements
Kalman filter corrects the model in two steps Forecast step Analysis step
Result is a (multivariate) Gaussian distribution of the unknown
Mean Covariance matrix
![Page 12: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/12.jpg)
Kalman filter equations
Correction factor ( ) for each standard concentration field
Kalman filter calculates a Gaussian distribution for the unknown variable ( )
The concentration NO can be found in a log-normal distribution
88
1,
,
iiki
KFk
kiemc
x
![Page 13: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/13.jpg)
Kalman filter equations
Second equation not linear in ( ), thus a linearization around
H: projection matrix, A: correlation matrix represents the uncertainty of the
measurements on time k
),0(~
ln)ln()ln(88
1,
1
,
kk
ki
ikikkk
kk
R
emHcHy
A
ki
Ν
0
kR
![Page 14: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/14.jpg)
Kalman filter equations
The linearization results in:
with:kkk
kk
Hy
A
~~
1
88
1
,~
)ln()ln(~
j
mk
jkj
mkkk
c
mH
cyy
![Page 15: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/15.jpg)
Kalman filter equations
Forecast step
represents the uncertainty of the model is the covariance matrix after the forecast step The temporal correlation matrix A is determined
with information from the measurements
),0(~T1
1
QNAAPP
A
kfk
kfk
Q
kP
![Page 16: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/16.jpg)
Kalman filter equations
Analysis step
Minimum variance gain
TT)~
()~
(
)~~(
KKRHKIPHKIP
HyK
kfkk
fkk
fkk
1TT )~~
(~ k
fk
fk RHPHHPK
![Page 17: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/17.jpg)
Kalman filter equations
Start of the iteration process:
Screening process: Before the analysis step is executed, the measurements are
screened. If difference between simulation and observation is too large,
that observation will be thrown away. In this application the criterion is twice the standard deviation
QP 00 ,0
![Page 18: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/18.jpg)
Results
![Page 19: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/19.jpg)
Results
For the whole domain on each hour an uncertainty interval for the concentration NO can be calculated
Annual mean of the widths of these uncertainty intervals
Population density on the whole domain Connection between population and
uncertainty
x
![Page 20: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/20.jpg)
Results
![Page 21: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/21.jpg)
Results
![Page 22: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/22.jpg)
Results
Connection between uncertainty and population
Kalman filter reduced the uncertainty Absolute uncertainty: 14.5 % Relative uncertainty: 16.1 %
gp
gp
n
iii
n
iii
uU
uU
1,rel
rel
1 abs,
abs
pop
pop
![Page 23: Kalman Filter in the Real Time URBIS model](https://reader035.vdocument.in/reader035/viewer/2022062308/56812d9e550346895d92bfb3/html5/thumbnails/23.jpg)
Extensions of the Kalman filter
Goal: Minimize the uncertainty connected with the population
Methods: Add extra monitoring stations to the system Apply Kalman filter with different time scale and
add stations with other time scales Analyse the values of the correction factors