Etude théorique de radars géologiques - Epublications - Université ...
Etude théorique de radars géologiques - Epublications - Université ...
Etude théorique de radars géologiques - Epublications - Université ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
Samuel BESSE : Étu<strong>de</strong> Théorique <strong>de</strong> Radars Géologiques Analyses <strong>de</strong> sols, d'antennes et interprétation <strong>de</strong>s signaux 41<br />
Phase d'initialisation 1 ère itération : étape 1 1 ère itération : étape 2 1 ère itération : étape 3<br />
Figure 33 : Représentation <strong>de</strong> l'algorithme <strong>de</strong> la métho<strong>de</strong> dite DiamondSquare4D sur un motif élémentaire.<br />
L'application d'une itération génère huit cubes sur lesquels on peut répéter l'algorithme.<br />
Comme dans le cas <strong>de</strong>s surfaces la valeur moyenne quadratique est liée à la valeur du<br />
paramètre h et la répartition <strong>de</strong>s pondérations suit une loi gaussienne. L'expression (20) se déduit<br />
d'une étu<strong>de</strong> statistique sur un grand nombre <strong>de</strong> réalisations (une centaine). L'étu<strong>de</strong> statistique révèle<br />
également que h est indépendante du nombre <strong>de</strong> motifs (NX.NY.NZ) et converge pour un<br />
nombre d'itérations Nit supérieur ou égal à quatre.<br />
h=0,31 0,45 ⋅h−1,3 ⋅h 2 1,4 ⋅h 3<br />
Il est souvent important <strong>de</strong> connaître la pondération maximale possible notamment dans les<br />
applications <strong>de</strong> la troisième partie. Par exemple, si l'on utilise la matrice pour modéliser un sol<br />
hétérogène en permittivité avec la métho<strong>de</strong> FDTD, il faut être certain <strong>de</strong> mailler suffisamment<br />
finement pour limiter la dispersion numérique (la taille d'une maille doit être inférieure à λ milieu/10).<br />
Pour cela, l'abaque présenté figure 34 permet <strong>de</strong> déterminer la valeur <strong>de</strong> pondération maximale<br />
<strong>théorique</strong> présente dans le matériau. Il est obtenu par un calcul numérique en substituant toutes les<br />
variables aléatoires intervenant dans l'algorithme par la borne supérieure <strong>de</strong> l'intervalle dans lequel<br />
elles sont choisies. Dans la pratique, le majorant <strong>théorique</strong> a une probabilité infime d'être atteint et<br />
cette probabilité est d'autant plus petite que le nombre d'itérations est important.<br />
Le passage <strong>de</strong> la matrice normalisée M0 à la matrice <strong>de</strong> permittivité s'obtient par (21) où les<br />
paramètres 〈 r 〉 et k laissés au choix <strong>de</strong> l'utilisateur désignent respectivement la permittivité<br />
moyenne et un facteur multiplicatif.<br />
r i , j , k = 〈 r 〉k⋅M 0 i , j , k (21)<br />
(20)