FR 3077383 A1 20190802 FR 1850782 A 20180131 L'invention concerne le domaine des dispositifs pour convertir l'énergie des vagues en énergie électrique. En particulier, l'invention concerne le domaine de la détermination de la force d'excitation par la houle incidente sur un moyen mobile, notamment pour la commande d'un système houlomoteur. Les ressources énergétiques renouvelables connaissent un fort intérêt depuis quelques années. Ces ressources sont propres, gratuites et inépuisables, autant d'atouts majeurs dans un monde rattrapé par la diminution inexorable des ressources fossiles disponibles et prenant conscience de la nécessité de préserver la planète. Parmi ces ressources, l'énergie houlomotrice, source relativement méconnue au milieu de celles largement médiatisées comme l'éolien ou le solaire, contribue à la diversification indispensable de l'exploitation des énergies renouvelables. Les dispositifs, couramment appelés appareils "houlomoteurs", sont particulièrement intéressants, car ils permettent de produire de l'électricité à partir de cette source d'énergie renouvelable (l'énergie potentielle et cinétique des vagues) sans émission de gaz à effet de serre. Ils sont particulièrement bien adaptés pour fournir de l'électricité à des sites insulaires isolés. Par exemple, les demandes de brevet FR 2876751, et WO 2009/081042 décrivent des appareils pour capter l'énergie produite par le flot marin. Ces dispositifs sont composés d'un support flottant sur lequel est disposé un pendule monté mobile par rapport au support flottant. Le mouvement relatif du pendule par rapport au support flottant est utilisé pour produire de l'énergie électrique au moyen d'une machine de conversion de l'énergie (par exemple une machine électrique). La machine de conversion fonctionne en tant que générateur et en tant que moteur. En effet, pour fournir un couple ou une force qui entraîne le mobile, on fournit une puissance à la machine de conversion afin de mettre le mobile en résonance avec les vagues (mode moteur). En outre, pour produire un couple ou une force qui résiste au mouvement du mobile, on récupère une puissance via la machine de conversion (mode générateur). Dans l'intention d'améliorer le rendement et donc la rentabilité des dispositifs pour convertir l'énergie des vagues en énergie électrique (systèmes houlomoteurs), il est intéressant d'estimer la force d'excitation exercée par la houle incidente sur le système houlomoteur. La force appliquée par la houle aux dispositifs pour convertir l'énergie des vagues en énergie électrique ou d'autres formes d'énergie exploitable, couramment appelés houlomoteurs, n'est pas directement mesurable en conditions de fonctionnement normales dans un contexte industriel. Pour certains systèmes houlomoteurs, il est possible de réaliser des essais spécifiques, consistant à bloquer la partie flottante de la machine et à mesurer 2 l'effort nécessaire pour contrer l'action de la vague (c'est-à-dire pour maintenir le flotteur immobile) à travers un ou plusieurs capteurs (de force ou de couple) placés sur la machine de conversion, appelée aussi système de prise de puissance (PTO), par exemple un générateur électrique assorti d'un dispositif permettant de transmettre le mouvement 5 oscillatoire du moyen mobile. Selon la terminologie courante, on appelle la quantité ainsi mesurée, force d'excitation de la houle (ou vague), en la distinguant notamment de la force (et de la houle) générée par le mouvement même du flotteur (force de radiation). Durant le fonctionnement normal du flotteur, en revanche, les mêmes capteurs d'effort mesureront seulement la force appliquée par le PTO, et non pas la force d'excitation de la houle 10 incidente. Il serait possible, en principe, de calculer l'ensemble des forces appliquées par la houle en intégrant des mesures de pression provenant de capteurs distribués sur toute la surface, mais il s'agit là d'une solution coûteuse et peu robuste qui peut difficilement être envisagée dans un contexte industriel. 15 Il y a peu de travaux scientifiques dans le domaine de l'estimation en temps-réel de la force (ou du couple) d'excitation de la vague. Par exemple, on peut citer le document Peter Kracht, Sebastian Perez-Becker, Jean-Baptiste Richard, and Boris Fischer. "Performance Improvement of a Point Absorber Wave Energy Converter by Application of an Observer- 20 Based Control: Results from Wave Tank Testing". Dans IEEE Transactions on Industry Applications, juillet 2015, qui décrit un algorithme d'estimation de la force de la houle, celui-ci est basé sur banc d'oscillateurs harmoniques indépendants et un observateur de Luenberger. Les résultats obtenus avec cette méthode montrent que les estimations sont retardées (déphasées) de manière significative par rapport à la vraie force d'excitation, ce 25 qui est gênant notamment pour une utilisation dans le cadre du contrôle du système houlomoteur. Plus généralement, cette approche considère une plage de fréquence de fréquence fixe (choisie a priori) pour le spectre de la houle. Pour que la méthode puisse marcher dans des conditions réalistes, où le spectre est variant dans le temps, il faut prendre en compte un nombre très élevé de fréquences, ce qui rend l'approche très lourde en termes 30 de calculs. Dans le document Bradley A Ling. "Real-Time Estimation and Prediction of Wave Excitation Forces for Wave Energy Control Applications",publié dans ASME 2015 34th International Conference on Ocean, Offshore and Arctic Engineering, une approche basée sur un filtre de Kalman étendu est proposée pour reconstruire le spectre en considérant la force d'excitation de la vague comme une (seule) sinusoïde dont les paramètres (amplitude, 35 fréquence et phase) varient dans le temps. Toutefois, la méthode ne peut être efficace que pour des vagues dans une bande de fréquences très étroite. 3 En outre, la demande de brevet FR2973448 (US9261070) décrit une méthode de commande pour convertisseurs ponctuels oscillants dont une des étapes est l'estimation (du spectre) de la force de la houle sur le flotteur (ou du couple sur la partie mobile du convertisseur), à partir d'un ensemble de sinusoïdes et d'un observateur de Luenberger. La 5 méthode s'apparente à celle publiée dans les documents précédemment cités, et présente a priori les mêmes désavantages. Une autre méthode est décrite dans la demande de brevet FR 3049989 (WO 2017/174244). Cette méthode consiste à déterminer en temps réel les efforts exercés par la 10 houle incidente sur le moyen mobile, afin d'adopter en conséquence les meilleurs réglages des stratégies de production de l'énergie électrique. La méthode est basée sur la construction d'un modèle de la force de radiation exercée sur le moyen mobile, et d'un modèle de la dynamique du système houlomoteur. Cette méthode utilise seulement des mesures de la cinématique du flotteur (position, vitesse et éventuellement accélération) et de 15 la force appliquée par la machine de conversion, mesures normalement disponibles sur un système houlomoteur, car utilisées pour son contrôle et pour sa supervision. Ainsi, grâce aux modèles utilisés, il est possible d'estimer la force exercée par la houle pour toute plage de fréquences des vagues, tout en gardant un temps de calcul adapté à une détermination et un contrôle en temps réel. 20 Cette méthode donne satisfaction pour les systèmes houlomoteurs « linéaires », c'est à dire, dont le comportement dynamique est décrit de manière suffisamment précise par des équations différentielles linéaires. Néanmoins, cette méthode ne permet pas de déterminer la force de la houle de manière précise pour les système houlomoteurs « non linéaires » (avec des équations différentielles non linéaires), notamment les systèmes houlomoteurs de type à 25 volet (« flap »), ou les systèmes houlomoteurs de type à convertisseur submergé à différentiel de pression de forme plate, ou les systèmes houlomoteurs de type machine à masse tournante tanguant en surface. Afin de pallier ces inconvénients, la présente invention consiste en un procédé de 30 détermination en temps réel des efforts exercés par la houle incidente sur un moyen mobile d'un système houlomoteur, afin d'adopter en conséquence les meilleurs réglages des stratégies de récupération de l'énergie électrique. La méthode est basée sur la construction d'un modèle de la force de radiation exercée sur le moyen mobile, d'un modèle de la force de traînée exercée sur le moyen mobile et d'un modèle non linéaire de la dynamique du 35 système houlomoteur. L'invention utilise seulement des mesures de la cinématique du flotteur (position, vitesse et éventuellement accélération) et de la force appliquée par la machine de conversion, mesures normalement disponibles sur un système houlomoteur, car 4 utilisées pour son contrôle et pour sa supervision. La détermination de la force d'excitation exercée par la houle incidente sur le moyen mobile utilise ces modèles, ces mesures, et un filtre de Kalman sans parfum. Il intègre notamment un modèle de la force de traînée visqueuse et rend possible, grâce au filtre de Kalman sans parfum, d'estimer la force 5 exercée par la houle incidente de manière précise pour tous les systèmes houlomoteurs, en particulier pour les systèmes houlomoteurs « non linéaires ». Le procédé selon l'invention est adapté pour toute plage de fréquences des vagues, tout en gardant un temps de calcul adapté à une détermination et un contrôle en temps réel. 10 Le procédé selon l'invention La présente invention concerne un procédé de détermination d'une force d'excitation exercée par la houle sur un moyen mobile d'un système houlomoteur, ledit système houlomoteur convertissant l'énergie de la houle en énergie au moyen dudit moyen mobile coopérant avec une machine de conversion, ledit moyen mobile effectuant un mouvement 15 par rapport à ladite machine de conversion sous l'action de la houle. Pour ce procédé, on réalise les étapes suivantes : a) on mesure la position et la vitesse dudit moyen mobile ; b) on mesure la force u exercée par ladite machine de conversion sur ledit moyen mobile ; 20 c) on construit un modèle de la force Arad de radiation exercée sur ledit moyen mobile, ledit modèle de la force de radiation reliant ladite force de radiation i, à la vitesse dudit moyen mobile ; d) on construit un modèle de force de traînée id exercée sur ledit moyen mobile, ledit modèle de force de traînée reliant ladite force de traînée id à la vitesse dudit moyen 25 mobile ; e) on construit un modèle de dynamique dudit système houlomoteur qui relie ladite force d'excitation -cw exercée par la houle incidente sur ledit moyen mobile, à ladite position dudit moyen mobile, à ladite vitesse du moyen mobile, à ladite force u exercée par ladite machine de conversion sur ledit moyen mobile, à ladite force de radiation Tract 30 exercée sur ledit moyen mobile et à ladite force de traînée id exercée sur ledit moyen mobile ; et f) on détermine ladite force d'excitation -cw exercée par la houle incidente sur ledit moyen mobile au moyen dudit modèle de dynamique, dudit modèle de la force de radiation, dudit modèle de la force de traînée 'Cd, desdites position et vitesse 35 mesurées et de ladite force u mesurée exercée par ladite machine de conversion sur 5 ledit moyen mobile et au moyen d'un filtre de Kalman sans parfum construit à partir d'un modèle de marche aléatoire de ladite force d'excitation exercée par la houle incidente sur ledit moyen mobile. 5 Selon un mode de réalisation, on construit ledit modèle de dynamique dudit système houlomoteur au moyen d'une équation du type : ied «t) = Thd(t)+T'd(t)+Td(t) + Tw(t) - u(t) avec ied le moment d'inertie total dudit moyen mobile (2), 8(0 l'angle de rotation dudit moyen mobile par rapport à la position d'équilibre, avec ..(t) l'accélération angulaire dudit moyen mobile et .(t) la vitesse angulaire dudit 10 moyen mobile, Thd(t) le moment de rappel hydrostatique, Trad(t) le moment de radiation, Td(t) le moment de trainée, Tw(t) le moment d'excitation de la vague, u(t) le moment exercé par ladite machine de conversion sur ledit moyen mobile. Conformément une mise en oeuvre, on détermine le moment de rappel hydrostatique Thd(t) au moyen d'une formule du type : Thd(t) = -K8(t) 15 où K est le coefficient de raideur hydrostatique. Selon un aspect, on construit ledit modèle de force de radiation par une équation du type : Trad(t) = -Ico«t) - Tr(t) avec 40 est le moment d'inertie ajouté à fréquence infiniment élevée et 20 Tr(t) = fot h(t - s)«s)ds = h(t) * .(t), avec h la réponse impulsionnelle qui relie la vitesse du moyen mobile à l'amortissement de radiation et .(t) la vitesse angulaire dudit moyen mobile. Avantageusement, on construit ledit modèle de force de trainée par une équation du type Td(t) = p8(op(t)1 25 où fi est le coefficient de trainée et .(t) la vitesse angulaire dudit moyen mobile. Alternativement, ledit modèle de la force de trainée Fd peut s'écrire par une équation de la forme : Fd(t) = R2(t)12(t)1 avec fl le coefficient de trainée, et 2 la vitesse du moyen mobile. 6 Conformément à une mise en oeuvre, on détermine ladite force d'excitation Ti,/ exercée par la houle incidente sur ledit moyen mobile en mettant en oeuvre les étapes suivantes : i) on initialise k=0, le vecteur d'état 2,(010)=m(0) et l'état de la matrice de covariance, P(010) = P0 ; 5 ii) à chaque instant k, on acquiert lesdites mesures de position et de vitesse dudit moyen mobile y(k) = [8(k) ((k)] et la mesure de ladite force exercée par ladite machine de conversion sur ledit moyen mobile u(k) ; et iii) à chaque instant k, on détermine la force d'excitation exercée par la houle incidente sur ledit moyen mobile (2) tw(k) au moyen des équations suivantes : 10 îw(k) = [0 l]xafklk) Px(k1k) = Px(klk - 1) - KPy(klk - 1)KT K = Pxy(klk - 1)Py(klk - 1)-1 2(k) = 2a(klk - 1) + K(3î(k)-3î(klk - 1)) 2n Pxy(klk - 1) =IVVic (2j(klk - 1)- xa(klk - 1)) (j77/(klk - 1)- y(kIk - 1)) T i=0 jkklk - 1) = Ej??1,3V1InyT(klk - 1), Py(klk - 1) = Ej?"20VVic (j77/(klk - 1)- y(kIk - 1)) (j77/(klk - 1)- y(kIk - 1)) +R 15 MkIk - 1) = Ca2j(klk - 1) 2(kIk - 1) = Ej??2'0VVirn2j(klk - 1), Px(kik - 1) =e2'0Wjc (2j(klic - 1)- xa(klk - 1))(2j(klk - 1)- xa(klk - 1))T + Q. îcj(klk - 1) = Aaxj(k - 1) + fa (xj(k - 1)) + Bau(k - 1), j = 0,1, ...,2n xo(k - 1) = £a (k - 11k - 1), xi(k - 1) = 2a(k - 11k - 1) + À Si(k - 1),i = 1,2, ..., n xi+n(k - 1) = £a (k- 11k - 1)-nE À Si (k - 1),i = 1,2, ..., n 20 où Si(k - 1) est la i-ème colonne de la racine carrée matricielle de Px(k - 11k - 1), soit Px(klk - 1) = S(k - 1)T S(k - 1) xk))1 5 xa(k) le vecteur d'état inconnu xa(k) = [ Tw( Aa,Ba,Ca,Da des matrices de la (k réalisation d'état, P la matrice de covariance du vecteur d'état, Q et R des matrices de calibration. 25 Selon un mode de réalisation, on commande ledit système houlomoteur en fonction de ladite force déterminée d'excitation Ti,/ exercée par la houle incidente sur ledit moyen mobile. 7 De manière avantageuse, ledit système houlomoteur est un système houlomoteur de type à convertisseur submergé à différentiel de pression de forme plate, de type machine à masse tournante tanguant en surface ou du type système houlomoteur à volet. 5 En outre, l'invention concerne un système houlomoteur. Il comprend des moyens de commande, notamment des moyens de commande informatique, mettant en oeuvre un procédé selon l'une des caractéristiques précédentes. Présentation succincte des figures 10 D'autres caractéristiques et avantages du procédé selon l'invention, apparaîtront à la lecture de la description ci-après d'exemples non limitatifs de réalisations, en se référant aux figures annexées et décrites ci-après. La figure 1 illustre un système houlomoteur de type à volet selon un mode de réalisation de l'invention. 15 La figure 2 illustre une estimation de la force exercée par la houle au moyen d'un procédé selon l'art antérieur pour une vague régulière. La figure 3 illustre une estimation de la force exercée par la houle au moyen d'un procédé selon un mode de réalisation de l'invention pour une vague régulière. La figure 4 illustre une estimation de la force exercée par la houle au moyen d'un 20 procédé selon un mode de réalisation de l'invention pour une vague irrégulière. Description détaillée de l'invention La présente invention concerne un procédé de détermination de la force d'excitation exercée par la houle incidente sur un moyen mobile d'un système houlomoteur, appelée 25 aussi force d'excitation de la houle ou de la vague. Un système houlomoteur est un système qui convertit l'énergie de la houle en énergie récupérable, en particulier en énergie électrique. Un système houlomoteur comporte généralement un moyen mobile, également appelé volet (« flap »), pendule ou flotteur, qui a un mouvement oscillatoire sous l'action de la houle. Le moyen mobile coopère avec une machine de conversion, appelée aussi système 30 de prise de puissance (PTO de l'anglais « power take off »), qui comporte dans la plupart des cas un générateur électrique couplé à un dispositif permettant d'adapter la transmission du mouvement oscillatoire, pour convertir le mouvement du moyen mobile en énergie récupérable. Dans certains cas, la machine de conversion peut agir en tant que moteur en générant un effort sur le moyen mobile. En effet, pour récupérer une puissance via la 35 machine de conversion, on produit un couple ou une force qui résiste au mouvement du mobile (mode générateur). En outre, si la machine de conversion le permet, on peut fournir 8 une puissance à la machine de conversion pour fournir un couple ou une force qui entraîne le mobile afin de l'aider à se mettre en résonance avec les vagues (mode moteur). Le procédé selon l'invention est adapté à tout type de système houlomoteur avec au moins un moyen mobile, par exemple ceux décrits dans la demande de brevet FR2973448 5 (US9261070). Le procédé de commande selon l'invention peut également être appliqué à un système houlomoteur appartenant à la catégorie de systèmes houlomoteurs avec colonnes d'eau oscillantes (OWC de l'anglais Oscillating Water Column). Toutefois, le procédé selon l'invention est particulièrement adapté à un système houlomoteur non linéaire, par exemple du type à volet (« flap »), de type à convertisseur submergé à différentiel de pression de 10 forme plate, ou de type machine à masse tournante tanguant en surface. La figure 1 illustre, de manière schématique, un exemple non limitatif de système houlomoteur 1 à volet. Le système houlomoteur 1 comporte un moyen mobile 2, qui a la forme d'un volet sensiblement vertical (la houle est représentée de manière schématique par deux traits courbes et la direction de la houle est représentée par la flèche H). Le volet 2 15 possède une partie immergée et une partie émergée. Le moyen mobile 2 coopère, au moyen d'une articulation 4, avec une machine de conversion basée sur un générateur électrique 3, qui est, dans ce cas, un générateur rotatif. L'articulation 4 permet un mouvement d'oscillation rotatif du volet 2 autour d'un axe horizontal par rapport à un support 5 qui est fixe par rapport au fond de la mer 6. 20 Dans la suite de la description et pour les revendications, les termes vagues, flots marins, et houle sont considérés comme équivalents. De plus, dans la description et pour les revendications, le terme force désigne un effort ou un couple. De la même manière, les termes position, vitesse et accélérations désignent à 25 la fois des valeurs « linéaires » ou « angulaires », les valeurs « linéaires » pouvant être associées aux efforts, et les valeurs « angulaires » pouvant être associées aux couples. Dans la suite de la description, seul le cas relatif aux couples sera illustré, mais le cas des efforts peut se déduire par transposition des équations dans un référentiel orthogonal. En outre, pour une meilleure compréhension, les différents modèles sont représentés 30 en une seule dimension, toutefois, le procédé selon l'invention est adapté à des modèles à plusieurs dimensions, pour des systèmes ayant plusieurs degrés de liberté dans leur mouvement. Le procédé selon l'invention comporte les étapes suivantes : 35 1) mesure de la position et de la vitesse du moyen mobile 2) mesure de la force exercée par la machine de conversion (ou PTO) 3) construction du modèle de force de radiation 9 4) construction d'un modèle de force de traînée 5) construction du modèle de dynamique 6) détermination de la force d'excitation de la houle incidente 7) (étape facultative) commande du système houlomoteur 5 Les étapes 3) et 4) peuvent être réalisées dans cet ordre, dans l'ordre inverse ou simultanément. Les étapes du procédé peuvent être mise en oeuvre par un équipement informatique (ou calculateur). Cet équipement peut comprendre des moyens de traitement de données, et 10 avantageusement des moyens de stockage de données. Les moyens de traitement de données sont configurés pour implémenter : - les mesures de position et de vitesse du moyen mobile, - la mesure de la force exercée par la machine de conversion, - la construction des modèles de force de radiation, de force de traînée et du modèle 15 dynamique, et - la détermination de la force d'excitation de la houle incidente. 1) Mesure de la position et de la vitesse du moyen mobile Lors de cette étape, la position et la vitesse du moyen mobile sont mesurées. La 20 position correspond au mouvement (par exemple distance ou angle) par rapport à la position d'équilibre du moyen mobile. Ces mesures peuvent être réalisées au moyen de capteurs, généralement implémentés sur un système houlomoteur pour son contrôle et/ou sa supervision. Selon une mise en oeuvre de l'invention, lors de cette étape, on peut également 25 mesurer l'accélération du moyen mobile, et l'utiliser pour estimer la vitesse, ou bien l'utiliser directement dans les modèles mis en oeuvre par le procédé selon l'invention. Par exemple, l'accélération peut être mesurée au moyen d'un accéléromètre placé sur le moyen mobile. 2) Mesure de la force exercée par la machine de conversion (PTO) 30 Lors de cette étape, on mesure la force (l'effort ou le cas échéant le couple) exercée par la machine de conversion PTO sur le moyen mobile. Cette mesure peut être réalisée au moyen d'un capteur, qui peut être un capteur d'effort ou un capteur de couple. Ce type de capteur est souvent installé, ou peut être facilement installé dans les systèmes houlomoteurs, pour leur contrôle et/ou leur supervision. Alternativement, la mesure peut être 35 remplacée par une estimation réalisée à partir de la consigne de force (ou de couple) envoyée au PTO. 10 Pour l'exemple de système houlomoteur illustré sur la figure 1, un capteur de couple peut être disposé au niveau de l'articulation 4 (ou du générateur électrique 3). 3) Construction d'un modèle de force de radiation 5 Lors de cette étape, on construit un modèle de la force de radiation exercée sur le moyen mobile. Selon la théorie linéaire des vagues (comme décrit par exemple dans le document Faines J, Kurniawan A. "Fundamental formulae for wave-energy conversion". R. Soc. open sci. 2: 140305, 2005, http://dx.doi.org/10.1098/rsos.140305), la force de radiation résulte de l'oscillation d'un corps immergé (donc dépend du mouvement du moyen mobile), 10 tandis que la force d'excitation, résultant de la présence même d'un corps dans l'eau, ne dépend pas du mouvement du corps immergé, mais de la vague incidente. En l'absence de houle incidente, la force de radiation amortit l'oscillation résiduelle du corps immergé jusqu'à l'arrêter. Il est important de noter que, bien que la théorie linéaire permette de relier la force d'excitation à l'élévation de la vague incidente à travers un modèle linéaire (dans le domaine 15 fréquentiel ou temporel), en pratique on ne peut pas s'en servir pour calculer la force d'excitation en ligne, même si c'était possible de mesurer l'élévation de la vague au centre de gravité du moyen mobile comme le requiert la théorie. En effet, la relation linéaire entre élévation de la vague et force d'excitation est non-causale, ce qui veut dire qu'on ne peut pas calculer la force d'excitation à un instant donné sans connaître l'élévation de la vague 20 dans les instants futurs (le calcul peut en revanche être fait hors ligne, une fois la vague passée). Dans un contexte de contrôle en temps-réel, la force d'excitation ne peut donc être considérée que comme une force exogène totalement inconnue agissant sur le flotteur. En revanche, toujours selon la théorie linéaire des vagues, la force de radiation est reliée au mouvement du flotteur, et plus précisément à son accélération et à sa vitesse par un modèle 25 linéaire (dans le domaine fréquentiel ou temporel) causal. Elle peut donc être calculée en ligne en utilisant les mesures d'accélération et vitesse courantes (et passées pour la vitesse). Selon une mise en oeuvre de l'invention, on construit ledit modèle de la force Tract (t) de radiation par une équation du type : Trad(t) = -ij(t) - -cr(t) 30 avec tr(t) = f(t) h(t - s)«s)ds = h(t) * .(t) la composante de ladite force T - rad de radiation qui dépend de la vitesse (courante et passée) du moyen mobile, que l'on peut appeler amortissement de radiation, l'accélération angulaire du moyen mobile, 11 'co le moment d'inertie ajouté à la fréquence infiniment élevée, que l'on peut obtenir par des codes de calcul BEM (« Boundary Element Method » signifiant méthode des éléments de frontières), tels que WAMIT (WAMIT, USA), ou Nemoh (Ecole Centrale de Nantes, France), à partir de la géométrie du moyen mobile, 5 s la vitesse angulaire du moyen mobile, h la réponse impulsionnelle qui relie la vitesse du moyen mobile à l'amortissement de radiation, dont les coefficients s'obtiennent à partir des paramètres hydrodynamiques du moyen mobile calculés par les mêmes codes de calcul BEM. 10 La construction de ce modèle permet une détermination à tout instant de la force de radiation, avec un temps de calcul limité. Ainsi, la détermination de la force exercée par la houle peut être déterminée à tout instant avec un temps de calcul court. 4) Construction d'un modèle de la force de traînée 15 Lors de cette étape, on construit un modèle de force de traînée exercée sur le moyen mobile. La force de traînée correspond aux frottements visqueux sur le moyen mobile du système houlomoteur. Il s'agit d'un modèle non linéaire, pouvant dépendre de la vitesse du moyen mobile. La force de traînée (due à des frottements visqueux) est souvent considérée négligeable pour les systèmes houlomoteurs de type convertisseurs ponctuels oscillants et 20 est généralement exclue de leur modélisation. Ce n'est pas le cas, toutefois, pour les systèmes houlomoteurs de type à volet (« flaps »), ni pour des machines telles que des convertisseurs submergés à différentiel de pression de forme plate ni pour des machines à masse tournante tanguant en surface, entre autres. 25 Selon un mode de réalisation de l'invention, pour un système houlomoteur ayant un mouvement rotatif, le modèle de la force de traînée id peut s'écrire par une équation de la forme : td(t) = (OP (t) I avec fi le coefficient de traînée, ce coefficient peut être déterminé par des essais expérimentaux du système houlomoteur ou par des simulations numériques de type CFD 30 de l'anglais Computational Fluid Dynamics, correspondant à Mécanique des fluides numérique) et s la vitesse angulaire du moyen mobile. 12 Selon une variante de réalisation de l'invention, pour un système houlomoteur ayant un mouvement de translation, le modèle de la force de traînée Fd peut s'écrire par une équation de la forme : Fd(t) = le i(t)Ii(t) avec fi le coefficient de traînée, ce coefficient peut être déterminé par des essais 5 expérimentaux du système houlomoteur ou par des simulations numériques de type CFD de l'anglais Computational Fluid Dynamics, correspondant à Mécanique des fluides numérique) et 2 la vitesse du moyen mobile. 10 La construction de ce modèle permet une détermination précise de la force de traînée. Ainsi la détermination de la force exercée par la houle sur le moyen mobile peut être précise, car elle ne néglige pas la force de traînée. 5) Construction du modèle de dynamique du système houlomoteur 15 Lors de cette étape, on construit un modèle de dynamique du système houlomoteur. On appelle modèle de dynamique un modèle qui relie la force d'excitation exercée par la houle incidente sur le moyen mobile, la force de radiation exercée sur le moyen mobile, la force de rappel hydrostatique exercée sur le moyen mobile, la force exercée par la machine de conversion sur le moyen mobile, à la position et la vitesse du moyen mobile. Ce type de 20 modèle permet d'obtenir des résultats représentatifs du comportement du système houlomoteur, si les mouvements ne sont pas trop grands. Avantageusement, le modèle de dynamique est obtenu par application du principe fondamental de la dynamique (deuxième loi de Newton) au niveau du moyen mobile. 25 Selon une mise en oeuvre de l'invention, pour lequel on considère les efforts, on peut construire le modèle de dynamique du système houlomoteur par une équation du type : /eqS(t) = thd(t)+-c'd(t)+-cd(t) + tw(t) - u(t) avec /eq le moment d'inertie du moyen mobile, l'accélération angulaire du moyen mobile, 30 -cw le couple d'excitation exercé par la houle incidente sur le moyen mobile, tradle couple de radiation exercé sur le moyen mobile, Thd le couple de rappel hydrostatique exercé sur le moyen mobile, rd le couple de traînée exercé sur le moyen mobile, et u le couple exercé par la machine de conversion sur le moyen mobile. 13 Ce modèle traduit un mouvement de rotation autour d'un axe horizontal (typique pour le système houlomoteur de la figure 1). Ce modèle est dérivé de la théorie linéaire des vagues. 5 Conformément à une première variante de réalisation, la force de rappel hydrostatique exercée sur le moyen mobile peut être approximée par une fonction linéaire de la position z définie par rapport à la position d'équilibre. Dans ce cas, la force de rappel hydrostatique peut s'écrire par une fonction du type : xhd(t) = -K8(t) avec ô la position angulaire du moyen mobile définie par rapport à sa 10 position d'équilibre et K le coefficient de raideur hydrostatique. Ainsi, la force de rappel hydrostatique peut être calculée à partir d'un modèle simple si l'on dispose de la mesure de la position ô. Cette fonction est particulièrement adaptée pour les petits déplacements ô. Selon une mise en oeuvre de l'invention, pour lequel on considère les efforts, on peut 15 construire le modèle de dynamique du système houlomoteur par une équation du type : M 2(t) = F'(t) + Fhd(t)+F'd(t) + F d(t) - Fu(t) avec M la masse du moyen mobile, 2 l'accélération du moyen mobile, F' la force d'excitation exercée par la houle incidente sur le moyen mobile, 20 Frad la force de radiation exercée sur le moyen mobile, Fhd la force de rappel hydrostatique exercée sur le moyen mobile, Fd la force de traînée exercée sur le moyen mobile, et Fu la force exercée par la machine de conversion sur le moyen mobile. Ce modèle traduit un mouvement de translation vertical (typique des flotteurs ayant un 25 mouvement de pillonement). Ce modèle est dérivé de la théorie linéaire des vagues. Ce modèle est le « miroir » du modèle pour lequel on considère les couples ; les différents termes du modèle sont de même nature. Selon un mode de réalisation de l'invention, la force de rappel hydrostatique et/ou le 30 couple de rappel hydrostatique peuvent être approximés par une fonction linéaire, ou par une fonction affine par morceaux. 6) Détermination de la force d'excitation exercée par la houle incidente Lors de cette étape, on détermine, en temps réel, la force d'excitation exercée par la 35 houle incidente sur le moyen mobile au moyen : 14 - des mesures de position et de vitesse (et éventuellement d'accélération) du moyen mobile déterminées à l'étape 1), - de la mesure de la force exercée par la machine de conversion PTO sur le moyen mobile déterminée à l'étape 2), 5 - du modèle de la force de radiation déterminé à l'étape 3), - du modèle de la force de traînée déterminé à l'étape 4), et - du modèle de dynamique du système houlomoteur déterminé à l'étape 5). Selon l'invention, on détermine la force d'excitation exercée par la houle incidente sur 10 le moyen mobile au moyen d'un observateur d'état basé sur un filtre de Kalman sans parfum (UKF pour « Unscented Kalman Filter ») construit à partir d'un modèle de marche aléatoire de la force d'excitation exercée par la houle incidente sur le moyen mobile. Le filtre de Kalman sans parfum permet de prendre en compte les non linéarités des modèles, en particulier du modèle de la force de traînée. 15 On rappelle qu'un observateur d'état, ou un estimateur d'état, est, en automatique et en théorie des systèmes, une extension d'un modèle représenté sous forme de représentation d'état. Lorsque l'état du système n'est pas mesurable, on construit un observateur qui permet de reconstruire l'état à partir d'un modèle. Un filtre UKF est basé sur la théorie de la transformation « non parfumée », qui permet 20 d'obtenir un estimateur pour un système non linéaire sans avoir besoin de le linéariser au préalable pour l'application du filtre. Le filtre UKF utilise une distribution statistique de l'état qui est propagée à travers les équations non linéaires. Un tel filtre présente l'avantage de permettre une stabilité, donc une robustesse de l'estimation. 25 Pour ce mode de réalisation, préalablement à cette étape, on connaît donc : - les mesures de position ô, et de vitesses du moyen mobile, - la mesure de la force u exercée par la machine de conversion PTO sur le moyen mobile, - le modèle de la force de radiationTrac(t) = - __ /' (t) - Tr(t) avec 30 'Fr (t) = fot h(t - s)«s)ds = h(t) * «t), - le modèle de la force de traînée Td(t) = p«op(01, et - le modèle de dynamique du système houlomoteur /eqe(t) = T hd(t)+Trad(t)+T d (t) + Tw(t) - U(t) . 35 Dans cette approche, le problème d'estimation de la force d'excitation de la houle est transformé en un problème classique d'estimation d'état (qui peut être résolu avec un filtre 15 de Kalman sans parfum), en exprimant la dynamique de la force d'excitation de la houle par un modèle de marche aléatoire. L'avantage principal de cette méthode est la prise en compte d'incertitudes permettant de prendre en compte les bruits de mesure et les dynamiques non modélisées. 5 En remplaçant dans l'équation qui décrit le mouvement du moyen mobile /eq SM = Thd (0+-c'd(t) +-cd (t) + -cw(t) - u(t) les expressions pour la force de rappel hydrostatique (avec le modèle linéaire), le modèle de la force de radiation, et le modèle de la force de traînée T'id (t) = -K8(t) Trad (t) = -Ico«t) - Tr(t)Td (0 = p«op(t)1 on obtient le modèle non linéaire suivant : 10 (I eq + I co)«t) + fl «OP MI + Ka (t) = Tw(t)+Tr(t) - u(t). Cette équation peut être mise sous forme d'état, en définissant rx1(t) = «0 lx2(t) = .(t) ce qui permet d'écrire le modèle précédent sous forme de représentation d'état r(t) = x2 (t) 1 ±i(t) =/eq + Lm ( lai(t) - fix2(01x2(01 - Tr (0 + Tw(t) - U(0) Ce système d'équation contient le terme intégral t Tr(t) = 1h(t - s)«s)ds = h(t) * «t) 0 que l'on peut considérer comme un système linéaire dont l'entrée est .(t) et la sortie 15 Tr(t). Avec la méthode de Prony, ce système peut alors être transformé en la représentation d'état équivalente : {±r(t) = Arxr(t) + BAt) Tr(t) = Crxr(t) + Dr «t) OU P.Cr (0 = ArXr(t) + BrX2 (t) (Tr(t) = Crxr(t) + DrX2 (t) où xr est un état interne (inaccessible) n'ayant pas de signification physique particulière et (Ar, Br, Cr, Dr) sont les matrices de la réalisation d'état. 20 En combinant les deux représentations d'état, on obtient 16 {±1(0 = x2(0 1 x1(t) , ( Kx1(t) - Drx2(t) - Crxr(t) - flx2(t)lx2 MI + Tw(t) - u(t)) led + I,0 jcr(t) = Brx2(t) + Arxr(t) ou, de manière équivalente {r±(t) = Acx(t) + fc(x(0) + Bc(rw(t) - u(0) Y(t) = Ccx(t) où x(t) = [x1(t) x2(t) xi:(01T [x1(01 l_x2 (01 et 0 1 0 0 K Dr Cr fi Ac =[ led + 40 led + 1,0 led + 1,01' fc(x(t)) = leq + Lmx2(01x2(01 [ 0 Br Ar 0 0 11 0 01 [ i Bc -[led + 1,0 ' Cc =0 1 0 0 A cause du terme fc(x(t)), le système est non linéaire. Le système a deux entrées, 5 Tw(t) et u(t), et deux sorties x1(t) = 8(0 et x2(t) = .(t). L'entrée Tw(t) n'est pas mesurable et est inconnue. Le problème à résoudre est de l'estimer à partir des quantités mesurées u(t) (moment appliqué par la machine de conversion PTO), 8(0 et .(t) (position et vitesse angulaire de la partie mobile). Si le mouvement principal du système houlomoteur était un mouvement de translation 10 et le problème était d'estimer la force d'excitation de la vague, les mêmes développements s'appliquent en remplaçant position et vitesse angulaires par des positions et vitesses et tous les moments par des forces. Pour réaliser l'estimation, dans un premier temps on discrétise le système d'état ci- 15 dessus (car les mesures sont échantillonnées et l'algorithme d'estimation est exécuté par un calculateur), avec la méthode d'Euler. Ce qui donne, pour une période d'échantillonnage donnée T., : fx(k + 1) = Adx(k) + fd(x(k)) + Bdtw(k) - Bdu(k) (Y(k) = Cdx(k) où Ad = I + Ts, fd(X(k)) = Tsfc(x(k)), Bd = TsBc, Cd = Cc et I est la matrice identité de dimensions appropriées. 20 Pour estimer le moment d'excitation Tw(k), on la considère comme un état, en introduisant un modèle mathématique qui relie Tw(k) et Tw(k + 1), en l'occurrence Y= 17 rw(k + 1) = rw(k) + Em(k) OU Em(k) décrit la variation de rw(k) et est considéré comme un nombre aléatoire. En d'autres termes, ce modèle suppose qu'à chaque instant k, le moment d'excitation s'éloigne d'un pas (d'une quantité) aléatoire de sa valeur précédente, et que ces pas sont distribués indépendamment et identiquement en taille. 5 On couple ce modèle de marche aléatoire de la force d'excitation à un modèle plus réaliste de la dynamique de la partie mobile f x(k + 1) = Adx(k) + fd(x(k)) + Bdrw(k) - Bdu(k) + Ex(k) y(k) = Cdx(k) + v(k) où Ex(k) représente les dynamiques non modélisées (frottements du PTO, non-linéarité hydrostatique, ...) et v(k) décrit les bruits de mesure sur la position et la vitesse du 10 flotteur. En combinant le modèle de marche aléatoire de la force d'excitation et dynamique non linéaire de la partie mobile, on obtient le système augmenté x(k + 1) = Adx(k)+fd(x(k)) + Bdrw(k) - Bdu(k) + Ex(k) { w(k + 1) = rw(k) + Em(k) y(k) = Cdx(k) + v(k) ou, de manière équivalente, 15 (xa(k + 1) = Aaxa(k)+fa(xa(k)) + Bau(k) + E(k) (y(k) = Caxa(k) + v(k) OÙ x(k) 1 , , [fd(x0 (k))1 E(k) - [Em(k)] xa(k) - [ rw(k) , fayca(k)) = et Aa=[Aod Bidl,Ba=[-oBal,Ca= [Cd 0] Ainsi, le problème d'estimer rw(k), devient un problème d'estimation d'état. 20 Une manière d'estimer le vecteur d'état inconnu xa(k), en prenant en compte des informations sur E(k) et v(k), est d'appliquer l'algorithme du filtre de Kalman (KF). Dans le procédé proposé, le filtre de Kalman sans parfum (UKF) est utilisé afin de traiter la non linéarité du système. L'UKF est généralement plus robuste et précis que le filtre de Kalman 25 étendu (EKF), qui traite la non linéarité en la linéarisant. De plus, dans ce cas de figure, la présence du terme de traînée, dont la dérivée n'est pas continue, rend l'EKF inapplicable. Comme l'EKF, l'UKF réalise l'estimation à travers deux étapes, prédiction d'état et correction des mesures, excepté que ces deux étapes sont précédées d'une étape préalable pour le calcul des « points sigma ». Les points sigma sont un ensemble d'échantillons 18 calculés de manière à pouvoir propager de manière exacte les informations de moyenne et de variance dans l'espace d'une fonction non linéaire. Conformément à une mise en oeuvre de l'invention, on peut adopter les hypothèses 5 suivantes : l'état initial x,(0) est un vecteur aléatoire de moyenne m(0) = E[x,(0)] et covariance P(0) = E[(xa(0) - m(0))(xa(0) - m(0))T] ; E(k) et v(k) sont des bruits gaussiens avec matrices de covariance Q et R, respectivement ; 10 et les notations suivantes : Xa(klk - 1) est l'estimation de xa(k) à partir des mesures jusqu'au temps k - 1, soit y(k - 1), y(k - 2), ... et u(k - 1), u(k - 2), ... .5ea(k lk) est l'estimation de xa(k) à partir des mesures jusqu'au temps k, soit y(k), y(k - 1), ... et u(k), u(k - 1), ... 15 Px(kik - 1) est la matrice de covariance de xa(k) à partir des mesures jusqu'au temps k - 1, soit y(k - 1), y(k - 2), ... et u(k - 1), u(k - 2), ... Px(kik) est la matrice de covariance de x(k) à partir des mesures jusqu'au temps k, soit y(k), y(k - 1), ... et u(k), u(k - 1), ... 20 Pour cette mise en oeuvre, les trois étapes de la méthode UKF peuvent être : 1. Calcul des points sigma Soient m - = w n-FÂ: vvoc = + (1 - a2 + y), 25 _ wc _ Â- - 2(+À.)' 1,2, ..., 2n où À = (a2 - 1)n est un paramètre de mise à l'échelle, n est la dimension de l'état xa(k), a est un paramètre qui détermine la dispersion des points sigma autour de x(k - 11k - 1) et auquel on donne généralement une valeur positive petite, par exemple 10-3, y est un paramètre utilisé pour incorporer la connaissance a priori de la distribution de 30 x : pour une distribution gaussienne, y = 2 est optimal. Au temps k - 1, on considère le choix suivant de points sigma (ensemble de points qui encode de manière exacte les informations de moyenne et variance) : xo(k - 1) = 2,(k - 11k - 1), xi(k - 1) = fca(k - 11k - 1) + À Si(k - 1), i = 1,2, ... , n xi+n(k - 1) = xa(k - 11k - 1) -n'N/E À S i(k - 1), i = 1,2, ... , n 19 où Si(k - 1) est la i-ème colonne de la racine carrée matricielle de Px(k - 11k - 1), soit Px(klk - 1) = S(k - 1)T S(k -1) 2. Mise à jour des prédictions 5 Chaque point sigma est propagé à travers le modèle non linéaire qui représente l'évolution des états : Xj(klk - 1) = Aaxi(k - 1) + fa (2j(k - 1)) + Bau(k - 1), j = 0,1, , 2n La moyenne et la covariance de 2a(klk - 1), la prédiction de xa(klk - 1), sont calculées comme X(kIk - 1) = VITiniXj(klk - 1), 10 Px(kik -1) = e2'0Wic - 1) - xa(klk - 1)) (2 j(klk - 1) - xa(klk - 1))T +Q. Les états prédits Xj(klk - 1) sont utilisés dans l'équation d'état de sortie, ce qui donne : MkIk - 1) = Ca21(klk - 1) La moyenne et la covariance de 3î(kik - 1) sont calculées comme 15 3î(Iclic - 1) _ lno 2, WirnMkik - 1), Py(klk - 1) = Ej2"2,3VVjc (3îj(klk - 1) - y(klk - 1)) (Mkik -1) - y(klk -1)) +R tandis que la covariance croisée entre 2a(klk - 1) et j(kik - 1) est : 2n Pxy(klk -1) _I (2j(kik - 1) - xa(klk - 1)) (j77/(kik - 1) - y(klk - 1)) T i=0 3. Mise à jour à partir des mesures 20 Comme dans le filtre de Kalman linéaire, l'estimation finale de l'état est obtenue en corrigeant la prédiction avec une rétroaction sur l'erreur de prédiction sur la sortie (mesurée) : 2(k) = 2a(klk - 1) + K(9(k) - 3î(kik - 1)) où le gain K est donné par K = Pxy(klk -1)Py(klk - 1)-1 25 La covariance a posteriori de l'estimation est mise à jour avec la formule suivante : Px(k1k) = Px(klk - 1) - KPy(klk - 1)KT Les règles suivantes sont utilisées pour choisir les matrices P0 (covariance de l'état initial) et R : - Si l'état initial .2a(k) au temps k = 0 est bien connu, c.-à-d. m(0) xa(0), alors Pc71 est grand ; 20 - S'il y a beaucoup de bruit dans les mesures y(k), alors R est petit. Il est plus complexe de choisir Q. On le prend généralement sous forme diagonale, comme il suit 0 i Qrni avec Qm » Qx. 5 Le procédé basé sur l'UKF et sur le modèle de marche aléatoire de la force de la vague selon cette mise en oeuvre peut être résumé comme il suit : - On initialise k = 0, le vecteur d'état 2,(010)=m(0) et l'état de la matrice de covariance, P(010) = P0 10 - À chaque instant k: on utilise : - les mesures de la position et de la vitesse du moyen mobile y(k) = [8(k) ((k)]et de la force exercée par le PTO sur le moyen mobile u(k) - les résultats des estimations du pas précédent 2a(k - 1Ik - 1), 15 P(k - 1Ik - 1) - les paramètres Q, R (matrices de covariance) on détermine la force d'excitation Tw exercée par la houle incidente sur le moyen mobile, noté pour ce mode de réalisation î-w(k), en mettant en oeuvre les étapes suivantes : 20 - on applique les trois étapes de l'algorithme du filtre de Kalman sans parfum pour obtenir 2a(klk), P(klk) décrites ci-dessus, ainsi on estime l'état complet avec sa matrice de covariance : 2(k) = 2a(klk - 1) + K(1'(k) -3î(k lk - 1)) K = Pxy(klk - 1)Py(klk - 1)-1 Px(k1k) = Px(klk - 1) - KPy(klk - 1)KT avec xa(k) le vecteur d'état inconnu xa(k) =[Tw((kk) x )1, Aa,Ba,Ca,Da des matrices de la réalisation d'état, P la matrice de covariance du vecteur 25 d'état, Q et R des matrices de calibration. - et on extrait la dernière composante du vecteur d'état estimé, c'est-à-dire la force d'excitation Tw exercée par la houle incidente sur le moyen mobile par une équation de la forme : 0(k) = [0 1] 2a(klk) 30 21 7) Commande du système houlomoteur Il s'agit d'une étape facultative. Lors de cette étape, on contrôle le système houlomoteur, de manière à prendre en compte la force d'excitation exercée par la houle incidente. Ainsi, il est possible de piloter le système houlomoteur, afin d'optimiser l'énergie 5 récupérée. La commande peut consister en un contrôle du moyen mobile du système houlomoteur, par exemple au moyen d'une machine de conversion électrique, pneumatique ou hydraulique, nommé système PTO (de l'anglais « power take-off »). Ce système PTO influence le mouvement du moyen mobile et permet de transférer l'énergie mécanique au 10 réseau électrique, pneumatique ou hydraulique. La commande prédictive par modèle (MPC) est un exemple de méthode de commande de systèmes houlomoteurs. Exemples comparatifs Les caractéristiques et avantages du procédé selon l'invention apparaîtront plus 15 clairement à la lecture des exemples comparatifs ci-après. Pour les exemples, on considère un système houlomoteur à volet, tel qu'illustré en figure 1, à petite échelle pour des tests en bassin. Les paramètres hydrodynamiques du système houlomoteur sont décrits dans le 20 tableau 1. Tableau 1 - paramètres hydrodynamiques du système houlomoteur Paramètres hydrodynamiques Moment d'inertie total /eq + Lm 0,2753 kg-m2 Coefficient de raideur hydrostatique K 118,25 N-m-rac11 Coefficient de traînée ig 2 N-m-racrls2 De plus, pour la représentation d'état du système houlomoteur, on utilise les matrices 25 suivantes : Ar =[ 1 -5,1899 -36,074341 Br = [11 Cr = [13,2612 0],D, = 0 ' Dans un premier temps, pour ce système houlomoteur on détermine la force d'excitation exercée par la houle sur le moyen mobile au moyen du procédé selon l'art antérieur décrit dans la demande de brevet FR 3049989, au moyen d'un filtre de Kalman 30 linaire après linéarisation du terme de traînée. La figure 2 illustre les courbes du couple 22 d'excitation Tw en Nm exercé par la houle sur le moyen mobile en fonction du temps T en s, pour une vague régulière. La courbe REF en trait continu correspond au couple de référence, réellement exercé, et la courbe EST en traits pointillés correspond au couple estimé par ce procédé de l'art antérieur. On remarque que le procédé selon l'art antérieur 5 avec le filtre de Kalman linéaire donne des résultats imprécis, avec une imprécision qui augmente avec l'amplitude de la houle. Dans un deuxième temps, pour ce système houlomoteur, on détermine la force d'excitation exercée par la houle sur le moyen mobile au moyen du procédé selon l'invention 10 (au moyen d'un filtre Kalman sans parfum). La figure 3 illustre les courbes du couple d'excitation Tw en Nm exercé par la houle sur le moyen mobile en fonction du temps T en s, pour une vague régulière. La courbe REF en trait continu correspond au couple de référence, réellement exercé, et la courbe EST en traits pointillés correspond au couple estimé par ce procédé selon l'invention. Les courbes REF et EST sont superposées, ce qui 15 montre que le procédé selon l'invention permet une détermination précise du couple d'excitation exercée par la houle sur le moyen mobile. Dans un troisième temps, on répète le deuxième test de la figure 2 pour une houle irrégulière. La figure 4 illustre les courbes du couple d'excitation Tw en Nm exercé par la 20 houle sur le moyen mobile en fonction du temps T en s, pour une vague irrégulière. La courbe REF en trait continu correspond au couple de référence, réellement exercé, et la courbe EST en traits pointillés correspond au couple estimé par ce procédé selon l'invention. Les courbes REF et EST sont superposées, ce qui montre que le procédé selon l'invention permet une détermination précise du couple d'excitation exercée par la houle sur le moyen 25 mobile, même pour une houle irrégulière. Le procédé selon l'invention est donc bien adapté pour déterminer de manière précise la force d'excitation exercée par la houle sur le moyen mobile, en particulier d'un REVENDICATIONS1) Procédé de détermination d'une force d'excitation exercée par la houle sur un moyen mobile (2) d'un système houlomoteur (1), ledit système houlomoteur (1) convertissant l'énergie de la houle en énergie au moyen dudit moyen mobile (2) coopérant avec une machine de conversion (3), ledit moyen mobile (2) effectuant un mouvement par rapport à ladite machine de conversion (3) sous l'action de la houle, caractérisé en ce qu'on réalise les étapes suivantes : a) on mesure la position et la vitesse dudit moyen mobile (2) ; b) on mesure la force u exercée par ladite machine de conversion (3) sur ledit moyen mobile (2) ; c) on construit un modèle de la force T -rad de radiation exercée sur ledit moyen mobile (2), ledit modèle de la force de radiation reliant ladite force de radiation Trad à la vitesse dudit moyen mobile (2) ; d) on construit un modèle de force de traînée id exercée sur ledit moyen mobile (2), ledit modèle de force de traînée reliant ladite force de traînée id à la vitesse dudit moyen mobile (2) ; e) on construit un modèle de dynamique dudit système houlomoteur qui relie ladite force d'excitation -cw exercée par la houle incidente sur ledit moyen mobile (2), à ladite position dudit moyen mobile (2), à ladite vitesse du moyen mobile (2), à ladite force u exercée par ladite machine de conversion (3) sur ledit moyen mobile (2), à ladite force de radiation T -rad exercée sur ledit moyen mobile (2) et à ladite force de traînée id exercée sur ledit moyen mobile (2) ; et f) on détermine ladite force d'excitation -cw exercée par la houle incidente sur ledit moyen mobile (2) au moyen dudit modèle de dynamique, dudit modèle de la force de radiation, dudit modèle de la force de traînée id, desdites position et vitesse mesurées et de ladite force u mesurée exercée par ladite machine de conversion (3) sur ledit moyen mobile (2) et au moyen d'un filtre de Kalman sans parfum construit à partir d'un modèle de marche aléatoire de ladite force d'excitation exercée par la houle incidente sur ledit moyen mobile (2). 2) Procédé selon la revendication 1, dans lequel on construit ledit modèle de dynamique dudit système houlomoteur au moyen d'une équation du type : /eq «t) = Tha(t)+Trad(t)+Ta(t) + Tw(t) - u(t) avec /eq le moment d'inertie total dudit moyen mobile (2), 8(0 l'angle de rotation dudit moyen mobile (2) par rapport à la position d'équilibre, avec ..(t) l'accélération angulaire dudit moyen mobile (2) et .(t) la 24 vitesse angulaire dudit moyen mobile (2), Thd(t) le moment de rappel hydrostatique, Trad(t) le moment de radiation, Td(t) le moment de trainée, Tw(t) le moment d'excitation de la vague, u(t) le moment exercé par ladite machine de conversion (3) sur ledit moyen mobile (2). 5 3) Procédé selon la revendication 2, dans lequel on détermine le moment de rappel hydrostatique Thd(t) au moyen d'une formule du type : Thd(t) = -K8(t) où K est le coefficient de raideur hydrostatique. 10 4) Procédé selon l'une des revendications précédentes, dans lequel on construit ledit modèle de force de radiation par une équation du type : Trad(t) = -100S(0 - 'rr(t) avec I, est le moment d'inertie ajouté à fréquence infiniment élevée et Tr(t) = fot h(t - s)«s)ds = h(t) * .(t), avec h la réponse impulsionnelle qui relie la 15 vitesse du moyen mobile à l'amortissement de radiation et .(t) la vitesse angulaire dudit moyen mobile (2). 5) Procédé selon l'une des revendications précédentes, dans lequel on construit ledit modèle de force de trainée par une équation du type : Td(o= p«op(t)1 20 où fi est le coefficient de trainée et .(t) la vitesse angulaire dudit moyen mobile (2). 6) Procédé selon l'une des revendications 1 à 4, dans lequel ledit modèle de la force de trainée Fd peut s'écrire par une équation de la forme : Fd(t) = fl2(t)12(t) avec fi le coefficient de trainée, 25 et 2 la vitesse du moyen mobile. 7) Procédé selon l'une des revendications précédentes, dans lequel on détermine ladite force d'excitation -cw exercée par la houle incidente sur ledit moyen mobile (2) en mettant en oeuvre les étapes suivantes : 25 i) on initialise k=0, le vecteur d'état 2,(010)=m(0) et l'état de la matrice de covariance, P(010) = P0 ; ii) à chaque instant k, on acquiert lesdites mesures de position et de vitesse dudit moyen mobile y(k) = [8(k) ((k)] et la mesure de ladite force exercée par 5 ladite machine de conversion sur ledit moyen mobile u(k) ; et iii) à chaque instant k, on détermine la force d'excitation exercée par la houle incidente sur ledit moyen mobile (2) tw(k) au moyen des équations suivantes : îw(k) = [0 1]2'(k1 k) Px(kik) = Px(klk - 1) - K Py(klk - 1)KT K = Pxy(klk - 1)Py(klk - 1)-1 10 2(k) = 2a(klk - 1) + K (j(k) - 3î(klk - 1)) 2n Pxy(klk - 1) =IVVic (klk - 1) - xa(klk - 1)) (j7Ti (kik - 1) - y (kik - 1)) T =0 jkklk - 1) = e1,31/11jrn99 (kik - 1), Py(klk - 1) _ Elno 1/5 (3i (kik - 1) - y(klic - 1)) (3î (kik - 1) - y (kik - 1)) +R MkIk - 1) = Ca2 (kik - 1) 2(kIk - 1) = 1'2'0Wri 2i (kik - 1), 15 Px(kik - 1) = e2'0Wic (2i (kik - 1) - xa(klk - 1)) (2 j(k lk - 1) - xa(klk - 1))T +Q. j(k lk - 1) = A axi (k - 1) + fa (k - 1)) + Bau(k - 1), j = 0,1, ... , 2n xo(k - 1) = £a (k - 11k - 1), xi(k - 1) = 2a(k - 11k - 1) + N./E À Si(k - 1), i = 1,2, ..., n xi+n(k - 1) = 2 a(k - 11k - 1) -nE À Si (k - 1), i = 1,2, ..., n où Si(k - 1) est la i-ème colonne de la racine carrée matricielle de Px(k - 1 ik - 1), soit Px(klk - 1) = S(k - 1)T S(k - 1) xk))15 A 20 xa(k) le vecteur d'état inconnu xa(k) = [ Tw( des matrices de la (kBcoCa,Da réalisation d'état, P la matrice de covariance du vecteur d'état, Q et R des matrices de calibration. 8) Procédé selon l'une des revendications précédentes, dans lequel on commande ledit 25 système houlomoteur en fonction de ladite force déterminée d'excitation Ti,/ exercée par la houle incidente sur ledit moyen mobile (2). 26 9) Procédé selon l'une des revendications précédentes, dans lequel ledit système houlomoteur est un système houlomoteur de type à convertisseur submergé à différentiel de pression de forme plate, de type machine à masse tournante tanguant en surface ou du type système houlomoteur à volet. 10) Système houlomoteur caractérisé en ce qu'il comprend des moyens de commande, notamment des moyens de commande informatique, mettant en oeuvre un procédé selon l'une des revendications 1 à 9. 5 10