statistical framework is proposed for the geotechnical inverse problems. ...
interest for solving mathematical inverse problems, not for data interpretation.
! " • • • # $ % # "& • ! • " • " " • # $ % $ • & %" % '("( )( * (+ " • ' ($ • ' ($ (++( !( • ! $ $ • $
, ) -. "( • )! ' • )! ' * /& ( 0+( (( &
* % ($+,)- % ($. ,/- $ $ 0 % ($1 ($ %$1 % ($* $ $ $ %$ $ %$ $ 1 % 2 $ % %$ $ % $ 1 $ $0 % ($*$ %$ $ $ $1 3* $ $* % * $ ( $ $ ( 0 $ ( ( $ $ ( 0 1 4 % ($ $ 0 %1 4 % ($ ($. $ $ 0 % ($ ( ( ($1 # (% $ %$ (% $ 0%$ $ 0 %14 % 15 $ ( 2 (%*%$ $, $ - , - ( $ $ 1 % 0 % 6 (% 0 % $ $1 4 * ($ $ ($13* ($ 0 % % %$$ 1
4 $ ($ * $ * $ $ % $ $ $ * ( $ 1 ($ ( (2 % 13* % 1 4 $$ 0 $ ( "5" % $ 1 1,)778* )779- % % $ $ % 1 1 ,)778- $ % $ 1 * $ ( $ $ ( (% $$6 0 $ 1 3* 0 $ 2 1 4 $ 0 % $ $ ,11 ($-1! % $ ($14$ ( ($ 0 % 1 $* $ $ 1 4$ 0 $ ( $ ( % $ 1 4 $ ($ (% % (( % ( (%* ( % % ,3 * )7:;-1 & $ %$ $ 0 %* (% % %$ ( $ 14 * $ ($ $ (%( * % $
%$1
4 6 $ $ ( (% 1! (% d $ $(% m1 ! +d $+m
=) >
=)(>
( 6 $( $ $1 4 $ $ $ % 1 4 $1 & % $ $ % $ % $ $ $ $ $ 1 %$( %* $ $ d = (m )
=/>
,1- $1 ( 1 4 ($ $ $ $ * * , $-*$$ $ $ )1 ! " $ $ ? @ A A1 4% ( % ($ (% $ $ ( , %- $1 4 ( % % $ ($1 # $ $ ($ , ($ $ 0 % $ -1 $ $ $ $ * % $ ( 1 # % $ %
$ ($ % $
,B %*)7:C*)7:9.D % B *)77;-14$ % % $ 1"% %
$ ( 14% $ $ ($* 1
" ,"- %$ (%$ ,( - ( 1 " $( 14 $$% 1 $ 1
$ 1 4 $ ( 1 4 $ $ * $ 1 / % % 1 %1 4 % * % 1 $ %* % ( (%
$ $ ( 1 4 (% 1 $ % ( $1 ( (%$ 14( $ $ ( $ ,E$ 1*)7:F-14 % $ 14 (%$ % ( 1 ( $ % ( % ( % 0 %$
$ % % , $ (- %
%* %$ % (,3 1 1 )7:7-1 4 $ ( $ % ( $ 1 (% ( 1
( *4 % % $ $ ( (1 ( ( %14$ ( % $
1 ( (%
%1 G 0 $ ( 1 & * % $ 1 3* 0 $ $( % $
$ (01
($ % ( % $* ( $ $ (14$ $ $% (% (% $ 1 % ( $ ( $ $ 0 ( $ ( % ( $( , *)7:7.4 *)7:;-1" (* $ $ $ % ((%0 =/>1 $ $ 0 $ % $ $ 0% % 1 # % 0 $ ($ %$ %1 ( $
%(
$14$ $ ( $ 6 1
0 $ 1 ,)7:7 $ ,)77;- $1 4 $
% %6 ($
% $ 1 #
% ($
$14 $ ( 6 +
% 1 4 ,)7:;- $ ($ $ % % $ $ ( $( ( ($1 4 % ($ $ % $ ( ($1 ($ $$ 6 C1 4 $ ($ C1( ( 14 % 1 $
% $1
# $ % # "&& ! ! 4 ( ($ ($14
$ ,m- ( % $ $1 4 ,d- $ ,g(.)- % $ 1 * $ $ * φ $( * τ = + σ φ 1 $ $ $ 9H)HHI *=/>( $ τ ) 9H τ = ) )HH φ
=C >
$ $ $ 9H*;9*)HH* )/9 I *=/>( $ τ ) 9H τ ) ;9 = τ ) )HH φ τ ) )/9
=C(>
$ $ $ 9H*9H*)HH* )HH I ,11 $ $ -*=/>( $ τ ) 9H τ ) 9H = τ ) )HH φ τ ) )HH "
=C >
%* $ % $ $
$1# 0 % $ (1 "
0 =/> =C> ( $ $1 3* $ (2 %1# % % $ ,4 * )7:;-1 $$ $ ($ $ (1 5 $ $ , ( (% ( $
( -(% $ 1 Dcorrected = ,m - + ∆ m
=8>
0=8>( ( (%%* $ ( $ $ ( (% $ (15 $ $ 0=8>
$ $1 Dobs = Dcorrected + ∆
=9>
( 0 =8>=9> Dobs = ,m - + ∆ + ∆ = ,m - + ∆
=F>
∆ ,$ ∆ ∆-1 4 * ( %$ ((% 81
! # #
$ 4 ( $ $ % + ,)- ,/- % 1 $* $ $ $ ( 1 4 $ $ $ $ $ ( $* $ $$*$ $$ * & J& ,& $ $$ J &( $ -1 =)77H> $1 4
% $ $ $ ( (% (1 % $ % ( $ $(% ( (%% $ 1*$ $ $ ( ∆ Dobs1 8 ( $ (% mM $ ( 914 * 0 =F> ( $ Dobs = ,M - + ∆
=;>
& % * M Dobs ( $ 2% ( $ ( 2%!,dobs*m-1 $( Dobs$ * $ $ 1E $ $* $ , $ - ( M ,11 ( M-1 " $$ ( * ,% $ - ($ $ $ (($ 13 * % ($*
0 + $ $ $ $ ,!,m--* %$, % *,M-K∆ ∆- $( ( $ $ $* $ ( $ $ ( $ $ ( $ ( L % ( (% ( (% * % $ ( (%, ( (%( $ $
( -, *)77)-
!(m M d obs ) =
! (d obs M m ) ! (m ) = !(d obs )
! (d obs M m ) !(m )
∫
Λ
! (d obs M m ) ! (m )m
=:>
"
%*$ ( % %$ $ $ $ 01 !
5 ( %$ ,M-L ,5 $ , -L-
!
3
% ,∆ ∆-L,5 ,dobs|m-* ( (% =;> ( (% ( ∆1
$ $ $ L-
0 ( ( $ * ( (% ( $ $ ( (
(% 0 =:>1 0 =:> % % % % ($14% 4 =)7:;> (% A $ $ ($( ( (% * % ( (%1A 4 ,)7;/- ,)7;8-*
0 % % $ ( % %1 " % * %,)7;9- % (/) % 1 " & % $* ($ % ( $ % % ,( (% $ $ $-1 0 =:> ($1 4$ % % $ ( $ (% % ( (%%!,m|dobs-1 % % $$ $( $ $ 1# $ $ ( (% $ ( (%( $ $1!$ (% Λ* m∈Λ Λ,Λ×Λ×N×Λ -*$ ( (
!( M d obs ) = ∫ ∫ ! ∫
− ∫
+ ! ∫ ! (m M d obs )
∫ =
− ∫
+ ! ∫ ! (d obs M m ) ! (m )
∫ !(d
M m ) ! (m )m
Λ
Λ
Λ
Λ −
Λ +
∫ ! ∫ Λ
Λ −
Λ
Λ +
Λ
obs
Λ
=7>
* 0 =7>$$ $ 1 ,dobs|m m $ ,,m--* $ (% $ ,m- $ 1 $ $$ $ $ 0 ,)779* )77F-1 4$ AA ($ %
1
# $* ( A(A$ $1 %! 4 $ % $ (% $( $ ( $ ($ $ ($$ 1!( (% Γ $ (% Λ* dobs ∈Γ Γ m∈Λ Λ1 4 $ $$ mO + Γ → Λ mO ,d obs - ( $ $ $ ,mtrue- dob 1 4 $ % $$$ 1 $ ($
+ E×E" E *
$ * ,O - $ &$ ' $%
F14% O = (O − ) , 0 -
=)H >
(O ) = O − ,"( -
=)H(>
H O − ≤ ε ,& $ (O ) = O ε ) − >
=)H >
' * , $ $ $-* ( (O ) ( (O ) = ){ (O , Dobs -* ) M = } )*+,- ./ +- .' $ 6 $ (*4 ( (O ) = ) {( (O )} = ) {) { (O M ) M = m}}
= ){ (O M )} = ) {) { (O M ) Dobs = d obs }}
=))>
4 $ $$6 (O ) 1 4 $ % ( (% $$6 ) { (O M ) Dobs = d obs } 0 =))>14 * $ $ O (d obs ) = $ ∫ (O ,d obs - , m ) ! (m d obs )
Λ
=)/>
$ $$6 B, -1 (
* $ ,)7::-* (% 0 =)H>* =)/> ( $ $ * $$ $ $$ ( (%(* %14 $ ( $$ $ ,11 m - ,)7::-1 4 $ $$ 0 , - $ * $$ "( , " - $ $$ " , "$ * %1 mO = ∫ m ! (m | d obs )m
=)C >
mO = mO O = = ! , M d obs ->
=)C(>
mO = $ ! (m | d obs )
=)C >
Λ
m
* $ , - $ , " -0 $ 1 $ $$ , "- ( $6 $1 "$ $ $ " $ $ ( $6 $ ( " 1 % $ ,1-* ( ( % $ $ ,1-1 % $$ * $6 $ % $$$1'( $6 0
$ " ' "$ ( $$$ ,)779*)77F-1 0 $"$ 5 ($ $ $ , $ * $ $ $$ -* % %1 4 $
$$ % % (1 M | Dobs = ∫ (m − mO MMSE )(m − mO MMSE ) ! ,m | d obs -m
Λ
=)8>
3* % $$ * $ %(% $ 1 4 $ $ % ( (% % ( (% $ $ $ $ 14 ( (% $ $ $ $ *(%Λ Λ2* 1 (m ∈ Λ P) = ∫ ! (m M d obs )m = Λ
∫ !(d ∫ !(d Λ
obs
| m ) ! (m )m
Λ
obs
| m ) ! (m )m
=)9>
'("( )( * (+ " " (*$ =7>*$ =)C >*$ =)C(>* =)8>* =)9> 0$ 13* % ( %( $( ($14 ($ ($ $ ( $ ' (1 "$ ' ( $ 1 " $ (' ( !( x ) =
)
) − − (x − µ ) (x − µ ) / ,/π -
=)F>
N,µ µ* -* µ µ$ $ 1 "$' ( $M ∆,11 MQN,mp* -*
∆QN,H* ∆--*!,m-QN,mp* - !,dobs|m-QN,,m-* ∆- ( (% %=:>( $ ! (m d obs ) =
)
! (d obs )(/π )
∆
⋅
[
]
) − − − (d obs − ,m - ) ∆ (d obs − ,m - ) + (m − m p ) (m − m p ) /
=);>
4 m "m( 14 % $ ' $ ( 1 ( ( ' ($1 4 ' ($ ( (% 6 1 2 3 ! & m "m %! ) ,m-
" $ $ * M ∆ ' 1 $ $$ $ ' ( , "
-* ( Dobs
M 2%' m p M N Q * D "m p " obs
" " " + ∆
=)7>
Mest
4 *Mest QN, mO * t-1 mO * " " $ 1! "
"
−
(
= " " " + ∆
)
−
=/)>
0 =/H> ($ − mO = m p + " (d obs − "m p )
=// >
−
M est = − " "
=//(>
( 0 %, " "
−
(
= " " " + ∆
(
−
= " ∆ " +
)
− −
-
)
−
" ∆
−
=/C>
& 0 =// > =/C>* ( ' "( (1 =)7:7>*4 ( ) "( $ ( ( (%
' $ $ $1
$ * 0 $ $ ,11 MN ∞ ' ,11∆N1
4 ()1E ( ' ($ "( Dobs
= "M + ∆ MQN,m1 -*∆ ∆QN,3*ε ∆-
= m p + m !
−
(d
− − m p ) * = ( + ∆ )
−
(
obs
! "# $'(
"# $%&
[
]
mO = " "
! %)&
−
[
m =
+&,
]
MN∆ ∆N
−
d obs
MNε∆ ∆N
−
MN∆ ∆Nε
obs
[
.#/-!
mO = " 4 "
%)& .#/-+& , .#/- ! %)&
]
−
[
] [d −
]
obs
− m p
[
−
= mp + +ε dobs − mp m −
[
MN∞∆ ∆N
" 4 d obs
= m p + − − m
[
−
MN∞∆ ∆N
−
∆
" d obs
[ ] = [ +ε ] d
%)&
)
− −
)&( *
= +ε dobs m
-!
−
= ∆ +
] [d
− −
−
= mp + +ε
obs
]
MNm∆ ∆N
]
MNmε∆ ∆N
MNm∆ ∆Nε
]
− mp
, 0 /! "# $%& $ - 1#ε2/ dobs m dobs m -m m (m m 3%+4 5 .
2#/# 6
.-2#/# 6
ε- #$
0$ $ 7( $/&$
--
8 9 6$ # 1 : 8∞9 6$ # ∞ :
'& 6
0 $"$ &0 % 0 =/C>* $ (
(
M est = − " " " + ∆
(
−
= − " ∆ " +
(
−
= " ∆ " +
)
)
− −
)
−
" −
" ∆ "
=/8>
− −
4 * % $* ( $ ( %( (%=/8> 4 (/14$ 0 $
, - ? ( @ ( $ $1" $ $ %1 (
ρ =
=/9>
% 56 ≤ ρ ≤ 6' ρ 7 $ * ρ 6 $ % * ρ 86 $ % 1" $ $ ( % (% % $
$( $1 4 (/1 "( Dobs
= "M + ∆ MQN,m1 -*∆ ∆QN,3*ε ∆−
;( $ !
(
Mest = − ( + ∆ ) = − ∆ +
(
−
= ∆ +
)
−
)
− −
−
∆
− −
)&( *
;( $+ 6
! %)&
MN∞∆ ∆N
+&,
MN∆ ∆N
-!
MN∆ ∆Nε
ε
%)&
MNε∆ ∆N
ε ε
MN∞∆ ∆N
MNm∆ ∆N
MNm∆ ∆Nε
MNmε ∆ ∆N
ε
ε ε
"# $%&
.#/-! %)& .#/-+& , .#/- ! %)& ,
$ $ 7( $/&$
--
8 9 6$ # 1 8∞9 6$ # ∞ '& 6
( ! $ $ % dobs ( mp1 $ ( mp 14$( $ 2 ( % $13*(2 $* $ * ( $ $ 0 0 1 mp % 2$ 0 $ , $ % P -1 4 % $ $ mp ((% $ 1 4 $ ( $ $ $ mO $ $ mtrue , * )7:7-1 $ $ % $ $mtrue,1'dobs "mtrue9δ-* − mO = m + " (d obs − "m )
= m + "
−
("m
true
−
+ δ − "m p )
=/F>
−
−
= " "m true + , − " "-m + " δ δ δ 6 $ ∆ ∆*
$14 ( −
(
−
( = " " = " ∆ " +
(
)
− −
= " " " + ∆
−
" ∆ "
)
−
"
=/;>
4 * $ $ $ 0 =/H> =//> ($ $ − mO = (m true + , − ( -m + " δ
=/: >
= − (
=/:(>
0 %* $ $ 0 $ $ " δ % % $ ( $6,11 ( $6$ -1 *$ ( %1"% $ $($ 6 %1 4 $ %*
1" '(,)7F:-* 14$ 2$ (%
$ $15 15 % 13* ($1" $ %* $ ( $ $ % 2$ %$$ $ $1 $$ $ % ( $ $ $ ,11 ($-1 ( %6$ $ (% $,,1-P-14 $ $ 1 D prd = "M est ( $ ( (%% 1 $ 0 =/H>* $ − ) {D prd }= O = "mO = "m + "" ,d obs − "m -
=/7 >
{D prd }= = " " = ", − " " - "
=/7(>
−
4 ( ,11 $ O - , -1 " 0 =/;>* , *)7:7: = ""
−
(
−
= " " ∆ " +
(
)
− −
= " " " " + ∆
" ∆
−
)
−
=CH>
4 * 0 =/:>* $ $ $ 0 =/7> ($ O = :d obs + , − : - "m p = :d + , − : - "m p + :δ = , − : - " " = , − : -, − ∆ -
=C) > =C)(>
4 $ $ ( 1 $
$ $ ( ( $*$ ( % $14 % $ ' ($* % $$ 64 (C1 4 (C1 $$ % $ ' ($ , * % ,dobs- m1 * $ $6 =);>0 $$6$ 1
{
}
− − mO MAP = $ (d obs − ,m - ) ∆ (d obs − ,m - ) + (m − m p ) (m − m ) m
=C/>
0 $"$ $$ % $ % $ ,1-*
$ $ %(%$ 1 * 0 =7>=)9> ($ % %6 %1 3* ($ %,11 ($ ,m- 6 ( ( (%-* % % (
$
(% $ , - 6 ($ ,4 *)7:;-1 mO ( "$ ( (%=C/>14$
(
$ (% 4 %P
$
∂ ,m - ≅ ,mO - + ,m − m ∂ m = ,mO MAP - + " ,m − mO = ∂ ∂ m
=CC>
% m = mO MAP 1 % %
$1 & =CC>* ( (% % =);> (
$ (% ' (
) − ! ,m | d obs - ≅ ⋅ − ,m − mO MAP - ,m − mO MAP - / ( (%(" " =/8>14
=C8> $
$ ( $
(
≅ " ∆ " +
)
− −
( ≅ − −
=C9> =CF>
* % (
$ %(%(" "1
(++( !(
+!* (!& ( ,' $ * $$ %TE 6 -
, ) )! $ 6 ( (% % 1 " ' $ $ ($ 1 4 % )! ' ($14 $ $ (%% 0 =C>1 -. "( 682 3 ) 1 $ )! ($* "9∆ QN,* σ- ∆QN,H* σ∆-1 4 2 ( (% ! $ (% =)7>1 4 ( (%(' $ (%=/H>1
"14
( 6 $ 7, -1 4 ( ($ ( $ 6 $ 1 4 $6 $ " $
$
% % ( ' ($1 $ % 2 ( (% * ( (% 6 '
$ 7,-* ,-* , - 7'777;8
7'6;' 4 ($ $$ 6 $ % % %1 % ($* ( (
( $$ "$ % %( % % $ 1 4 0 )! ($ )H1 ( $$ ($ 01 $ % $ 0
(% %$ %$, $ 7, -)H,(--1 4 * $ $ % $ % $ ($ % $1 * /& 5 $ * $ ($ $ $1%=C(>( ( 1"$ $ $ (( $ 14( $ $ % $ 1 "$ $ ( $
9H = CH /9 H = H , 9- / H ∆ = H / 4$ 0 =C> ( (%=/;> H17;C7 )C1/:H7 (= H1HHHC H1://8
=C; >
H17:8C ;1HCH8 (= H1HHH/ H17)79
=C;(>
H17:9F ;1C7/8 (= H1HHH/ H17H)C
=C; >
4 =C(> ( =C > =C > ( =C >13* 0 $ % 0 % ( $ 0 % $ 1 4 $ % (
( ( % 1 % $
( $ $ ((1 (&0+( ,&4!E$ $ (( & ,&$ $ $
)
4$$ % ( $ %$ 1# $ (* % %$ 13*$ $ 0 % ($ $ $ ( % %$1 4
$ ($ 1 %$ ($* ( $ % $ $*
$ (%( 14 $ ($ %$1 4 %$ % $1 4 $ ( $ % $ (% $ $ % %1 "$ ' $ 1 4 ( $ $ % % $ 14 ' ($ % % $ $ $$ 0 14 % 0 $ (% $ $ 14 % $ $ ($(
$$ % 14 % $ ($ $ $* $ 1 * % $ $ $ ( ( $1 $ $1 $ $ %* $ ( $?@ ($1
, $$ ) ' !( V J VQ,µ*Σ-*JQ,H*Σ-13 $ $ = + 14 V W 2%' $ (% Σ µ Aµ AΣ
Σ A AΣ A + Σ
$$ / ' !(
V B 2% ' $ $
(% µ Σ µ Σ
Σ Σ
4 $ V B
