US20170003215A1 - Method for quantitative estimation of fouling of the spacers plates in a steam generator - Google Patents

Method for quantitative estimation of fouling of the spacers plates in a steam generator Download PDF

Info

Publication number
US20170003215A1
US20170003215A1 US15107323 US201415107323A US2017003215A1 US 20170003215 A1 US20170003215 A1 US 20170003215A1 US 15107323 US15107323 US 15107323 US 201415107323 A US201415107323 A US 201415107323A US 2017003215 A1 US2017003215 A1 US 2017003215A1
Authority
US
Grant status
Application
Patent type
Prior art keywords
fouling
indicators
vectors
θ
database
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.)
Pending
Application number
US15107323
Inventor
Bruno RENARD
Nicolas Paul
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

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N17/00Investigating resistance of materials to the weather, to corrosion, or to light
    • G01N17/008Monitoring fouling
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B23/00Testing or monitoring of control systems or parts thereof
    • G05B23/02Electric testing or monitoring
    • G05B23/0205Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
    • G05B23/0218Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
    • G05B23/0221Preprocessing measurements, e.g. data collection rate adjustment; Standardization of measurements; Time series or signal analysis, e.g. frequency analysis or wavelets; Trustworthiness of measurements; Indexes therefor; Measurements using easily measured parameters to estimate parameters difficult to measure; Virtual sensor creation; De-noising; Sensor fusion; Unconventional preprocessing inherently present in specific fault detection methods like PCA-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/18Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • GPHYSICS
    • G21NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21CNUCLEAR REACTORS
    • G21C17/00Monitoring; Testing ; Maintaining
    • G21C17/017Inspection or maintenance of pipe-lines or tubes in nuclear installations
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F01MACHINES OR ENGINES IN GENERAL; ENGINE PLANTS IN GENERAL; STEAM ENGINES
    • F01KSTEAM ENGINE PLANTS; STEAM ACCUMULATORS; ENGINE PLANTS NOT OTHERWISE PROVIDED FOR; ENGINES USING SPECIAL WORKING FLUIDS OR CYCLES
    • F01K7/00Steam engine plants characterised by the use of specific types of engine; Plants or engines characterised by their use of special steam systems, cycles or processes; Control means specially adapted for such systems, cycles or processes; Use of withdrawn or exhaust steam for feed-water heating
    • F01K7/16Steam engine plants characterised by the use of specific types of engine; Plants or engines characterised by their use of special steam systems, cycles or processes; Control means specially adapted for such systems, cycles or processes; Use of withdrawn or exhaust steam for feed-water heating the engines being only of turbine type

Abstract

The invention relates to a method for evaluating fouling of the passages of a spacer plate (10) of a shell-and-tube heat exchanger (11), wherein said passages (12 a, 12 b) are provided along the tubes (11) to allow fluid to flow through the spacer plate (10), wherein, for each of a plurality of at least one passage (12 a, 12 b): at least one measurement of a parameter dependent on fouling or the presence of magnetite is made in the vicinity of the passage, using an eddy current probe; this measurement is used to derive at least one indicator of the fouling of said passage, characterized in that said fouling is evaluated by comparing a set of one or more fouling indicator vectors of dimension of at least two, built from the fouling indicators thus obtained with a plurality of fouling indicator vector sets contained in a database, wherein each of said indicator vector sets is associated with a quantitative descriptor of fouling.

