Brevet FR2742876A1 — publié le 19970627

La présente invention se rapporte à un procédé de détermination du taux de précipitation et du facteur de réflectivité d'un hydrométéore par radar à impulsions à longueur d'onde atténuée par les hydrométéores et à un radar météorologique à double polarisation mettant en oeuvre ce procédé. En météorologie, on a besoin d'observer les chutes de pluies et autres hydrométéores notamment pour des problèmes de prévision à court terme. A cette fin, on a depuis de nombreuses années pensé à utiliser le radar pour déterminer quantitativement le taux de précipitation qui est un paramètre caractéristique de ces phénomènes. Pour cela on utilise des procédés basés sur la mesure de la réflectivité de la pluie ou des hydrométéores illuminés par un faisceau radar. Cependant pour la plupart de ces procédés, I'atténuation par les hydrométéores constitue un biais pour les mesures et il faut donc fonctionner à des longueurs d'onde aussi peu atténuées que possible, donc relativement élevées. Ceci entraine des difficultés de mise en oeuvre et un coût important. Un autre procédé plus attractif a été proposé celui de l'atténuation sur deux longueurs d'onde et une seule polarisation qui a l'avantage d'être indépendant de la calibration du radar. Mais ce procédé conduit à certaines difficultés techniques pour obtenir des mesures ayant la précision requise. En effet, on rencontre des problèmes relatifs à la coïncidence des faisceaux aux deux longueurs d'onde et à la simultanéité des mesures. De toutes les études faites, il ressort l'intérêt qu'il y aurait à utiliser des radars à plus faible longueur d'onde permettant de choisir des matériels plus compacts et moins coûteux. Cependant on se heurtait alors à des problèmes de correction d'atténuation très difficiles à résoudre, notamment en présence de fortes précipitations. L'invention a pour objet un procédé remédiant à ces inconvénients malgré l'utilisation de longueurs d'onde faibles entraînant une atténuation relativement élevée des ondes électromagnétiques sur le trajet aller et retour à travers le milieu absorbant constitué par les hydrométéores (pluie, grêle, nuages). L'invention est basée sur l'utilisation d'un radar à double polarisation pour mesurer les réflectivités selon les deux composantes de polarisation H et V en corrigeant les effets de l'atténuation sur les valeurs obtenues grâce aux mesures différentielles effectuées. Selon l'invention, il est donc prévu un procédé de détermination du taux de précipitation et du facteur de réflectivité d'un hydrométéore par radar à impulsions à longueur d'onde atténuée par ledit hydrométéore, ledit radar étant un radar à double polarisation fournissant une double série de mesures du facteur de réflectivité atténué (ZH" (r)ZvJ (r)) pour les différentes portes en distance selon chacune des polarisations, ledit procédé étant caractérisé en ce qu'il comprend les étapes suivantes A) déterminer, dans chaque porte en distance du radar, une première estimation du taux de précipitation (rus) à partir du facteur de réflectivité mesuré pour une seule polarisation par une combinaison des relations liant le facteur de réflectivité au taux de précipitation et le coefficient d'atténuation sur une seule polarisation au taux de précipitation; B) obtenir à partir de ladite première estimation du taux de précipitation une valeur estimée du facteur de réflectivité différentielle non atténué (ZDRS); C) déduire de ladite valeur estimée du facteur de réflectivité différentielle non atténué (ZDRs) et de la mesure radar du facteur de réflectivité différentielle atténué (ZDR) une valeur estimée de l'atténuation polarimétrique différentielle (ad);; et D) chercher à obtenir l'égalité entre l'intégrale (i (r, R))du taux de précipitation estimé dans l'étape A sur une seule polarisation du radar et l'intégrale ('2 (r, R)) du taux de précipitation déduite de ladite atténuation polarimétrique différentielle (ADP), par ajustement itératif d'un des paramètres (a, p) de la relation liant le facteur de réflectivité au taux de précipitation, de manière à pouvoir extraire les valeurs corrigées du facteur de réflectivité non atténué (ZHs (r)) et du taux de précipitation (Rs (r)) Le grand avantage de ce procédé est qu'il est indépendant de la calibration du radar et de l'atténuation par d'éventuels nuages non détectés car il est basé sur une mesure différentielle. Selon un autre aspect de l'invention, il est également prévu un radar météorologique à double polarisation fonctionnant à impulsions avec une longueur d'onde atténuée par les hydrométéores, pour la détermination du taux de précipitation et du facteur de réflectivité, ledit radar comportant: - une antenne à double polarisation associée à des moyens pour extraire les modes correspondant à deux voies de polarisation orthogonales H et 3 - un émetteur délivrant des impulsions aux deux voies de ladite antenne; - au moins une voie de réception pour traiter les signaux issus des deux voies de polarisation orthogonales; ledit radar étant caractérisé en ce qu'il comporte un dispositif de traitement numérique du signal prévu pour mettre en oeuvre le procédé décrit cidessus. L'invention sera mieux comprise et d'autres caractéristiques et avantages apparaîtront à l'aide de la description ci-après et des dessins joints ou: - la figure 1 est un schéma d'un radar météorologique selon l'invention; et - la figure 2 est un schéma de principe d'un dispositif de traitement numérique du signal pour la mise en oeuvre du procédé selon l'invention. Sur la figure 1 est représenté le schéma d'un radar météorologique bipolarisation. Ce radar comporte une antenne à double polarisation comprenant une source primaire 1 à double polarisation illuminant un réflecteur (non représenté). Cette source est associée à un extracteur de mode 2 (orthomode) à deux voies de polarisation H et V. Cette source primaire est alimentée par un émetteur à impulsions (6) fournissant des impulsions à fréquence de répétition adaptée à la portée désirée (par exemple Fr = i kHz pour une portée de 150 km). La puissance d'émission est distribuée aux deux voies de l'antenne par l'intermédiaire d'un té magique pour obtenir une polarisation orientée à 45" par rapport à l'horizontale ou par un coupleur 3dB, 5, pour obtenir une polarisation circulaire. Ces deux modes d'émission sont a priori équivalents puisqu'on veut émettre et recevoir deux signaux à polarisations orthogonales et les traiter en amplitude. La polarisation circulaire a un avantage secondaire en ce qui concerne le taux d'ondes stationnaires réfléchi par l'antenne vers l'émetteur. On prévoit deux voies de réception identiques traitant les signaux prélevés par deux circulateurs de puissance 3 et 4. Chaque voie comprend de manière classique un limiteur 7, 8, un amplificateur 9, 10, un mélangeur 11, 12 pour passer en fréquence intermédiaire, un amplificateur à fréquence intermédiaire 15, 16 puis un amplificateur logarithmique 17, 18 pour obtenir une grande gamme dynamique et un codeur analogique-numérique 19, 20. Les signaux codés sont dirigés vers un dispositif de traitement numérique selon l'invention 21 dont on verra le détail ultérieurement et qui restitue en sortie le facteur de réflectivité Z et le taux de précipitation R corrigés. Les mélangeurs 11 et 12 reçoivent par ailleurs le signal d'un oscillateur local 13 contrôlé par un circuit de commande automatique de fréquence 14. Pour certaines applications météorologiques où on veut obtenir des informations complémentaires sur la nature des hydrométéores, il est intéressant de mesurer la phase des signaux reçus par des voies de réception Doppler supplémentaires. Ces voies comprennent un circuit de démodulation vectorielle 22, 25 suivi de deux codeurs analogiquenumérique 23, 24 d'une part, 26, 27 d'autre part pour restituer les signaux codés en quadrature IHX QH pour la voie H et Iv, Qv pour la voie V. Le principe du procédé selon l'invention est le suivant. A partir d'une seule polarisation, par exemple la voie H, la relation entre le facteur de réflectivité atténué mesuré dans la pluie (ou autre hydrométéore) ZHa(r) à la distance r et le coefficient d'atténuation a peut s'écrire: où ZO(r) est le facteur de réflectivité non atténué, r est la largeur d'une porte en distance, a, est le coefficient d'atténuation propre à la porte en distance de rang i, exprimé en dB par unité de longueur, où Z11a(r) et Z0 (r) s'expriment en dBZ, c'est-à-dire = 10 log Zen mm6 m-3 et où r = nAr. Le coefficient d'atténuation a est lié au taux de précipitation R par la relation connue: a = k, R (2) où kHetyH sont des paramètres pour la polarisation H qui dépendent faiblement de la température et de la distribution de taille des gouttes (DSD pour "Drop Size Distribution") et qui peuvent être calculés pour chaque longueur d'onde. En combinant les relations (1 ) et (2), on obtient: Zo(r)=Z11a(r)+2rk11 11(r,R) (3) avec: lî représente l'intégrale du taux de précipitation le long du faisceau du radar pour la polarisation H choisie, c'est-à-dire qu'il traduit le cumul de l'atténuation dans les n-1 portes en distance précédant la porte n. Or il existe une relation connue entre le facteur de réflectivité Z dû à la pluie et le taux de précipitation R: Z=aRss (5) où aetfi sont des paramètres. On peut ainsi écrire la relation (4) sous la forme équivalente: Les paramètres aetfi dépendent fortement de la DSD. En utilisant des valeurs approchées pour les quatre paramètres k, y, &alpha; et ss dans les relations (3) et (5), on peut calculer des estimés corrigés Z'(r) à partir de la série de n valeurs mesurées Za(r) grâce une méthode itérative proposée par P.H. Hildebrand en 1978 dans "Iterative correction for attenuation of 5cm radar in rain" J.Appl.Meteor., 17, pages 508-514.Comme signalé dans cet article, des erreurs dans la DSD et la température supposées peuvent dégrader sérieusement ces estimés d'atténuation mais ces dégradations sont faibles par rapport à celles résultant d'erreurs de calibration du radar. Les paramètres concernés par l'erreur de calibration en radar étant ceux de la relation entre Z et R, l'incertitude essentielle sur I (r, R) résulte de l'incertitude sur les paramètres a et ss. Pour lever ces incertitudes, on fait donc appel à des mesures différentielles utilisant les deux polarisations d'un radar. La réflectivité différentielle ZDR est définie par: ZDR ZH - ZV (6) où Z s'exprime en dBZ et où les indices H et V correspondent aux deux polarisations du radar, par exemple horizontale et verticale. Pour les longueurs d'onde atténuées, en utilisant la relation (1) pour les polarisation H et V et en soustrayant les termes V des termes H, on obtient: ZDRa (r) = ZDRs(r)-2 ADP (r) (7) OÙ ZDRa est la mesure radar polarimétrique atténuée, ZDRs = ZOH - ZOV est le coefficient de réflectivité différentielle non atténué, dû seulement à la forme des gouttes, et A DP est l'atténuation polarimétrique différentielle pour un trajet de l'onde électromagnétique, soit: Avec une DSD supposée, le coefficient de réflectivité différentielle ZDRS dû aux formes des gouttes peut être calculé comme une fonction du taux de précipitation:: ZDRJ=f(R) (9) Les coefficients d'atténuation a,11 et au de la relation (8) sont donnés par la relation (2) pour les polarisations H et V. Pour la DSD supposée, les paramètres k11, kv, 711 et &gamma;v peuvent être calculés. Le calcul montre que la différence aH-a;, ne dépend pas de manière significative de la DSD. Pour une longueur d'onde donnée, ky etkg, sont nettement différents mais 711 ~ YlZ ~ 7. En combinant les relations (2) et (8), on a donc: ADP(r)=#r(kH-kv)I2(r,R)) (10) avec: On peut remarquer que l'on peut écrire la relation (11) comme la relation (4'), c'est-à-dire en fonction de Z. I(r,R) représente l'intégrale du taux de précipitation déduite de l'atténuation différentielle ADP,kH-kV étant supposé connu avec précision. Comme elle résulte d'une mesure différentielle, A DP ne dépend pas de la calibration du radar. Cela étant, dans l'estimation de I, et '2' seuls les paramètres aetfi sont variables en raison de la variabilité naturelle de la DSD dans la pluie et des erreurs de calibration. Dans le procédé selon l'invention, il est proposé de maintenir p à une valeur constante, égale à une valeur moyenne climatologique qui correspond à la longueur d'onde et au type de chute de pluie observé, et d'ajuster le paramètre a en utilisant l'atténuation différentielle ADP. Cette hypothèse sur p a été confirmée par de nombreuses mesures disponibles dans la littérature. Le principe est donc de déterminer I, à partir de la série de mesures sur une seule des polarisations en choisissant une valeur arbitraire intermédiaire de a, d'estimer I, à partir des mesures différentielles du radar et d'un calcul utilisant I, pour estimer le facteur de réflectivité différentielle non atténué et de comparer I, et I2 pour modifier éventuellement a jusqu'à ce que I1 = I2, ce qui permet alors d'extraire les estimations corrigées des coefficients de réflectivité non atténués Z, (r) et des taux de précipitation g(r). Un dispositif de traitement numérique permettant de mettre en oeuvre ce procédé est représenté sur la figure 2. La partie gauche de la figure concerne le traitement d'une seule polarisation, la polarisation H, cependant que la partie droite est relative au traitement différentiel. Les signaux logarithmiques fournis par le codeur 19 de la voie H sont retraduits en signaux linéaires ZHa(r) dans l'étape 30. Ces signaux sont moyennés sur N impulsions du radar (étape 31) puis on calcul le taux de précipitation estimé Ra à partir de paramètres a et ss prédéterminés, arbitrairement pour a, selon la relation (5) (étape 32), ceci pour chaque porte en distance i. L'étape 33 consiste à calculer les coefficients d'atténuation al à partir de la relation (2). On en déduit alors la valeur estimée T'(r,R) selon les relations (2) et (4) (étape 34). Par ailleurs, par les étapes 36 et 37, on calcule le facteur de réflectivité atténué ZHa(r) déduit de la série de mesures sur la seule polarisation H et, par l'étape 35, on estime le terme d'atténuation AH de la relation (3). En appliquant cette relation (3) par l'étape 38 consistant en une sommation, on obtient un estimé ZH,(r) du coefficient de réflectivité non atténué pour la polarisation H puis, par l'étape 47, un estimé du taux de précipitation corrigé Rs(r) en appliquant la relation (5) avec les paramètres a et ss déjà choisis. A partir de cette valeur estimée, on peut calculer par l'étape 48 une valeur estimée du facteur de réflectivité différentielle non atténué ZDR5(r) selon la relation (9). D'autre part, en partant des signaux codés des voies H et V fournis par 30 et 40, on effectue (étape 41) le rapport ZHa(r)lZvu,(r) qui donne la mesure radar polarimétrique atténuée ZDRa(r) moyennée ensuite sur N impulsions dans l'étape 42. La différence (étape 43) entre cette valeur moyennée et la valeur estimée du facteur de réflectivité différentielle non atténué fournie par 48 constitue l'estimation de l'atténuation polarimétrique différentielle A,(r). L'étape 44 fournit à partir des valeurs Ar et kH - kV la valeur estimée de 72 (r,R) selon la relation (10). Dans l'étape 45, on compare alors 7, et 12. Si I, i 12, alors l'étape 46 consiste à modifier le paramètre a d'un incrément Aa dans le sens convenable et à appliquer la nouvelle valeur a' aux étapes 32 et 47 pour recommencer les calculs. On renouvelle ces opérations jusqu'à ce que l'on obtienne 7, = 7. ce qui correspond à la bonne valeur donnée à a.On peut alors (étape 49) extraire les valeurs corrigées du facteur de réflectivité non atténué Z,(r) = Z111 (r) et du taux de précipitation R,(r) = Dans l'exposé du principe du procédé selon l'invention, on a supposé pour la relation (9) que ZDRs, de même que k11 et k,, pour la relation (10), étaient calculés pour la seule DSD supposée. Ceci est acceptable car, pour des conditions climatiques données, les DSD moyennées ont des paramètres stables. On peut donc choisir une DSD particulière à partir de considérations climatologiques. D'autre part, il est préférable de choisir des portes en distance à la distance r où le facteur de réflectivité est faible car alors ZDRS est faible et l'erreur d'estimation selon la relation (9) est, elle aussi, faible. Ce choix d'une porte particulière est d'autant plus justifié et facilité qu'on peut remplacer une porte à la distance r par une autre autour de r avec des facteurs de réflectivité différents sans que les valeurs '1 12 et a soient modifiées de façon significative. On peut d'ailleurs calculer a à partir de deux ou trois portes contiguës et utiliser la moyenne pour réduire les effets de dispersion de la DSD. On peut aussi remarquer qu'on peut subdiviser une direction radiale observée en plusieurs segments auxquels le procédé selon l'invention peut être appliqué successivement. Le procédé décrit ci-dessus est très peu sensible aux erreurs d'échantillonnage aléatoires sur le coefficient de réflectivité mesuré par le radar sur une polarisation ou sur les valeurs différentielles. Un des principaux avantages du procédé selon l'invention a déjà été mentionné : il s'agit de son indépendance vis à vis de la calibration du radar et de l'atténuation due à des nuages non détectés. On peut ajouter qu'une autre utilisation possible est la détection qualitative de la grêle. En effet quand de la grêle est présente dans une cellule de pluie, comme les grains de grêle ont une forme presque sphérique, ils produisent un facteur de réflectivité élevé pour les deux polarisations avec Z11 = Z,, et donc un facteur de réflectivité différentielle sensiblement nul. Par contre, on a une forte atténuation indépendante de la polarisation d'où une contribution nulle à l'atténuation différentielle AD?. Ainsi la présence de grêle modifie I, mais pas de manière correspondante I2, contrairement au cas d'une pluie homogène. On constate alors une modification notable de a qui peut être utilisée comme indicateur de grêle. Il faut également remarquer que, si la description ci-dessus a été faite dans le cadre de l'utilisation de l'atténuation polarimétrique différentielle ADP, on pourrait aussi faire fonctionner le procédé selon l'invention sur le déphasage différentiel KDP = KH - K,, exprimé en degré.km de propagation. Bien entendu, les exemples décrits ne sont nullement limitatifs de l'invention. En particulier, le radar météorologique décrit pourrait utiliser une seule voie de réception au lieu des deux voies décrites, les signaux des deux polarisations étant commutés en entrée pour utiliser la voie en partage dans le temps et "décommutés" en sortie. Le rythme de commutation doit simplement être choisi pour que les signaux reçus sur les polarisations H et V ne soient pas décorrélés pendant la durée de la commutation. D'autre part, peut remarquer que les signaux I et O disponibles sur les sorties des voies Doppler permettent de calculer l'amplitude de signal dans chacune des voies H et Y, amplitude qui pourrait être utilisée à la place des signaux issus des codeurs 19 et 20 pour éviter l'utilisation des voies 17, 19 et 18, 20. L'inconvénient serait cependant de ne pas disposer de la dynamique résultant de l'emploi d'amplificateurs logarithmiques. Enfin, dans le dispositif de traitement selon l'invention de la figure 2, on pourrait penser à continuer à utiliser les signaux logarithmiques au lieu de repasser en signaux linéaires. REVENDICATIONS 1. Procédé de détermination du taux de précipitation et du facteur de réflectivité d'un hydrométéore par radar à impulsions à longueur d'onde atténuée par ledit hydrométéore, ledit radar étant un radar à double polarisation fournissant une double série de mesures du facteur de réflectivité atténué (ZHa (r), ZvO (r)) pour les différentes portes en distance selon chacune des polarisations, ledit procédé étant caractérisé en ce qu'il comprend les étapes suivantes A) déterminer, dans chaque porte en distance du radar, une première estimation du taux de précipitation (Rs) à partir du facteur de réflectivité mesuré pour une seule polarisation par une combinaison des relations liant le facteur de réflectivité au taux de précipitation et le coefficient d'atténuation sur une seule polarisation au taux de précipitation; B) obtenir à partir de ladite première estimation du taux de précipitation une valeur estimée du facteur de réflectivité différentielle non atténué (Z); C) déduire de ladite valeur estimée du facteur de réflectivité différentielle non atténué (ZDRS) et de la mesure radar du facteur de réflectivité différentielle atténué (ZDRa) une valeur estimée de l'atténuation polarimétrique différentielle (ad);; et D) chercher à obtenir l'égalité entre l'intégrale (I, (r, R))du taux de précipitation estimé dans l'étape A sur une seule polarisation du radar et l'intégrale (I. (r, R)) du taux de précipitation déduite de ladite atténuation polarimétrique différentielle (AD?), par ajustement itératif d'un des paramètres (a, p) de la relation liant le facteur de réflectivité au taux de précipitation, de manière à pouvoir extraire les valeurs corrigées du facteur de réflectivité non atténué (ZH5 (r)) et du taux de précipitation (Rs (r)). 2. Procédé selon la revendication 1, caractérisé en ce que l'étape A) comprend les étapes suivantes: a) calculer sur une première desdites séries de mesures (ZHa (r)) des valeurs estimées du taux de précipitation (R') pour les différentes portes en distance (1 à n-1) à partir de paramètres (a, p) fixés arbitrairement; b) déduire desdites valeurs estimées du taux de précipitation (R) une première valeur estimée (II (r, R)) caractéristique de l'atténuation due aux précipitations, ladite première valeur estimée étant l'intégrale du taux de précipitation estimé sur une seule polarisation;; c) déterminer à partir de ladite première valeur estimée (I, (r, R)) et de la première série de mesures (ZHa (r)) une valeur estimée du facteur de réflectivité non atténué (ZH5 (r)) et du taux de précipitation corrigé (R1). 3. Procédé selon l'une des revendications 1 ou 2, caractérisé en ce que ladite étape B) comprend l'étape suivante: d) calculer, à partir de ladite valeur estimée du taux de précipitation corrigé (rus) et d'une distribution donnée des gouttes (DSD), une valeur estimée du facteur de réflectivité différentielle non atténué (ZDRS (r))- 4.Procédé selon l'une des revendications 1 à 3, caractérisé en ce que ladite étape C) comprend les étapes suivantes: e) calculer, à partir de ladite double série de mesures du facteur de réflectivité atténué (ZHa (r), Zva (r)), la mesure du facteur de réflectivité différentielle atténué (ZDRa (r)); f) déduire de la mesure du facteur de réflectivité différentielle atténué (ZDRa (r)) et de la valeur estimée du facteur de réflectivité différentielle non atténué (ZDRs (r)) obtenue lors de l'étape d) une valeur estimée de l'atténuation polarimétrique différentielle (ADP). 5. Procédé selon l'une quelconque des revendications 1 à 4, caractérisé en ce que ladite étape D) comprend les étapes suivantes: g) déduire de ladite valeur estimée de l'atténuation polarimétrique différentielle (AD?) une seconde valeur estimée (I.(r,R)) caractéristique de l'atténuation due aux précipitations, qui est l'intégrale du taux de précipitation, h) comparer lesdites première (I, (r, R)) et seconde (1, (r, R)) valeurs estimées caractéristiques de l'atténuation due aux précipitations, i) si ladite comparaison selon l'étape h) indique que desdites première et seconde valeurs estimées sont différentes, modifier au moins un desdits paramètres (a, p) de l'étape a) et reprendre les étapes a) à h);; j) si ladite comparaison selon l'étape h) indique que lesdites première et seconde valeurs estimées sont égales, extraire lesdites valeurs estimées du facteur de réflectivité non atténué (ZH5(r)) et du taux de précipitation corrigé (R,). 6. Procédé selon l'une quelconque des revendications 2 à 5, caractérisé en ce que l'étape a) consiste à: a. 1) moyenner lesdites mesures (ZHa (r)) de ladite première série sur N impulsions du radar; a.2) calculer pour chaque porte en distance (i) des valeurs estimées du taux de précipitation (R) à partir de paramètres (a, p) fixés arbitrairement selon la relation: 7.Procédé selon l'une quelconque des revendications 2 à 6, caractérisé en ce que l'étape b) consiste à: b.1) calculer à partir des valeurs estimées du taux de précipitation (R,) les coefficients d'atténuation (a,) dans chaque porte en distance selon la relation: a, / kH = R# YH où k11 et 711 sont des paramètres prédéterminés pour une température et une distribution de taille de goutte données; et b.2) déduire de ces coefficients d'atténuation (a,) ladite première valeur estimée (t(r, R)) selon la relation où n est le nombre de portes en distance pour une distance r du radar. 8. Procédé selon l'une quelconque des revendications 2 à 7, caractérisé en ce que l'étape c) consiste à: c.l) calculer, à partir de ladite première valeur estimée (I,(r,R)), une valeur d'atténuation A telle que: A=2 #r kH I1(r,R) où Ar est la largeur d'une porte en distance; c.2) moyenner lesdites mesures (Z11a(r)) de ladite première série sur N impulsions du radar; c.3) faire la somme de ladite atténuation A et de la valeur moyennée du facteur de réflectivité atténué ZHa(r) pour obtenir une valeur estimée Z111 (r) du facteur de réflectivité non atténué; et c.4) calculer une valeur estimée R, du taux de précipitation corrigé selon la relation: 9.Procédé selon l'une quelconque des revendications 4 à 8, caractérisé en ce que l'étape e) consiste à: e.1) faire le rapport des mesures des séries de ladite double série (ZNa(r) Zwa (r)), pour obtenir les valeurs du facteur de réflectivité différentielle atténue (ZDRa(r)) du radar; e.2) moyenner lesdites valeurs du facteur de réflectivité différentielle atténué sur N impulsions du radar. 10. Procédé selon la revendication 5, caractérisé en ce que l'étape f) consiste à effectuer la différence entre ladite valeur estimée du facteur de réflectivité différentielle non atténué (ZDRs(r)) obtenue selon l'étape d) et ladite valeur moyennée du facteur de réflectivité différentielle atténué (ZDRa (r)) pour obtenir une valeur estimée de l'atténuation polarimétrique différentielle (ADp(r)) sur un trajet radar. 11. Procédé selon l'une quelconque des revendications 5 à 10, caractérisé en ce que l'étape g) consiste à calculer, à partir de ladite valeur estimée de l'atténuation polarimétrique différentielle (ADp(r)), ladite seconde valeur estimée 12 (r, R) selon la relation: ADP(r) I2(r,R) = #r (kH - kV) où k, est un paramètre prédéterminé de la même manière que k11 pour la seconde série de mesures (Z;, (r)). 12. Radar météorologique à double polarisation fonctionnant à impulsions avec une longueur d'onde atténuée par les hydrométéores, pour la détermination du taux de précipitation et du facteur de réflectivité, ledit radar comportant: - une antenne (1) à double polarisation associée à des moyens (2, 3, 4) pour extraire les modes correspondant à deux voies de polarisation orthogonales H et V - un émetteur (5, 6) délivrant des impulsions aux deux voies de ladite antenne -au moins une voie de réception (7, 9, 11, 13, 15, 17, 19; 10,12, 13, 16, 18, 20) pour traiter les signaux issus des deux voies de polarisation orthogonales; ledit radar étant caractérisé en ce qu'il comporte un dispositif de traitement numérique du signal (21) prévu pour mettre en oeuvre le procédé selon l'une quelconque des revendications 1 à 11. 13. Radar selon la revendication 12, caractérisé en ce qu'il comprend à la sortie desdites voies de réception des amplificateurs logarithmiques (17, 18) pour augmenter la gamme dynamique et en ce que ledit dispositif de traitement numérique (21) comporte en entrée des moyens de conversion (30, 36, 40) pour retrouver les valeurs des signaux reçus. 14. Radar selon l'une des revendications 12 ou 13, caractérisé en ce qu'il comprend en outre deux voies de réception Doppler (22 à 24 ; 25 à 27) associées respectivement aux sorties des voies de réception des deux polarisations. 15. Radar selon l'une quelconque des revendications 12 à 14, caractérisé en ce qu'il comporte une seule voie de réception associée à des moyens de commutation en entrée et en sortie commandés en synchronisme pour traiter alternativement les deux polarisations, le rythme de commutation étant tel que les signaux reçus selon les deux polarisations ne sont pas décorrélés pendant la durée de la commutation.