savel'eva geostats
Post on 05-Apr-2018
215 Views
Preview:
TRANSCRIPT
-
7/31/2019 Savel'eva Geostats
1/55
..
()
-
7/31/2019 Savel'eva Geostats
2/55
;
(
);
.
-
7/31/2019 Savel'eva Geostats
3/55
, (Z*(x0)).
().
-
7/31/2019 Savel'eva Geostats
4/55
(Z(:)), S;
xS Z(x),
();
Z(x) ; ,
Z(xi): z(xi)=zi, i=1,,N.
-
7/31/2019 Savel'eva Geostats
5/55
min)]()(*[ 00 xZxZVar(...))),(),...,(),(()(* 210 NxZxZxZFxZ
,
:
0 )}()(*{ xZxZE
(... - ,
Z(x)
-
7/31/2019 Savel'eva Geostats
6/55
()
.
. ,
h.
,
.
-
7/31/2019 Savel'eva Geostats
7/55
C(x,h)
C(x,h)= E{(Z(x)-m(x))(Z(x+h)-m(x+h))}
, :
12
2
( , ) ( ) ( )
( ) ( )
Var Z Z
E Z Z
x h x x h
x x h
-
7/31/2019 Savel'eva Geostats
8/55
BLUE:
Best
Linear
Unbiased
Estimator
)()(* j
n
j
j xZxZ
1
0
min)()(* xZxZVar
0 )(*)( xZxZE
-
7/31/2019 Savel'eva Geostats
9/55
:
1. E[Z(x+h)-Z(x)] = 0, x2. h:
2(h) = Var[Z(x+h)-Z(x)],
0111
mxxZxZxExZxZE
xn
i
i
xn
i
ii )()()()()()(*
)()(
11
)(
)(
xn
i
ix
-
7/31/2019 Savel'eva Geostats
10/55
( )
01
( )
1
1
n x
j ij ij
n x
j
j
i=1,,n(x)
i
iiOK 0
2
-
7/31/2019 Savel'eva Geostats
11/55
;
;
; ;
2{ *} { }
OKVar Z Var Z
2{ }OK Var Z
-
7/31/2019 Savel'eva Geostats
12/55
Rhifcc
RhifR
h
R
hcc
h
,
,2
1
2
3
)(
0
3
0
-
7/31/2019 Savel'eva Geostats
13/55
-
7/31/2019 Savel'eva Geostats
14/55
-
7/31/2019 Savel'eva Geostats
15/55
()
2 : ;
.
.
-
7/31/2019 Savel'eva Geostats
16/55
)(
1
2)()(
)(2
1)(
h
hxxh
h
N
i
ii ZZN
N(h) h.
hb
w
4h
()
(bw)
h
h
-
7/31/2019 Savel'eva Geostats
17/55
-
7/31/2019 Savel'eva Geostats
18/55
- B,
)()(
T
hBhh
hh
cossin
sincos0
01
m ax
m in
aa
-
7/31/2019 Savel'eva Geostats
19/55
nugget
sill
Range 1 Range 2
(Sill)
(Nugget) (Range)
-
7/31/2019 Savel'eva Geostats
20/55
( )*
0
1
( ) ( )n x
j j
j
Z x Z x
( )
0
1
( )
1
1
n x
j ij i
j
n x
j
j
i=1,,n(x)
i
iiOK 0
2
-
7/31/2019 Savel'eva Geostats
21/55
;
(
);
(-);
( );
(
);
.
-
7/31/2019 Savel'eva Geostats
22/55
)()()()( xyxaxaxm 10
)(
)()()(*xn
j
jj xZxxZ1
)(
)(
)(
)()()(
)(
)(,,)()()()()()(
xn
i
ii
xn
i
i
xn
ijRjijRi
xyxyx
x
xnjxxCxyxxxxCx
1
1
110
1
1
-
7/31/2019 Savel'eva Geostats
23/55
1,400,000 2, 400 .
(C) .
-
7/31/2019 Savel'eva Geostats
24/55
= 0.874, RMSE = 3.13 = 0.984, RMSE = 1.19
-
7/31/2019 Savel'eva Geostats
25/55
-
7/31/2019 Savel'eva Geostats
26/55
-
7/31/2019 Savel'eva Geostats
27/55
90Sr286
137Cs680
-
7/31/2019 Savel'eva Geostats
28/55
1: 137Cs
( )
137Cs
( )
( )
-
7/31/2019 Savel'eva Geostats
29/55
137Cs 10Ci/km2
-
7/31/2019 Savel'eva Geostats
30/55
2: 90Sr
-
7/31/2019 Savel'eva Geostats
31/55
-( )
*
1 1
( ) ( )n xK
i i
i
Z x Z x
K
)(
1
)(
1
* )()()()()(xn
j
jj
xn
i
ii
VU
xVxbxUxaxU
2 U(x) V(x)
)(
1
)(
1
)(
1
)(
1
)(
1
)(
1
0
1
)(,...,1),,(),(),(
)(,...,1),,(),(),(
xn
i
i
xn
i
i
V
xn
i
xn
i
jUVVjiVijiUVi
U
xn
i
xn
i
jUUjiVUijiUi
V
U
U V
U V
b
a
xnjxxxxbxxa
xnjxxxxbxxa
)(),()(),()()()(
1
)(
1
2xxxxbxxxax
U
xn
i iUVi
xn
i iUUiOCK
VU
-
-
-
7/31/2019 Savel'eva Geostats
32/55
90Sr
(- )
90Sr
( -)
-
7/31/2019 Savel'eva Geostats
33/55
-
3
-
7/31/2019 Savel'eva Geostats
34/55
3:
( )
:
Kkz
zxI kjkj ,...,,otherwise
)z(xif,);( 1
0
1
0
0
I0
1
1
0
zk
Z
-
7/31/2019 Savel'eva Geostats
35/55
:
);())(|;(*))(|;(* kikk zxInzxFnzxI
-
7/31/2019 Savel'eva Geostats
36/55
-
7/31/2019 Savel'eva Geostats
37/55
-
7/31/2019 Savel'eva Geostats
38/55
:
629
-
7/31/2019 Savel'eva Geostats
39/55
)(|)(,,)(Pr))(|,,,;,,,(
nzxZzxZ
nzzzxxxF
NN
NN
11
2121
-
7/31/2019 Savel'eva Geostats
40/55
))(|;())(|;())(|;(
))(|,,;,,(
nzxFNnzxFNnzxF
nzzxxF
NNNN
NN
1111
11
21
-
7/31/2019 Savel'eva Geostats
41/55
24
50.2
14.9?
24
50.2
14.9
24
50.2
14.935.2
?30.1
39.121.1
-
7/31/2019 Savel'eva Geostats
42/55
: 5.32 110 1.46
2158 77 1.33
-
7/31/2019 Savel'eva Geostats
43/55
-
7/31/2019 Savel'eva Geostats
44/55
-
7/31/2019 Savel'eva Geostats
45/55
-
7/31/2019 Savel'eva Geostats
46/55
-
7/31/2019 Savel'eva Geostats
47/55
-
7/31/2019 Savel'eva Geostats
48/55
+ :
-
7/31/2019 Savel'eva Geostats
49/55
+
3D
-
7/31/2019 Savel'eva Geostats
50/55
90Sr: 5
1, 5, 15 cm
90Sr: 20
10, 15, 30 m
-
7/31/2019 Savel'eva Geostats
51/55
0.2 Ci/km2 : 1-2 cm0.5 i/km2 :
( 10 )
-
7/31/2019 Savel'eva Geostats
52/55
-
7/31/2019 Savel'eva Geostats
53/55
;
.
-
7/31/2019 Savel'eva Geostats
54/55
;
; ;
;
.
-
7/31/2019 Savel'eva Geostats
55/55
( );
;
- ;
;
.
top related