Description

    TECHNICAL FIELD AND BACKGROUND OF THE INVENTION
  • The present invention generally relates to the field of inspecting tubes of heat exchanger with tubes. More specifically, the invention relates to a method for evaluating fouling of a spacer plate of a heat exchanger with tubes, said passages being made along the tubes and used for circulation of a fluid in said heat exchanger through said plate.
  • A steam generator generally consists of a bundle of tubes in which flows the hot fluid, and around which flows the fluid to be heated. For example, in the case of a steam generator of a nuclear power plant of the EPR type, the steam generators are heat exchangers which use the energy of the primary circuit from the nuclear reaction for transforming the water of the secondary circuit into the steam which will supply the turbine and thus produce electricity.
  • The steam generator brings the secondary fluid from a liquid water condition to the steam condition just at the limit of saturation, by using the heat of the primary water. The latter flows in tubes around which flows the secondary water. The outlet of the steam generator is the point with the highest temperature and pressure of the secondary circuit.
  • The exchange surface, physically separating both circuits, thus consists of a tubular bundle, consisting of 3,500 to 5,600 tubes, depending on the version, in which flows the primary water brought to a high temperature (320° C.) and to a high pressure (155 bars).
  • These tubes of the steam generator are maintained by spacer plates generally positioned perpendicularly to the tubes which cross them.
  • In order to let through the fluid which vaporizes, the passages of these spacer plates are foliated, i.e. their shape has lobes around the tubes. As the water passes from the liquid condition to the steam condition, it deposits all the materials which it contained. If material deposits are made in the lobes, they reduce the free passage: this is fouling, which is therefore the gradual obturation by deposits of holes intended for letting through the water/steam mixture.
  • FIG. 1 schematically illustrates a top view of a foliated passage in a spacer plate 10, in which passes a tube 11. The lobes 12 a and 12 b allow the water to cross the spacer plate 10 along the tube 11, thereby allowing flow of the water in the steam generator. A deposit 13 is visible at the lobe 12 b, fouling said lobe 12 b. The deposit may be located on the side of the tube and/or on the side of the plate.
  • Fouling leads to modifications of the flow of the water in the steam generator, and thus promotes the occurrence of excessive vibrations of the tubes, as well as to inducing significant mechanical forces on the internal structures of the steam generators. This degradation therefore has both effects on the safety and on the performances of the facilities. It is therefore indispensable to be well aware of the nature and of the time-dependent change in this degradation.
  • It is therefore sought to estimate the fouling level of these passages. This fouling level corresponds to the ratio between the blocked surface area of these passages over the total surface area of the latter. More generally this consists in quantitatively evaluating this fouling level.
  • Presently, the only non-destructive examination system which is capable of accessing the totality of the tubes/spacer plate intersections of the steam generators is the eddy current axial probe (SAX probe). Eddy currents appear in a conductive material when the magnetic flux is varied near the material. Thus a multi-frequency eddy current probe is thus circulated in a tube of said exchanger and with the latter a measurement signal is measured, depending on the environment in which the probe is found, from which it is possible to extract information as to the anomalies in the heat exchanger.
  • A variation of the magnetic induction, notably by a coil in which circulates an alternating current, generates eddy currents, for which the induced variation in the magnetic field is detected. Typically, the voltage difference generated by the variation of the impedance of the coil is measured.
  • Utilization of the measurement signals of this probe with eddy currents does not induce any extension of the standstill of the steam generator, since this eddy current probe is already used during shutdowns, notably for inspecting the integrity of the tubes of the steam generator.
  • This eddy current probe, initially intended for detecting damaging of the tubes, is also sensitive to fouling. Further, the interpretation of this signal is presently achieved manually by specialized operators, which takes a very long time, of the order of about one week of processing for analyzing a single steam generator. Further, the intervention of an operator for noting down the measurements from a piece of analysis software often gives rise to a bias which is difficult to quantify.
  • The evaluation of the fouled aspect of a foliated passage by an operator from the measurement signal is further not very reliable, as it is generally carried out empirically upon examination of the received signal.
  • The article “Tube Support Plate Blockage Evaluation with Televisual Examination and eddy Current Analysis” of L. Châtelier et al., AIP Conference proceedings, Vol. 1096, Jul. 25, 2008, pages 766, 773, describes the determination of a fouling level by means of calculation of an indicator called SAX ratio rsax, which is a scalar indicator giving the possibility of measuring the fouling level from the amplitude difference of the signals of both sides of the spacer plate. This ratio rsax is defined as being the ratio between the amplitude difference between the upper and lower edges of the spacer plate at the passage of the tube through the spacer plate, and the maximum of the two:
  • r sax = Y 1 - Y 2 max ( Y 1 , Y 2 )
  • A correlation is demonstrated between the values assumed by the SAX ratio and the fouling level determined by remote viewing examination of the passages of the spacer plates. However, the correlation is limited to fouling levels of less than 50%, and because of the lack of accuracy between the value of the SAX ratio and the fouling level, the obtained accuracy is not satisfactory, so that it is only possible to obtain a broad estimation of the fouling, per fouling classes (0-15%, 15-25%, 25-50%).
  • Further, this correlation between the SAX ratio and the estimated fouling level by remote viewing examination depends on the type of steam generator. Further, in the case of perturbations which is common since the eddy current probe reacts to all the defects, the SAX ratio integrates all these perturbations and is then not representative of the fouling. For example, a fouling failure in the vicinity of the edge of the spacer plate has an influence on the signal of the eddy current probe, therefore on the SAX ratio, since the latter does not allow discrimination of the causes of the perturbation.
  • Document EP 2 474 826 A1 proposes a method for evaluating the fouling, wherein the eddy current measurements are noted, and the signals corresponding to the passages of the spacer plates are identified and then an average value is determined, which is used as an evaluation signal. Characteristics extracted from these signals are then used as fouling indicators and a fouling level is inferred therefrom. For example, the distance between the extreme points of a Lissajous representation may be used for determining the fouling level, by means of a predetermined calibration curve relating said distance to a fouling level.
  • In another example, the comparison is carried out by means of a predetermined equation. The equation is described as having been obtained by data read-out means on a device for which the fouling characteristics are known. The fouling level corresponds to the result of the equation assuming as variables the fouling indicators.
  • These methods therefore imply predetermination of a model which is expressed by a calibration curve or by an equation. They therefore assume a simple relationship between the extracted characteristics and the fouling level, and assume that a fouling level is expressed by a same signal shape. Now, it is found that this is not the case, and that a same fouling level may give varied signals, and that consequently it is not possible to simply model the relationship between the measurement signal and the fouling level. These methods therefore do not allow proper appreciation of the fouling level.
  • Document EP 2 584 254 A2, in a similar context to the preceding one, proposes a method aiming at predicting the time-dependent change in the fouling by determining for each passage through the spacer plate, the fouling level. For this purpose, a model for viewing the spacer plates illustrating the fouling of the passages and their time-dependent change is provided. This time-dependent change is determined by means of a fouling curve which is determined by measurements via visual inspections and via eddy currents.
  • In order to determine the evaluation fouling curve, a relationship is established between the actual fouling level obtained by visual inspection of a passage and a fouling indicator from the signal of the corresponding eddy current probe, for example a distance on a Lissajous representation as earlier. Subsequent measurements with eddy currents give the possibility of readjusting the estimation of the fouling rate by estimating the fouling level by comparing the fouling indicator determined from the signal of the eddy current probe with the fouling evaluation curve.
  • The method proposed by document EP 2 584 254 A2 is therefore based on the same assumption as that of document EP 2 474 826 A1, i.e. the existence of a representative model of a simple relationship between the measurement signal and the fouling level. Now, experiment has shown that this is not the case, and therefore the proposed methods do not allow correct estimation of the fouling level.
  • PRESENTATION OF THE INVENTION
  • A general object of the invention is to overcome all or part of the defects in the methods for evaluating fouling of foliated passages around tubes in the spacer plates of the state of the art, by proposing a comparison of a vector consisting of several fouling indicators with other vectors consisting of several fouling indicators.
  • A method for evaluating the fouling of passages of a spacer plate of a heat exchanger with tubes is notably proposed, said passages being made along the tubes for crossing the spacer plate with a fluid, wherein, for each of at least one passage:
  • at least one measurement of a parameter depending on the fouling or on the presence of magnetite is conducted in the vicinity of the passage by means of an eddy current probe,
  • at least one fouling indicator of said passage is derived from this measurement,
  • characterized in that the fouling is evaluated by comparing a set of one or several vectors of fouling indicators with a dimension of at least two, built from the thereby obtained fouling indicators, with a plurality of sets of vectors of fouling indicators contained in a database, each of said sets of vectors of fouling indicators of the database being associated with a quantitative fouling descriptor.
  • This method is advantageously completed by the following characteristics, taken alone or in any of their technically possible combinations:
    • the sets of vectors of fouling indicators are represented by distributions of vectors of fouling indicators of passages of a spacer plate portion and the quantitative descriptor associated with each distribution is an average fouling level of the passages of said spacer plate portion, said database dealing with at least N portions of spacer plates of different heat exchangers, N≧2, and including N distributions of vectors of indicators each associated with an average fouling level of the passages of said spacer plate portion;
    • the distributions of vectors of fouling indicators are associated with spatial information, so as to correspond to representative images of the spatial distribution of the fouling values;
    • the method comprises steps according to which
      • the distribution of vectors of indicators Ptest(θ) of the inspected plate portion is determined,
      • a similarity measurement dn, between the distribution of vectors of indicators Ptest(θ) of the inspected spacer plate portion and each of the distributions of vectors of indicators Pn(θ) of the database is calculated,
      • K distributions of vectors of indicators Pn(θ) of the database are selected, for which the similarity measurements dn, with the distribution of indicators Ptest(θ) of the inspected spacer plate portion are the greatest,
      • the fouling is determined from fouling levels associated with said selected K distributions of vectors of indicators Pn(θ) of the database;
    • the determination of the fouling comprises a step according to which:
      • an average of the fouling levels associated with said selected K distributions of vectors of indicators Pn(θ) of the database is calculated, each fouling level being weighted by the similarity measurements between the distribution of vectors of indicators Pn(θ) of the database with which it is associated and the distribution of vectors of indicators Ptest(θ) of the inspected spacer plate portion;
    • the method further comprises a determination of an evaluation of the uncertainty on the thereby determined fouling, on the basis of the measurement of the similarity between the thereby selected K distributions of vectors of indicators Pn(θ) of the database and the distribution of vectors of indicators Ptest(θ) of the inspected spacer plate portion and/or of the variability of the quantitative fouling descriptors associated with the distributions of vectors of indicators of the selected spacer plate portion;
    • the calculation of the measurement of similarity dn, between the distribution of vectors of indicators Ptest(θ) of the inspected spacer plate portion and each of the distributions of vectors of indicators Pn(θ) of the database comprises an estimation of the distributions by means of a probability law model
      Figure US20170003215A1-20170105-P00001
      , preferably a Gaussian law, a Parzen modeling or a weighted average of probability laws;
    • a set of vectors of fouling indicators is a vector of fouling indicators of a tube and the quantitative descriptor associated with said vector is a fouling level of said tube, said database dealing with at least M tubes of different heat exchangers, M≧2, said database including M vectors of fouling indicators of a passage, each associated with a fouling level of said passage;
    • the method comprises the steps according to which:
      • the vectors of indicators θ of the inspected tube is determined,
      • the a posteriori fouling c distribution, (p(c|θ)), is calculated for the vectors of indicators θ of the inspected tube from the vectors of indicators of the database,
      • the fouling is determined by the sum of the a posteriori fouling distribution p(c|θ) weighted by the fouling;
    • the calculation of the a posteriori law comprises an estimation of the a priori law p(c) and of the likelihood p(θ|c);
    • the a priori law is determined by a ratio between:
      • the number Mk of tubes of the database having a fouling level c comprised in an interval [ck; ck+1], and
      • the total number of tubes in the database;
    • the likelihood law is approached on c comprised in an interval [ck: ck+1] by a probability law
      Figure US20170003215A1-20170105-P00001
      , preferably a Gaussian law, a Parzen modeling or a weighted average of laws;
    • the set of vectors of fouling indicators of the database are packets each having a center or an average and grouping on the basis of a similarity measurement dealing with the values of said vectors of fouling indicators, the vectors of fouling indicators which are the closest, in the sense of the similarity measurement, to said center or to said average, a quantitative fouling descriptor being associated with each of said packets, and wherein, for a set of one or several vectors of fouling indicators of the inspected tube or plate portion:
      • each of the vectors of indicators of the set of vectors of fouling indicators of the inspected tube or plate portion are each compared with the respective centers or averages of the packets of the database by a similarity measurement,
      • m packets of vectors of fouling indicators are selected on the basis of this comparison,
      • the fouling level of the inspected tube or plate portion is determined from quantitative descriptors associated with m packets of selected vectors of indicators.
  • Preferably, the fouling level of the inspected tube or plate portion is determined from an average of the quantitative descriptors of each packet weighted from the calculated similarity measurements.
  • The invention also relates to a computer program product comprising program code instructions for executing the steps of the method according to the invention when said program is executed on a computer.
  • PRESENTATION OF THE FIGURES
  • Other features, objects and advantages of the invention will become apparent from the following description, which is purely illustrative and non-limiting, and which should be read with reference to the appended drawings wherein:
  • FIG. 1, having already been commented, schematically illustrates, in a top view a foliated passage in a spacer plate, in which passes a tube, according to a common configuration of a steam generator;
  • FIG. 2 schematically illustrates the steps of the method according to a first alternative of the invention;
  • FIG. 3 schematically illustrates the steps of the method according to a second alternative of the invention.
  • DETAILED DESCRIPTION
  • The method in a way known to one skilled in the art begins by, in the vicinity of passages, conducting at least one measurement of a parameter depending on the fouling or on the presence of magnetite, typically by means of an eddy current probe, the measurement of which is representative of the impedance variations which fouling may cause, for example by magnetite.
  • Next from this measurement at least one fouling indicator of said passage is derived. The description below gives a non-limiting example of deriving such a fouling indicator.
  • After extraction from the measurement signal of a signal corresponding to the passage of the downstream edge of the spacer plate 10 by the probe, and from a signal corresponding to the passage of the upstream edge of the spacer plate 10 by the probe, it is then proceeded with the determination from the measurement signal of a lower edge signal corresponding to the passing of the downstream edge of the spacer plate 10 by the probe, and of an upper edge signal corresponding to the passage of the upstream edge of the spacer plate 10 by the probe.
  • The eddy current probe typically acquires at least partly the measurement signal in a differential mode, and the measurement signal is a multi-frequency signal consisting of at least two signals at different frequencies.
  • Preferably, only the signals corresponding to the differential mode (z1 and z3) are used since they are more sensitive to the passage of the spacer plate 10. These signals are acquired at different frequencies and the lower edge signal is determined as a linear combination of at least two signals at different frequencies of said measurement signal, in this case z1 and z3.
  • This linear combination involves a complex coefficient ∝ optimized for minimizing the signal power along the tube 11 outside the spacer plate areas 10.
  • Thus, the lower edge signal zinf is determined from signals obtained in a differential mode on the frequencies f3 and f1, so that

  • z inf [n]=z 3inf [n]−∝·z 1inf [n],

  • with ∝=argmin∥z 3 [n]−∝×z 1 [n]∥ 2,
  • for the indices n corresponding to the signal outside the spacer plate areas 10, and z3inf corresponding to the response of the probe in a differential mode on frequency f3 during the passage of the downstream edge, i.e. lower edge, of the spacer plate 10 by the probe, and z1inf corresponding to the response of the probe in a differential mode on frequency f1 during the passage of the downstream edge, i.e. lower edge, of the spacer plate 10 by the probe.
  • One proceeds in the same way with the upper edge signal, with preferably the same coefficient ∝, so that zsup [n]=z3sup[n]−∝z1sup[n], with Z3sup corresponding to the response of the probe in a differential mode on frequency f3 during the passing of the upstream edge, i.e. upper edge, of the spacer plate 10 by the probe, and z1sup corresponding to the response of the probe in a differential mode on frequency f1 during passing of the upstream edge, i.e. upper edge, of the spacer plate 10 by the probe.
  • Thus two complex signals are obtained. The lower edge signal zinf is written as:

  • z inf [n]=x inf [n]+i. y inf [n]
  • with xinf and yinf the respectively real and imaginary components of the lower edge signal and i the imaginary unit such that i2=−1. Also, the upper edge signal zsup is written as:

  • z sup [n]=x sup [n]+i. y sup [n]
  • with xsup and ysup being the respectively real and imaginary components of the upper edge signal and i being the imaginary unit such that i2=−1.
  • An adequate processing of these signals therefore remains to be applied in order to evaluate the fouling of the passage of the spacer plate 10. This processing is applied on the lower edge signal, which is a complex signal. Indeed, the fouling of the foliated passages, i.e. the lobes 12 a, 12 b in the spacer plates 10 occurs at the lower edge of the spacer plates 10, upstream from the passage for the fluid flow passing through the spacer plate 10. It is therefore from the lower edge signal that it is possible to estimate the fouling level.
  • More specifically, the lower edge signal is de-convoluted by the complex impulse response of the probe.
  • In fact, in an ideal case of a perfect SAX probe, the signal should only contain a sequence of complex pulses, corresponding to the passage through a spacer plate edge 10, to the encounter of a deposit, and the study of the sole lower edge signal should be sufficient for quantifying the fouling.
  • However, in practice, the response of the SAX probe to an impedance variation is not perfect. This is called the impulse response of the probe. Therefore it is necessary to restore the lower edge signal in order to again find the response of the probe representative of the fouling condition of the foliated passage in the spacer plate 10.
  • For this purpose an estimation of the impulse response of the probe is determined, preferably corresponding to the passage of a specific edge of the spacer plate 10 by the probe in the tube 11, for example from the upper edge signal. It is then sought to de-convolute the lower edge signal zinf[n] with a signal h[n] corresponding to the impulse response of the probe upon passing the spacer plate.
  • To do this, it is possible to use a filter. Such a filter is called a deconvolution filter or further a restoration filter. The deconvolution filter is calculated from the estimation of the impulse response, and it applies a deconvolution of the lower edge signal by means of said deconvolution filter. The deconvolution filter may be an approximation of the reciprocal of the impulse response of the probe. It may also be a Wiener filter and the deconvolution may thus be a Wiener deconvolution, which is a preferential embodiment of the described method. Other deconvolution methods exist and may be used.
  • For example, it is possible to search for the de-convoluted lower edge signal zinf,id[n] which at best corresponds to the lower edge signal zinf[n] which was observed:

  • z inf,id [n]=argminz[n]{ J 1(z inf [n]
    Figure US20170003215A1-20170105-P00002
    h[n])+λ×J 2(z[n])}
  • with J1 being the criterion suitable for the data (for example a standard L2, a standard L2 squared, a standard L1, . . . ) and J2 a criterion expressing a characteristic a priori known on the signal which one seeks to rebuild (for example a standard L2, a standard L2 squared, a standard L1, a function of the deviations between neighboring samples z[n]−z[n−1]). The term λ gives the possibility of assigning more or less importance to the a priori on the solution (J2) relatively to the suitability to the data (J1). This criterion may also be written in the frequency domain.
  • Therefore there are several alternative deconvolution criteria J1 and J2 which may be used, and for each alternative, several resolution methods, for example by filtering or by optimization methods.
  • In the case when the deconvolution filter is a Wiener filter, the frequency response of the Wiener filter is of the form:
  • G [ f ] = H * [ f ] H [ f ] 2 + B [ f ] S [ f ]
  • with the exponent * designating the complex conjugation, H[f] the Fourier transform of the impulse response of the probe, S[f] the spectrum power density of the signal to be estimated and B[f] the spectrum power density of the noise. Zero-padding, i.e. adding zeros within the signals, may be carried out during the calculation of the discrete Fourier transforms in order to increase the frequency resolution.
  • The impulse response h[n] of the probe may be estimated from the response of the probe upon passing the upstream edge of the spacer plate 10 by the probe, i.e. by means of the upper edge signal, according to the formula:

  • h[n]=−z sup [−n]
  • For example, from processing operations carried out for extracting the useful portions of the measurement signal, indices iinf and isup of the measurement signal respectively corresponding to the passages of the lower and upper edges of the spacer plate 10 are known. For a sampling frequency Fe=1,000 Hz, a velocity of the probe v=0.5 m.s−1 and a spacer plate length 10 of 30 mm, one has 60 signal samples corresponding to the spacer plate 10, and an impulse response of about 20 samples. It is then possible to select for the range of values of the upper edge signal zsup[n], the 60 samples according to the center of the spacer plate 10 determined to be about 0.5×(iinf+isup), i.e. a margin of 20 samples on each side of the impulse response. These figures are of course indicated as a non-limiting example of the use of the upper edge signal zsup[n] for estimating the impulse response of the probe.
  • Several approaches are possible for estimating the noise-to-signal ratio corresponding to the ratio of the power spectral density of the noise B[f] and of the spectral power density S[f] of the signal to be estimated. One of these approaches consists of approximating this ratio with a constant. Indeed, the signal to be estimated corresponds to an ideal lower edge signal which would have a sequence of pulses corresponding to the complex impedance variations encountered by the probe in the vicinity of the lower edge of the spacer plate 10. Consequently, the power spectral density S[f] of this signal may be considered as a constant. The power spectral density of the noise B[f] may be determined on the portions of the signal between the spacer plates 10. The latter may be assimilated to white noise, and therefore this power spectral density of the noise B[f] may be considered as a constant. Thus, the ratio of the power spectral densities of noise and of the signal to be estimated may be considered as a constant. This constant may be adjusted empirically, for example by assuming:
  • B [ f ] S [ f ] = 10 σ 2 ,
  • with σ2 the power of the noise, calculated on an area outside the plates.
  • Once the deconvolution filter is determined, it is then possible to proceed with deconvolution of the lower edge signal by means of said deconvolution filter. The deconvolution filter g is then applied to the lower edge signal zinf in order to obtain a complex de-convoluted lower edge signal zinf id introduced by the impulse response of the probe:

  • z inf id =z inf * g
  • In practice, this operation may be carried out in the frequency domain:

  • z inf id =TF {Z inf [f]×G[f]},
  • with Zinf[f] the Fourier transform of the lower edge signal zinf, G[f] being the Fourier transform of the deconvolution filter g, and TF−1 indicating the inverse Fourier transform.
  • In order to avoid too substantially amplifying certain frequencies only corresponding to measurement noise, filtering with a low-pass filter is applied to the de-convoluted lower edge signal, the cutoff frequency of said low-pass filter being determined by means of a standard deviation of a Gaussian function forming an approximation of the real portion of a lower edge signal pulse corresponding to the passing of the edge of the spacer plate 10.
  • Indeed, the real or imaginary portion of the lower edge signal impulse corresponding to the passing of the edge of the spacer plate 10 has shapes very close to Gaussian functions or to their derivatives. For example, the impulse 0 in the real portion of the lower edge signal corresponding to the passing of the lower edge of the spacer plate 10 in a configuration without fouling may be assimilated to a Gaussian function, and a linear combination of derivatives of the Gaussian function may be assimilated to pulses in the imaginary portion of the lower edge signal corresponding to the passing of the lower edge of the spacer plate 10 in a fouled configuration.
  • If σ is the standard deviation of this Gaussian function, generally of the order of 3 or 4 samples, the Fourier transforms of the signals to be de-convoluted no longer contain any energy beyond a maximum frequency fmax:
  • f max = 3 2 πσ .
  • Therefore it is possible to select this maximum frequency fmax as a cutoff frequency of the low-pass filter.
  • Once the lower edge signal is de-convoluted and thus filtered, there remains analysis of the latter for evaluating the fouling. Following the process described above, during which at least one measurement of a parameter depending on the fouling or on the presence of magnetite, of fouling indicators of said passage which are derived from this measurement are then available.
  • Diverse types of indicators may be used. For example, if one designates by y+ (respectively y−) the positive values (respectively negative values) assumed by the imaginary portion of the signal obtained in the vicinity of the lower edge of the plate after deconvolution, and if the various following quantities are defined as:
  • EY+/EY−:energy of y+ and of y−
  • PY+/PY−:power of y+ and y−
  • MY+/MY−: maximum value of y+ and of |y−|
  • ΓY+Y−: standard deviation of the values assumed by y+ and by y− It is also possible to take the minimum and maximum values of the quantities below, for example for each pair of physical quantities, XY+/XY−, with X corresponding to E, P, M or Γ, and a minimum indicator and a maximum indicator may be defined:
  • Xmin=min {XY+, XY−}
  • Xmax=max {XY+, XY−}.
  • The fouling is then evaluated by comparing a set of one or several vectors of fouling indicators, of the dimension of at least two, built from the thereby obtained indicators, with a plurality of sets of vectors of fouling indicators contained in a database, each of said sets being associated with a quantitative fouling descriptor. A vector of fouling indicators is preferably of a dimension of at least two, i.e. it is not preferably a scalar.
  • Each vector of indicators is of a dimension of at least two, which means that it comprises at least two indicators as components. For example, it is possible to build an vector of indicators comprising as components:
      • the energy EY+ of the positive values assumed by the imaginary portion of the signal obtained in the vicinity of the lower edge of the plate after deconvolution, and
      • the energy EY− of the negative values assumed by the imaginary portion of the signal obtained in the vicinity of the lower edge of the plate after deconvolution.
        The vector of indicators is then written as (EY+;Ey−). Other vectors of indicators may be used, combining two or more indicators,
  • The fouling indicators contained in a database are typically fouling levels from remote viewing examination (RVE). The design of the steam generators actually allows inspection of their upper spacer plate via an automated camera. On each photograph, one of the foliages of the tube/plate intersection is observed. The fouling level of the observed foliated passage is evaluated by measuring the section reduction at the lower edge. By considering that the obstruction phenomenon is homogeneous on each of the foliages of the tube/plate intersection, the fouling level of the passage of the latter is obtained.
  • This method, unlike the examination by an eddy current probe; gives the possibility of having a quantitative indication of the fouling level of the plate, which forms a quantitative fouling descriptor. But it is only applicable to the upper spacer plate, except for certain types of steam generators, also allowing the passing of the image-shooting apparatus over a few tubes of the intermediate plate. The RVE therefore does not give the possibility of obtaining the profile of the fouling on the whole of the steam generator. On the other hand, on an inspected plate, the peripheral tubes remain inaccessible.
  • These remote viewing examinations have however been carried out for a long time, and the results have been stored in databases, so that a large amount of data exists which may be utilized for inferring therefrom the quantitative evaluation of the fouling by combining the quantitative evaluations of the remote viewing examinations passed with an inspection of the tubes with eddy currents.
  • It should be noted that the steps of the following methods are applied with at least one computer, one central or processing unit, an analog electronic circuit, a digital electronic circuit, a microprocessor, and/or software means.
  • Estimation on the Whole of a Plate Portion
  • In a first alternative of the method, it is sought here to estimate directly the average fouling level per spacer plate portion, without requiring an evaluation of the fouling of each of its tubes.
  • The sets of fouling indicators of the database are therefore represented by distributions of vectors of fouling indicators of passages of a spacer plate portion and the quantitative descriptor associated with each distribution is an average fouling level of each of the passages of said spacer plate portion, said database dealing with at least N portions of spacer plates of different heat exchangers, N≧2, and including N distributions of vectors of indicators each associated with an average fouling level of the passages of said spacer plate portion.
  • In the example hereafter, the spacer plate portions are half-spacer-plates, corresponding to the portion of the spacer plates present in the cold or hot branch of the heat exchanger, here a steam generator. It is therefore assumed that one has N half-plates for which the fouling level is determined from the remote viewing examination Cn and the values of the vectors of indicators θ for each of their tubes.
  • Approach by Similarity Measurement Between Distributions of Vectors of Indicators
  • A first approach lies on similarity measurements between distributions of vectors of indicators. The principle of this approach is to recognize in the available database, the distributions of vectors of indicators the most similar to the one of the half-plate which is intended to be evaluated. The Ptest(θ), the distribution of the vectors of indicators of the inspected half-plate, and Pn,(θ) (n ∈ [[1, N]]) the distributions of the vectors of indicators of N half-plates available in the database are considered.
  • Thus, after having determined the distribution of the vectors of indicators Ptest(θ) of the inspected plate portion (step S20), a similarity measurement dn, is calculated between the distribution of vectors of indicators Ptest(θ) of the inspected spacer plate portion and each of the distributions of vectors of indicators Pn(θ) of the database (step S21).
  • The similarity measurement between distributions may for example be evaluated by means of a distance function which will be noted as D. The distance dn is thus calculated between the distribution of vectors of indicators of the plate to be evaluated Ptest(θ) and each of the distributions Pn(θ) of the database:

  • d n =D(P test , P n)
  • There exist several similarity measurements which may be used. It is notably possible to use for example the Kullback-Leibler divergence, the Bhattacharyya distance or further the Hellinger distance. The latter in particular provides the advantage of giving a result limited between 0 and 1, thus proving to be interpretable in absolute terms. Its formula is given by the following equation:

  • D H(P, Q)=1/√{square root over (2)}√{square root over (Σi(√{square root over (p i)}−√{square root over (q i)})2)}
  • with p and q the vectors of indicators of the distributions of vectors of indicators P and Q, respectively.
  • It is then possible to select (step S22), from among the N thereby calculated distances, the K smallest ones (K ∈ [[1, NJ]]), therefore corresponding to the most similar distributions. The K distributions of indicators Pn,(θ) of the database are thus selected, the similarity measurements of which with the distribution of indicators Ptest(θ) of the inspected spacer plate portion are the greatest, are therefore selected. When the similarity measurement is a distance, this therefore amounts to taking the smallest K distributions.
  • The fouling (step S23) is then determined from the fouling levels associated with said selected K distributions of vectors of indicators Pn,(θ) of the database.
  • It is possible to calculate an average of the fouling levels associated with said selected K distributions of vectors of indicators Pn,(θ) of the database, each fouling level being weighted by the similarity measurement between the distribution of vectors of indicators Pn,(θ) of the database with which it is associated and the distribution of vectors of indicators Ptest(θ) of the inspected spacer plate portion.
  • The average of the foulings of these K half-plates are then weighted with their respective distances with the distribution of vectors of indicators of the half-plate to be evaluated, so as to give more weight to the most similar ones:
  • c ^ test = Σ k = 1 K d k - 1 c k Σ k = 1 K d k - 1
  • The calculation of the similarity measurement dn, between the distribution of vectors of indicators Ptest(θ) of the inspected spacer plate portion and each of the distributions of vectors of indicators Pn(θ) of the database may comprise an estimation of the distributions by means of a model of a probability law
    Figure US20170003215A1-20170105-P00001
    , preferably a Gaussian law, a Parzen modeling or a weighted average of probability laws.
  • Indeed, measuring distances between distributions however requires estimation of the latter. They may be then approached with a probability law
    Figure US20170003215A1-20170105-P00001
    model, a function of parameters cu to be determined, built according to the observations which are available:

  • P(θ)≈
    Figure US20170003215A1-20170105-P00001
    (θ, ω)
  • The question of the selection of
    Figure US20170003215A1-20170105-P00001
    is open: the latter may for example deal with a multidimensional Gaussian law:

  • P(θ)≈N(θ, μ, Σ)
  • wherein μ represents the average vector of indicators of the half-plate, and Σ is its variance-covariance matrix (both of these elements forming the parameters of the distribution to be calculated).
  • Also it is possible to turn towards Parzen modeling. Its principle is to place a kernel function, for example a Gaussian on each of the observations of the statistical population for which one seeks to estimate the probability density. The sum of all these Gaussians gives the Parzen likelihood.
  • In the relevant case, the observations are the vectors of indicators extracted from each of the tubes. The distribution of the vectors of indicators of the half-plate n including Mn, tubes of vectors of indicators θm (m ∈ [[1, Mn]]) of dimension d is then given by:
  • P n ( θ ) = 1 Σ n ( 2 π ) d 2 h d m = 1 M n exp ( - 1 2 h 2 ( θ - θ m ) T Σ n - 1 ( θ - θ m ) )
  • The parameter h represents the variance of the Gaussian (i.e. its width) to be applied on each of the observations. Usually the same value is used on each of the observations. For example it is possible to assume that h=0.5. The variance-covariance matrix Σn is calculated on the whole of the vectors of indicators belonging to the half-plate n. This inverse matrix may be approached by the matrix formed with its sole inverse diagonal coefficients, according to
  • Σ n - 1 = ( σ i , j ) - 1 ( 1 σ 1 , 1 0 0 1 σ d , d )
  • with i, j ∈ [[1, d]].
  • Moreover, the similarity measurement between distributions gives the possibility of introducing an uncertainty score of the result: small distances between a half-plate to be evaluated and certain half-plates of the database predict strong similarity with fouling examples already observed, and therefore high confidence in the result. Conversely, an atypical half-plate differing from the available history, will see its distribution of vectors of indicators moved away from all the others, expressing greater uncertainty on the result of the estimation.
  • The method may thus comprise a determination of the uncertainty evaluation on the thereby determined fouling, on the basis of the similarity measurement between the K distributions of vectors of indicators Pn(θ) of the database which were thereby selected and the distribution of vectors of indicators Ptest(θ) of the inspected spacer plate portion, and/or the variability of the quantitative descriptors associated with the distributions of vectors of indicators of a spacer plate portion, i.e. the half-plates.
  • Approach by Vector Quantification
  • Another approach lies on vector quantification. The principle of vector quantification is to partition a large number of data (vectors of a given space) into a restricted number of packets (or “cluster”) in the sense of a similarity measurement (generally a distance). Thus the space of the vectors of indicators θ is separated into K packets, K being a parameter of the algorithm determined beforehand, each including a center or an average. Each vector of indicators θ belongs to the cluster for which the center or the average is the closest.
  • Thus, the sets of vectors of fouling indicators of the database are then packets having each a center or an average, and grouping on the basis of a similarity measurement dealing with the vectors of fouling indicators, the vectors of fouling indicators which are the closest to said center or said average, in the sense of the similarity measurement, a quantitative fouling descriptor being associated with each of said packets, and in which, for a set of one or several vectors of fouling indicators of the inspected tube or plate portion:
      • each of the vectors of indicators of the set of vectors of fouling indicators of the inspected tube or plate portion is compared with the respective centers or averages of the packets of the database,
      • m packets of vectors of fouling indicators are selected on the basis of this comparison,
      • the fouling level of the inspected tube or plate portion is determined from quantitative descriptors associated with the m packets of selected vectors of indicators.
  • The partitioning of the space of the indicators may be obtained in diverse ways. It may be set arbitrarily a priori, or may use an algorithm which gives the packets or clusters, being more naturally disengaged from the set of vectors of indicators of the database (such as the so-called “K-means” algorithm). Once this partitioning has been carried out, for the half-plate to be inspected the number of tubes making it up is determined, belonging to each partition of the space of the vectors of indicators. The vector r=(r1, . . . , rk) results from this wherein the rk component corresponds to the proportion of tubes of the half-plate belonging to the cluster k (one therefore has Σkrk=1).
  • The estimation is then directly carried out on these vectors r, by for example searching for the N half-plates of the database for which the vectors rn, (n ∈ [[1, N]]) are the closest to the vector of the half-plate to be estimated in the sense of a given similarity measurement, typically a distance. Finally it is possible to calculate the average fouling by the average of the fouling cn, of the selected half-plates weighted by the calculated distances dn:
  • c = Σ n = 1 N d n - 1 c n Σ n = 1 N d n - 1
  • It should be noted that the distributions of vectors of fouling indicators may be associated with spatial information, such as the position of tubes in the spacer plate, so as to correspond to the representative images of the spatial distribution of the fouling values. In this case, the fouling is estimated by means of an image recognition process, which may resume the principles stated above, in order to estimate the fouling by recognizing in the database, the image(s) the closest to those obtained for the inspected heat exchanger.
  • Tube by Tube Estimation
  • According to another alternative, a set of vectors of fouling indicators is a vector of fouling indicators of a tube and the quantitative descriptor associated with said vector is a fouling level of said tube, said database dealing with at least M tubes from different heat exchangers, M≧2, said database including M vectors of fouling indicators of one passage, each associated with a fouling level of said passage of said tube.
  • In the example hereafter, and as earlier, the spacer plate portions are spacer half-plates, corresponding to the portion of the spacer plates present in the cold or hot branch of the heat exchanger, here a steam generator. For a given tube/plate intersection, the fouling level obtained for example by remote viewing examination is noted as c and θ is the vector of the qualitative indicators calculated earlier, after having determined the vectors of indicators θ of the inspected tube. Therefore M pairs {θ, c} are available in the database.
  • Probabilistic Approach
  • According to this approach, after having determined the vectors of indicators of the tube (S30), the a posteriori fouling level c distribution p(c|θ) for the vectors of indicators θ is calculated from vectors of the database (step S31), and the fouling is determined by the sum of the a posteriori fouling distribution weighted with fouling levels (step S32).
  • Indeed, this so-called a posteriori least squares approach consists of minimizing the average quadratic estimation error defined by
  • ε ^ = c ( c est - c ) 2 p ( c θ )
  • wherein p(c|θ) designates the fouling distribution of the vector of indicators θ (this is the a posteriori law). The estimator cest of the fouling minimizing the previous equation is given by the equation:
  • c est = c cp ( c θ )
  • The a posteriori law may be given by Bayes theorem:

  • p(c|θ)∝p(θ|c)p(c)
  • wherein p(c) designates the a priori law and p(θ|c) is the likelihood of the indicators in the Bayes theory. The calculation of the a posteriori law may thus comprise an estimation of the a priori law p(c) and of the likelihood p(θ|c).
  • In order to express both of these probabilities, the interval of the possible fouling levels (from 0 to 100%) may be sampled in several consecutive windows Ck=[ck; ck+1] and p(c) and p(θ|c) calculated in each of the latter. Thus, the a priori law is obtained according to a ratio between:
  • the number Mk of tubes of the database having a fouling level c comprised in the interval [ck; ck+1], and
  • the total number of tubes in the database:
  • p ( c C k ) = N k Card ( database )
  • wherein Nk represents the number of tubes of the database having fouling c comprised in the interval Ck, and Card(database) is the total number of tubes in the database.
  • The likelihood law is approached on c comprised on an interval [ck; ck+1] by a probability law
    Figure US20170003215A1-20170105-P00001
    , preferably a Gaussian law, a Parzen modeling, or a weighted average of laws, depending on parameters cu to be determined, like what was indicated above:

  • p((θ)|c)≈
    Figure US20170003215A1-20170105-P00001
    (θ, ω)
  • For example, within the scope of Parzen modeling:
  • ( θ c ε C k ) = 1 Σ n ( 2 π ) d 2 h d j = 1 N k exp ( - 1 2 h 2 ( θ - θ j ) T Σ k - 1 ( θ - θ j ) )
  • Like in the case of the estimation per half-plate, the value h of the width of the Gaussians (for example here again h=0.5) should be set. In order to counter possible conditioning problems of the variance-covariance matrix of the vectors of indicators of the tubes belonging to the class Ck, Σk, is also possible here to approach its inverse matrix with the matrix formed of its sole inverse diagonal coefficients, according to
  • Σ k - 1 = ( σ i , j ) - 1 ( 1 σ 1 , 1 0 0 1 σ d , d )
  • with i, j ∈ [[1, d]].
  • This calculation of the likelihood gives the possibility of rewriting the equation of the a posteriori law and of obtaining the expression of the a posteriori probability law in the interval Ck:
  • p ( c ε C k θ ) = p ( θ c ε C k ) p ( c ε C k ) Σ i p ( θ c ε C i ) p ( c ε C i ) = ( θ , ω k ) p ( c ε C k ) Σ i ( θ , ω i ) p ( c ε C i )
  • The estimation cest of the fouling (step S33) for the vectors of indicators θ of the inspected tube is inferred from the previous equation:
  • c est - c cp ( c θ )
  • and is given by
  • c est = c c i p ( c ε C i θ )
  • wherein <ci> symbolizes the average fouling of the tubes belonging to the interval Ci.
  • Approach by Vector Quantification
  • The use of vector quantification techniques has the purpose of evaluating the fouling level of a tube according to the position of its vectors of indicators θ in the space of the vectors of indicators, relatively to the examples of vectors which make up the database. Schematically, the principle is to give to the inspected tube for which the fouling is to be evaluated, a fouling level similar to those of the tubes of the base which are close to its vectors of indicators θ.
  • As earlier, the sets of vectors indicating fouling of the database are packets each having a center or an average and grouping on the basis of a similarity measurement dealing with the values of said fouling indicators, the fouling indicators for which the values are the closest to said center or to said average, a quantitative descriptor of fouling being associated with each of said packets.
  • More specifically, the set of the vectors of indicators of the database is thus partitioned into K sets, for example with the so-called “K-means” algorithm. Once this step is carried out, each of the vectors of indicators of the set of fouling indicators are each compared with the respective centers or averages of the packets of the database, by means of a similarity measurement, for example by calculating the distances dk of the vectors of indicators θ of the inspected tube to each of the centers of the clusters (k ∈ [[1, K]]). Several types of distances may be used, from the customary Euclidean distance to similarity measurements taking into account the distribution of the data. From the latter, mention will be made of the Mahalanobis distance, which involves the variance-covariance matrix)Σk of the partition k and its center μk. The Mahalanobis distance between a vectors of indicators θ and the set of the data of the cluster k is then written as:

  • d k=√{square root over ((θ−Ξk)TΣk −1(θ−μk))}
  • Then there remains the estimation of the value of the fouling level of the tube according to these calculated distances. A solution is to proceed with an averaging of the average foulings <ck> of each packet weighted from similarity measurements, for example the reciprocal of the distances such as the distances dk calculated previously:
  • c est = Σ k = 1 K d k - 1 c k Σ k = 1 K d k - 1
  • It should be noted that the average may be calculated on the basis of the fouling levels of the totality of the packets, or on a selection of m of them on the basis of a comparison in the sense of a similarity measurement, as previously. It is therefore considered that the selection of the m packets may comprise all the packets, the closest packets, or all the packets except for certain of them which are set aside because of anomalies.
  • It will be noted that unlike the probabilistic approach, this estimation alternative by vector quantification does not make any a priori on the set of available data. Indeed, at any moment, the a priori probability law of the fouling, p(c), is not involved. Now, the latter may prove to be strongly biased if certain ranges of fouling values are over-represented (or under-represented). This method is thus less dependent on the representativity of the database.

