PROBLEMES INVERSES POUR L'ENVIRONNEMENT

Description

PROBLEMES INVERSES POUR L'ENVIRONNEMENT 15 5. Methodes variationnelles : controle optimal, etat adjoint Nous allons etudier dans cette partie le 4D-VAR, qui est une generalisation du 3D-VAR pour des observations distribuees dans le temps. On tient donc compte desormais d'un modele d'evolution (dans le temps) qui permettra de comparer l'etat du systeme avec les observations a l'instant approprie. La fenetre d'assimilation est un intervalle de temps donne, l'analyse est realisee a l'instant initial, et on suppose les observations distribuees sur m instants (ti)0≤i≤m dans l'intervalle. On notera y(ti) les observations, x(ti) l'etat du systeme, et xt(ti) l'etat vrai du systeme, a l'instant ti. La matrice de covariance des erreurs d'obser- vations a l'instant ti est notee Ri. L'operateur d'observation correspondant est note Hi. Par contre, la matrice de covariance d'erreur d'ebauche est toujours notee B, car elle n'est definie qu'a l'instant initial, l'ebauche xb etant une estimation a priori de l'analyse, donc a l'instant initial. 5.1. 4D-VAR : assimilation variationnelle en dimension 4. Dans sa formu- lation generale, le principe du 4D-VAR est de considerer la minimisation de la fonction cout suivante : (12) J(x0) = (x0?xb) TB?1(x0?xb)+ m∑ i=0 (yi ?Hi(x(ti))) T R?1i (yi ?Hi

• matrice de covariance des erreurs d'obser- vations


• condition initiale


• instant initial


• modele adjoint


• fac¸on retrograde en temps


• developpement de taylor des operateurs d'observation


• adjoint du modele lineaire tangent

Sujets

• Condition initiale

Informations

• Publié par: pefav
• Nombre de lectures: 25
• Langue: Français

Extrait

` PROBLEMES INVERSES POUR L’ENVIRONNEMENT

15

5.ontres:copti ˆ olee ´ atam,liotnatjdM ´ hoetvsedairanoitllen Nousallons´etudierdanscettepartiele4D-VAR,quiestuneg´ene´ralisationdu 3D-VARpourdesobservationsdistribue´esdansletemps.Ontientdonccompte de´sormaisd’unmod`eled’´evolution(dansletemps)quipermettradecomparer l’´etatdusyste`meaveclesobservationsa`l’instantappropri´e. Lafenˆetred’assimilationestunintervalledetempsdonne´,l’analyseestre´alis´ee`a l’instantinitial,etonsupposelesobservationsdistribu´eessurminstants (t i ) 0≤i≤m dans l’intervalle. On noteray(t i ) les observations,x(t i me,ettate´’l)e`tsysudx t (t i ) l’´etatvraidusyste`me,a`l’instantt i . La matrice de covariance des erreurs d’obser- vations`al’instantt i tontsee´eR i e´totera´eopL’.vrtaoicnru’dboesdantestnorrespon H i ot´eursneehcuabe´ojuottsed’ceanrid’urreerocraertnP.decevacoam,lriatB, carellen’estde´finiequ’a`l’instantinitial,l’e´bauchex b oianrpoiir´etantuneestimat del’analyse,donca`l’instantinitial. 5.1.4D-VAR : assimilation variationnelle en dimension 4.Dans sa formu- lationge´n´erale,leprincipedu4D-VARestdeconsid´ererlaminimisationdela fonctioncoˆutsuivante: m X T−1 T −1 (12)J(x() = t)))R( (x(t))). 0 (x 0 −x b )B(x 0 −x b () + y i −H i (x i i y i −H i i i=0 Laminimisationdecettefonctioncoˆutestre´alise´esouscontrainte.Ils’agitd’une contraintefortedemod`ele,puisquel’´ecrituredeJ(xd´)leesdvadsenuerpx(t i ), qui elles-meˆmesde´pendentdelaconditioninitialex 0 . D’un point de vuecontinusoppustuepno,noilotu´’veuqdeanimledyod`eelemerqu s’e´critdelafac¸onsuivante: dx (13) =F(x), x(0) =x 0 . dt Maisonpeut´egalemente´crireleprobl`emedefa¸concs`ridtee, en supposant que la suitedes´etatsdumod`elex(t i anivsuon¸cfaladercti’se´):et (14)x(t) =M(x), i0→t i 0 o`uM 0→t i r´laoleseld`oue(temrtnattnavep)eearetruomts’lpoe´ssreedapnitsed’lant initial`al’instantt i . Le4D-VARestdoncunprobl`emecomplique´deminimisationnonline´airesous contrainte.Onpeutlesimplifier`al’aidedesdeuxhypothe`sessuivantes: •ausaC:eil´tledeom`delod`elessuitedemmmocenuece´’erirnpiotseuev’´utol interm´ediaires,permettantdepasserd’uninstantausuivant.Onpeutsupposer queM 0 e),etsil’idesttsna0ta`ul-iˆmmedentsspadeerinl’itne(e´tmrepatte on noteM i `dlepereemttnadtlar´esolvantedumotnaseprdse’ieltanst i−1 a` t i , alorsx i =M i x i−1 e´ucaprrec,rrnedonc,et x i =M i M i−1 . . . M 1 x 0 . •ixorppAilnoitamnaegtn:e´naeriteoncoˆutplafonctia-dneruqeuˆtueerte dratiqueensupposant,enplusdelaline´arisationdesope´rateursd’observation H i eoldqeumle,`eMeupettˆilerae´n´sirnO.epourraalorsremplcareelom`dlee parsonapproximationline´airetangente(ousad´eriv´ee),etMsera alors le mode`leline´airetangent.Ils’agitdumeˆmeproce´d´equepourlesop´erateurs d’observationH´erag´enntlemeestrpiharcoi(v.)eCnest´cderpe´eesth`esypottteh validesilap´erioded’assimilationn’estpastroplongue.