Brevet FR2950796A1 — publié le 20110408

CARACTERISATION DES CYCLES RESPIRATOIRES. La présente invention concerne un procédé de détection des caractéristiques des cycles respiratoires d'un patient portant un masque d'assistance respiratoire. Elle concerne également un produit programme d'ordinateur mettant en oeuvre le procédé et un appareil d'assistance respiratoire. Dans le cadre de l'assistance respiratoire à des malades souffrant d'insuffisance respiratoire, il est important de pouvoir surveiller les cycles de respiration de façon à détecter les anomalies apparaissant dans les interactions patient-ventilateur, les apnées et autres anomalies et, ainsi, permettre une assistance respiratoire seulement quand cela est nécessaire et avec les bons paramètres. La possibilité de la détection de cycles respiratoires à partir de mesures non-invasives remonte aux travaux de Giannouli, K. Webster, D. Roberts & M. Younes « Response of ventilator-dependent patients to different levels of pressure support and proportional assist », American Journal of Respiratory and Critical Care Medicine, 159, 1716-1725, 1999. La première utilisation d'une procédure automatique remonte aux travaux de L. Achour, C. Letellier, A. Cuvelier, E. Verin & J.-F. Muir, « Asynchrony and cyclic variability in pressure support noninvasive ventilation », Computer in Biology and Medicine, 37, 1308-1320, 2007. Ces travaux ne concernaient que l'identification des efforts inefficaces et la procédure automatique d'identification n'avait été appliquée qu'à des séances de ventilation proposées à des patients éveillés. L'utilisation de la procédure décrite dans ce dernier document à des séances de ventilation nocturnes s'est rapidement révélée un échec : jusqu'à 50% des cycles n'étaient pas correctement identifiés et les apnées du sommeil, souvent fréquentes chez les insuffisants respiratoires chroniques, introduisaient de nombreux artefacts. Par ailleurs, certains ventilateurs offrent un bilan global des cycles non-déclenchés mais aucun ne proposent un bilan détaillé des autres types d'asynchronismes (déphasage patient-ventilateur, double déclenchement, etc.). Or l'analyse de la qualité mécanique de la ventilation non invasive peut constituer une aide au diagnostic médical et guider le pneumologue dans le réglage du ventilateur. Il apparait donc un besoin d'identification des cycles respiratoires qui soit suffisamment fiable pour des séances de ventilation nocturne, et permettant une identification en temps réel des asynchronismes, que ce soit sur des patients éveillés ou en état de sommeil. Pour résoudre un ou plusieurs des inconvénients cités précédemment, dans un aspect de l'invention, un procédé de détection de caractéristiques des cycles respiratoires d'un patient portant un masque d'assistance respiratoire, comporte : • acquisition de données numériques de débit respiratoire échantillonnées selon une fréquence d'acquisition prédéterminée, définissant ainsi une courbe de débit fonction du temps ; • définition d'une taille d'une fenêtre temporelle glissante de débit ; • détermination d'une pseudo-dérivée par : • pour chaque fenêtre glissante : • recherche du maximum et du minimum de débit à l'intérieur de la fenêtre glissante ; • détermination de la pseudo-dérivée à partir du maximum et du minimum trouvés, la pseudo-dérivée étant positive si le maximum arrive après le minimum, et étant négative autrement ; • déplacement de la fenêtre glissante vers l'échantillon de débit suivant ; • détection d'un plateau de débit lorsque la pseudo-dérivée est, en valeur absolue, inférieure à une valeur prédéterminée. Des caractéristiques ou des modes de réalisation particuliers sont : • la taille de la fenêtre temporelle glissante est déterminée en fonction d'une durée moyenne de cycle respiratoire et de la fréquence d'acquisition ; • la durée moyenne de cycle respiratoire est calculée à partir d'un maximum de densité spectrale extrait d'un périodogramme du débit généré par une méthode de Welch ; • le procédé comporte en outre : • recherche des points de croisement de la pseudo-dérivée correspondant au passage d'une valeur négative à une valeur positive de la pseudo-dérivée pour une première catégorie de points de croisement et au passage d'une valeur positive à une valeur négative de la pseudo-dérivée pour une seconde catégorie de points de croisement ; • détermination de la durée d'un cycle respiratoire sur la base de la durée entre deux points de croisement de même catégorie consécutifs. • il comporte en outre : • détermination de la dérivée du débit ; • recherche de la position du maximum de débit dans l'intervalle temporelle compris entre un point de croisement de seconde catégorie et le point de croisement de première catégorie suivant ; • définition du point milieu entre le point de croisement de première catégorie précédant ledit point de croisement de seconde catégorie et le maximum ; • détermination du maximum de la dérivée du débit et de sa position dans l'intervalle compris entre le point milieu et ledit maximum de débit ; • détermination du début du cycle respiratoire comme : • la position du maximum de la dérivée du débit si la valeur de celui-ci est inférieure à une valeur prédéterminée, • la position du point le plus tardif, et précédant la position du maximum de la dérivée, de la dérivée du débit inférieure à ladite valeur prédéterminée si la valeur de celui-ci est supérieure à ladite valeur prédéterminée. • la valeur prédéterminée est égale à 0,4 I/s2. 30 • la détermination de la dérivée du débit comporte : • détermination d'une fréquence de coupure d'un filtre passe-bas du débit, la fréquence de coupure étant déterminée comme le deuxième minimum de la réponse en fréquence d'un filtre à réponse impulsionnelle finie à coefficients Ci, lesdits coefficients correspondant à la première colonne d'une matrice A définie comme A= [VTVr1VT où v est une matrice à deux colonnes et 2M + 1 lignes, 2M+1 étant le nombre de points d'une fenêtre glissante sur la courbe de débit, VT étant la transposée de la ùM 1` ù(M -1) 1 matrice V et V = • détermination de la dérivée au centre de ladite fenêtre glissante comme la somme de la valeur du débit en chaque point de la 2M+1 fenêtre glissante pondérée par les coefficients Ci, soit 1CiXi ; i=1 répétition de la détermination de la dérivée au centre de ladite fenêtre 15 glissante en décalant d'un point ladite fenêtre glissante sur la courbe de débit jusqu'au dernier point de la courbe. Le procédé utilise avantageusement des opérations de traitement du signal relativement sophistiquées pour extraire des régularités et des irrégularités significatives d'un signal de débit très bruité par nature. 20 Le procédé ainsi décrit permet avantageusement d'obtenir une détection automatique de la plupart des asynchronismes courants, en temps réel et sans ajout de capteur additionnel. De plus ce procédé est fiable quel que soit l'état du patient (éveillé, endormi, etc.). En effet, il permet avantageusement la détection de chaque cycle 25 respiratoire ainsi que d'évènements de ventilation et/ou de respiration et cette détection est fiable aussi bien sur des sujets éveillés que sur des sujets endormis. Il se base avantageusement sur des mesures non invasives de débit. Et l'analyse cycle-à-cycle de l'occurrence des évènements permet avantageusement d'étudier de possibles causalités et fournit ainsi avantageusement une aide au diagnostic médical. Selon un deuxième aspect de l'invention, un produit programme d'ordinateur comprend des instructions de code de programme enregistrées sur un support lisible par un ordinateur, pour mettre en oeuvre les étapes du procédé lorsque le programme fonctionne sur un ordinateur. Selon un troisième aspect de l'invention, un appareil d'assistance respiratoire comprend : • un capteur de débit respiratoire fournissant des données numériques de débit respiratoire échantillonnées selon une fréquence d'acquisition prédéterminée, définissant ainsi une courbe de débit fonction du temps ; ledit capteur fournissant lesdites données à • un calculateur définissant une taille d'une fenêtre temporelle glissante de débit et déterminant une pseudo-dérivée par : • pour chaque fenêtre glissante : • recherche du maximum et du minimum de débit à l'intérieur de la fenêtre glissante ; • détermination de la pseudo-dérivée à partir du maximum et du minimum trouvés, la pseudo-dérivée étant positive si le maximum arrive après le minimum, et étant négative autrement ; • déplacement de la fenêtre glissante vers l'échantillon de débit suivant ; le calculateur détectant un plateau de débit lorsque la pseudo- dérivée est, en valeur absolue, inférieure à une valeur prédéterminée. L'invention sera mieux comprise à la lecture de la description qui suit, faite uniquement à titre d'exemple, et en référence aux figures en annexe dans lesquelles : - la figure 1 est une vue schématique d'un appareil respiratoire selon un mode de réalisation de l'invention ; - la figure 2 est un ordinogramme d'un procédé de détection de caractéristiques de cycle respiratoire utilisant l'appareil de la figure 1; - la figure 3 est un ordinogramme de calcul de la taille d'une fenêtre dans le procédé de la figure 2 ; - la figure 4 est une vue schématique d'une fonction de débit et de sa pseudo-dérivée issue du procédé de la figure 2 ; - la figure 5 est un ordinogramme de recherche de durée grossière de cycle respiratoire ; - la figure 6 est une courbe d'une pseudo-dérivée d'une courbe de débit ; - la figure 7 est un ordinogramme d'un procédé de recherche de début de cycle respiratoire ; - la figure 8 est une série de courbes permettant de déterminer le début de cycle respiratoire ; et - la figure 9 est un ordinogramme de calcul de la dérivée d'une courbe de débit. En référence à la figure 1, un patient 1 porte un masque respiratoire 3 connecté à un ventilateur, 5 formant un circuit de ventilation. Le circuit de ventilation comporte un capteur 7 de débit de la respiration du patient 1. Ce capteur est, par exemple, le capteur traditionnellement installé dans le ventilateur 5 ou bien est un capteur additionnel. Le capteur 7 n'interfère pas sur la respiration du patient et ne nécessite aucune intervention médicale pour être installée autre que la pose du masque respiratoire lui-même qui, d'ailleurs, est réalisée souvent par le patient lui-même. Ce capteur 7 transmet des informations de débit à un calculateur 9. Il est à noter que, soit au niveau du capteur, soit au niveau du calculateur, l'information analogique de débit est échantillonnée à une fréquence d'échantillonnage prédéterminée facq et numérisée. La fréquence d'échantillonnage facq est choisie pour être supérieure ou égale à 100Hz, par exemple 200Hz soit un échantillon toutes les 0,005 secondes. Le calculateur 9 reçoit ainsi un flux de données numériques de débit pouvant se représenter sous la forme d'une courbe de débit fonction du temps 11 visualisable sur un écran de visualisation 13. Le calculateur 9 comporte également des moyens de stockage 15 sous forme de mémoires ou de disques magnétiques de stockage dans lesquels les données de débit sont stockées sous forme de données numériques horodatées. Par exemple, les données sont stockées sous forme d'entiers signés, chaque donnée étant associée à son heure d'acquisition. Le calculateur 9 est adapté pour déterminer des caractéristiques du débit respiratoire à partir de ces données et, en particulier, la durée et le début de chaque cycle respiratoire. Le procédé de détection des cycles ventilatoires comporte, figure 2 : • acquisition 20 par le capteur de données numériques de débit respiratoire échantillonnées selon la fréquence d'acquisition prédéterminée facq, définissant ainsi une courbe de débit fonction du temps ; • définition 22 d'une taille d'une fenêtre temporelle glissante de débit ; • détermination 24 d'une pseudo-dérivée par : pour chaque fenêtre glissante : • recherche 26 du maximum et du minimum de débit à l'intérieur de la fenêtre glissante ; • détermination 28 de la pseudo-dérivée à partir du maximum et du minimum trouvés, la pseudo-dérivée étant positive si le maximum arrive après le minimum, et étant négative autrement ; • déplacement 30 de la fenêtre glissante vers l'échantillon de débit suivant ; • détection 32 d'un plateau du débit lorsque la pseudo-dérivée est, en valeur absolue, inférieure à une valeur prédéterminée. En reprenant en détail chaque étape de ce procédé, suite à l'acquisition des signaux de débit, le procédé détermine 22 la taille, ou largeur, d'une fenêtre temporelle glissante F, c'est-à-dire le nombre d'échantillons de la fenêtre. Pour cela, figure 3, une estimation spectrale, ou périodogramme, du signal de débit est calculé 34 par la méthode de Welch décrite dans « The Use of Fast Fourier Transform for the Estimation of Power Spectra : A Method Based on Time Averaging Over Short, Modified Periodograms », IEEE Transactions on Audio and Electroacoustics, 15 : 70-3, 1967. Dans cette méthode, une taille de fenêtre utilisée est de 214, la fenêtre étant une fenêtre HANNING, par exemple. Puis un maximum de la densité spectrale est déterminé 36 dans une région définie du spectre, correspondant approximativement à une plage de fréquence de respiration communément rencontrée, soit 9-30 cycles/min. Ce maximum permet d'identifier la fréquence respiratoire moyenne du débit étudié. La taille de la fenêtre temporelle F est alors déterminé 38 à partir du temps total moyen de la durée du cycle respiratoire Ttot, de la fréquence d'acquisition Facq et d'une constante prédéterminée K selon la formule F = Ttot H facq (1) dans laquelle la constante K est généralement choisie égale à environ 2,5. Avant traitement, les données d'acquisition sont filtrées par un filtre passe-bas de façon à éliminer les bruits de mesure. Le filtre passe-bas est, par exemple, constitué d'un premier filtre de type Butterworth passe-bas de 31ème ordre avec une fréquence de coupure basse, par exemple autour de 5Hz, ce filtre étant disposé en série avec un filtre identique mais de manière anti-causale afin d'éliminer le déphasage. A partir du signal filtré, la pseudo-dérivée s est calculée 24, figure 2, en recherchant 26 le maximum et le minimum du signal de débit à l'intérieur de la fenêtre glissante de taille F puis en définissant 28 la pseudo-dérivée s à partir de la différence entre le minimum et le maximum trouvés à l'intérieur de la fenêtre glissante. La pseudo-dérivée s est positive si le maximum arrive après le minimum et est négative si le maximum arrive avant le minimum. La procédure est alors répétée 30 en déplaçant la fenêtre glissante d'un point au point suivant et jusqu'au dernier point du signal de débit. Considérant la pseudo-dérivée s, les intervalles inférieures, en valeur absolue, à une valeur prédéterminée, par exemple 0,08, sont détectés 32 comme des plateaux de débit. La figure 4 est une illustration d'une courbe de débit 40 et de la courbe de sa pseudo-dérivée s 42 sur laquelle une zone 44 est définie comme un plateau de débit car la pseudo-dérivée y est voisine de 0 pendant une durée de presque mille échantillons. Durant ces plateaux, pour la suite des recherches de caractéristique exposée ci-après, la courbe de pseudo-dérivée est remplacée par un signal constant 46 égal à la dernière valeur de la pseudo-dérivée s trouvée avant l'intervalle considéré. Il est à noter que, au vu de la définition de la pseudo-dérivée, les 15 plateaux correspondent à des intervalles pendant lesquels il n'y a pas, ou peu, de variation de débit. A partir de cette pseudo-dérivée s, il est possible de rechercher la durée grossière d'un cycle respiratoire. Pour cela, figure 5, on recherche 50 des points de croisement de la 20 pseudo-dérivée correspondant au passage d'une valeur négative à une valeur positive de la pseudo-dérivée pour une première catégorie in(k) de points de croisement et au passage d'une valeur positive à une valeur négative de la pseudo-dérivée pour une seconde catégorie inn(k) de points de croisement. La durée d'un cycle respiratoire est déterminée 52 sur la base de la 25 durée entre deux points de croisement de même catégorie consécutifs. Par exemple, figure 6, les points in(k), in(k+1) et in(k+2) correspondent au passage de la pseudo-dérivée s d'une valeur négative à une valeur positive : ce sont des points de la première catégorie. Les points inn(k), inn(k+1) et inn(k+2) correspondent au passage de la pseudo-dérivée s 30 d'une valeur positive à une valeur négative : ce sont des points de la seconde catégorie. La durée séparant in(k) et in(k+1) définit une durée de cycle respiratoire. Il en est de même pour la durée séparant in(k+1) et in(k+2). Il est également possible d'utiliser plutôt les points de seconde catégorie. On remarque que le choix de la catégorie peut avoir une influence sur la durée d'un cycle déterminé pris isolément. Cependant, si les durées sont regardées dans une approche statistique, les différences s'annulent et il est possible de prendre l'une ou l'autre catégorie. Cependant, en recherchant précisément le début des cycles respiratoires, il est possible d'obtenir une durée précise de cycle respiratoire comme étant la durée entre deux débuts de cycle. Pour cela, figure 7, la dérivée du débit est calculée 70 préalablement. Puis la position du maximum de débit max(k) dans l'intervalle temporelle compris entre un point de croisement de seconde catégorie inn(k) et le point de croisement de première catégorie suivant in(k+1) est recherchée 72. On recherche 74 la position du point milieu in2(k) entre le point de croisement de première catégorie in(k) précédant le point de croisement de seconde catégorie inn(k) et le maximum max(k), soit in2(k)=(max(k)-in(k))/2. Le maximum maxd(k) de la dérivée du débit et sa position sont recherché 76 dans l'intervalle compris entre le point milieu et ledit maximum de débit, soit [in2k, max(k)]. Le début du cycle respiratoire est alors défini 78 comme : • la position du maximum de la dérivée du débit si la valeur de celui-ci est inférieure à une valeur prédéterminée, par exemple 0,41/s2 ; la position du point le plus tardif, et précédant la position du maximum de la dérivée, de la dérivée du débit inférieure à la valeur prédéterminée si la valeur de celui-ci est supérieure à la valeur prédéterminée. La figure 8 illustre les différents éléments du calcul. La première courbe 80 est la courbe de la pseudo-dérivée, la deuxième courbe 82 la courbe de débit et la troisième courbe 84 la courbe de la dérivée du débit. La première partie 86 de chaque courbe illustre le calcul lorsque le maximum de la dérivée du débit est inférieur à la valeur prédéterminée et la seconde partie 88 illustre le second cas. La dérivée du débit est calculé, par exemple, à partir du signal de débit filtré, en considérant, figure 9, la détermination 90 d'une fréquence de coupure d'un filtre passe-bas du débit, la fréquence de coupure étant déterminée comme le deuxième minimum de la réponse en fréquence d'un filtre à réponse impulsionnelle finie à coefficients Ci, les coefficients Ci correspondant à la première colonne d'une matrice A définie comme A= [VTV] 1VT (2) où V est une matrice à deux colonnes et 2M + 1 lignes, 2M+1 étant le nombre de points d'une fenêtre glissante sur la courbe de débit, VT étant la transposée de la matrice V et ûM 1` û(M -1) 1 (3) V= La dérivée de au centre de ladite fenêtre glissante est calculée 92 10 comme la somme de la valeur du débit en chaque point de la fenêtre glissante pondérée par les coefficients Ci, soit 2M +1 dc= 1CiXi (4) =1 Ce calcul de la dérivée au centre de ladite fenêtre glissante est répété 94 en décalant d'un point ladite fenêtre glissante sur la courbe de débit 15 jusqu'au dernier point de la courbe. Le procédé décrit ci-dessus dans différentes variantes peut être mis en oeuvre sous la forme d'un produit programme d'ordinateur comportant des instructions de programme telles qu'un ordinateur exécutant ces instructions met en oeuvre le procédé. 20 Dans une mise en oeuvre particulière, l'appareil d'assistance respiratoire comprend le capteur de débit et le calculateur dans un ensemble intégré de sorte que l'appareil soit facilement manipulable. De plus, les automatismes d'assistance respiratoire sont adaptés pour intégrer les informations de début et de durée de cycle de façon à adapter l'assistance aux besoins du patient, par exemple, en adaptant les temps de forçage de l'apport en oxygène aux rythmes respiratoires du patient. L'invention a été illustrée et décrite en détail dans les dessins et la description précédente. Celle-ci doit être considérée comme illustrative et donnée à titre d'exemple et non comme limitant l'invention a cette seule description. De nombreuses variantes de réalisation sont possibles. En particulier, le procédé peut être adapté en fonction des besoins. Par exemple, si la caractéristique importante est le début de cycle respiratoire, les calculs peuvent être organisés pour optimiser le calcul de cette caractéristique. Dans les revendications, les mots « comprenant » et « comportant » n'excluent pas d'autres éléments et l'article indéfini « un/une » n'exclue pas une pluralité. Le traitement des données de débit ainsi décrit répond donc au besoin d'identification des cycles respiratoires en étant suffisamment fiable pour des séances de ventilation nocturne, et en permettant une identification en temps réel des asynchronismes, que ce soit sur des patients éveillés ou en état de sommeil. REVENDICATIONS1. Procédé de détection de caractéristiques des cycles respiratoires d'un patient portant un masque d'assistance respiratoire, ledit procédé comportant : • acquisition (20) de données numériques de débit respiratoire échantillonnées selon une fréquence d'acquisition prédéterminée, définissant ainsi une courbe de débit fonction du temps ; • définition (22) d'une taille d'une fenêtre temporelle glissante de débit ; • détermination (24) d'une pseudo-dérivée par : pour chaque fenêtre glissante : • recherche (26) du maximum et du minimum de débit à l'intérieur de la fenêtre glissante ; • détermination (28) de la pseudo-dérivée à partir du maximum et du minimum trouvés, la pseudo-dérivée étant positive si le maximum arrive après le minimum, et étant négative autrement ; • déplacement (30) de la fenêtre glissante vers l'échantillon de débit suivant ; • détection (32) d'un plateau de débit lorsque la pseudo-dérivée est, en valeur absolue, inférieure à une valeur prédéterminée. 2. Procédé selon la revendication 1, caractérisé en ce que la taille de la fenêtre temporelle glissante est déterminée en fonction d'une durée moyenne de cycle respiratoire et de la fréquence d'acquisition. 3. Procédé selon la revendication 2, caractérisé en ce que la durée moyenne de cycle respiratoire est calculée à partir d'un maximum de densité spectrale extrait d'un périodogramme du débit généré par une méthode de Welch. 4. Procédé selon la revendication 1, caractérisé en ce qu'il comporte en outre : 30• recherche (50) des points de croisement de la pseudo-dérivée correspondant au passage d'une valeur négative à une valeur positive de la pseudo-dérivée pour une première catégorie de points de croisement et au passage d'une valeur positive à une valeur négative de la pseudo-dérivée pour une seconde catégorie de points de croisement ; • détermination (52) de la durée d'un cycle respiratoire sur la base de la durée entre deux points de croisement de même catégorie consécutifs. 5. Procédé selon la revendication 4, caractérisé en ce qu'il comporte en outre : • détermination (70) de la dérivée du débit ; • recherche (72) de la position du maximum de débit dans l'intervalle temporelle compris entre un point de croisement de seconde catégorie et le point de croisement de première catégorie suivant ; • définition (74) du point milieu entre le point de croisement de première catégorie précédant ledit point de croisement de seconde catégorie et le maximum ; • détermination (76) du maximum de la dérivée du débit et de sa position dans l'intervalle compris entre le point milieu et ledit maximum de débit ; • détermination (78) du début du cycle respiratoire comme : • la position du maximum de la dérivée du débit si la valeur de celui-ci est inférieure à une valeur prédéterminée, • la position du point le plus tardif, et précédant la position du maximum de la dérivée, de la dérivée du débit inférieure à ladite valeur prédéterminée si la valeur de celui-ci est supérieure à ladite valeur prédéterminée. 6. Procédé selon la revendication 5, caractérisé en ce que ladite valeur prédéterminée est égale à 0,4 I/s2. 7. Procédé selon la revendication 5, caractérisé en ce que la détermination de la dérivée du débit comporte :30 • détermination (90) d'une fréquence de coupure d'un filtre passe-bas du débit, ladite fréquence de coupure étant déterminée comme le deuxième minimum de la réponse en fréquence d'un filtre à réponse impulsionnelle finie à coefficients C;, lesdits coefficients correspondant à la première colonne d'une matrice A définie comme A = [VTV]1VT où v est une matrice à deux colonnes et 2M + 1 lignes, 2M+1 étant le nombre de points d'une fenêtre glissante sur la courbe de débit, VT étant la ùM 1' ù(M -1) 1 transposée de la matrice V et V = Mù1 1 M 1J • détermination (92) de la dérivée au centre de ladite fenêtre glissante 10 comme la somme de la valeur du débit en chaque point de la fenêtre 2M+1 glissante pondérée par les coefficients C;, soit 1CiXi ; i=1 • répétition (94) de la détermination de la dérivée au centre de ladite fenêtre glissante en décalant d'un point ladite fenêtre glissante sur la courbe de débit jusqu'au dernier point de la courbe. 15 8. Produit programme d'ordinateur comprenant des instructions de code de programme enregistrées sur un support lisible par un ordinateur, pour mettre en oeuvre les étapes du procédé selon l'une quelconque des revendications 1 à 7 lorsque ledit programme fonctionne sur un ordinateur. 20 9. Appareil d'assistance respiratoire comprenant : • un capteur (7) de débit respiratoire fournissant des données numériques de débit respiratoire échantillonnées selon une fréquence d'acquisition prédéterminée, définissant ainsi une courbe de débit fonction du temps ; 25 ledit capteur fournissant lesdites données à 510• un calculateur (9) définissant une taille d'une fenêtre temporelle glissante de débit et déterminant une pseudo-dérivée par : • pour chaque fenêtre glissante : • recherche du maximum et du minimum de débit à l'intérieur de la fenêtre glissante ; • détermination de la pseudo-dérivée à partir du maximum et du minimum trouvés, la pseudo-dérivée étant positive si le maximum arrive après le minimum, et étant négative autrement ; • déplacement de la fenêtre glissante vers l'échantillon de débit suivant ; • ledit calculateur détectant un plateau de débit lorsque la pseudo-dérivée est, en valeur absolue, inférieure à une valeur prédéterminée.