Claims (15)

  1. 1. A method for evaluating the fouling of passages of a spacer plate of a heat exchanger with tubes, said passages being made along the tubes for crossing the spacer plate by a fluid, wherein, for at least one passage:
    at least one measurement of a parameter depending on the fouling or on the presence of magnetite is conducted in the vicinity of the passage by means of an eddy current probe,
    at least one fouling indicator of said passage is derived from this measurement, wherein the fouling is evaluated by comparing a set of one or several vectors of fouling indicators with a dimension of at least two, built from the thereby obtained fouling indicators, with a plurality of sets of vectors of fouling indicators contained in a database, each of said sets of vectors of indicators being associated with a quantitative fouling descriptor.
  2. 2. The method according to claim 1, wherein the sets of vectors of fouling indicators are represented by distributions of vectors of fouling indicators of passages of a spacer plate portion and the quantitative descriptor associated with each distribution is an average fouling level of the passages of said spacer plate portion, said database dealing with at least N portions of spacer plates of different heat exchangers, N≧2 and including N distributions of vectors of indicators each associated with an average fouling level of the passages of said spacer plate portion.
  3. 3. The method according to claim 2, wherein:
    the distribution of vectors of indicators Ptest(θ) of the inspected plate portion is determined,
    a similarity measurement dn, between the distribution of vectors of indicators Ptest(θ) of the inspected spacer plate portion and each of the distributions of vectors of indicators Pn(θ) of the database is calculated,
    the K distributions of vectors of indicators Pn(θ) of the database are selected, for which the similarity measurements dn with the distribution of vectors of indicators Ptest(θ) of the inspected spacer plate portion are the greatest,
    the fouling is determined from fouling levels associated with said selected K distributions of vectors of indicators Pn(θ) of the database.
  4. 4. The method according to claim 3, wherein the determination of the fouling comprises a step according to which:
    an average of the fouling levels associated with said selected K distributions of vectors of indicators Pn(θ) of the database is calculated, each fouling level being weighted by the similarity measurements between the distribution of vectors of indicators Pn(θ) of the database with which said fouling level is associated and the distribution of vectors of indicators Ptest(θ) of the inspected spacer plate portion.
  5. 5. The method according to claim 3, further comprising a determination of an evaluation of the uncertainty on the thereby determined fouling, on the basis of the measurement of the similarity between the thereby selected K distributions of vectors of indicators Pn(θ) of the database and the distribution of vectors of indicators Ptest(θ) of the inspected spacer plate portion and/or of the variability of the quantitative descriptors associated with the selected distributions of vectors of indicators of a spacer plate portion.
  6. 6. The method according to claim 3, wherein the calculation of the similarity measurement dn between the distribution of vectors of indicators Ptest(θ) of the inspected spacer plate portion and each of the distributions of vectors of indicators Pn(θ) of the database comprises an estimation of the distributions by means of a probability law model
    Figure US20170003215A1-20170105-P00003
    , preferably a Gaussian law, a Parzen modeling or a weighted average of probability laws.
  7. 7. The method according to claim 2, wherein the distributions of vectors of fouling indicators are associated with spatial information, so as to correspond to representative images of the spatial distribution of the fouling values.
  8. 8. The method according to claim 1, wherein a set of vectors of fouling indicators is a vector of fouling indicators of a tube and the quantitative descriptor associated with said vector is a fouling level of said tube, said database dealing with at least M tubes of different steam generators, M≧2, said database including M vectors of fouling indicators of a passage, each associated with a fouling level of said passage.
  9. 9. The method according to claim 8, wherein:
    the vectors of indicators θ of the inspected tube is determined,
    the a posteriori fouling level distribution p(c|θ) is calculated for the vectors of indicators θ, from vectors of indicators of the database,
    the fouling is determined by the sum of the a posteriori fouling distribution p(c|θ) weighted by the fouling levels c.
  10. 10. The method according to claim 9, wherein the calculation of the a posteriori law comprises an estimation of the a priori law p(c) and of the likelihood p(θ|c).
  11. 11. The method according to claim 10, wherein the a priori law is determined by a ratio between:
    the number Mk of tubes of the database having a fouling level c comprised in an interval [ck; ck+1], and
    the total number of tubes in the database.
  12. 12. The method according to claim 9, wherein the likelihood law is approached on c comprised in an interval [ck; ck+1] by a probability law
    Figure US20170003215A1-20170105-P00001
    , preferably a Gaussian law, a Parzen modeling, or a weighted average of laws.
  13. 13. The method according to claim 1, wherein the set of vectors of fouling indicators of the database are packets of vectors of indicators each having a center or an average and grouping, on the basis of a similarity measurement dealing with said vectors of fouling indicators, the vectors of fouling indicators which are the closest, in the sense of the similarity measurement, to said center or to said average, a quantitative fouling descriptor being associated with each of said packets, and wherein, for a set of one or several vectors of fouling indicators of the inspected tube or plate portion:
    each of the vectors of indicators of the vectors of fouling indicators of the inspected tube or plate portion is compared with the respective centers or averages of the packets of the database by a similarity measurement,
    m packets of vectors of fouling indicators are selected on the basis of this comparison,
    the fouling level of the inspected tube or plate portion is determined from quantitative descriptors associated with the m packets of selected vectors of indicators.
  14. 14. The method according to claim 13, wherein the fouling level of the inspected tube or plate portion is determined from an average of the quantitative descriptors of each packet weighted from the calculated similarity measurements.
  15. 15. A computer program product comprising program code instructions for executing the steps of the method according to claim 1 when said program is executed on a computer.
