WO2014009366A1 - Procédé de détection du colmatage d'un échangeur thermique - Google Patents

Procédé de détection du colmatage d'un échangeur thermique Download PDF

Info

Publication number
WO2014009366A1
WO2014009366A1 PCT/EP2013/064480 EP2013064480W WO2014009366A1 WO 2014009366 A1 WO2014009366 A1 WO 2014009366A1 EP 2013064480 W EP2013064480 W EP 2013064480W WO 2014009366 A1 WO2014009366 A1 WO 2014009366A1
Authority
WO
WIPO (PCT)
Prior art keywords
detection signal
signal
peaks
probe
measurement signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Ceased
Application number
PCT/EP2013/064480
Other languages
English (en)
Inventor
Nicolas Paul
Mabrouka EL-GUEDRI
Didier BOLDO
Alexandre Girard
Lorenzo Audibert
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Electricite de France SA
Original Assignee
Electricite de France SA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Electricite de France SA filed Critical Electricite de France SA
Publication of WO2014009366A1 publication Critical patent/WO2014009366A1/fr
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/72—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables
    • G01N27/82—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws
    • G01N27/90—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws using eddy currents
    • G01N27/9046—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating magnetic variables for investigating the presence of flaws using eddy currents by analysing electrical signals
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F22—STEAM GENERATION
    • F22B—METHODS OF STEAM GENERATION; STEAM BOILERS
    • F22B37/00—Component parts or details of steam boilers
    • F22B37/002—Component parts or details of steam boilers specially adapted for nuclear steam generators, e.g. maintenance, repairing or inspecting equipment not otherwise provided for
    • F22B37/003—Maintenance, repairing or inspecting equipment positioned in or via the headers
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21C—NUCLEAR REACTORS
    • G21C17/00—Monitoring; Testing ; Maintaining
    • G21C17/017—Inspection or maintenance of pipe-lines or tubes in nuclear installations
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21D—NUCLEAR POWER PLANT
    • G21D1/00—Details of nuclear power plant
    • G21D1/006—Details of nuclear power plant primary side of steam generators
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Definitions

  • the present invention relates generally to the field of tube inspection of a tube heat exchanger. More specifically, it relates to a method for detecting various degradations, for example clogging, of a tube heat exchanger in which a multifrequency eddy current probe is circulated in a tube of said exchanger and is measured therewith. a measurement signal depending on the environment in which the probe is located.
  • a steam generator is generally composed of a bundle of tubes in which circulates the hot fluid, and around which circulates the fluid to
  • the steam generators are heat exchangers that use the energy of the primary circuit resulting from the nuclear reaction to transform the circuit water secondary steam that will power the turbine and thus produce electricity.
  • the steam generator brings the secondary fluid from a state of liquid water to a vapor state just at the saturation limit, using the heat of the primary water. This circulates in tubes around which the secondary water circulates.
  • the output of the steam generator is the highest point in temperature and pressure of the secondary circuit.
  • the exchange surface physically separating the two circuits, thus consists of a tubular bundle, composed of 3500 to 5600 tubes, according to the model, in which the primary water carried at high temperature (320 ° C.) circulates and high pressure (155bars).
  • spacers arranged generally perpendicular to the tubes which pass through the passages of the spacer plates.
  • the passages of these spacer plates are folished, that is to say whose shape has lobes around the tubes.
  • the water changes from the liquid state to the vapor state, it deposits all the contents it contained. If the deposits of matter are in the lobes, they reduce the free passage: it is the clogging, which is therefore the gradual filling, by deposits, holes for the passage of the water / steam mixture.
  • FIG. 1 schematically illustrates a folate passage in a spacer plate 10, through which a tube 11 passes.
  • the lobes 12a and 12b allow the water to pass through the spacer plate 10 along the tube 11, thus allowing the flow of the water in the steam generator.
  • a deposit 13 is visible at the level of the lobe 12b, clogging said lobe 12b.
  • SAX probe eddy current axial probe
  • the eddy current probe is already used during slice shutdowns, in particular to inspect the integrity of the tubes of the steam generator.
  • This eddy current probe initially intended for the detection of damage to the tubes, is also sensitive to clogging. Moreover, the interpretation of this signal is currently carried out manually by specialized operators, which is very long, of the order of about a week of treatment for the analysis of a single steam generator. In addition, the intervention of an operator to take measurements from an analysis software often gives rise to a bias that is difficult to quantify.
  • the measurement signal is not calibrated and is noisy, so that its operation can be difficult.
  • This analysis is based on the construction, from one-dimensional data, of a two-dimensional image representative of the surface of a tube, taking advantage of the rotation of the probe and interpolating. A Sobel filter is then applied to the resulting image to identify intensity changes in it, and changes above a threshold are
  • the purpose of the invention is to overcome at least one of these defects, all of them preferentially, by proposing a method that automates the extraction of the portions of the signal emitted by a SAX probe corresponding to the passage of the spacer plates of a tube heat exchanger.
  • a PWR nuclear steam generator preferentially a PWR nuclear steam generator, without making use of considerations on the geometry of the steam generator, while providing a calibrated and denoised signal.
  • a method for detecting spacer plates of a heat exchanger in which a multifrequency eddy current probe is circulated in a tube of said exchanger and a measurement signal according to the invention is measured therewith. the environment in which the probe is located.
  • Different components of the measurement signal are combined to determine a detection signal, the combination thus performed being chosen to minimize the energy of the detection signal,
  • the combination of the different components of the measurement signal is an adaptive linear combination, the values of the weighting coefficients used for this combination being determined, for a sample of the detection signal, by minimizing the variance of the detection signal on a window observing the measurement signal around this sample; different components of the measurement signal combined to determine the detection signal are real and imaginary components of the measurement signal;
  • different components of the measurement signal combined to determine the detection signal are components of different frequencies
  • the detection signal is compared with a detection threshold, said threshold being a function of the minimum value that the standard deviation of the detection signal takes on a part of the corresponding signal. a number of samples smaller than the number of theoretical samples between two consecutive spacer plates;
  • the detected peaks are limited to the peaks corresponding to a local maximum of the detection signal over a range of the detection signal corresponding to a number of samples on each side of the local maximum less than or equal to the minimum of the following three values:
  • At least one subset of peaks minimizing the difference between the number of samples of the detection signal between two successive peaks of said subset and a median value of the number of samples of the detection signal between two successive peaks;
  • the measurement signal corresponds to the measurement of an eddy current detector probe acquiring said measurement signal during the return of the probe in a back-and-forth displacement of said probe in a tube of a heat exchanger;
  • the method comprises a preliminary calibration of the measurement signal, comprising:
  • the method further comprises a subsequent verification step in which the positions of the peaks of the detection signal corresponding to the passage of the spacer plates are compared with known data on the geometry of the heat exchanger;
  • the method comprises a final verification step in which the opposite of the detection signal is taken, the steps of detecting and selecting the positions of the peaks of the detection signal corresponding to the passage of the spacer plates are applied and verifying that the deviations between the positions of the peaks on
  • the detection signal and on its opposite correspond to the edge-to-edge distance of a spacer plate.
  • the invention also relates to a computer program product comprising program code instructions for performing process processing steps according to one of the preceding claims, when said
  • FIG. 2 is a block diagram illustrating the calibration of the measurement signal in the context of the invention
  • FIG. 3 is a schematic diagram of the method according to the invention.
  • a multifrequency eddy current probe is circulated in a tube 11 of said exchanger and a measurement signal is measured therewith which is a function of the environment in which the probe is located.
  • the measurement signal is therefore acquired by the eddy current axial probe during its displacement in the tube 11 of the steam generator.
  • This measurement signal is multifrequency, and is in complex form, with real and imaginary components.
  • the probe makes a round trip in the tube, the signals corresponding to the phase of displacement where the probe is pulled by the mechanism that moves it, typically during its return, will preferably be selected. Indeed, it results in a better mechanical stability of the probe, and consequently a higher speed
  • the selection of the part of the signal of the probe corresponding to the return of the probe is based on the detection of a significant drop in the voltage measured at the moment when the probe switches to traverse the tube in the direction of the return.
  • the selection of a measurement signal limited to one direction of the probe allows
  • a preliminary calibration step is preferably implemented.
  • step S0 is carried out in the same manner regardless of the acquisition frequency considered, with, with reference to FIG. 2, for each
  • step S20 identification (step S20) of the parts of the measurement signal corresponding to known defects of the tube of the heat exchanger
  • step S21 estimation (step S21) of the phase and the amplitude of the parts of the measurement signal corresponding to the known defects, determination (step S22) of a transformation to be applied to the measurement signal from said parts of the measurement signal corresponding to the known defects for calibrating the measurement signal,
  • step S23 application of the transformation to the measurement signal.
  • a calibration defect is characterized by two parameters A 0 and ⁇ 0 which are respectively the amplitude and the phase of the vector max-min (vector i od 'maximum amplitude on the representation according to Lissajous). These are the characteristics that the calibration defect must have after calibration.
  • the geometrical configuration of the tube or the heat exchanger is not used to detect the zone of the fault, but only to the characteristics of the signal (phase and amplitude). This provides an independent calibration of the probe handling or signal sampling problems, which may bias the judgment of the position of the part of the signal corresponding to the defect.
  • the measurement signal Once the measurement signal has been calibrated, it remains to determine the parts of the multifrequency measurement signal corresponding to the passage at the level of the spacer plates 10 of the eddy current detector probe acquiring said measurement signal in the tube of the heat exchanger .
  • the measurement signal is very noisy and can not be used directly to detect the spacer plates 10.
  • a detection signal denoted d [n] is constructed by linear combination of components of the frequency measurement signal. different to minimize the energy of the detection signal.
  • the detection signal then takes a minimum value between the spacer plates 10, while taking a large value when passing the spacer plates 10.
  • the signals available at the input of the algorithm are the three complex calibrated signals z A [n], z ⁇ n] and z 3 [n] corresponding to the return of the probe in the tube. Only the signals corresponding to the differential mode (z1 and z3), more sensitive to the passage of the probe by a plate, are used in this part.
  • x ( [n] and [n] respectively the real and imaginary components of the signal 3 ⁇ 4 [n].
  • a high-pass filter for example a Butterworth filter
  • z 1f and z 3f The signals resulting from the filtering of the signals z 1 and z 3 are denoted z 1f and z 3f .
  • the real and imaginary parts of these signals are therefore noted:
  • zi f [n] xi f [n] + jy 1f [n]
  • This first filtering operation makes it possible to construct signals having a baseline of constant average and zero.
  • each distance can be converted into a number of signal samples according to the has formula:
  • an observation window corresponding to a length of M samples smaller than the number of theoretical samples 25 may be obtained between two successive spacer plates ⁇ ⁇ , ⁇ , for example
  • This search is performed for example on the signal X 3f .x 3f having a high variance when passing the probe at a spacer plate 10. It is sufficient to find such a window to find the portion of the signal length M where the variance of x 3f is minimal. Noting m 0 the first index of 30 this window we have:
  • the combination of the different components of the measurement signal is an adaptive linear combination, the values of the weighting coefficients used for this combination being determined, for a sample of the detection signal, by minimizing the variance of the detection signal on the window of observation of the measurement signal around this sample.
  • the different components of the measurement signal combined to determine the detection signal are real and imaginary components of the measurement signal.
  • different components of the measurement signal combined to determine the detection signal are different frequency components.
  • the coefficients ⁇ , ⁇ and ⁇ are updated for each sample of signal. More precisely, the triplet ⁇ [ ⁇ ], ⁇ [ ⁇ ], ⁇ [ ⁇ ] ⁇ chosen at an instant n is the one that minimizes the power of the reconstruction error on a signal window centered on the nth sample: [k] - "if [k] - ⁇ 3 f [k] - yyi / [k]
  • a detection signal d [n] is thus obtained by adaptive linear combination of real and imaginary components of the measurement signal.
  • the combination can be done between all the components (real and imaginary parts) of the different frequencies;
  • the combination coefficients are not constant along the tube, but vary in order to adapt to the tube. Thus, the variations of the state of the tube are managed.
  • the detection signal Once the detection signal has been constructed, it remains to detect peaks of the detection signal likely to correspond to the passage of spacer plates 10, compared with a detection threshold. Indeed, the signal of
  • the detection thus created contains different peaks, which correspond to the passage of the probe by the different plates, and not only the spacer plates 10, and other structural elements.
  • the detection threshold for detecting the peaks of the detection signal is a function of the minimum value that the standard deviation ⁇ of the detection signal has on a part of the signal corresponding to a number of samples less than the number of samples. theoretical samples between two consecutive spacer plates 10.
  • the standard deviation ⁇ can be determined as follows:
  • the detection threshold must then be a compromise between the risk of false detections (a threshold that is too low results in the detection of numerous peaks that do not correspond to the passage of a spacer plate 10by the probe) and the risk of not detecting a spacer plate. 10. It may for example be between 5 and 15 ⁇ , preferably 10 ⁇ .
  • n1 and n1 are two indices of successive peaks, we say that ni is a "secondary peak" of n1 if d [n1]> d [n2] and if the number of samples between n1 and n1 is much smaller than the theoretical minimum distance between two plate peaks.
  • the peaks detected are limited to the peaks corresponding to a local maximum of the detection signal over a range of the detection signal corresponding to a number of samples on each side of the maximum local less than or equal to at least the following three values:
  • the secondary peaks being suppressed, it then remains to select peaks of the detection signal corresponding to the passage of the spacer plates 10 by a determination of the peaks of the signal having a regular spacing.
  • the number of theoretical samples between two plates can be calculated from a knowledge of the distance between two spacer plates. Ideally, if the probe has a known and constant velocity v, then the number of samples between two plates is exactly proportional to the distance between plates:
  • the speed of the probe is not exactly constant or precisely known and can vary along the tube.
  • the difference between two peaks corresponding to spacer plates 10 is therefore not exactly equal to N P e, p e and it may be different from one pair of spacer plates 10 to another. This average gap between spacer plates 10 can however be estimated.
  • the peaks of the detection signal corresponding to the passage of the spacer plates 10 are selected by selecting at least one subset of peaks minimizing the difference between the number of samples of the detection signal between two peaks. successive of said subset and a median value of the number of samples of the detection signal between two peaks
  • a tube of a heat exchanger extends along the latter through two parts, generally called cold leg and hot leg, depending on the direction of flow of the fluids carrying out the heat exchange, separated by a plate-free structure spacer 10 designated as the hanger. Therefore, such a structure breaks the regularity of spacing of the spacer plates. It suffices to call the following method for one branch, then for the other, to select all the peaks corresponding to the spacer plates 10.
  • N pe eg, the number of k Average samples between two spacer plates 10.
  • the spacer plates are the majority plates along the tube, the majority of the values taken by ⁇ [k] is located around N pe , eg, real.
  • An estimate of N pe , pe, real is thus provided by the median value of the set of values taken by A [k]:
  • Npe.pe.réi median ⁇ A [k] ⁇
  • the indices of the samples of the signal corresponding to the passage of the spacer plates 10 are then obtained.
  • the first verification criterion concerns the length of the hanger: for each tube, the number of samples between the last plate of the hot branch and the last plate of the cold branch is compared with the known distance according to the plans. A margin is allowed, to take into account the uncertainty on the speed of the probe.
  • the second validation criterion makes it possible to check the length of all the detected plates.
  • the detection algorithm having made it possible to estimate the indices of the spacer plates 10 uses the positive peaks of the detection signal. However the passage of the probe by a spacer plate 10 generates on the signal at least one positive peak and one negative peak (the two peaks correspond to the passages of the probe by the plate edges).
  • a new estimate of the indexes of the spacer plates 10 can thus be obtained by detecting this time the negative peaks of the detection signal, without taking into account the first peaks detected.
  • the method may thus comprise a final verification step in which the opposite of the detection signal is taken, the detection S2 and S3 selection steps are applied to the positions of the peaks of the detection signal corresponding to the passage of the spacer plates 10 and verifies that the differences between the positions of the peaks on the detection signal and on its opposite correspond to the distance edge-to-edge of a spacer plate 10.
  • the invention also relates to a computer program product comprising program code instructions recorded on a medium usable in a computer for performing process processing steps according to one of the preceding claims, when said program is executed. on a computer.
  • the measurement signal is transmitted from the eddy current probe to a memory for storage therein for processing.
  • This processing of the measurement signal to which the present invention relates is implemented by a processing unit provided with a computer, typically a computer provided with display and communication means, by which it acquires the measurement signal and transmits the results of the implementation of the process.

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Plasma & Fusion (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Analytical Chemistry (AREA)
  • Thermal Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mechanical Engineering (AREA)
  • Electrochemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)

Description

Procédé de détection du colmatage d'un échangeur thermique
DOMAINE TECHNIQUE GENERAL
La présente invention se rapporte d'une façon générale au domaine de 5 l'inspection des tubes d'un échangeur de chaleur à tubes. Plus précisément, elle concerne un procédé de détection de différentes dégradations, par exemple le colmatage, d'un échangeur thermique à tubes, dans lequel on fait circuler dans un tube dudit échangeur une sonde à courant de Foucault multifréquence et on mesure avec celle-ci un signal de mesure fonction de l'environnement dans lequel î o la sonde se trouve.
CONTEXTE DE L'INVENTION
Un générateur de vapeur est généralement composé d'un faisceau de tubes dans lesquels circule le fluide chaud, et autour desquels circule le fluide à
15 chauffer. Par exemple, dans le cas d'un générateur de vapeur d'une centrale nucléaire de type REP, les générateurs de vapeur sont des échangeurs de chaleur qui utilisent l'énergie du circuit primaire issue de la réaction nucléaire pour transformer l'eau du circuit secondaire en vapeur qui alimentera la turbine et ainsi à produire de l'électricité.
20 Le générateur de vapeur amène le fluide secondaire d'un état d'eau liquide à l'état de vapeur juste en limite de saturation, en utilisant la chaleur de l'eau primaire. Celle-ci circule dans des tubes autour desquels circule l'eau secondaire. La sortie du générateur de vapeur est le point le plus haut en température et pression du circuit secondaire.
25 La surface d'échange, séparant physiquement les deux circuits, est ainsi constituée d'un faisceau tubulaire, composé de 3500 à 5600 tubes, selon le modèle, dans lesquels circule l'eau primaire portée à haute température (320° C) et haute pression (155bars).
Ces tubes du générateur de vapeur sont maintenus par des plaques
30 entretoises disposés généralement perpendiculaires aux tubes qui les traversent au niveau de passages des plaques entretoises.
Afin de laisser passer le fluide qui se vaporise, les passages de ces plaques entretoises sont foliés, c'est-à-dire dont la forme présente des lobes autour des tubes. Comme l'eau passe de l'état liquide à l'état vapeur, elle dépose toutes les matières qu'elle contenait. Si les dépôts de matière se font dans les lobes, ils diminuent le passage libre : c'est le colmatage, qui est donc l'obturation progressive, par des dépôts, des trous destinés au passage du mélange eau/vapeur.
5 La figure 1 illustre schématiquement un passage folié dans une plaque entretoise 10, dans lequel passe un tube 11. Les lobes 12a et 12b permettent à l'eau de traverser la plaque entretoise 10 le long du tube 11 , permettant ainsi la circulation de l'eau dans le générateur de vapeur. Un dépôt 13 est visible au niveau du lobe 12b, colmatant ledit lobe 12b.
î o Le colmatage conduit à des modifications de l'écoulement de l'eau dans le générateur de vapeur, et ainsi à favoriser l'apparition de vibrations excessives des tubes, ainsi qu'induire des efforts mécaniques importants sur les structures internes des générateur de vapeur. Cette dégradation a donc des effets à la fois sur la sûreté et sur les performances des installations. Il est donc indispensable de
15 bien connaître la nature et l'évolution de cette dégradation.
Actuellement, le seul système d'examen non destructif qui soit capable d'accéder à la totalité des tubes des générateurs de vapeur est la sonde axiale à courant de Foucault (sonde SAX). Les courants de Foucault apparaissent dans un matériau conducteur lorsque l'on fait varier le flux magnétique à proximité. On
20 fait ainsi circuler dans un tube dudit échangeur une sonde à courant de Foucault multifréquence et on mesure avec celle-ci un signal de mesure fonction de l'environnement dans lequel la sonde se trouve, duquel on peut extraire des informations quant à des anomalies dans l'échangeur thermique.
Une variation de l'induction magnétique, typiquement par une bobine dans
25 laquelle circule un courant alternatif, engendre des courants de Foucault, dont la variation induite du champ magnétique est détectée. Typiquement, on mesure la différence de tension engendrée par la variation d'impédance de la bobine.
L'exploitation des signaux de mesure de cette sonde à courant de Foucault n'induit pas d'allongement de l'arrêt du générateur de vapeur, puisque cette
30 sonde à courant de Foucault est déjà utilisée lors des arrêts de tranche, notamment pour inspecter l'intégrité des tubes du générateur de vapeur.
Cette sonde à courant de Foucault, initialement destinée à la détection d'endommagement des tubes, est également sensible au colmatage. De plus, l'interprétation de ce signal est actuellement réalisée manuellement par des opérateurs spécialisés, ce qui est très long, de l'ordre d'une semaine de traitement environ pour l'analyse d'un seul générateur de vapeur. De plus, l'intervention d'un opérateur pour relever des mesures à partir d'un logiciel d'analyse donne souvent lieu à un biais difficile à quantifier.
5 Un problème important consiste à extraire le signal utile, c'est-à-dire les parties du signal de mesure multifréquence correspondant au passage au niveau des plaques entretoises 10 de la sonde à courant de Foucault acquérant ledit signal de mesure. En effet, la variabilité de la vitesse de la sonde à courants de Foucault, ainsi que des imprécisions quant à la géométrie réelle du générateur de î o vapeur, rend peu fiable l'extraction de ces parties sur la seule base d'une estimation de leur position par des considérations spatiales.
En outre, le signal de mesure n'est pas calibré et est bruité, de sorte que son exploitation peut se révéler difficile.
Le document "Automated analysis of rotating probe multi-frequency eddy
15 current data from steam generator tubes", par P. Xiang et al., International Journal of Applied Electromagnetics and Mechanics, vol.12, p.151 -164, propose un algorithme pour l'analyse automatisée de données de sonde de Foucault multifréquence dans le cadre de tubes de générateurs de vapeur de centrale nucléaires.
20 Cette analyse se base sur la construction, à partir de données en une dimension, d'une image en deux dimensions représentative de la surface d'un tube, en tirant parti de la rotation de la sonde et en interpolant. Un filtre de Sobel est ensuite appliqué sur l'image résultante pour identifier les changements d'intensité dans celle-ci, et les changements supérieurs à un seuil sont
25 sélectionnés.
Une telle approche nécessite une sonde rotative, et, basé sur des images, ne permet pas une détection précise parties du signal correspondant au passage de la sonde au niveau des plaques entretoises, en raison notamment de la présence de bruit, qui peut être important.
30 Le document "Estimation of mixing parameters for cancellation of discretized eddy current signais using time and frequency domain techniques", par C. K. Sword et al., Journal of Nondestructive Evaluation, vol.5, n° 1 , pages 27-35, décrit un procédé dans lequel les composantes d'un signal d'une sonde à courant de Foucault sont combinées dans le but de supprimer des signaux de support de plaque indésirables. Le procédé en question ne permet donc pas d'identifier les parties du signal correspondant au passage de la sonde au niveau des plaques entretoises.
5 PRESENTATION DE L'INVENTION
L'invention a pour but de pallier au moins un de ces défauts, préférentiellement tous, en proposant un procédé automatisant l'extraction des portions du signal émis par une sonde SAX correspondant au passage des plaques entretoises d'un échangeur de chaleur à tubes, préférentiellement un générateur î o de vapeur de centrale nucléaire REP, sans faire appel à des considérations sur la géométrie du générateur de vapeur, tout en procurant un signal étalonné et débruité.
A cet effet, on propose un procédé de détection des plaques entretoises d'un échangeur thermique, dans lequel on fait circuler dans un tube dudit 15 échangeur une sonde à courant de Foucault multifréquence et on mesure avec celle-ci un signal de mesure fonction de l'environnement dans lequel la sonde se trouve.
Pour identifier les parties de ce signal correspondant au passage de la sonde au niveau des plaques entretoises,
20 - on combine différentes composantes du signal de mesure pour déterminer un signal de détection, la combinaison ainsi effectuée étant choisie pour minimiser l'énergie du signal de détection,
- on détecte sur ce signal des pics de détection,
- on sélectionne parmi ces pics ceux qui présentent un espacement 25 régulier, les pics ainsi sélectionnés correspondant au passage de la sonde au niveau des plaques entretoises.
L'invention est avantageusement complétée par les caractéristiques suivantes, prises seules ou en une quelconque de leur combinaison techniquement possible:
30 - la combinaison des différentes composantes du signal de mesure est une combinaison linéaire adaptative, les valeurs des coefficients de pondération utilisés pour cette combinaison étant déterminés, pour un échantillon du signal de détection, par minimisation de la variance du signal de détection sur une fenêtre d'observation du signal de mesure autour de cet échantillon ; - différentes composantes du signal de mesure combinées pour déterminer le signal de détection sont des composantes réelles et imaginaires du signal de mesure ;
- différentes composantes du signal de mesure combinées pour déterminer le 5 signal de détection sont des composantes de fréquences différentes ;
- pour détecter des pics de détection sur le signal de détection, on compare le signal de détection à un seuil de détection, ledit seuil étant fonction de la valeur minimum que prend l'écart-type du signal de détection sur une partie du signal correspondant à un nombre d'échantillons inférieur au nombre d'échantillons î o théorique entre deux plaques entretoises consécutives ;
- préalablement à l'étape de sélection des pics, on restreint les pics détectés aux pics correspondant à un maximum local du signal de détection sur une plage du signal de détection correspondant à un nombre d'échantillons de part et d'autre du maximum local inférieur ou égal au minimum des trois valeurs suivantes :
15 - 0,5 fois un nombre d'échantillons théorique entre deux plaques entretoises ,
0,8 fois un nombre d'échantillons théorique entre la plaque entretoise de la branche chaude la plus proche de la branche froide dudit échangeur thermique et la plaque entretoise de la branche froide la plus proche de la branche chaude dudit échangeur thermique,
20 - 0,8 fois un nombre d'échantillons théorique entre une première des plaques entretoises et une autre plaque précédente dans le sens de circulation de la sonde.
- pour sélectionner les pics du signal de détection présentant un espacement régulier et correspondant au passage des plaques entretoises, on détermine au
25 moins un sous-ensemble de pics minimisant la différence entre le nombre d'échantillons du signal de détection entre deux pics successifs dudit sous- ensemble et une valeur médiane du nombre d'échantillons du signal de détection entre deux pics successifs ;
- le signal de mesure correspond à la mesure d'une sonde détectrice de courants 30 de Foucault acquérant ledit signal de mesure lors du retour de la sonde dans un déplacement aller-retour de ladite sonde dans un tube d'un échangeur thermique ;
- le procédé comporte un étalonnage préalable du signal de mesure, comprenant :
- l'identification des parties du signal de mesure correspondant à des défauts connus du tube de l'échangeur thermique, - l'estimation de la phase et de l'amplitude des parties du signal de mesure correspondant aux défauts connus,
- la détermination d'une transformation à appliquer au signal de mesure à partir desdites parties du signal de mesure correspondant aux défauts connus pour
5 étalonner le signal de mesure,
- l'application de la transformation au signal de mesure ;
- le procédé comporte en outre une étape subséquente de vérification dans laquelle les positions des pics du signal de détection correspondant au passage des plaques entretoises sont comparées avec des données connues sur la géométrie de î o l'échangeur thermique;
- le procédé comprend une étape finale de vérification dans laquelle on prend l'opposé du signal de détection, on applique les étapes de détection et de sélection des positions des pics du signal de détection correspondant au passage des plaques entretoises et on vérifie que les écarts entre les positions des pics sur
15 le signal de détection et sur son opposé correspondent à la distance bord à bord d'une plaque entretoise.
L'invention concerne également un produit programme d'ordinateur comprenant des instructions de code de programme pour l'exécution d'étapes de traitement du procédé selon l'une des revendications précédentes, lorsque ledit
20 programme est exécuté sur un ordinateur.
PRESENTATION DES FIGURES
D'autres caractéristiques, buts et avantages de l'invention ressorti ront de la description qui suit, qui est purement illustrative et non limitative, et qui doit 25 être lue en regard des dessins annexés parmi lesquels:
- la figure 1 , déjà commentée, illustre schématiquement un passage folié dans une plaque entretoise, dans lequel passe un tube, selon une configuration courante d'un générateur de vapeur ;
- la figure 2 est un schéma de principe illustrant l'étalonnage du signal de 30 mesure dans le cadre de l'invention ;
- la figure 3 est un schéma de principe du procédé selon l'invention.
DESCRIPTION DETAILLEE On fait circuler dans un tube 11 dudit échangeur une sonde à courant de Foucault multifréquence et on mesure avec celle-ci un signal de mesure fonction de l'environnement dans lequel la sonde se trouve. Le signal de mesure est donc acquis par la sonde axiale à courant de Foucault lors de son déplacement dans le 5 tube 11 du générateur de vapeur. Ce signal de mesure est multifréquence, et se présente sous forme complexe, avec des composantes réelles et imaginaires. Dans le cadre de cette description, nous considérerons un signal avec des composantes complexes selon trois fréquences issues de trois canaux d'acquisition, dont deux en mode différentiels notés z^n] et z3[n] et un en mode absolu noté zA[n],. Les î o fréquences sont par exemple comprises dans la plage 100-600 kHz.
La sonde effectuant un aller- retour dans le tube, on choisira de préférence les signaux correspondant à la phase de déplacement où la sonde est tirée par le mécanisme qui la déplace, typiquement lors de son retour. En effet, il en résulte une meilleure stabilité mécanique de la sonde, et par conséquent une vitesse plus
15 régulière de celle-ci.
La sélection de la partie du signal de la sonde correspondant au retour de la sonde est fondée sur la détection d'une chute importante de la tension mesurée au moment où la sonde bascule pour parcourir le tube dans le sens du retour. La sélection d'un signal de mesure limité à un seul sens de la sonde permet
20 également de limiter la quantité de données à sauvegarder et à traiter.
Compte tenu de la sensibilité des acquisitions des signaux des courants de Foucault aux conditions d'acquisition, au basculement de la sonde et afin de standardiser les signaux analysés, on met de préférence en œuvre une étape préalable d'étalonnage.
25 II s'agit d'étalonner les signaux par rapport à des signaux de référence, appelés signaux d'étalonnage. Ces derniers correspondent à des défauts implantés sur le tube d'étalonnage, et dont on connaît les caractéristiques (amplitude et phase). L'étalonnage (étape S0) s'effectue de la même manière quelle que soit la fréquence d'acquisition considérée, avec, en référence à la figure 2, pour chaque
30 tube contrôlé :
- identification (étape S20) des parties du signal de mesure correspondant à des défauts connus du tube de l'échangeur thermique,
- estimation (étape S21 ) de la phase et de l'amplitude des parties du signal de mesure correspondant aux défauts connus, - détermination (étape S22) d 'une transformation à appliquer au signal de mesure à partir desdites parties du signal de mesure correspondant aux défauts connus pour étalonner le signal de mesure,
- application (étape S23) de la transformation au signal de mesure.
5 Soient zbrut(f)=x+jy (où j est le nombre imaginaire) issu d'une acquisition à la fréquence f, et
Figure imgf000010_0001
l'impédance complexe mesurée pour un défaut d 'étalonnage donné à la même fréquence.
Un défaut d 'étalonnage est caractérisé par deux paramètres A0 et θ0 qui sont respectivement l'amplitude et la phase du vecteur max-min (vecteur î o d 'amplitude maximale sur la représentation selon Lissajous). Ce sont les caractéristiques que le défaut d 'étalonnage doit présenter après étalonnage.
Avant étalonnage, ces deux paramètres valent respectivement A et Θ. L'étalonnage du signal z consiste alors en une rotation d 'angle δφ=θ0-θ et une homothétie de paramètre h=A0/A tel que décrit par l'équation suivante :
15 z = h∑o e}5^
De fait, il n'est pas fait appel à la configuration géométrique du tube ou de l'échangeur thermique pour détecter la zone du défaut, mais uniquement aux caractéristiques du signal (phase et amplitude). On dispose ainsi d'un étalonnage indépendant des problèmes de manipulation de la sonde ou d 'échantillonnage du 20 signal, susceptible de fausser l'appréciation de la position de la partie du signal correspondant au défaut.
Une fois le signal de mesure étalonné, il reste à déterminer les parties du signal de mesure multifréquence correspondant au passage au niveau des plaques entretoises 10 de la sonde détectrice de courants de Foucault acquérant ledit 25 signal de mesure dans le tube de l'échangeur thermique.
Le signal de mesure est très bruité et ne peut être utilisé directement pour détecter les plaques entretoises 10. Afin de faciliter cette détection, un signal de détection, noté d [n] , est construit par combinaison linéaire de composantes du signal de mesure de fréquences différentes pour minimiser l'énergie du signal de 30 détection. Ainsi, le signal de détection prend alors une valeur minimale entre les plaques entretoises 10, tout en prenant une valeur importante au passage des plaques entretoises 10.
Les signaux disponibles en entrée de l'algorithme sont les trois signaux complexes étalonnés zA[n], z^n] et z3[n] correspondant au retour de la sonde dans le tube. Seuls les signaux correspondant au mode différentiel (z1 et z3), plus sensibles aux passages de la sonde par une plaque, sont utilisés dans cette partie. On notera x([n] et [n] les composantes respectivement réelles et imaginaires du signal ¾[n] .
5 Dans un premier temps l'ensemble des mesures disponibles est filtré par un filtre passe haut pour supprimer les faibles variations des signaux dues au mouvement horizontal de la sonde dans le tube. Pour chaque signal disponible on applique un filtre passe-haut (par exemple un filtre de Butterworth), par exemple de fréquence réduite de coupure 0.01 . Les signaux résultant du filtrage des î o signaux zi et z3 sont notés z1f et z3f. Les parties réelles et imaginaires de ces signaux sont donc notées :
zif[n]= xif[n]+j y1f[n]
Figure imgf000011_0001
Cette première opération de filtrage permet de construire des signaux 15 présentant une ligne de base de moyenne constante et nulle. Pour ensuite construire un signal de détection qui soit le plus faible possible entre les plaques, on commence par rechercher une fenêtre d 'observation exprimée en nombre d 'échantillons de signal qui limite l'impact du passage d 'une plaque entretoise 10.
Pour ce faire, on dispose d 'information sur la vitesse moyenne v de la 20 sonde, ainsi que sur la fréquence d'échantillonnage fe. Connaissant la vitesse de la sonde et sa fréquence d 'échantillonnage, chaque distance peut être convertie en un nombre d'échantillons de signal selon l
Figure imgf000011_0002
a formule :
M — distance x
J V échantillons
Figure imgf000011_0003
V
On peut par exemple chercher une fenêtre d 'observation correspondant à une longueur de M échantillons inférieure au nombre d 'échantillons théorique 25 entre deux plaques entretoises 10 successives Νρθ,ρθ, par exemple
M = Npe'pe
2
Cette recherche s'effectue par exemple sur le signal X3f.x3f présentant une forte variance lors du passage de la sonde au niveau d 'une plaque entretoise 10. I l suffit pour trouver une telle fenêtre de rechercher la portion du signal de longueur M où la variance de x3f est minimale. En notant m0 le premier indice de 30 cette fenêtre on a :
Figure imgf000012_0001
Cependant, les portions de signal correspondant au passage par les plaques entretoises 10 ne sont pas connues, puisque c'est justement ce que l'on cherche à déterminer. Ainsi, s'il devait se trouver une plaque entretoise 10 dans la fenêtre d'observation, cela pourrait grandement impacter la suite du procédé.
5 Afin de limiter l'impact du passage d'une plaque entretoise 10 dans le calcul, une fenêtre suffisamment grande est de préférence choisie. Par exemple, la fenêtre d'observation peut correspondre à un nombre d'échantillons théorique du signal de mesure représentant la distance entre des plaques entretoises 10 consécutives le long du tube 11 , c'est-à-dire de taille Μ=Νρθ,ρθ.
î o Pour limiter l'impact du passage par une plaque entretoise 10, on peut également choisir l'ensemble des indices du signal n tels que | x3f [n] | est inférieur à un certain seuil. Ce seuil peut être par exemple l'écart-type du signal x3f.
On combine alors différentes composantes du signal de mesure pour déterminer un signal de détection, la combinaison ainsi effectuée étant choisie 15 pour minimiser l'énergie du signal de détection. La combinaison des différentes composantes du signal de mesure est une combinaison linéaire adaptative, les valeurs des coefficients de pondération utilisés pour cette combinaison étant déterminés, pour un échantillon du signal de détection, par minimisation de la variance du signal de détection sur la fenêtre d'observation du signal de mesure 20 autour de cet échantillon.
De préférence, les différentes composantes du signal de mesure combinées pour déterminer le signal de détection sont des composantes réelles et imaginaires du signal de mesure. De même, différentes composantes du signal de mesure combinées pour déterminer le signal de détection sont des composantes 25 de fréquences différentes.
On cherche ainsi des coefficients de pondération α, β et γ tel que x3f[n] * a.x1f[n] + .y3f[n] + y.y1f[n]
afin de construire le signal de détection d[n] défini par
] = X3f[n]- a.x1f[n] - .y3f[n] - y.y1f[n]
30 De plus, pour gérer la forte non-stationnarité des composantes du signal de mesure, les coefficients α, β et γ sont mis à jour pour chaque échantillon de signal. Plus précisément le triplet {α[η], β[η], γ[η]} choisi à un instant n est celui qui minimise la puissance de l'erreur de reconstruction sur une fenêtre de signal centrée sur le n-ème échantillon : [k] - « if [k] - βγ3 f [k] - yyi / [k] )
Figure imgf000013_0001
L'optimisation de {α[η], β[η], γ[η]} s'effectue à l'aide de l'algorithme des 5 moindres carrés. On notera que l'on peut changer le nombre de coefficients.
On obtient ainsi un signal de détection d[n] choisi par combinaison linéaire adaptative de composantes réelles et imaginaires du signal de mesure.
Une telle construction du signal de détection présente plusieurs avantages:
- une estimation automatique du coefficient optimal de combinaison des î o fréquences (i.e. recherche du coefficient minimisant l'énergie du signal hors plaque) ;
- la combinaison peut se faire entre toutes les composantes (parties réelles et imaginaires) des différentes fréquences ;
- les coefficients de combinaison ne sont pas constants le long du tube, 15 mais varient afin de s'adapter au tube. On gère ainsi les variations de l'état du tube.
Une fois le signal de détection construit, il reste à détecter des pics du signal de détection susceptibles de correspondre au passage de plaques entretoises 10, par comparaison avec un seuil de détection. En effet, le signal de
20 détection ainsi créé contient différent pics, qui correspondent au passage de la sonde par les différentes plaques, et pas seulement les plaques entretoises 10, et d'autres éléments de structures.
Ces pics de signal peuvent être identifiés en ne gardant que les pics de signal supérieurs à un seuil de détection s.
25 De préférence, le seuil de détection pour détecter les pics du signal de détection est fonction de la valeur minimale que prend l'écart-type σ du signal de détection sur une partie du signal correspondant à un nombre d'échantillons inférieur au nombre d'échantillons théorique entre deux plaques entretoises 10 consécutives.
30 L'écart-type σ peut être déterminé de la façon suivante :
Figure imgf000014_0001
Le seuil de détection doit alors être un compromis entre le risque de fausse détections (un seuil trop faible entraîne la détection de nombreux pics ne correspondant pas au passage d'une plaque entretoise 10par la sonde) et le risque de ne pas détecter une plaque entretoise 10. Il peut par exemple être compris entre 5 et 15 σ, de préférence 10 σ.
Il se peut que le passage de la sonde par une plaque entretoise se traduise par plusieurs passages du signal au-dessus du seuil de détection. Pour ne garder qu'un seul pic par plaque, pour ensuite repérer les pics correspondant aux plaques entretoises 10, le test suivant est effectué sur les couples de pics successifs: si n1 et ni sont deux indices de pics successifs, on dira que ni est un « pic secondaire » de n1 si d[n1 ]>d[n2] et si le nombre d'échantillons entre n1 et ni est fortement inférieur à la distance minimale théorique entre deux pics de plaques.
Ainsi, préalablement à l'étape de sélection des pics, on restreint les pics détectés aux pics correspondant à un maximum local du signal de détection sur une plage du signal de détection correspondant à un nombre d'échantillons de part et d'autre du maximum local inférieur ou égal au minimum des trois valeurs suivantes :
0,5 fois un nombre d'échantillons théorique entre deux plaques entretoises 0,8 fois un nombre d'échantillons théorique entre la plaque entretoise de la branche chaude la plus proche de la branche froide dudit échangeur thermique et la plaque entretoise de la branche froide la plus proche de la branche chaude dudit échangeur thermique,
0,8 fois un nombre d'échantillons théorique entre une première des plaques entretoises et une autre plaque précédente dans le sens de circulation de la sonde.
Les pics secondaires étant supprimés, il reste alors à sélectionner des pics du signal de détection correspondant au passage des plaques entretoises 10 par une détermination des pics du signal présentant un espacement régulier.
Le nombre d'échantillons théorique entre deux plaques peut être calculé à partir d'une connaissance de la distance entre deux plaques entretoises. Idéalement, si la sonde a une vitesse connue et constante v, alors le nombre d'échantillons entre deux plaques est exactement proportionnel à la distance entre plaque:
distance entre plaques x fe
Npe,pe = ~
Dans ce cas idéal, il suffirait pour identifier les indices des échantillons correspondant de ne garder que les pics dont les indices sont exactement espacés 5 de NPe,pe. En pratique, la vitesse de la sonde n'est pas exactement constante ni précisément connue et peut varier le long du tube. L'écart entre deux pics correspondant à des plaques entretoises 10 n'est donc pas exactement égal à NPe,pe et il peut être différent d'un couple de plaques entretoises 10 à un autre. Cet écart moyen entre plaques entretoises 10 peut cependant être estimé.
î o Pour ce faire, on sélectionne les pics du signal de détection correspondant au passage des plaques entretoises 10 par la sélection d'au moins un sous- ensemble de pics minimisant la différence entre le nombre d'échantillons du signal de détection entre deux pics successifs dudit sous-ensemble et une valeur médiane du nombre d'échantillons du signal de détection entre deux pics
15 successifs.
Un tube d'un échangeur thermique s'étend le long de celui-ci à travers deux parties, généralement nommées branche froide et branche chaude, en fonction du sens d'écoulement des fluides réalisant l'échange thermique, séparées par une structure sans plaque entretoise 10 désignée comme le cintre. Par 20 conséquent, une telle structure rompt la régularité de l'espacement des plaques entretoises. Il suffit d'appeler le procédé suivant pour une branche, puis pour l'autre, pour sélectionner l'ensemble des pics correspondant aux plaques entretoises 10.
On appelle ik l'indice du k-ème pic détecté, A[k]=ik+1-ik l'écart entre les pic 25 successifs k et k+1 , et Npe,pe,réei le nombre d'échantillons moyen entre deux plaques entretoises 10. Les plaques entretoises étant les plaques majoritaires le long du tube, la majorité des valeurs prises par Δ [k] est située autour de Npe,pe,réei. Une estimation de Npe,pe,réei est donc fournie par la valeur médiane de l'ensemble des valeurs prises par A[k] :
Npe.pe.réei = mediane{A[k]}
30 II reste alors à trouver, parmi l'ensemble I des indices de pics, le sous- ensemble lpic d'indices correspondant aux nombre de plaques entretoises NbPe, tels que les écarts entre deux indices successifs soient toujours à peu près égale à
N pe,pe,réel ·
Figure imgf000016_0001
On obtient alors les indices des échantillons du signal correspondant au passage des plaques entretoises 10.
Une fois estimé l'ensemble d'indices correspondant au passage de la sonde par les plaques entretoises 10, on peut ensuite procéder à une étape subséquente de vérification dans laquelle les positions des pics du signal de détection correspondant au passage des plaques entretoises 10 sont comparées avec des données connues sur la géométrie de l'échangeur thermique. Deux critères sont utilisés : la longueur des plaques entretoises 10 et la longueur du cintre.
Le premier critère de vérification concerne la longueur du cintre : pour chaque tube, le nombre d'échantillon entre la dernière plaque de la branche chaude et la dernière plaque de la branche froide est comparé à la distance connue selon les plans. Une marge est autorisée, pour tenir compte de l'incertitude sur la vitesse de la sonde.
Le deuxième critère de validation permet de vérifier la longueur de l'ensemble des plaques détectées. L'algorithme de détection ayant permis d'estimer les indices des plaques entretoises 10 utilise les pics positifs du signal de détection. Cependant le passage de la sonde par une plaque entretoise 10 génère sur le signal au moins un pic positif et un pic négatif (les deux pics correspondent aux passages de la sonde par les bords de plaques).
Une nouvelle estimation des indices des plaques entretoises 10 peut donc être obtenue en détectant cette fois les pics négatifs du signal de détection, sans prendre en compte les premiers pics détectés. En pratique, il suffit de rappeler l'algorithme de détection des pics en prenant l'opposé du signal de détection (d[n] devient -d[n]). Le procédé peut ainsi comprendre une étape finale de vérification dans laquelle on prend l'opposé du signal de détection, on applique les étapes de détection S2 et de sélection S3 des positions des pics du signal de détection correspondant au passage des plaques entretoises 10 et on vérifie que les écarts entre les positions des pics sur le signal de détection et sur son opposé correspondent à la distance bord à bord d'une plaque entretoise 10.
On dispose donc ensuite, pour chaque branche, de deux ensembles d'indices estimés. Si ces deux ensembles d'indices sont corrects, l'écart entre les deux estimations d'une même plaque entretoise 10 doit être de l'ordre de grandeur de la taille d'une plaque entretoise 10. En tenant compte des imprécisions sur la vitesse de la sonde et des effets de bords de plaques, une marge est autorisée pour valider l'estimation.
Une fois l'identification des parties de ce signal correspondant au passage de la sonde au niveau des plaques entretoises 10, il suffit d'analyser lesdites parties du signal, par exemple par une comparaison de la variation de l'intensité ou de la puissance du signal de détection ou de mesure avec un modèle de référence et/ou un modèle calibré, pour déterminer une anomalie au niveau du passage d'une plaque entretoise 10 susceptible de correspondre par exemple à un colmatage.
L'invention concerne également un produit programme d'ordinateur comprenant des instructions de code de programme enregistrées sur un support utilisable dans un ordinateur pour l'exécution d'étapes de traitement du procédé selon l'une des revendications précédentes, lorsque ledit programme est exécuté sur un ordinateur.
De fait, le signal de mesure est transmis depuis la sonde à courant de Foucault à une mémoire pour y être stocké en vue de son traitement. Ce traitement du signal de mesure sur lequel porte la présente invention est mis en œuvre par une unité de traitement munie d'un calculateur, typiquement un ordinateur muni de moyen d'affichage et de communication, par lesquels il acquiert le signal de mesure et transmet les résultats de la mise en œuvre du procédé.

Claims

Revendications
Procédé de détection de dégradations d'un échangeur thermique à tubes, dans lequel on fait circuler dans un tube (1 1 ) dudit échangeur une sonde à courant de Foucault multifréquence et on mesure avec celle-ci un signal de mesure fonction de l'environnement dans lequel la sonde se trouve, caractérisé en ce que pour identifier les parties de ce signal correspondant au passage de la sonde au niveau de plaques entretoises (10),
- on combine (S1 ) différentes composantes du signal de mesure pour déterminer un signal de détection, la combinaison ainsi effectuée étant choisie pour minimiser l'énergie du signal de détection,
- on détecte (S2) sur ce signal des pics de détection,
- on sélectionne (S3) parmi ces pics ceux qui présentent un espacement régulier, les pics ainsi sélectionnés correspondant au passage de la sonde au niveau des plaques entretoises (10).
Procédé selon la revendication 1 , dans lequel la combinaison des différentes composantes du signal de mesure est une combinaison linéaire adaptative, les valeurs des coefficients de pondération utilisés pour cette combinaison étant déterminés, pour un échantillon du signal de détection, par minimisation de la variance du signal de détection sur une fenêtre d'observation du signal de mesure autour de cet échantillon.
Procédé selon l'une des revendications précédentes, dans lequel différentes composantes du signal de mesure combinées pour déterminer le signal de détection sont des composantes réelles et imaginaires du signal de mesure.
Procédé selon l'une des revendications précédentes, dans lequel différentes composantes du signal de mesure combinées pour déterminer le signal de détection sont des composantes de fréquences différentes.
Procédé selon l'une des revendications précédentes, dans lequel, pour détecter des pics de détection sur le signal de détection, on compare le signal de détection à un seuil de détection, ledit seuil étant fonction de la valeur minimum que prend l'écart-type du signal de détection sur une partie du signal correspondant à un nombre d'échantillons inférieur au nombre d'échantillons théorique entre deux plaques entretoises (10) consécutives.
Procédé selon l'une des revendications précédentes, dans lequel préalablement à l'étape (S3) de sélection des pics, on restreint les pics détectés aux pics correspondant à un maximum local du signal de détection sur une plage du signal de détection correspondant à un nombre d'échantillons de part et d'autre du maximum local inférieur ou égal au minimum des trois valeurs suivantes :
0,5 fois un nombre d'échantillons théorique entre deux plaques entretoises,
0,8 fois un nombre d'échantillons théorique entre la plaque entretoise de la branche chaude la plus proche de la branche froide dudit échangeur thermique et la plaque entretoise de la branche froide la plus proche de la branche chaude dudit échangeur thermique,
0,8 fois un nombre d'échantillons théorique entre une première des plaques entretoises et une autre plaque précédente dans le sens de circulation de la sonde.
Procédé selon l'une des revendications précédentes, dans lequel pour sélectionner les pics du signal de détection présentant un espacement régulier et correspondant au passage des plaques entretoises (10), on détermine au moins un sous-ensemble de pics minimisant la différence entre le nombre d'échantillons du signal de détection entre deux pics successifs dudit sous-ensemble et une valeur médiane du nombre d'échantillons du signal de détection entre deux pics successifs.
Procédé selon l'une des revendications précédentes, dans lequel le signal de mesure correspond à la mesure d'une sonde détectrice de courants de Foucault acquérant ledit signal de mesure lors du retour de la sonde dans un déplacement aller-retour de ladite sonde dans le tube (1 1 ) d'un échangeur thermique.
9. Procédé selon l'une des revendications précédentes, comprenant un étalonnage préalable (S0) du signal de mesure, comprenant :
- l'indentification (S20) des parties du signal de mesure correspondant à des défauts connus du tube (11 ) de l'échangeur thermique,
- l'estimation (S21 ) de la phase et de l'amplitude des parties du signal de mesure correspondant aux défauts connus,
- la détermination (S22) d'une transformation à appliquer au signal de mesure à partir desdites parties du signal de mesure correspondant aux défauts connus pour étalonner le signal de mesure,
- l'application (S23) de la transformation au signal de mesure.
10. Procédé selon l'une quelconque des revendications précédentes, comprenant en outre une étape subséquente de vérification dans laquelle les positions des pics du signal de détection correspondant au passage des plaques entretoises (10) sont comparées avec des données connues sur la géométrie de l'échangeur thermique.
1 1 . Procédé selon l'une quelconque des revendications précédentes, comprenant une étape finale de vérification (S4) dans laquelle on prend l'opposé du signal de détection, on applique les étapes de détection (S2) et de sélection (S3) des positions des pics du signal de détection correspondant au passage des plaques entretoises (10) et on vérifie que les écarts entre les positions des pics sur le signal de détection et sur son opposé correspondent à la distance bord à bord d'une plaque entretoise (10).
12. Produit programme d'ordinateur comprenant des instructions de code de programme enregistrées sur un support utilisable dans un ordinateur pour l'exécution d'étapes de traitement du procédé selon l'une des revendications précédentes, lorsque ledit programme est exécuté sur un ordinateur.
PCT/EP2013/064480 2012-07-09 2013-07-09 Procédé de détection du colmatage d'un échangeur thermique Ceased WO2014009366A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
FR1256584A FR2993090B1 (fr) 2012-07-09 2012-07-09 Procede de detection du colmatage d'un echangeur thermique
FR1256584 2012-07-09

Publications (1)

Publication Number Publication Date
WO2014009366A1 true WO2014009366A1 (fr) 2014-01-16

Family

ID=46826802

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/EP2013/064480 Ceased WO2014009366A1 (fr) 2012-07-09 2013-07-09 Procédé de détection du colmatage d'un échangeur thermique

Country Status (2)

Country Link
FR (1) FR2993090B1 (fr)
WO (1) WO2014009366A1 (fr)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106501356A (zh) * 2016-12-29 2017-03-15 核动力运行研究所 一种内穿式涡流检测探头对中性校验方法

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2999776B1 (fr) 2012-12-19 2015-01-09 Electricite De France Procede d'evaluation du colmatage d'un echangeur thermique
FR3074885B1 (fr) * 2017-12-13 2019-12-27 Electricite De France Procede d'evaluation du colmatage d'un echangeur thermique

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2950434A1 (fr) * 2009-09-24 2011-03-25 Electricite De France Perfectionnements a la detection de depots comportant au moins un materiau ferromagnetique sur ou a proximite de la paroi externe d'un tube

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2950434A1 (fr) * 2009-09-24 2011-03-25 Electricite De France Perfectionnements a la detection de depots comportant au moins un materiau ferromagnetique sur ou a proximite de la paroi externe d'un tube

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
C K SWORD ET AL: "Estimation of Mixing Parameters for Cancellation of Discretized Eddy Current Signals Using Time and Frequency Domain Techniques", JOURNAL OF NONDESTRUCTIVE EVALUATION, vol. 5, no. 1, 1 March 1985 (1985-03-01), pages 27 - 35, XP055055499, ISSN: 0195-9298, DOI: 10.1007/BF00568761 *
C. K. SWORD ET AL.: "Estimation of mixing parameters for cancellation of discretized eddy current signals using time and frequency domain techniques", JOURNAL OF NONDESTRUCTIVE EVALUATION, vol. 5, no. 1, pages 27 - 35, XP055055499, DOI: doi:10.1007/BF00568761
INOUE ET AL: "Vibrations of probe used for the defect detection of helical heating tubes in a fast breeder reactor", NUCLEAR ENGINEERING AND DESIGN, AMSTERDAM, NL, vol. 237, no. 8, 24 March 2007 (2007-03-24), pages 858 - 867, XP022001442, ISSN: 0029-5493, DOI: 10.1016/J.NUCENGDES.2006.11.001 *
P. BASART JOHN ET AL: "Distinguishing signals from noise in eddy current inspection of steam generator tubes", SPIE PROCEEDINGS, vol. 2944, 15 November 1996 (1996-11-15), pages 201 - 209, XP055076294, DOI: 10.1117/12.259059 *
P. XIANG ET AL.: "Automated analysis of rotating probe multi-frequency eddy current data from steam generator tubes", INTERNATIONAL JOURNAL OF APPLIED ELECTROMAGNETICS AND MECHANICS, vol. 12, pages 151 - 164, XP008158200
XIANG ET AL: "Automated analysis of rotating probe multi-frequency eddy current data from steam generator tubes", INTERNATIONAL JOURNAL OF APPLIED ELECTROMAGNETICS AND MECHANICS, IOS PRESS, AMSTERDAM, NL, vol. 12, 1 January 2000 (2000-01-01), pages 151 - 164, XP008158200, ISSN: 1383-5416 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106501356A (zh) * 2016-12-29 2017-03-15 核动力运行研究所 一种内穿式涡流检测探头对中性校验方法
CN106501356B (zh) * 2016-12-29 2019-08-13 核动力运行研究所 一种内穿式涡流检测探头对中性校验方法

Also Published As

Publication number Publication date
FR2993090B1 (fr) 2014-08-08
FR2993090A1 (fr) 2014-01-10

Similar Documents

Publication Publication Date Title
EP2936135B1 (fr) Procédé d'évaluation du colmatage d'un échangeur thermique
JP6727229B2 (ja) 粒子の分析方法
JP6091493B2 (ja) 試料に存在する成分を決定するための分光装置と分光法
CA2968487C (fr) Procede de detection et de caracterisation par ultrasons de defauts dans un materiau heterogene
FR2917500A1 (fr) Compensation automatique d'eloignement pour inspection par courants de foucault pulses
FR3057903A1 (fr) Outils a courant de foucault a champ lointain
FR2953964A1 (fr) Procede de traitement d'images obtenues par tomographie ou tomosynthese a faible nombre de projections
FR3085481A1 (fr) Procede de detection et de caracterisation par ultrasons de defauts dans un materiau heterogene
Maher et al. Determination of ideal offset for spatially offset Raman spectroscopy
WO2014009366A1 (fr) Procédé de détection du colmatage d'un échangeur thermique
EP3087444B9 (fr) Procédé d'estimation quantitative du colmatage des plaques d'un générateur de vapeur
EP2827143A1 (fr) Procédé et dispositif d'imagerie par ultrasons avec prediction des artefacts induits entre modes de reconstruction
FR2696550A1 (fr) Procédé de traitement de signaux recueillis par un capteur ponctuel absolu à courants de Foucault.
FR3071611A1 (fr) Procede et dispositif de controle non destructif d'une piece aeronautique
WO1999004253A1 (fr) Sonde a courants de foucault pour le controle non destructif de la paroi d'un tube et procede de traitement des signaux de la sonde
EP4361626A1 (fr) Procédé de contrôle non destructif d'une pluralité d'éléments
EP4362046B1 (fr) Procédé de détection de défauts dans un échangeur thermique à tubes
EP4118616B1 (fr) Procédé et système de contrôle non destructif d'une pièce aéronautique
EP4451060B1 (fr) Procédé d'identification d'un défaut affectant un motif de contrôle porté par un composant microélectronique ; procédé de contrôle, système instrumental et produit programme d ordinateur associés
US20250060318A1 (en) Methods and systems for evaluating fiber qualities
FR2951561A1 (fr) Procede de detection de la position d'un front d'onde correspondant a un evenement dans un signal recu par un detecteur
FR3140442A1 (fr) Méthode d’évaluation d’une ligne de transmission par analyse automatique d’un réflectogramme

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 13739176

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 13739176

Country of ref document: EP

Kind code of ref document: A1