DEFORMACIJSKA ANALIZA PO POSTOPKU FREDERICTON DEFORMATION ANALYSIS: THE FREDERICTON APPROACH Anja Vrečko, TomažAmbrožič UDK: 528.33 IZVLEČEK V članku je opisan postopek Fredericton, ki je eden izmed postopkov deformacijske analize. Za opis razlik med geodetskimi opazovanji ali koordinatami točk v geodetski mreži v več terminskih izmerah predpostavimo več mogočih deformacijskih modelov. Na podlagi statističnih testov in razpoložljivih informacij o dogajanju na obravnavanem območju izberemo najboljši deformacijski model. V članku je najprej podano teoretično ozadje postopka, nato je postopek uporabljen na primeru generiranih opazovanj v dveh terminskih izmerah. Rezultati postopka Fredericton na obravnavanem primeru niso bistveno odstopali od rezultatov postopkov Delft, Karlsruhe in Hannover. KLJUČNE BESEDE Klasifikacija prispevka po COBISS-u: 1.01 ABSTRACT In this article, the Fredericton approach to deformation analysis is presented. It is possible to use several deformation models to determine the differences between the geodetic observations or between the coordinates of points in geodetic network in more epochs. The most appropriate deformation model has been chosen based on statistical testing and available information about dynamics at the area of interest. First, a theoretical background of the approach is described. Then it is applied to the generated observations in two epochs. In the present example, the results of the Fredericton approach differ only slightly from the results obtained with the Delft, Karlsruhe in Hannover approaches. KEY WORDS deformacijska analiza, postopek Fredericton, računski primer deformation analysis, Fredericton approach, numerical example 1 UVOD Postopek Fredericton so razvili A. Chrzanowski, Y. Q. Chen in J. Secord na oddelku za geodezijo univerze New Brunswick v mestu Fredericton v Kanadi. Avtorji so želeli zagotoviti splošno uporabnost postopka. Pri njegovem razvoju so postavili naslednje zahteve (Chen, 1983): • Postopek mora biti uporaben za absolutne in relativne mreže. Absolutne so mreže, pri katerih nekatere točke, ki jih imenujemo osnovne (referenčne), ležijo zunaj območja deformacije, se ne premikajo in določajo koordinatni sistem. Pri relativnih mrežah pa vse točke ležijo na območju deformacij ali na objektu oziroma nismo prepričani glede njihovega premikanja (Mozetič in sod., 2006). >Î5 S S- S jg SÜ o ■ to • Postopek mora biti uporaben za katerikoli tip deformacij, tako za določanje premika posamezne točke kot pri analizi deformacij bloka točk. • Postopek mora biti uporaben za obravnavo 1D- , 2D- in 3D- geodetskih mrež. • Postopek mora omogočati vključitev geodetskih in fizikalno-mehanskih opazovanj (nakloni, napetosti ...) v izračun deformacijskih parametrov. • Postopek mora biti uporaben za poljubno geometrijo geodetske mreže, tudi za primere nepopolnih mrež z defektom konfiguracije. • Postopek mora biti uporaben ne glede na določitev geodetskega datuma (vpeta mreža z minimalni številom vezi, prosta mreža) v koraku izravnave. Osnovna ideja postopka je na podlagi statističnih značilnosti polja premikov izbrati najustreznejši deformacijski model z metodo najmanjših kvadratov določiti deformacijske parametre, ki jih (razen parametra rotacije) obravnavamo kot neodvisne od določitve geodetskega datuma (Marjetič, 2011; Stanek in Turk, 1998). Pri tem imamo dve možnosti (Chrzanowski in sod., 1982): • V izračun vključimo opazovanja. S poljem premikov tako podamo razliko med opazovanji v prvi in drugih terminskih izmerah (za vsako naslednjo terminsko izmero dobimo svoje polje premikov) za vsako točko obravnavanega objekta oziroma območja. Ta način poimenujmo pristop z opazovanji. • Opazovanja predhodno izravnamo in izračunamo koordinate točk. S poljem premikov je v tem primeru podana razlika med položajem točke v prvi in drugih terminskih izmerah. Ta način poimenujmo koordinatni pristop. vi r-C-"i s o Prednosti pristopa z opazovanji sta, da se s tvorjenjem razlik med opazovanji eliminirajo morebitni sistematični pogreški ter da so razlike med opazovanji v različnih terminskih izmerah neodvisne od geodetskega datuma. Prednost koordinatnega pristopa je, da je formulacija deformacijskega modela enostavna in pregledna. Grafična predstavitev razlik koordinat je bolj nazorna kakor grafična predstavitev razlik opazovanj, kar bistveno pripomore k boljši izbiri mogočih deformacijskih modelov. Druga pomembna prednost koordinatnega pristopa je, da omogoča odkrivanje in izločitev grobo pogrešenih opazovanj pred obravnavo deformacijskega modela. Potek postopka Fredericton lahko razdelimo v pet korakov, ki so v nadaljevanju tudi podrobneje opisani (Secord, 1985; primerjaj tudi Savšek-Safic, 2002): 1. izravnava opazovanj v geodetski mreži po metodi najmanjših kvadratov (MNK) ter izločanje grobo pogrešenih opazovanj, ta korak se izvede, samo kadar uporabimo koordinatni pristop; 2. določitev mogočih deformacijskih modelov; 3. določitev deformacijskih parametrov za izbrane deformacijske modele; 4. pregled deformacijskih modelov in izbira najboljšega; 5. grafični prikaz izbranega modela. 2 TEORETIČNO OZADJE 2.1 Izravnava opazovanj v geodetski mreži po metodi najmanjših kvadratov in izločanje grobo pogrešenih opazovanj Najprej vsako terminsko izmero posebej izravnamo po metodi najmanjših kvadratov (MNK), tako da minimiziramo uteženo vsoto kvadratov popravkov opazovanj. Postopek je podrobno opisan v Grigillo in Stopar (2003). V nadaljevanju opisovanja postopka Fredericton uporabimo spodaj naštete količine iz izravnave posamezne terminske izmere: • vektor izravnanih koord in at x = x + x0 (A + A0), • pripadajočo matriko kofaktorjev koordinatnih nezn ank Q^ (QAA), • referenčno varianco a posteriori og2 (cr02), • pripadajoče prostostne stopnje f (r), • matriko normalnih enačb N (N), • datumsko matriko H (H). Oznake, ki jih uporabljamo, so identične oznakam, uporabljenim v Chen (1983), v oklepajih pa podajamo pripadajoče oznake iz članka Grigillo in Stopar (2003). 2.2 Določitev mogočih deformacijskih modelov Cilj tega koraka je določiti/identificirati deformacijske modele, s katerimi bi lahko opisali premike v geodetski mreži in torej na obravnavanem objektu. Deformacijski model je povezava med opazovanimi količinami in deformacijskimi parametri (Chen, 1983). Mogoče deformacijske modele izberemo na podlagi predhodnih podatkov in/ali izkušenj o dogajanju na območju. V primeru premikov površja so uporabne informacije o geološki sestavi tal, prelomih in premikih tektonskih plošč itd. Pri umetnih objektih se pri določanju mogočih deformacijskih modelov lahko opremo na informacije o lastnostih materiala, izpostavljenosti naravnim pojavom (npr. vetru, soncu, mrazu) ter delovanju drugih mehanskih sil na objekt (npr. posedanje terena v okolici objekta). Pomagamo si tudi z grafično predstavitvijo vektorja premikov d med terminskima izmerama, ki je lahko obremenjen z defektom geodetskega datuma. Ta vpliv reduciramo z iterativnim dodeljevanjem uteži vplivom posamezne točke na določitev datuma, tako da imajo točke z večjimi premiki manjši vpliv. Osnovna ideja iterativnega dodeljevanja uteži IWP (angl. iterative weighted projection) je vektor d transformirati v datum, ki ga bolj določajo (imajo večjo utež) točke z manjšo koordinatno razliko med dvema terminskima izmerama. Teoretično s tem postopkom določimo datum, ki minimizira | ^ prvo normo končne projekcije vektorja d (Chen, 1983). Postopek je izpeljanka transformacije | g S, kjer sta v matriki uteži mogoči le vrednosti 0 ali 1 (Marjetič in Stopar, 2007). To je iterativni | ^ II to postopek, začetne vrednosti priredimo kot: || „g d0 = X - x„ W0 = I (1) I I vi S kjer sta: d0 - vektor razlik izravnanih koordinat iz dveh terminskih izmer, W0 - matrika uteži, ki določi, kolikšen naj bo vpliv posamezne točke pri določitvi geodetskega datuma. Število v indeksu označuje, za kateri korak iteracije gre. V/tem koraku iteracije izračunamo: d. = (I - H(HTW._1H)-1HTW._1)d-1. (2) = di3g fe) (3) Na podlagi matrike W, določene v zadnjem koraku, transformiramo tudi matriko kofaktorjev, ki jo potrebujemo za izris elips: Q = (I - H(HTWH)_1HTW)Q..(I - H(HTWH)-1HTW)T (4) vi r-C-"i S o to kjer je: Qdd = (Qx x + ) - matrikako faktorjev razlik izravnanih koordinat iz dveh terminskih izmer. Iteracije končamo, ko je razlika med vektorjema d.. in d._; dovolj majhna. V praksi izračun izvedemo v od štirih do šestih korakih. Pri tem se lahko zgodi, da bodo nekatere komponente v vektorju premikov zelo blizu 0, in torej pripadajoče uteži v matriki W zelo velike, kar bi povzročilo numerično nestabilnost. Zato premikom, manjšim od neke vrednosti (npr. |d, |<10-6m = 8), dodelimo utež 0 v matriki W (Chen, 1983). Problem bi lahko rešili tudi tako, da bi uteži takih 1 premikov izračunali kot Wj = diag dj\+5 (Chen in sod., 1990). Polje premikov nato smiselno razdelimo na bloke (podobmočja), tako da lahko znotraj posameznega bloka predpostavimo homogene deformacije. Polje premikov modeliramo z aproksimacijskim polinomom, kjer upoštevamo le linearne člene (Secord, 1985): dy = a0 + a1y + a2x, dx = b0 + b1x + b2x. (5) (6) (7) (8) Koeficienti v zgornjih enačbah imajo fizikalni pomen. Enačbi lahko pretvorimo v zapis: dy = t + s y + (s + co)x, s y yy yx y dx = t + s x + (s - cc)y, x xx yx ^ kjer so: ty, tx - komponenti translacije togega telesa (bloka) v smereh koordinatnih osi, syy, sx - normalni deformaciji togega telesa (bloka), syx - strižna deformacija togega telesa (bloka) in ca - parameter rotacije togega telesa (bloka). Enačbi (7) in (8) uporabimo za opise premikov posamezne točke v izbranih deformacijskih modelih. Nekatere od šestih deformacijskih parametrov lahko izpustimo iz nekaterih deformacijskih modelov. Za vse točke v posameznem bloku zapišemo enačbi (7) in (8) v matrični obliki: d = Be, dy 1 y x 0 0 x dx 0 0 y 1 x -y (9) (10) "yy £yx tX w J kjer sta: B - matrika deformacijskega modela, ki povezuje deformacijske parametre z razlikami koordinat iz dveh terminskih izmer, e - vektor deformacijskih parametrov. 2.3 Določitev deformacijskih parametrov za izbrane deformacijske modele 2.3.1 Izračun za več terminskih izmer Matematični model sloni na hipotezi, da je pričakovana vrednost opazovanj (koordinat) v /-ti izmeri enaka seštevku pričakovane vrednosti opazovanj (koordinat) v prvi terminski izmeri in premika, ki ga opisuje deformacijski model. Ničelno (H0) in alternativno (Ht) hipotezo matematično formiramo (Chen, 1983): - pri pristopu z opazovanji: H0 :E [l,.] = E [l, ] + A,B,.e H,:E [l, ]* E [l, ]+ A,. B,.e - pri koordinatnem pristopu: H0:E[l, ]= A,. E[x, ] = A,. E[ x, ] + A,. B, e H0:E[x,. ]= E[x, ] + B,.e H,:E[x, ]* E[x, ] + B,.e kjer so: E [ ] - pričakovana vrednost, l - vektor opazovanj, A - matrika koeficientov enačb, ki povezujejo koordinate z opazovanji. (11) (12) (13) (14) ■is S vi S j5 to * £ o to Oba pristopa se razlikujeta le v matriki A, ki je pri koordinatnem pristopu enotska in jo lahko izpustimo. Zato je v nadaljevanju podan enak potek postopka za oba primera. Na podlagi ničelne hipoteze (11) ali (13) lahko tvorimo matematični model: X + Sj = % , x + 8i = % + A Be (15) (16) kjer so: x - vektor opazovanj (pristop z opazovanji) ali ocenjenih koordinat (koordinatni pristop) posamezne izmere, 5 - vektor popravkov opazovanj (pristop z opazovanji) ali popravkov koordinat (koordinatni pristop), ki se določijo pri ocenjevanju deformacijskih parametrov, 4 - vektor konstant. Cilj tega koraka je torej določiti vrednosti v vektorju e. To je mogoče, kadar je velikost (dimenzija) vektorja x enaka v vseh terminskih izmerah. To pomeni, da morajo biti v vseh terminskih izmerah opazovana ista opazovanja in v izmero vključene iste točke. V praksi se to ne zgodi vedno. Ker je ena izmed zahtev za postopek, da ga je mogoče uporabiti tudi ob morebitnem defektu konfiguracije mreže, v takem primeru vektor x ustrezno dopolnimo s približno vrednostjo manjkajočega opazovanja (koordinate), v matriki uteži pa na ustrezna mesta postavimo ničle ter tako eliminiramo vpliv tega opazovanja (koordinate) na rezultate. Sistem (15) in (16) je rešljiv, če je premikov točk v mreži (število enačb) vsaj toliko, kolikor je deformacijskih parametrov (število neznank). Če je premikov točk več kot deformacijskih parametrov, slednje ocenimo z metodo najmanjših kvadratov. Dobimo sistem normalnih enačb (izpeljava je zaradi preglednosti v prilogi A): rt I- ¿=i k TA;TP; rt ¿=2 k TAjTPjAjBj 0- rt ¿=1 k IB, TA¿TP¿X¿ (17) d r-C- s TS o kjer je: P. - matrika uteži /-te terminske izmere. i Za koordinatni pristop to matriko določimo kot: P. = N. = Qx x 1 za pristop z opazovanji kot: P = P,,, (18) kjer je: Pu - matrika uteži opazovanj. Število normalnih enačb je enako številu neznanih parametrov in sistem (17) je enolično rešljiv. Iz enačb eliminiramo 4 tako, da ga izrazimo iz prve enačbe in vstavimo v drugo, nato izrazimo neznan vektor e : i=1 r fc e = £ B^V^Bi - £ Bi%TPi l£ Pi I £ PiAiBi i=2 i=2 Vi=l / i=2 fc fc / fc x _ fc £BiTAiTPiXi-£BiTAiTPi(£Pi) £PiXi i=l A . (20) (21) Pripadajočo kovariančno matriko vektorja e izračunamo kot Zee = ôl Q, ee uO Vee k ■•ee — ôo kjer je: iv rv / \ £ B^A^AA - £ B^Vp l£ I £ P^B; -i=i i=i \i=i / i=i (22) q2 _ - referenčna varianca a posterioriiz izravnave k-izmer. Pri koordinatnem pristopu pomeni težavo še singularnost matrike P;, saj v nadaljnjem izračunu potrebujemo inverzijo matrike Pr Uporabimo tako imenovano generalizirano inverzijo (Ek P()-. Generalizirana inverzija ni enolično določena, vendar so deformacijski parametri neodvisni od načina določitve generalizirane inverzije. Če je defekt ranga matrike 2kwP. posledica le defekta geodetskega datuma (in ne tudi defekta konfiguracije mreže), je izračun generalizirane inverzije preprost: i tI (23) k /k . + hh \i=i / \i=i / 2.3.1 Izračun za dve terminski izmeri » Izračun deformacijskih parametrov se poenostavi, kadar primerjamo le dve terminski izmeri. Odštejemo enačbi v matematičnem modelu (15) in (16) in dobimo: d + vd = ABe kjer je: d = x2 - Xj. Kovariančno matriko vektorja premikov izračunamo kot: (24) § a I 5 S „s 1 ^ o Edd=<^Qdd= (Qx,i,+Q i, i,)- v! r-C-"i s o to (25) Deformatijrko earametrn določimo po poenostavljeni enačbi (poenostavitevje narcdi pre^ednosti izjp eljana v prilogi B): e = BTAT ( + P^A^B^BTA(p + P2)d , = 60[bJAJ (P! + P2)A2E12 =O02Qee , , kjer so: - kovariančna matrika deformacijskih parametrov, Qee - matrika kofaktorjev deformacijskih parametrov, (26) (27) Pi +P, paralelna vsota matrik uteži iz posamezne terminske izmere (definirana je v izpeljavi v prilogi B). Za pristop z opazovanji to matriko določimo kot: P, ± P2 = QUi-1 ± Q„2-1 =(QUi + Q„J za koordinatni pristop kot: p, ± P2 = NJCN, + N2>- N2 (28) (29) Deformacijske parametre smo torej pri koordinatnem pristopu izračunali iz premikov, izračunanih neposredno iz izravnanih koordinat. Čeprav so ti premiki odvisni od določitve geodetskega datuma, lahko ocenjene deformacijske parametre, če ne pride do spremembe orientacije mreže (pojasnjeno na koncu poglavja), obravnavamo kot neodvisne od določitve geodetskega datuma. Vektor premikov d, ki je izračunan iz izravnanih koordinat, lahko namreč zapišemo kot linearno kombinacijo vektorja dh, ki je odvisen od zagotovitve datuma v izravnavi in »pravega« vektorja premika dp(Chen, 1983): d = d + d,. p h Vektor dh je linearna kombinacija stolpcev matrike H. Stolpci matrike H predstavljajo bazo ničelnega prostora matrike N. Vektor dh torej leži v ničelnem prostoru matrike N. Pri določevanju deformacijskih parametrov, tj. določanju vrednosti elementov vektorja e, vektor premikov množimo z matriko P, ± P2 = N^N, + N2)- N2. Ničelni prostor matrike P, ± P2 je enak kot ničelni prostor matrike N. Za vektor dh torej velja, da leži tudi v ničelnem prostoru matrike P, ± P2 in se pri preslikavanju s to matriko preslika v 0: (P, ± P2)dh = 0. (31) To lastnost uporabimo v enačbi za določitev deformacijskih parametrov: e = (BT(Pi ± P2)B)" BT(Pi ± P2)d = = (Bt(Pi ± P2)B)"1Bt(P1 ± P2) (dp + dh) = = (BT(Pi ± P2)B)"1BT(P1 ±P2)dp + + (Bt(Pi ± P2)B)_1Bt (Pi ± P2)dh = (33) (34) = (BT(Pi ± P2)B)_1Bt(P1 ± P2)dp = e (35) Vidimo torej, da je rezultat (izračunani deformacijski parametri) enak ne glede na to, kako je določen geodetski datum. Pri vpeti mreži je pri odpravi vpliva določitve datuma pomembno, da je parameter rotacije m vključen v deformacijski model. Če se namreč med terminskima izmerama premakne katera od točk, med katerima je definirana smer, se zasuka celotna mreža. V deformacijskem modelu prevzame parameter m nase ta vpliv in tako ostanejo drugi deformacijski parametri neodvisni od določitve geodetskega datuma. Element rotacije v model vključimo tako, da ga dodamo v vektor neznank e ter v matriko B dodamo ustrezne odvode po neznankah (Chen, 1983). 2.4 Pregled deformacijskih modelov in izbira najboljšega V fazi identifikacije kot mogoče podamo več modelov. Nato za vsakega določimo deformacijske parametre. V tej fazi je treba modele pregledati in izbrati najboljšega. To naredimo na podlagi statističnih testov. Računamo globalni test modela in stopnje zaupanja za posamezni deformacijski parameter. Globalno primernost modela testiramo z drugo normo popravkov iz matematičnega modela. Modela ne moremo zavrniti, če velja (Chen, 1983): r2_£L_vk ItiSiP;«; < F (a, dfe, df) (36) (T2 (T2 dfe kjer so: df = £k=f - število prostostnih stopenj pri izračunu referenčne variance a posteriori v vseh terminskih izmerah, f - število prostostnih stopenj za izračun a posteriori variance v posamezni terminski izmeri, dfe = rang (P) - ne = rang (P) - u + d - število prostostnih stopenj pri izračunu referenčne variance a posteriori pri določevanju deformacijskih parametrov, P = diag (Pj, ..., P;) - matrika uteži, ne - število deformacijskih parametrov v vektorju e, u - število neznank v matematičnem modelu (17), ■is S S- vi S j5 * o ■ to d - skupni defekt datuma iz vseh terminskih izmer1, F(a, dfe, df) - kritična vrednost po Fisherjevi porazdelitvi, 1 - a - stopnja zaupanja testa. Stopnjo zaupanja za posamezni deformacijski parameter i določimo tako, da izračunamo njegovo standardizirano formo ter poiščemo tisto vrednost stopnje zaupanja 1 - a, pri kateri je izpolnjen pogoj: e2 t2 >F(a,\, df), (37) kjer sta: e i - element vektorja e, izračunanega v enačbi (26), ki pripada obravnavanemu deformacijskemu parametru i, qe - kofaktor obravnavanega deformacijskega parametra, tj. i-ti diagonalni element matrike Qe, izračunane v enačbi (27). Najboljši model izberemo glede na naslednja merila: - manjše število deformacijskih parametrov, - preprostejši deformacijski model, - razlogi, ki imajo geofizikalno ozadje, - minimalni vsota kvadratov popravkov, - stopnja zaupanja parametrov. 2.5 Grafični prikaz izbranega deformacijskega modela Izbrani deformacijski model prikažemo še grafično. Prikažemo razporeditev točk v bloke in določene premike vsakega od blokov. 3 RAČUNSKI PRIMER Postopek Fredericton smo uporabili na podatkih iz literature (Mihailovic in Aleksic, 1994). Na istih podatkih so že bili preizkušeni postopki Hannover (Ambrožič, 2001), Karlsruhe (Ambrožič, 2004) in Delft (Marjetič in sod., 2012). Ponovitev izračuna je mogoča na podlagi vhodnih podatkov, ki so zapisani v Ambrožič (2001). Simulirani premiki so zaradi lažje primerjave z rezultati deformacijske analize po različnih postopkih zapisani v preglednici 3. Rezultate postopka Fredericton podajamo po zgoraj navedenih petih korakih. 1 ^s S 2 I to H0 Deformacijski Vrednosti Stopnja Globalni test df parametri parametrov [mm] zaupanja [%] [a = 0.05] dfe Model 0: 91.43 > 1.80 60 ni premikov H0 zavrnjena 17 dy1 - 19.6 100 Model 1: dxi dy2 - 38.0 - 38.7 100 100 premiki točk: dx2 49.0 100 0.12 < 2.04 60 1, 2, 3 in 7 dy3 dx3 dy7 dxl 20.6 - 44.3 23.6 42.9 100 100 100 100 H0 ni zavrnjena 9 dy1 - 18.0 99.9 Model 2: dx1 dy2 - 37.3 - 36.6 100 100 premiki točk: dx2 51.3 100 1, 2, 3, 4 in 7 dy3 22.1 100 0.09 < 2.17 60 dx3 - 41.0 100 H0 ni zavrnjena 7 dy4 1.1 22.6 dx4 3.2 41.2 dy7 24.6 100 dx7 44.4 100 d/2 - 21.7 100 Model 3: dx2 d/i 87.5 27.7 100 100 premiki točk: dx3 - 4.2 48.8 11.32 > 2.04 60 2, 3, 4 in 7 d/. - 0.6 12.3 H0 zavrnjena 9 dx4 30.1 100 d/7 30.4 100 dxi 67.7 100 vi r-C- S o to Preglednica 1: Vrednosti deformacijskih parametrov, njihove stopnje zaupanja in vrednosti testne statistike za globalni test modela za štiri izbrane deformacijske modele, po koordinatnem pristopu. Za predlagana deformacijska modela 0 in 3 je vrednost testne statistike presegla kritično vrednost. Zato lahko v teh dveh primerih zavrnemo hipotezo H0, ki pravi, da lahko polje premikov pojasnimo s predlaganim deformacijskim modelom. »Najboljši« deformacijski model smo zato izbirali med modeloma 1 in 2, ki se razlikujeta v tem, da slednji kot premaknjeno obravnava tudi točko 4. Glede na nizki stopnji zaupanja za premika te točke v smereh koordinatnih osi smo kot »najboljšega«izbrali model 1, torej premike točk 1, 2, 3 in 7. Za primerjavo rezultatov smo za izbran deformacijski model določili vrednosti parametrov tudi s pristopom z opazovanji. Rezultati so podobni kot pri koordinatnem pristopu, premiki se razlikujejo za največ 0,1 mm, nekoliko manjši sta tudi stopnji zaupanja v komponento y premika točk 1 in 3. Rezultati so podani v preglednici 2. Ho Deformacijski parametri Vrednosti parametrov [mm] Stopnja zaupanja [%] Globalni test [a = 0.05] dy1 - 19.7 97.4 Model 1: dx1 dy2 - 38.1 - 38.8 100 100 premiki točk: dx2 49.0 100 5.79 < 12.59 1, 2, 3 in 7 dy3 20.6 99.6 H0 ni zavrnjena dx3 dy7 - 44.4 23.6 100 100 dx7 42.9 100 Preglednica 2: Ocenjene vrednosti parametrov in njihove stopnje zaupanja ter globalni test modela po pristopu z opazovanji. 4 PRIMERJAVA Z REZULTATI DRUGIH METOD Obravnavan računski primer smo primerjali tudi z drugimi postopki deformacijske analize. Rezultati so zbrani v preglednici 3, ki je povzeta po Marjetič in sod. (2012). Dodani so rezultati postopka Fredericton. Iz Preglednice 3 je razvidno, da so premiki, določeni po vseh postopkih, zelo podobni. Razlike med njimi so reda velikosti nekaj desetink milimetra za posamezno koordinatno komponento, medtem ko od pravih (simuliranih vrednosti) odstopajo za približno 5 milimetrov, največje pa je odstopanje komponente y premika na točki 2, ki v absolutnem smislu znaša 8,8 milimetra. vi s „s ■is s S- o ■ to Točka 1 2 3 4 5 6 7 o dy r [mm] - 20.0 - 30.0 25.0 0.0 0.0 0.0 25.0 is dx [mm] - 34.6 52.0 - 43.3 0.0 0.0 0.0 43.3 6 So d [mm] 40.0 60.0 50.0 0.0 0.0 0.0 50.0 F n 210 330 150 - - - 30 d y y [mm] - 19.4 - 38.1 21.4 0.7 - 0.8 0.0 24.0 dx [mm] - 37.5 49.5 - 43.5 1.0 - 2.3 1.3 42.9 a d [mm] 42.2 62.5 48.5 1.2 2.4 1.3 49.2 F [°] 207 322 154 35 199 0 29 Premik da da da ne ne ne da d y y [mm] - 19.7 - 38.8 20.6 - - - 23.6 ü dx [mm] - 38.0 49.0 - 44.4 - - - 42.9 Ti d [mm] 42.8 62.5 48.9 - - - 49.0 F [°] 207 322 155 - - - 29 Premik da da da ne ne ne da dy [mm] - 19.6 - 38.7 20.6 - 4.0 - 6.4 3.3 23.6 £ o dx [mm] - 38.0 49.0 - 44.3 5.1 - 7.1 - 10.6 42.9 Ö Ö cS d [mm] 42.8 62.4 48.9 6.5 10.0 11.1 49.0 X F [°] 207 322 155 322 222 163 29 Premik da da da ne ne ne da dy r [mm] - 19.6 - 38.7 20.6 - - - 23.6 ö o dx [mm] - 38.0 49.0 - 44.3 - - - 42.9 •o o d [mm] 42.8 62.5 48.9 - - - 48.9 £ F [°] 207 322 155 - - - 29 Premik da da da ne ne ne da Preglednica 3: Rezultati deformacijske analize po postopkih Delft, Karlsruhe, Hannover in Fredericton. vi S TS o Sü to 5 SKLEP Prednost postopka Fredericton pred drugimi vidimo v tem, da poskuša premike smiselno povezati v sklenjene bloke in njihove deformacije utemeljiti z geofizikalnim dogajanjem na obravnavanem območju ali objektu. Postopka Fredericton ni mogoče popolnoma avtomatizirati, saj je v fazi identificiranja mogočih deformacijskih modelov potrebna strokovnjakova odločitev, katere točke je smiselno povezati v bloke in katere deformacijske parametre je, glede na polje premikov, smiselno določevati za posamezen blok. Pri tem je treba vedno obravnavati polje premikov, ki smo ga predhodno transformirali z iterativnim dodeljevanjem uteži (poglavje 2.2). Korak identifikacije mogočih deformacijskih modelov je pomemben, saj uporabnih rezultatov ne moremo dobiti, dokler ne identificiramo dovolj dobrega modela (v nasprotnem primeru bo globalni test modela zavrnjen). Zanimivo pa bi bilo preskusiti, ali bi lahko ta korak zadovoljivo avtomatizirali z uporabo katerega od modernejših postopkov gručenja vektorjev premikov, na primer z nevronskimi mrežami. Rezultati štirih postopkov, ki smo jih primerjali v poglavju 4, se le nekoliko razlikujejo med seboj. Pri tem je treba poudariti, da gre za test na simuliranem primeru, kjer so bili simulirani premiki glede na natančnosti določitve položajev točk veliki. Predvidevamo, da bi večje razlike med rezultati postopkov nastale pri obravnavi manjših premikov. Za boljše razumevanje vpliva posebnosti posameznega postopka na rezultate bi bilo nujno izvesti primerjavo rezultatov še za več različnih geodetskih mrež v več terminskih izmerah. Pri tem bi bilo smiselno opazovati, kako na različnost rezultatov med metodami vplivajo geometrija mrež, defekti konfiguracij med terminskimi izmerami, polje premikov in število nadštevilnih opazovanj. Taka primerjava bi lahko bila osnova za odločitev, katerega od postopkov uporabiti v konkretni situaciji. Literatura in viri: Ambrožič, T. (2001). Deformacijska analiza, po postopku Hannover. Geodetski vestnik, 45(1-2), 38-53. Ambrožič, T. (2004). Deformacijska analiza, po postopku Karlsruhe. Geodetski vestnik, 48(3), 315-331. Chen, Y. Q. (1983). Analysis of Deformation Surveys - A Generalized Approach. Doktorska, disertacija. Fredericton: University of New Brunswick, Department of Geodesy and Geomatics Engineering; pridobljeno 2. 6.2010 s spletne strani: http://gge.unb.ca/Pubs/TR94. Chen, Y. Q., Chrzanowski, A., Secord J. M. (1990). A Strategy for the Analysis of the Stability of Reference Points in Deformation Surveys. CISM Journal ACSGC, 44(2), 141-149. Chrzanowski, A., Chen, Y. Q., Secord J. M. (1982). A Generalized Approach to the Geometrical Analysis of Deformation Surveyes. 3rd international symposium on deformation measurments by geodetic methods, Budimpešta, 25.-27. avgust 1982. Grigillo, D., Stopar, B. (2003). Metode odkrivanja grobih pogreškov v geodetskih opazovanjih. Geodetski vestnik, 47(4), 387-403. Marjetič, A., Stopar, B. (2007). Geodetski datum in S-transformacija. Geodetski vestnik, 51(3), 549-564. Marjetič, A., Zemljak, M., Ambrožič, T. (2012). Deformacijska analiza, po postopku Delft Geodetski vestnik, 56(1), 9-26. Marjetič, A. (2011). Statistična analiza značilnih premikov v geodetskih mrežah. Doktorska disertacija. Ljubljana: § Fakulteta za gradbeništvo in geodezijo. | I Mihailovic, K, Aleksic, I. (1994). Deformaciona analiza geodetskih mreža. Beograd: Univerzitet u Beogradu. ^ Graifevinski fakultet, Institut za geodeziju. % , § -v Mozetič, B., Kogoj, D., Ambrožič, T. (2006). Uporabnost izbranih metod deformacijske analize na praktičnih primerih ^ geodetskih mrež. Geodetski vestnik, 50(4), 620-631. | ^ Rao, C. R., Mitra, S. K. (1971). Generalized Inverse of Matrices and its Applications. New York, London, Sydney, g Cn Toronto: John Wileyand Sons. Savšek -Safic, S. (2002). Optimalna metoda, določevanja, stabilnih točk v deformacijski analizi. Doktorska, disertacija. Ljubljana: Fakulteta za gradbeništvo in geodezijo. Secord, M. J. (1985). Implementation of a. Generalized Method for the Analysis of Deformation Surveys. Magistrska naloga.. Fredericton: University of New Brunswick, Department of Geodesy and Geomatics Engineering; pridobljeno 2. 6.2010 s spletne strani: http://gge.unb.ca/Pubs/TR117. f Stanek, M., Turk, G. (1998). Osnove mehanike trdnih teles. Ljubljana: Fakulteta za gradbeništvo in geodezijo. ^ to 3 t-O I I ris ^ C PRILOGA A: 3 v! I ^ H "i | •V I S ^ Š 1 1 •iS to ■S Izpeljava modela s pogojem najmanjših kvadratov popravkov Imamo matematični model: X + 6i = % x + 8( = % +A,.B,.e, i = 2, ..., k Poiskati moramo tako rešitev, da bo utežena vsota kvadratov popravkov minimalna. Ta pogoj zapišemo kot: k $ = STPi s1 + 6,- F,-6,- = min i=2 = s:[p181 + ^8[pi8i =mi i=2 k = G - x1)TP1(? - x1) + ^(š + AiBie - xi)TPi(? + A^e - Xi) = min Člene zmnožimo in seštejemo istovrstne člene. Pri tem upoštevamo, da so matrike Pi diagonalne in velja Pi T = Pi ter pravilo za množenje transponiranih matrik (AB)T = BTAT. $ = - ^TP1x1 - x1TP1^ + x1TP1x1 + k + + ?TPiAiBie - ?TPixi + e TBiTAiTPi^ + e TBiTAiTPiAiBie - i=2 - e TBiTAiTPixi -xiTPi^-xiTPiAiBie + xiTPixi) = min. $ = ^P^- 2 ?1P1x1 + x11P1x1 + k + ;TPi; + 2 ;TPiAiBie - 2 ;TPixi + e TBiTAiTPiAiBie - 2 e TBiTAiTPixi + i=2 + xiTPixi) = min. £ A Ker moramo poiskati minimum, enačbo odvajamo po obeh neznankah, po s in po e, ter oba odvoda enačimo z 0: k — = 2?^- 2P1x1+^(2Pi? + 2PiAiBie -2PiXi) = 0 dt d® = 2( ?TPiAiBi)T + 2BjTAjTPjAjBje - 2 BiTAiTPixi) = 0 ■ * i s to t Obe enačbi delimo z 2 in seštejemo istovrstne člene: k k k ^(Pi? ) + ^(PiAiBie ) -^(PiXi) = 0 i=l i=2 i=l k k k ^(BiTA[Pi?) + ^(BiTA7PiAiBie ) - ^(B/A^P^) = 0 1=2 1=2 1=2 Enačbi zapišemo v matrični obliki in izpeljava je končana: Li=2 k k y p^B; t i ¿ = 1 i=2 k BiTAiTPi V BiTAiTPiAiB ( i ¿=2 1 = 1 k ^BiTAiTPixi vi S is ■is s i o ■ to PRILOGA B: Izpeljava enačbe za izračun deformacijskih parametrov za dve terminski izmeri Izpeljava je bila narejena na podlagi poglavij o paralelni vsoti matrik avtorjev Rao in Mitra (1971). V izpeljavi bomo uporabili naslednje zveze: A ± B M A(A + B)"B = A(A + B)"B+A(A + B)"A-A(A + B)-A = = A (A + B)"(A + B) -A(A + B)"A = A -A(A + B)"A vi r-C- S •v o to Enačba za izračun deformacijskih parametrov iz poljubne št. terminskih izmer je (enačba 21): r fc e = ^BiTPiBi-^BiTAiTPii^Pij ^AA ¿=2 ¿=2 M=l / ¿=2 fc fc / fc \ " fc ¿=l ¿=l Zaradi preglednosti izpeljave obravnavamo matriko A kot enotsko oziroma kot del matrike B. Enačbo zapišemo za i = 1, 2: e = [B2TP2B2 -B2tP2(P1 + P2)"P2B2]_1[B2tP2x2 - B2TP2(P1 + P2)"(P1x1 + P2x2)] Vidimo, da so vsi členi množeni z B2T, v prvem oklepaju pa tudi z B2. Izpostavimo: e = B2T(P2-P2(PI + P2)-P2)B2 Pl ±p2 [B2t(P2x2 - P2(PX + P2)"(P1x1 + P2x2))] e= [B2T(P1 ± P2)B2]" [B2T(P2x2- P2(P1 + P2)"P1x1-P2(P1 + P2)-P2x2)] e= [B2T(P1 ± P2)B2]" B2T( (P2 - P2(Pi + P2)-P2)X2 - P2(Pi + P^)-PiXi ) e = [B2t(P1 ± P2)B2]"1B2t(P1 ± P2) (x2 - Xl) e = [B2T(P1 ± P2)B2f B2T(P1 ± P2)d LLI ^ O oc H Z P O z < z N N Prispelo v objavo: 14. marec 2013 oŠ Sprejeto: 20. junij 2013 ^ S asist. Anja Vrečko, univ. dipl. inž. geod. ^ FGG - Oddelek za geodezijo, Jamova 2, SI-1000 Ljubljana 'S e-pošta: anja.vrecko@fgg.uni-lj.si izr. prof. dr. Tomaž Ambrožič, univ. dipl. inž. geod., univ dipl. inž. rud. ■I FGG - Oddelek za geodezijo, Jamova 2, SI-1000 Ljubljana e-pošta: tomaz.ambrozic@fgg.uni-lj.si