US15107323 2013-12-23 2014-12-23 Method for quantitative estimation of fouling of the spacers plates in a steam generator Pending US20170003215A1 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
FR1363489A FR3015757A1 (en) 2013-12-23 2013-12-23 Method for quantitative estimation of sealing plates of a steam generator
FR1363489 2013-12-23
PCT/EP2014/079163 WO2015097221A3 (en) 2013-12-23 2014-12-23 Method for quantitative estimation of fouling of the spacers plates in a steam generator

Publications (1)

Publication Number Publication Date
US20170003215A1 true true US20170003215A1 (en) 2017-01-05

Family

ID=50543185

Family Applications (1)

Application Number Title Priority Date Filing Date
US15107323 Pending US20170003215A1 (en) 2013-12-23 2014-12-23 Method for quantitative estimation of fouling of the spacers plates in a steam generator

Country Status (9)

Country Link
US (1) US20170003215A1 (en)
EP (1) EP3087444B9 (en)
JP (1) JP2017502293A (en)
KR (1) KR20160102049A (en)
CN (1) CN106030428A (en)
CA (1) CA2934949A1 (en)
ES (1) ES2655520T3 (en)
FR (1) FR3015757A1 (en)
WO (1) WO2015097221A3 (en)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030059837A1 (en) * 2000-01-07 2003-03-27 Levinson Douglas A. Method and system for planning, performing, and assessing high-throughput screening of multicomponent chemical compositions and solid forms of compounds
US20040039968A1 (en) * 2000-09-11 2004-02-26 Kimmo Hatonen System, device and method for automatic anomally detection
US20050011278A1 (en) * 2003-07-18 2005-01-20 Brown Gregory C. Process diagnostics
US20060111844A1 (en) * 2004-11-12 2006-05-25 Chandler Larry S Implementing two dimensional segment inversions with inversion-conforming data sets processing being rendered to include generalized composite weight factors in the processing of error-affected multivariate data samples
US20070276270A1 (en) * 2006-05-24 2007-11-29 Bao Tran Mesh network stroke monitoring appliance
US7428514B2 (en) * 2004-01-12 2008-09-23 Honda Research Institute Europe Gmbh System and method for estimation of a distribution algorithm
US20120179402A1 (en) * 2011-01-06 2012-07-12 Mitsubishi Heavy Industries, Ltd. Deposit measurement apparatus, deposit measurement method, and computer-readable storage medium storing deposit measurement program
US20130041636A1 (en) * 2011-08-12 2013-02-14 Anna CRAIG Methods and systems for comparing vertical axis turbine arrays and providing configurations thereof
US20130101153A1 (en) * 2011-10-20 2013-04-25 Mitsubishi Heavy Industries, Ltd. Insertion-hole blockage-rate evaluation system, insertion-hole blockage-rate evaluation method, and insertion-hole blockage-rate evaluation program
US20140039834A1 (en) * 2012-08-01 2014-02-06 Hitachi Power Solutions Co., Ltd. Method and apparatus for monitoring equipment conditions
US20150363925A1 (en) * 2013-01-22 2015-12-17 Hitachi, Ltd. Anomaly Diagnosis Method and Apparatus

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7866211B2 (en) * 2004-07-16 2011-01-11 Rosemount Inc. Fouling and corrosion detector for process control industries
US7405558B2 (en) * 2006-07-10 2008-07-29 Areva Np, Inc. Steam generator mapping with reflections of eddy current signal
CN102564747B (en) * 2011-12-13 2014-12-10 中国农业大学 Comprehensive evaluation method and test system for clogging characteristics of emitters of drip irrigation system

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030059837A1 (en) * 2000-01-07 2003-03-27 Levinson Douglas A. Method and system for planning, performing, and assessing high-throughput screening of multicomponent chemical compositions and solid forms of compounds
US20050089923A9 (en) * 2000-01-07 2005-04-28 Levinson Douglas A. Method and system for planning, performing, and assessing high-throughput screening of multicomponent chemical compositions and solid forms of compounds
US20040039968A1 (en) * 2000-09-11 2004-02-26 Kimmo Hatonen System, device and method for automatic anomally detection
US7519860B2 (en) * 2000-09-11 2009-04-14 Nokia Corporation System, device and method for automatic anomaly detection
US20050011278A1 (en) * 2003-07-18 2005-01-20 Brown Gregory C. Process diagnostics
US7290450B2 (en) * 2003-07-18 2007-11-06 Rosemount Inc. Process diagnostics
US7428514B2 (en) * 2004-01-12 2008-09-23 Honda Research Institute Europe Gmbh System and method for estimation of a distribution algorithm
US20060111844A1 (en) * 2004-11-12 2006-05-25 Chandler Larry S Implementing two dimensional segment inversions with inversion-conforming data sets processing being rendered to include generalized composite weight factors in the processing of error-affected multivariate data samples
US20070276270A1 (en) * 2006-05-24 2007-11-29 Bao Tran Mesh network stroke monitoring appliance
US20120179402A1 (en) * 2011-01-06 2012-07-12 Mitsubishi Heavy Industries, Ltd. Deposit measurement apparatus, deposit measurement method, and computer-readable storage medium storing deposit measurement program
US20130041636A1 (en) * 2011-08-12 2013-02-14 Anna CRAIG Methods and systems for comparing vertical axis turbine arrays and providing configurations thereof
US20130101153A1 (en) * 2011-10-20 2013-04-25 Mitsubishi Heavy Industries, Ltd. Insertion-hole blockage-rate evaluation system, insertion-hole blockage-rate evaluation method, and insertion-hole blockage-rate evaluation program
US20140039834A1 (en) * 2012-08-01 2014-02-06 Hitachi Power Solutions Co., Ltd. Method and apparatus for monitoring equipment conditions
US20150363925A1 (en) * 2013-01-22 2015-12-17 Hitachi, Ltd. Anomaly Diagnosis Method and Apparatus

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
John C, Sloan "Finite safety models for high-assurance systems", A Dissertation Submitted to the Faculty of The College of Engineering and Computer Science in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy, Florida Atlantic University Boca Raton, Florida August 2010 *
K.G. Murti and et al, "Classification of Spectral Patterns Obtained from Eustachian Tube Sonometry", IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. BME-29, NO. 6, JUNE 1982 *
M. Halkidi and et al, "On Clustering Validation Techniques", Journal of Intelligent Information Systems, 17:2/3, 107–145, 2001, © 2001 Kluwer Academic Publishers. Manufactured in The Netherlands. *

Also Published As

Publication number Publication date Type
EP3087444A2 (en) 2016-11-02 application
EP3087444B1 (en) 2017-10-11 grant
CA2934949A1 (en) 2015-07-02 application
ES2655520T3 (en) 2018-02-20 grant
KR20160102049A (en) 2016-08-26 application
WO2015097221A3 (en) 2015-08-13 application
WO2015097221A2 (en) 2015-07-02 application
JP2017502293A (en) 2017-01-19 application
CN106030428A (en) 2016-10-12 application
FR3015757A1 (en) 2015-06-26 application
EP3087444B9 (en) 2018-06-27 grant

Similar Documents

Publication Publication Date Title
Fassois et al. Time-series methods for fault detection and identification in vibrating structures
Hirata et al. Galaxy–galaxy weak lensing in the Sloan Digital Sky Survey: intrinsic alignments and shear calibration errors
Marinoni et al. The VIMOS VLT Deep Survey-Evolution of the non-linear galaxy bias up to z= 1.5
Young et al. Automated spectral analysis II: application of wavelet shrinkage for characterization of non‐parameterized signals
Dennis et al. On the limitations of Taylor's hypothesis in constructing long structures in a turbulent boundary layer
Nishikawa et al. Concrete crack detection by multiple sequential image filtering
Chen et al. Feature extraction and selection for defect classification of pulsed eddy current NDT
Steyn et al. The detection of mixed layer depth and entrainment zone thickness from lidar backscatter profiles
Men’Shchikov et al. A multi-scale, multi-wavelength source extraction method: getsources
Timmins et al. A method for automatic estimation of instantaneous local uncertainty in particle image velocimetry measurements
Uttaro et al. Factors influencing the Mantel-Haenszel procedure the detection of differential item functioning
Takehara et al. A Kalman tracker for super-resolution PIV
Páez et al. A general framework for estimation and inference of geographically weighted regression models: 1. Location-specific kernel bandwidths and a test for locational heterogeneity
Ishikawa et al. A novel algorithm for particle tracking velocimetry using the velocity gradient tensor
Choi et al. Morphological analysis and classification of types of surface corrosion damage by digital image processing
Connolly et al. Calibration of the cloud particle imager probes using calibration beads and ice crystal analogs: The depth of field
Fadda et al. Density estimation with non–parametric methods
Demyanov et al. Neural network residual kriging application for climatic data
Liu et al. Fast leak detection and location of gas pipelines based on an adaptive particle filter
Nenes et al. A theoretical analysis of cloud condensation nucleus (CCN) instruments
Murray et al. An application of Gappy POD
Poelma et al. Turbulence statistics from optical whole-field measurements in particle-laden turbulence
US5424640A (en) Method for removal of random noise in eddy-current testing system
US3312274A (en) Calorimeter for measuring fouling resistance of a surface condenser tube
Dervilis et al. On robust regression analysis as a means of exploring environmental and operational conditions for SHM data

Legal Events

Date Code Title Description
AS Assignment

Owner name: ELECTRICITE DE FRANCE, FRANCE

Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:RENARD, BRUNO;PAUL, NICOLAS;SIGNING DATES FROM 20160818 TO 20160822;REEL/FRAME:040024/0070