[go: up one dir, main page]

WO2008002192A1 - Procédé de création de modèles de graduation multidimensionnels - Google Patents

Procédé de création de modèles de graduation multidimensionnels Download PDF

Info

Publication number
WO2008002192A1
WO2008002192A1 PCT/RU2007/000317 RU2007000317W WO2008002192A1 WO 2008002192 A1 WO2008002192 A1 WO 2008002192A1 RU 2007000317 W RU2007000317 W RU 2007000317W WO 2008002192 A1 WO2008002192 A1 WO 2008002192A1
Authority
WO
WIPO (PCT)
Prior art keywords
samples
calibration
properties
results
changes
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/RU2007/000317
Other languages
English (en)
French (fr)
Inventor
Konstantin Anatolyevich Zharinov
Pavel Aleksandrovich Luzanov
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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to EA200900029A priority Critical patent/EA012950B1/ru
Priority to CN200780023164.4A priority patent/CN101473197B/zh
Priority to CA2655745A priority patent/CA2655745C/en
Priority to EP20070794046 priority patent/EP2040042B1/en
Priority to US12/305,122 priority patent/US20090204352A1/en
Publication of WO2008002192A1 publication Critical patent/WO2008002192A1/ru
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01JMEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
    • G01J3/00Spectrometry; Spectrophotometry; Monochromators; Measuring colours
    • G01J3/28Investigating the spectrum
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/27Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
    • G01N21/274Calibration, base line adjustment, drift correction

Definitions

  • the claimed invention relates to analytical instrumentation, in particular, to methods for creating calibration models for various types of measuring instruments.
  • determining the properties of a sample using direct measurement methods does not provide the necessary analysis speed or leads to destruction of the sample, and in some cases it may not be possible to directly determine the required properties.
  • An example is analytical chemistry, where the traditional direct method of chemical analysis for determining the concentrations of the components that make up a sample is based on chemical reactions, which leads to the destruction of the sample and, in addition, requires considerable time.
  • indirect methods are widely used, where the analyzed properties of the samples are determined by measuring other properties of the samples, depending on the analyzed, which, however, unlike the analyzed properties, can be easily measured directly for a short period of time and without destroying the sample. Often, the directly measured properties of the samples are called “primary properties”, and the analyzed properties of the samples are called “primary properties *.
  • One of the most effective indirect research methods is spectroscopic analysis, in which, on the basis of measuring the optical absorption, reflection, or scattering spectra ("primary" properties), the required "secondary" properties of the samples (for example, chemical composition) are determined.
  • calibration set In order for the constructed calibration model to provide the specified accuracy in determining the parameters characterizing the analyzed secondary properties of an arbitrary sample, it is necessary to analyze a large number of samples representative of those samples that will be analyzed on the device in the future (calibration set).
  • the selection of calibration calibration samples is regulated by standards for various indirect analysis methods, for example, the standard for qualitative analysis using NIR spectroscopy [1].
  • the values of the parameters characterizing the analyzed secondary properties of the samples of the calibration set are preliminarily determined using standard methods called reference. When conducting reference analyzes, special attention should be paid to the accuracy of the analysis results, since the accuracy of the reference analysis limits the accuracy of the entire calibration.
  • the analyzed properties of the samples of the calibration set can be previously determined by standard chemical methods using chemical reactions. Samples of the calibration set are selected so that their secondary properties cover the possible range of changes in these properties when analyzing unknown samples and are evenly distributed over this range. Fulfillment of these conditions increases stability calibration model, due to which small noise changes in the measured primary properties do not lead to statistically significant changes in the results of the analysis of secondary properties.
  • the created calibration model should be subjected to a standard verification procedure called validation [1], as a result of which statistical parameters characterizing the quality of the calibration model, as well as samples falling out of the calibration, are determined.
  • validation [1] a standard verification procedure
  • a statistical analysis of emissions is used, for example, according to the distance of Mahalanobis [1], using the measurement data of the primary properties of the complete calibration set.
  • Mahalanobis [1] using the measurement data of the primary properties of the complete calibration set.
  • such samples should be excluded from the calibration set.
  • the measured primary properties of the sample are also analyzed by emission statistics. This problem is similar to the tasks solved in a qualitative analysis, where, based on measuring the primary properties of the sample (spectral data) and comparing them with library data, a conclusion is drawn about the set of components in the sample. It should be noted that the measurement conditions of the analyzed sample and calibration samples should be identical.
  • the results of measuring the primary properties of the samples can be significantly affected by various properties, for example, changes in external conditions or technical characteristics of the measuring device.
  • the technical characteristics of the device may change over time due to aging or during repair work, as well as the replacement of individual typical structural elements of the device.
  • this process is lengthy, time consuming and expensive. Therefore, several methods were proposed for creating calibration models that are resistant to changing conditions. measurements, technical characteristics of the device and other properties that affect the measurement results, which allows you not to repeat the complex process of constructing a calibration model on each individual device after it has been repaired, replaced with some typical structural elements or changed device parameters over time.
  • This method takes into account the variability of spectral instruments by including in the calibration model the spectra of samples recorded when changing certain parameters of the instrument in order to cover the entire possible range of such changes.
  • This idea is almost identical to the method described in the patent [2]. For example, a wavelength shift is introduced by a certain amount for all points at which measurements are made, or the magnitude of this shift is proportional to the position of a particular point on the wave number axis. It is also allowed to account for such changes using mathematical processing of spectra.
  • this method requires a large number of measurements, which leads to an increase in labor costs, time and material resources at the stage of creating the calibration model.
  • the considered algorithm is not without other shortcomings, namely, a decrease in the accuracy of the model, because the model initially assumes a variation in the results of measuring primary properties, and the difficulty of assessing the appropriateness of using a calibration model to analyze an unknown sample, which can lead to errors in the analysis of secondary properties.
  • spectrometers are characterized by determining whether their spectra (spectral characteristics) correspond to a predetermined limited number of clusters. Cluster membership is determined based on the similarity of spectral features and performance characteristics. The spectral features used for classification can be attributed to known parameters of the device or represent some abstract characteristic features obtained using computational methods. The calibration model for each cluster compensates for instrument changes; it will be simpler and more accurate.
  • the proposed algorithm includes four stages: measurement, classification, calibration, determination of anomalous samples.
  • the spectra of certain standards are measured, which can be used, for example, to classify displacements observed along the axis of wavelengths or intensities.
  • the type of measured standards is determined by what effects need to be considered in each case.
  • instrument changes can be classified on the basis of spectral data into the following types: changes in signal intensity, bandwidth, wavelength, or a combination of these phenomena.
  • the source in question gives a fairly detailed classification of them. It should be noted that in some In cases, real samples can be used as standards, modeling the measurement features using this device.
  • the obtained spectra must be classified.
  • characteristic spectral features are distinguished using mathematical transformations of the obtained data that improve certain aspects useful for interpretation, for example, analysis based on the main components, calculation of the Mahalanobis distance. Also, these features can be obtained on the basis of a priori knowledge of the system (noise characteristics, detector linearity, etc.). Subsequently, classification is performed, i.e., the spectra are assigned to specific clusters, and a classification model is developed by finding some law that allows one to determine which cluster corresponds to the spectrum of the measured sample. At the third stage, a calibration model is constructed for each cluster, taking into account the instrumental changes characteristic of this cluster.
  • the spectra within a particular cluster have a high degree of internal constancy and have similar features that determine not all possible changes in the parameters of the device, but only one or more taken into account within the cluster. This approach allows us to simplify the procedure for creating graduations for a single cluster. Another advantage is the ability to determine anomalous spectra. Since each cluster has its own set of classification requirements, then if the measured spectrum does not satisfy them, an attempt is made to choose another cluster whose parameters allow the analysis of the specified spectrum. If this is not possible, this spectrum is subjected to mathematical processing to match one of the clusters. If compliance cannot be achieved, the spectrum is considered abnormal.
  • the advantages of this method include the simplicity of the calibration model and the rather high accuracy of determining the secondary properties for a particular cluster, as well as the possibility of determining anomalous spectra, which makes it possible to assess the legitimacy of the use of the calibration model for the analysis of one or another unknown sample, and, accordingly, reduce the probability of errors in the analysis secondary properties.
  • the shortcomings one can single out the need to perform a large number of operations related to the definition of clusters, finding qualification algorithms and constructing calibration models for each cluster, as well as the narrow focus of the proposed method on spectrometric instruments.
  • the relationships are determined for converting the results of measurements carried out on the original instrument to a form equivalent to the measurements on the instrument in its current state, in which there are some changes in the characteristics in the process operation. Subsequently, using the relations found, the data for the calibration set samples measured on the initial device are transformed to a form equivalent to the measurement results of the calibration set samples on the device in the current state, and a calibration model is created from these converted data.
  • the calibration model created using this method takes into account changes in the device characteristics that occurred during operation, provides a fairly high accuracy in determining the analyzed properties, and there is no need to repeat measurements on the instrument of all samples of the calibration set, which reduces the duration and complexity of the calibration process.
  • a significant disadvantage of this invention is that the created calibration model takes into account only the current state of the device and its characteristics, as well as other properties that affect the measurement results. Thus, after a certain time, these properties can change again due to various reasons, in particular aging, repair of the device, replacement of individual structural elements, changes in operating conditions. As a result, the analysis error by the created calibration model may increase, and as a result, there will be a need to build a new calibration, that is, repeat all the operations described above in accordance with this method even if there is a multiple change in the same characteristics.
  • a measurement is carried out on the instrument, for which a calibration model is determined, of a large number of parameters characterizing the primary properties of the sample, after which mathematical processing of the results is carried out and a number of calibration constants are determined in the relations between the values of the quantities characterizing the primary and secondary properties calibration samples.
  • the peculiarity of the method lies in the fact that during the measurement process for one or more samples from the calibration set, at least one of the parameters of the measuring device is intentionally changed, in addition, a change in external conditions can be introduced.
  • the magnitude of the changes in these parameters during the construction of the calibration model should be of the same order or greater than the expected magnitude of the changes in these parameters between different devices during operation.
  • Changes in instrument parameters or other measurement conditions can also be introduced not in the process of real measurements, but using mathematical transformations.
  • spectrometers based on a monochromator it is proposed to measure one or more samples used for calibration, intentionally producing a wavelength shift in the range of values that can be expected during operation. This can be achieved by making physical changes to the monochromator (depending on the design, a shift or a change in the angle of inclination of the grating or filter, a change in the position of the holes on the grating, a change in the slope of the incident radiation is used) or by changing the constants used to calculate the wavelengths of the monochromator.
  • results of the analysis of the secondary properties of the sample depend little on the measurement conditions and the technical parameters of the measuring device.
  • the deliberate introduction of the scatter of data from the measurement results of samples of the calibration set increases the stability of the model, and the range of its applicability becomes wider.
  • the objective of the present invention is to create a multidimensional calibration model that provides high accuracy in determining the analyzed properties and is resistant to changes in properties that affect the measurement results of the device, and even in the case of repeated changes of such factors, in particular, to possible linear and nonlinear differences in the technical parameters of the device arising from aging, repair or replacement of certain typical structural elements, changes in operating conditions, as well as other factors, this reduces the complexity and duration of the procedure for constructing a similar model, which is subsequently used to determine one or more secondary properties of an unknown sample from the results of measurements of many primary properties of this sample, and not necessarily spectral ones.
  • the solution of this problem is achieved by the method of creating multidimensional calibration models, consisting of a sequence of actions united by a single inventive concept.
  • the method includes selection of a calibration set of samples with known secondary properties determined by reference methods; measuring on the instrument the primary properties of each of the samples of the calibration set; artificially making changes to at least one of the properties that affect the measurement results of the device; measurement on the device in a state of primary properties altered in this way of at least one sample; creating a calibration model that is resistant to these changes by finding, using the methods of multivariate regression analysis using the data obtained on the device after making changes, calibration ratios that allow you to determine the analyzed secondary properties of an unknown sample from the results of measurements of many primary properties of this sample, characterized in that before making changes, a set of samples is formed to calculate the correcting ratios, the primary properties of each sample are measured and from this set on the device before and after making changes, and comparing using the methods of multivariate regression analysis the results of measurements of the primary properties of the samples of the set obtained on the device before making changes
  • All measurement results of calibration samples including those that complement the calibration set, the results of measuring the primary properties of samples of the calibration set, converted to a form equivalent to the measurement results on the device, when at least one of the properties affecting the results of its measurements is changed, stored in the computer of the calibrated instrument. This ensures the independence of the created model from measurements on the device in the initial state, and allows you to evaluate the validity of applying calibration to analyze an unknown sample on the device in a state where one or more properties affecting the measurement results have changed, and also determine the outliers using the statistics of emission prediction, which guarantees high accuracy in determining the analyzed secondary properties.
  • the use of emission prediction statistics and the exclusion from the model of dropping out samples of the calibration set before determining calibration ratios increases the stability of the created calibration model.
  • the normalization procedure is the choice of one or another method of mathematical preprocessing.
  • the selection criterion is the accuracy of the analysis of the secondary properties of the samples, which provides a graduated instrument, which was created using this type of mathematical preprocessing.
  • the quantitative parameters of the calibration model validation procedure are used (for example, standard validation error) [1].
  • the essence of the invention lies in the fact that the proposed combination of features allows you to create a multidimensional calibration model on a calibrated measuring device that is resistant to possible variations in properties that affect the measurement results of the device, in particular changes in the technical parameters of the device resulting from aging, repair or replacement of individual typical elements design, as well as the characteristics of the operating conditions, and even in the case of repeated changes of these factors, while zhaetsya the complexity and duration of procedures for constructing such a model compared with existing known methods.
  • the created calibration model makes it possible to accurately predict the secondary properties of unknown samples from the results of measurements of many primary, not necessarily spectroscopic, properties, and the calibration model is built on the basis of data on the primary properties of the samples from the calibration set measured on the calibrated instrument and the same data, transformed to look, as if measurements were taken on the device, when some changes were made to at least one of the properties that affect x on the measurement results of the device, for example, corresponding to the maximum expected during operation changes for devices of this type.
  • Correlation ratios that allow you to convert the measurement results of calibration samples to a form equivalent to the measurement results on the device in an altered state are determined from measurements of a representative set of samples on the device before and after changes are made, and the set for calculating correction ratios consists of a much smaller number of samples than the calibration . Samples from the set for calculating corrective relations provide significant differences in the measurement results over the entire range of primary properties on the device, both in the initial and in the changed state. Using a set of samples to calculate corrective ratios allows you to determine the nonlinear relationship between the measurement results of the same samples on the device before and after making changes by correlation analysis using regression methods.
  • the inventive method for creating multidimensional calibration models can be used for any devices where the analyzed properties of the sample are determined on the basis of multiple measurements of other properties, in particular, spectrophotometry for various types of spectrometers measuring the absorption of light radiation by a sample at many different wavelengths. Data characterizing the result of such measurements is called the spectrum.
  • spectrophotometry for various types of spectrometers measuring the absorption of light radiation by a sample at many different wavelengths.
  • Data characterizing the result of such measurements is called the spectrum.
  • a spectrometer like any analytical instrument, requires preliminary calibration, which determines the relationship between the analyzed properties of the sample and spectral characteristics. Note that often with the analyzed properties of the sample, the results of measurements are not directly compared with spectral data that have already passed the normalization procedure (preliminary mathematical preprocessing). Thus, for example, spectral smoothing, baseline subtraction, or differentiation can be performed.
  • the type of preliminary mathematical processing is selected based on the criterion of maximum accuracy in determining the analyzed secondary properties and minimizing the influence of side factors associated with parasitic scattering and sample preparation features.
  • the same mathematical pretreatment is applied to all spectra of samples of the calibration set; after the normalization procedure, the transformed spectral data have pronounced characteristic features even with minor changes in the analyzed properties [6].
  • Various statistical characteristics are used as criteria for assessing the quality of the calibration model, for example, standard calibration error (SEC), standard validation error (SEV), and standard cross-validation error (SECV) [1].
  • SEC standard calibration error
  • SEV standard validation error
  • SECV standard cross-validation error
  • the most common type of mathematical preprocessing in spectral analysis is to find the weighted average values of spectral data [6], which reduces by 1 the number of degrees of freedom in the calibration model. In this preprocessing find the spectrum averaged over the entire calibration set and subtract it from each spectrum of calibration samples. Similarly, find the average values of the reference data. Then, when analyzing an unknown sample, before applying the constructed calibration model, the spectrum averaged over the calibration set is also subtracted from the measured spectrum.
  • a set of samples is selected that is representative of those samples that will be further analyzed.
  • the following criteria are used [1]: a) the samples should contain all the chemical components that are planned to be analyzed; b) the range of changes in the concentration of the analyzed components in the samples of the calibration set should exceed the range of changes in the analyzed unknown samples; c) the magnitude of the change in the concentration of chemical components from sample to sample should be evenly distributed over the entire range of changes; d) the number of samples should ensure that, using statistical methods, mathematical relationships between spectroscopic data and the concentration of individual chemical components are found. Samples falling out of the calibration set are determined using a statistical analysis of emissions, for example, by calculating the Mahalanobis distance [1], which is defined as:
  • analyzed properties of r-calibration samples according to the range in which they vary, i.e. when the composition of the sample is significantly different from the rest of the samples in the calibration set. Samples should also be excluded from the calibration set, the values of the analyzed properties of which, determined using the constructed model, are significantly different from the values that the reference method gives. These samples are determined from student discrepancies calculated by the formula: Here ⁇ j - unlike obtained using the calibration model of the chemical component concentration or the analyzed property from the referential value for i-th calibration sample, SEC - standard error of calibration [1], Dj 2 - Mahalanobis distance for the i-th calibration sample. Student discrepancies should be evenly distributed according to normal law. The discrepancy is compared with the student coefficient, for confidence probability
  • the analyzed properties of the calibration samples should be known in advance or determined using reference methods.
  • the obtained numerical values are considered real values of the analyzed properties, therefore, the accuracy of the reference analysis determines the accuracy of the entire calibration, and quite high requirements are imposed on this procedure. It is possible to increase the accuracy of the reference analysis by conducting repeated repetitions of the analyzes with subsequent averaging of the results.
  • the spectra of the calibration set samples are recorded on the instrument, after which a multidimensional calibration model can be created.
  • multivariate mathematical methods are used, such as principal component analysis (PCA), fractional least squares (PLS) method, etc.
  • PCA principal component analysis
  • PLS fractional least squares
  • the analyzed properties for example, chemical composition. Accuracy such measurements will be sufficiently high as long as the properties that influence the measurement results of the device remain constant. But, as already mentioned, during operation, during repair work or the replacement of individual typical structural elements, for example, the technical characteristics of the device can change, which in turn can lead to a decrease in the accuracy of predictions and the need to build a new calibration model.
  • the inventive method allows you to create a multidimensional calibration model for a calibrated device that is resistant to possible changes in one or more factors affecting the results of its measurements, without measuring the full set of calibration samples on the device, when the specified factor or factors have changed, and this is true even in the case of multiple changes.
  • spectral data measured on a calibrated instrument are used, supplemented by spectral data obtained by converting the initial data to a form equivalent to measurements on the instrument in a state when at least one of the properties affecting the measurement results is changed, for example, in such a way to cover the expected range of similar changes for devices of this type during operation.
  • the scope and stability of the calibration model are analyzed based on the totality of these spectral data. Since all the data is stored in the computer of the calibrated instrument, this allows the evaluation of the outliers data when analyzing unknown samples, for example, using Mahalanobis statistics.
  • a specially selected set of samples is used, hereinafter referred to as a set for calculating correction ratios, the number of samples in which is much smaller than in the full calibration set, at
  • their properties may not be known; the only important thing is that this set of samples provides significant variations in the measured spectral data, allowing one to find the expressions and I.
  • the spectrum of each sample of the set is measured on the device in a state corresponding to the measurement of the spectra of calibration samples, and in the same device, when at least one of the properties that affect the measurement results has been made some changes (“changed” device).
  • the spectral data can be subjected to normalization, which consists in carrying out the same mathematical transformations for all measured spectra. This ensures the identification of obvious differences in the spectral data measured on the instrument before and after the introduction of changes, which in turn contributes to a more accurate definition of expressions for transforming spectral data.
  • the calibration model for the calibrated device is created from the totality of the initial and converted spectral data of the calibration samples using standard mathematical methods of multivariate regression analysis (MLR, PCA, PLS, etc. [6]), and falling samples should be excluded from the calibration, which guarantees stability of the created model, after which it can be used to determine the properties of an unknown sample.
  • MLR multivariate regression analysis
  • the obtained spectral data can be subjected to a normalization procedure that takes into account the characteristics of transmitting instruments and using the principles of Fourier spectroscopy [6].
  • a normalization procedure that takes into account the characteristics of transmitting instruments and using the principles of Fourier spectroscopy [6].
  • the type of mathematical processing is determined in such a way as to ensure a minimum error in determining the analyzed properties. If the spectral data are subjected to the normalization procedure, then the calibration ratios are determined based on a comparison of these data with the known ones that also underwent normalization, properties of samples of the calibration set determined by reference methods.
  • SEC standard calibration error
  • SECV cross-validation standard error
  • the standard validation error characterizes the deviation error between the reference and the values predicted by the calibration equation for the samples of the additional set that were not included in the calibration.
  • d v is the total number of reference values of the analyzed indicator for all spectra of the additional set, y. - reference values of the analyzed indicator for the i-th spectrum of the additional set, y.- predicted values of the analyzed indicator for the i-th spectrum of the additional set.
  • Rf 1 are the values of the spectral data measured on the device when changes are made to its design (i-th wavelength, j-th sample from the set for calculating corrective relations), i ?, '"are similar spectral data measured on the device Before making changes, the spectral data can undergo a normalization procedure, and the mathematical processing used must be the same for the spectra measured on the instrument both before and after making the changes. Drains are regression coefficients.
  • the multidimensional calibration model created in accordance with the claimed invention makes it possible to determine the properties of unknown samples with high accuracy and is much less sensitive to the possible change in the design of the device taken into account using the proposed algorithm, in this case, the replacement of the beam splitter.
  • any other changes in properties that affect the measurement results of the device can be taken into account, the example given here is chosen as a clear illustration of the claimed invention due to the fact that such changes are one of the most characteristic for spectrometers of this type, and can lead to a significant decrease in the accuracy of analysis according to the originally constructed calibration model of the device, if they are not taken into account when creating it.
  • Table 2 Verification of the initial calibration on a device with different beam splitter m and.
  • Table 3 Calibration verification, adjusted taking into account the replacement of the beam splitter, on a device with different beam splitters.
  • the spectral data obtained during the measurement on a graduated instrument of a calibration set of samples were adjusted according to the corrections found earlier to the form corresponding to the instrument for which the beam splitter was replaced.
  • a calibration model was created from the totality of the initial and adjusted spectral data of the calibration samples of the device.
  • several different beam splitters were installed on a calibrated instrument, and in each case, the spectra of the validation set of samples were recorded.
  • the basic statistical parameters of the initial calibration model constructed for this device and a calibration model resistant to the effect of the beam splitter replacement were estimated on the basis of correction ratios calculated earlier for another device.
  • the data obtained are shown in tables 4 and 5, respectively.
  • the results of tables 4 and 5 confirm the possibility of using corrective relations defined for any one device and taking into account the influence of changes in its parameters or other properties that affect the measurement results on any other device in this series when constructing a calibration model that is resistant to these changes.
  • This feature allows, after creating, in accordance with the proposed method, a calibration model that compensates for differences in a particular property or properties, for example, the technical parameters of the device, changes that occur due to aging, repair or replacement of individual typical structural elements, for any one instrument, in the future use the found correcting relations for all other instruments of this type in the construction of calibration models her, taking into account with the help of these correlation ratios the corresponding changes in the affecting properties.
  • the scope of the proposed method is not limited to spectroscopy.
  • the claimed invention can be used for various devices that determine some properties of the sample based on multiple measurements of other properties.

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Mathematical Physics (AREA)
  • Biochemistry (AREA)
  • Theoretical Computer Science (AREA)
  • Engineering & Computer Science (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Indication And Recording Devices For Special Purposes And Tariff Metering Devices (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Description

СПОСОБ СОЗДАНИЯ МНОГОМЕРНЫХ ГРАДУИРОВОЧНЫХ МОДЕЛЕЙ
Заявляемое изобретение относится к аналитическому приборостроению, в частности, к способам создания градуировочных моделей для различного вида измерительных приборов.
Во многих отраслях промышленности и научных исследований определение свойств образца при помощи прямых методов измерения не обеспечивает необходимой скорости анализа или приводит к разрушению образца, а в некоторых случаях непосредственное определение требуемых свойств вообще может оказаться невозможным. В качестве примера можно привести аналитическую химию, где традиционный прямой метод химического анализа для определения концентраций образующих образец компонентов основан на проведении химических реакций, что приводит к разрушению образца, а, кроме того, требует значительных временных затрат.
В подобных случаях широкое распространение получили непрямые методы, где анализируемые свойства образцов определяют путем измерения других свойств образцов, зависящих от анализируемых, которые, однако, в отличие от анализируемых свойств, могут быть легко измерены напрямую в течение короткого промежутка времени и без разрушения образца. Часто при этом непосредственно измеряемые свойства образцов называют «пepвичными cвoйcтвaми», а анализируемые свойства образцов «втopичными свойствами*. Одним из наиболее эффективных непрямых методов исследования является спектроскопический анализ, при котором на основе измерения оптических спектров поглощения, отражения или рассеяния («пepвичныx» свойств) определяют требуемые «втopичныe» свойства образцов (например, химический состав).
После измерения множества первичных свойств необходимо найти математические соотношения, связывающие результаты измерений со значениями величин, характеризующих анализируемые вторичные свойства.
Эти соотношения между измеряемыми на приборе первичными характеристиками и анализируемыми вторичными свойствами образца называются градуировочными моделями или градуировками, в результате их нахождения определяется ряд градуировочных констант. Процедура создания градуировочных моделей достаточно трудоемкая, длительная и дорогостоящая. Особенно этот процесс усложняется для случая многомерного анализа, когда для определения количественных характеристик вторичных свойств образца используется набор результатов измерения большого числа параметров, характеризующих первичные свойства. Например, в случае спектроскопического анализа для определения концентрации различных компонентов проводят измерения большого количества спектральных данных (величины поглощения, отражения или рассеяния) для разных значений волновых чисел (длин волн, частот).
Чтобы построенная градуировочная модель обеспечивала заданную точность определения параметров, характеризующих анализируемые вторичные свойства произвольного образца, необходимо провести анализ большого числа образцов, представительных к тем образцам, которые будут анализироваться на приборе в дальнейшем (градуировочный набор). Подбор образцов градуировочного набора регламентируется стандартами на различные непрямые методы анализа, например, стандартом на качественный анализ при помощи БИК спектроскопии [1]. Значения параметров, характеризующих анализируемые вторичные свойства образцов градуировочного набора, предварительно определяются при помощи стандартных методов, называемых референтными. При проведении референтных анализов следует обратить особое внимание на точность результатов анализов, так как точность референтного анализа лимитирует точность всей градуировки. Например, в рассмотренном ранее случае спектроскопического анализа химического состава образцов, анализируемые свойства образцов градуировочного набора могут быть предварительно определены стандартными химическими методами с использованием химических реакций. Образцы градуировочного набора выбираются таким образом, чтобы их вторичные свойства охватывали возможный диапазон изменений этих свойств при анализе неизвестных образцов и были равномерно распределены по этому диапазону. Выполнение данных условий повышает устойчивость градуировочной модели, благодаря чему небольшие шумовые изменения в измеряемых первичных свойствах не приводят к статистически значимым изменениям в результатах анализа вторичных свойств.
Созданная градуировочная модель должна быть подвергнута стандартной процедуре проверки, называемой валидацией [1], в результате которой определяются статистические параметры, характеризующие качество градуировочной модели, а также выпадающие из градуировки образцы. Для поиска таких образцов применяется статистический анализ выбросов, например, по расстоянию Махаланобиса [1], использующего данные измерений первичных свойств полного градуировочного набора. Для повышения устойчивости градуировки такие образцы необходимо исключить из градуировочного набора.
Чтобы оценить применимость градуировочной модели для описания неизвестного образца, измеренные первичные свойства образца также подвергают анализу по статистике определения выбросов. Эта проблема схожа с задачами, решаемыми при качественном анализе, где на основании измерения первичных свойств образца (спектральных данных) и сопоставления их с библиотечными данными делается вывод о наборе компонентов в образце. Следует заметить, что условия измерений анализируемого образца и градуировочных образцов должны быть идентичными.
На результаты измерения первичных свойств образцов, а как следствие, и на точность определения параметров, характеризующих анализируемые вторичные свойства образцов, могут существенно влиять различные свойства, например, изменения внешних условий или технических характеристик измерительного прибора. Необходимо отметить, что технические характеристики прибора могут изменяться со временем вследствие старения или при проведении ремонтных работ, а также замены отдельных типовых элементов конструкции прибора. В таком случае может потребоваться создание новой градуировочной модели, учитывающей особенности текущего состояния прибора. Как отмечалось ранее этот процесс длительный, трудоемкий и дорогостоящий. Поэтому было предложено несколько способов создания градуировочных моделей, устойчивых к изменению условий измерений, технических характеристик прибора и прочих влияющих на результаты измерений свойств, что позволяет не повторять сложный процесс построения градуировочной модели на каждом отдельном приборе после выполнения его ремонта, замены отдельных типовых элементов конструкции или изменения параметров прибора с течением времени.
В одном из подобных способов [3] рассматривается построение градуировочной модели приборов для спектроскопического анализа химического состава, которая мало восприимчива к изменению характеристик прибора. Данный подход применили T. Sсесiпа и др., которые основывались на том, что искусственное добавление в градуировочную модель спектров образцов, измеренных при изменении конкретного влияющего параметра (например, температура или давление) свыше ожидаемого диапазона изменения при эксплуатации улучшает анализ с использованием подобной модели тех образцов, для которых изменение данного параметра характерно. Применяя указанную концепцию по отношению к изменению технических характеристик измеряющего прибора, был разработан алгоритм, позволяющий учесть подобные изменения. Указанный метод принимает во внимание изменчивость спектральных приборов за счет включения в градуировочную модель спектров образцов, зарегистрированных при изменении определенных параметров прибора, чтобы перекрыть весь возможный диапазон подобных изменений. Данная идея практически идентична методу, описанному в патенте [2]. Например, вводится сдвиг длины волны на определенную величину для всех точек, в которых производятся измерения, или величина данного сдвига пропорциональна положению конкретной точки на оси волновых чисел. Допускается также учет подобных изменений с помощью математической обработки спектров.
С другой стороны предлагается возможность модификации рассмотренного алгоритма, заключающаяся в создании градуировочных моделей на нескольких приборах, предпочтительно при этом, чтобы технические характеристики этих приборов охватывали возможный диапазон изменения указанных характеристик на всех приборах данной серии. Затем полученные данные объединяются в единый градуировочный набор и дополнительно искусственно вводятся изменения характеристик прибора посредством математических манипуляций над спектральными данными. Такой подход позволяет в еще большей степени учесть изменчивость параметров прибора, делая создаваемую модель более устойчивой.
Однако, данный метод требует большого количества проводимых измерений, что приводит к увеличению затрат труда, времени и материальных ресурсов на стадии создания градуировочной модели. Рассмотренный алгоритм не лишен и других недостатков, а именно снижение точности модели, т.к. в модели изначально предполагается вариация результатов измерения первичных свойств, и сложность оценки правомерности применения градуировочной модели для анализа того или иного неизвестного образца, что может привести к ошибкам в анализе вторичных свойств.
Еще один способ, посвященный созданию градуировочных моделей, компенсирующих изменения параметров прибора, предложен в [4]. Согласно данному методу спектрометры характеризуются при помощи определения соответствия их спектров (спектральных характеристик) предварительно определенному ограниченному числу кластеров. Принадлежность кластеру определяется на основании подобия спектральных особенностей и рабочих характеристик. Спектральные особенности, используемые для классификации, могут быть отнесены к известным параметрам прибора или представлять собой некоторые абстрактные характерные черты, получаемые с помощью вычислительных методов. Градуировочная модель для каждого кластера компенсирует приборные изменения, она будет более простой и точной.
Предложенный алгоритм включает в себя четыре стадии: измерение, классификация, градуировка, определение аномальных образцов. На первой стадии производится измерение спектров определенных стандартов, которые могут быть использованы, например, для классификации смещений, наблюдаемых по оси длин волн или интенсивности. Вид измеряемых стандартов определяется тем, какие эффекты необходимо учесть в каждом конкретном случае. В общем случае приборные изменения могут быть классифицированы на основе спектральных данных на следующие виды: изменения интенсивности сигнала, ширины полосы, длины волны или комбинация указанных явлений. В рассматриваемом источнике дается довольно подробная их классификация. Нужно отметить, что в некоторых случаях в качестве стандартов могут выступать реальные образцы, моделирующие особенности измерений с помощью данного прибора. На второй стадии полученные спектры необходимо классифицировать. Для этого выделяют характерные спектральные особенности, используя математические преобразования полученных данных, улучшающие определенные аспекты, полезные для интерпретации, например, анализ на основе главных компонент, вычисление расстояния Махаланобиса. Также данные особенности можно получить на основе априорного знания о системе (шумовые характеристики, линейность детектора и т.п.). В дальнейшем производится классификация, то есть отнесение полученных спектров определенным кластерам, и разрабатывается классификационная модель путем нахождения некоторого закона, позволяющего определить, какому кластеру соответствует спектр измеряемого образца. На третьей стадии для каждого кластера строится градуировочная модель, учитывающая приборные изменения, характерные для этого кластера. Спектры внутри определенного кластера имеют высокую степень внутреннего постоянства и обладают схожими особенностями, определяющими не все возможные изменения параметров прибора, а лишь одно или несколько, учитываемых внутри кластера. Подобный подход позволяет упростить процедуру создания градуировки для единичного кластера. Другое его преимущество - возможность определения аномальных спектров. Поскольку каждый кластер имеет свой собственный набор классификационных требований, то если измеренный спектр им не удовлетворяет, производится попытка подобрать другой кластер, параметры которого позволяют провести анализ указанного спектра. Если это оказывается невозможным, данный спектр подвергают математической обработке, чтобы он соответствовал одному из кластеров. Если соответствия достигнуть не удается, спектр считается аномальным.
По сути, это модификация метода, описанного в предыдущих патентах с введением изменений в характеристики приборов при построении градуировочной модели. Отличие здесь в том, что рассматривают не весь возможный диапазон изменений, а разбивают его на несколько поддиапазонов (кластеров), в каждом из которых, например, изменяется только один параметр и в ограниченном диапазоне значений. Строится градуировочная модель, на основании которой вторичные свойства определяются с достаточно высокой точностью для всех приборов, первичные свойства которых характеризуются данным кластером. Для приборов из другого кластера (изменяется другой параметр или тот же параметр, но в другом диапазоне значений) строится независимая градуировочная модель или проводят перенос градуировки по известным методам.
К достоинствам данного метода относятся простота градуировочной модели и достаточно высокая точность определения вторичных свойств для конкретного кластера, а также возможность определения аномальных спектров, что позволяет оценить правомерность применения градуировочной модели для анализа того или иного неизвестного образца, и, соответственно, снизить вероятность ошибок в анализе вторичных свойств. Из недостатков можно выделить необходимость выполнения большого количества операций, связанных с определением кластеров, нахождением квалификационных алгоритмов и построением градуировочных моделей для каждого кластера, а также узкая ориентированность предложенного метода на спектрометрические приборы.
Известен способ переноса градуировочных моделей между приборами [5], который также может быть использован для переградуировки одного прибора, чтобы учесть изменения в характеристиках прибора, возникшие в процессе эксплуатации вследствие различных причин, в частности старения. Согласно данному методу, для учета подобных изменений, выбирается некоторый набор образцов, например, из тех, которые использовались при создании градуировочной модели. Первичные свойства каждого образца из этого набора измеряются на приборе, характеристики которого изменились с течением времени, и производится сопоставление при помощи регрессионных методов результатов измерений первичных свойств образцов набора с результатами измерения первичных свойств тех же самых образцов, полученными на приборе в исходном состоянии, для которого создавалась градуировка. Таким образом, определяются соотношения для преобразования результатов измерений, проведенных на исходном приборе, к виду, эквивалентному измерениям на приборе в его текущем состоянии, при котором произошли некоторые изменения в характеристиках в процессе эксплуатации. В дальнейшем при помощи найденных соотношений данные для образцов градуировочного набора, измеренных на исходном приборе, преобразуются к виду, эквивалентному результатам измерений образцов градуировочного набора на приборе в текущем состоянии, и по этим преобразованным данным создается градуировочная модель.
Созданная при помощи данного способа градуировочная модель учитывает произошедшие в процессе эксплуатации изменения в характеристиках прибора, обеспечивает достаточно высокую точность определения анализируемых свойств, при этом нет необходимости повторять измерения на приборе всех образцов градуировочного набора, что снижает длительность и трудоемкость процесса градуировки. При измерении неизвестного образца допускается возможность проведения анализа выпадающих данных при помощи статистики предсказания выбросов, что позволяет оценить правомерность применения градуировочной модели. Для учета влияния мешающих факторов, таких как различия в пробоподготовке и в состоянии образцов, предложено использовать процедуру нормализации с использованием различных видов математической предобработки результатов измерений и данных по вторичным свойствам (референтных данным).
Существенный недостаток данного изобретения в том, что созданная градуировочная модель учитывает лишь текущее состояние прибора и его характеристик, а также иных свойств, оказывающих влияние на результаты измерений. Таким образом, по прошествии определенного времени эти свойства могут измениться снова вследствие различных причин, в частности, старения, выполнения ремонта прибора, замены отдельных конструктивных элементов, изменения условий эксплуатации. В результате погрешность анализа по созданной градуировочной модели может увеличиться, и в итоге возникнет необходимость построения новой градуировки, то есть повторения всех операций описанных выше в соответствии с данным методом даже в том случае, если происходит многократное изменение одних и тех же характеристик.
Известен еще один способ [2] создания многомерных градуировочных моделей, маловосприимчивых к изменению параметров прибора, на котором проводятся измерения, а также к изменению внешних условий, при которых происходят измерения, и к изменению прочих свойств образца. По совокупности существенных признаков данный способ является наиболее близким к заявляемому изобретению и выбран в качестве прототипа. В рассматриваемом методе для создания градуировочной модели выбирается необходимое количество градуировочных образцов с известными вторичными свойствами, таким образом, чтобы прочие свойства образцов изменялись в максимально возможном диапазоне ожидаемых изменений. Затем для каждого из выбранных градуировочных образцов проводят измерение на приборе, для которого определяется градуировочная модель, большого числа параметров, характеризующих первичные свойства образца, после чего проводится математическая обработка полученных результатов и определяется ряд градуировочных констант в соотношениях между значениями величин, характеризующих первичные и вторичные свойства градуировочных образцов. Особенность метода заключается в том, что в процессе измерений для одного или нескольких образцов из градуировочного набора намеренно производят изменения, как минимум, одного из параметров измерительного прибора, дополнительно при этом может вводиться изменение внешних условий. Величина изменений указанных параметров при построении градуировочной модели должна быть того же порядка или большей, чем ожидаемая величина изменений этих параметров между различными приборами в процессе эксплуатации. Изменения параметров прибора или прочих условий измерений могут также вводиться не в процессе проведения реальных измерений, а с помощью математических преобразований. Например, в случае спектрометров на основе монохроматора предлагается провести измерение одного или нескольких образцов, использующихся для градуировки, намеренно произведя сдвиг длины волны в диапазоне значений, которые могут ожидаться при эксплуатации. Этого можно достигнуть, произведя физические изменения монохроматора (в зависимости от конструкции используется сдвиг или изменение угла наклона решетки или фильтра, изменение положения отверстий на решетке, изменение наклона падающего излучения) или с помощью изменения констант, используемых для вычисления длин волн монохроматора. Результаты анализа вторичных свойств образца, полученные с использованием многомерной градуировочной модели, созданной в соответствии с данным способом, мало зависят от условий измерений и технических параметров измерительного прибора. Намеренное введение разброса данных результатов измерения образцов градуировочного набора повышает устойчивость модели, а область ее применимости становится более широкой.
Данный метод не лишен и недостатков. Прежде всего, необходимо отметить снижение точности результатов анализа с использованием подобной градуировочной модели, что обусловлено изначально предполагаемой вариацией результатов измерения первичных свойств. Это же обстоятельство затрудняет оценку примененимости градуировки для анализа того или иного неизвестного образца. Например, если градуировочная модель для определения процентного содержания различных химических компонентов в образце будет использована для анализа неизвестного образца, состоящего из других компонентов и в значительной степени отличающегося от образцов градуировочного набора, это приведет к ошибкам определения вторичных свойств. Кроме того, построение математической модели, малочувствительной к воздействию всех факторов, влияющих на результаты измерений, не всегда возможно, вследствие большого количества подобных факторов. Поэтому процесс создания градуировочной модели, учитывающей влияние большого количества дополнительных факторов, требует проведения огромного объема измерений градуировочных образцов при различных условиях, что еще более увеличивает трудоемкость и длительность процедуры создания градуировки.
Задачей настоящего изобретения является создание многомерной градуировочной модели, которая обеспечивает высокую точность определения анализируемых свойств и является устойчивой к изменениям свойств, влияющих на результаты измерений прибора, причем даже в случае многократного изменения подобных факторов, в частности, к возможным линейным и нелинейным отличиям технических параметров прибора, возникающим вследствие старения, выполнения ремонта или замены отдельных типовых элементов конструкции, изменению условий эксплуатации, а также прочих факторов, при этом снижается трудоемкость и длительность процедуры построения подобной модели, впоследствии применяемой для определения одного или нескольких вторичных свойств неизвестного образца по результатам измерения множества первичных свойств этого образца, причем не обязательно спектральных.
Решение поставленной задачи достигается способом создания многомерных градуировочных моделей, состоящим из последовательности действий, объединенных единым изобретательским замыслом. Способ включает в себя подбор градуировочного набора образцов с известными вторичными свойствами, определенными референтными методами; измерение на приборе первичных свойств каждого из образцов градуировочного набора; искусственное внесение изменений, по крайней мере, в одно из свойств, оказывающих влияние на результаты измерений прибора; измерение на приборе в измененном таким образом состоянии первичных свойств как минимум, одного образца; создание градуировочной модели, устойчивой к указанным изменениям, путем нахождения при помощи методов многомерного регрессионного анализа с использованием данных, полученных на приборе после внесения изменений, соотношений градуировки, позволяющих определить анализируемые вторичные свойства неизвестного образца по результатам измерений множества первичных свойств этого образца, отличающийся тем, что до внесения изменений формируют набор образцов для расчета корректирующих соотношений, измеряют первичные свойства каждого образца из этого набора на приборе до и после внесения изменений, и, сопоставляя при помощи методов многомерного регрессионного анализа результаты измерений первичных свойств образцов набора, полученные на приборе до внесения изменений, с результатами измерений первичных свойств тех же самых образцов, полученными на приборе в состоянии, когда внесены изменения, определяют корректирующие соотношения, причем найденные для какого-либо одного прибора корректирующие соотношения, учитывающие влияние внесенных изменений, используются на любом другом приборе данного типа при построении устойчивой к указанным изменениям градуировочной модели; осуществляют преобразование результатов измерения первичных свойств образцов rрадуировочного набора при помощи полученных корректирующих соотношений к виду, соответствующему измерениям на приборе после внесения изменений; скорректированными таким образом под измененное состояние прибора результатами измерения первичных свойств образцов rрадуировочного набора дополняют результаты измерения первичных свойств градуировочных образцов, полученные на приборе до внесения изменений, причем проверку многомерной градуировочной модели, созданной по совокупности исходных и скорректированных результатов измерения первичных свойств образцов rрадуировочного набора, производят, осуществляя процедуру валидации и используя определяемые на основании нее количественные параметры валидации градуировки.
Определение корректирующих соотношений на основе измерения сформированного в соответствии с данным методом представительного набора образцов позволяет учесть нелинейные отличия в характеристиках прибора, обусловленные изменением свойств, влияющих на результаты его измерений, т.к. для нахождения соотношений преобразования результатов измерений используется несколько зависимостей изменения первичных свойств. При этом нет необходимости проводить огромный объем измерений образцов rрадуировочного набора при различных условиях, что значительно снижает трудоемкость и длительность процесса создания градуировочной модели, устойчивой к изменениям факторов, влияющих на результаты измерений прибора, в соответствии с предложенным алгоритмом.
Изменения вносятся таким образом, чтобы охватить их возможный диапазон в процессе эксплуатации прибора. Поэтому построенная градуировочная модель будет устойчива даже к многократным изменениям влияющих на результаты анализа свойств, если их воздействие было учтено предложенной коррекцией градуировки на стадии ее создания. Этот факт также является существенным достоинством данного метода.
Все результаты измерений градуировочных образцов, в том числе и дополняющие градуировочный набор результаты измерения первичных свойств образцов rрадуировочного набора, преобразованные к виду, эквивалентному результатам измерений на приборе, когда внесены изменения по крайней мере в одно из влияющих на результаты его измерений свойств, сохраняются в компьютере данного градуируемого прибора. Это обеспечивает независимость созданной модели от измерений на приборе в исходном состоянии, и позволяет оценить правомерность применения градуировки для анализа неизвестного образца на приборе в состоянии, когда изменились одно или несколько влияющих на результаты измерений свойств, а также определить выпадающие данные при помощи статистики предсказания выбросов, что гарантирует высокую точность определения анализируемых вторичных свойств.
Для оценки точности анализа вторичных свойств неизвестных образцов по рассчитанной градуировочной модели на заключительном этапе проводится ее проверка (валидация). Данная процедура осуществляется путем сравнения вторичных свойств образцов, определенных по данной градуировке на основании результатов измерений прибора, с прямыми результатами измерения вторичных свойств, полученными с использованием референтных методов.
Использование статистики предсказания выбросов и исключение из модели выпадающих образцов градуировочного набора перед определением соотношений градуировки повышает устойчивость создаваемой градуировочной модели. Для преобразования результатов измерений и референтных данных к наиболее оптимальному виду можно использовать процедуру нормализации. Это позволяет минимизировать погрешность определения анализируемых вторичных свойств и учесть технические особенности измерительного прибора, а также различие пробоподготовки и состояния исследуемого образца. Процедура нормализации представляет собой выбор того или иного метода математической предобработки. Критерием выбора является точность анализа вторичных свойств образцов, которую обеспечивает прибор с градуировкой, при создании которой использовался данный вид математической предобработки. В качестве основных количественных критериев используются количественные параметры процедуры валидации градуировочной модели (например, стандартная ошибка валидации) [1].
Найденные корректирующие соотношения, учитывающие определенный вариант изменений свойств, влияющих на результаты измерений прибора, предлагается использовать с целью компенсации таких изменений при построении градуировки и для любого другого прибора данного типа, используя заявляемый алгоритм.
Сущность изобретения заключается в том, что предложенная совокупность признаков позволяет создать на градуируемом измерительном приборе многомерную градуировочную модель, устойчивую к возможным вариациям свойств, влияющих на результаты измерений прибора, в частности изменениям технических параметров прибора, возникающим вследствие старения, выполнения ремонта или замены отдельных типовых элементов конструкции, а также особенностям условий эксплуатации, причем даже в случае неоднократного изменения указанных факторов, при этом снижается трудоемкость и длительность процедуры построения подобной модели по сравнению с существующими известными методами. Созданная градуировочная модель предоставляет возможность с высокой точностью предсказывать вторичные свойства неизвестных образцов по результатам измерения множества первичных, не обязательно спектроскопических, свойств, причем градуировочная модель строится на основе данных по первичным свойствам образцов из градуировочного набора, измеренных на градуируемом приборе, и этих же данных, преобразованных к виду, как если бы измерения проводились на приборе, когда внесены некоторые изменения, по крайней мере, в одно из свойств, влияющих на результаты измерений прибора, например, соответствующие максимальным ожидаемым в процессе эксплуатации изменениям для приборов данного типа. Корректирующие соотношения, позволяющие преобразовать результаты измерений градуировочных образцов к виду, эквивалентному результатам измерений на приборе в измененном состоянии, определяются из измерений представительного набора образцов на приборе до и после внесения изменений, причем набор для расчета корректирующих соотношений состоит из гораздо меньшего количества образцов, чем градуировочный. Образцы из набора для расчета корректирующих соотношений обеспечивают существенные различия в результатах измерений во всем диапазоне первичных свойств на приборе, как в исходном, так и в измененном состоянии. Использование набора образцов для расчета корректирующих соотношений позволяет определить нелинейную связь между результатами измерений одних и тех же образцов на приборе до и после внесения изменений путем корреляционного анализа с использованием регрессионных методов. Найденные корректирующие соотношения, учитывающие определенный вариант изменений по крайней мере одного из факторов, влияющих на результаты измерений прибора, предлагается использовать при построении градуировочной модели с целью компенсации таких изменений и для любого другого прибора данного типа. Таким образом, существует возможность построения градуировки для любого другого прибора, устойчивой к отличиям технических параметров градуируемого прибора, возникающим вследствие старения, выполнения ремонта или замены отдельных типовых элементов конструкции, изменениям условий эксплуатации, а также прочих свойств, влияющих на результаты измерений прибора, используя корректирующие соотношения, рассчитанные только для одного прибора на основе предложенного алгоритма, что в еще большей степени упрощает процесс создания градуировочной модели. Так как при построении градуировочной модели используется совокупность результатов измерений первичных свойств образцов градуировочного набора, полученных на приборе в исходном состоянии, а также приведенных к виду прибора, когда по крайней мере в одно из влияющих на результаты измерений свойств внесены некоторые изменения, можно оценить правомерность применимости построенной градуировки для анализа неизвестного образца, используя статистические методы анализа выпадающих данных, например, по расстоянию Махаланобиса. Выполнение процедуры нормализации с использованием различных видов математической предобработки результатов измерений и данных по вторичным свойствам (референтных данных), позволяет обеспечить минимальную погрешность определения анализируемых вторичных свойств, определив, таким образом, оптимальные для данных условий соотношения градуировки. Заявляемое изобретение поясняется фиг.1 , где представлено схематическое изображение заявляемого способа в виде потоковой диаграммы. Заявляемый способ создания многомерных rрадуировочных моделей может быть использован для любых приборов, где анализируемые свойства образца определяются на основе многократного измерения других свойств, в частности, в спектрофотометрии для различного вида спектрометров, измеряющих поглощение светового излучения образцом на множестве различных длин волн. Данные, характеризующие результат подобных измерений, называются спектром. Рассмотрим применение заявляемого способа на примере спектрометров для анализа химического состава образца, заметив, однако, еще раз, что область применения заявляемого способа не ограничивается спектроскопией.
Спектрометр, как и любой аналитический прибор, требует предварительной градуировки, определяющей взаимосвязь между анализируемыми свойствами образца и спектральными характеристиками. Отметим, что часто с анализируемыми свойствами образца сопоставляются не непосредственно результаты измерений, а уже прошедшие процедуру нормализации (предварительную математическую предобработку) спектральные данные. Так, например, может проводиться сглаживание спектров, вычитание базовой линии или дифференцирование. Вид предварительной математической обработки выбирается исходя из критерия максимальной точности определения анализируемых вторичных свойств и сведения к минимуму влияния побочных факторов, связанных с паразитным рассеянием и особенностями пробоподготовки. Одна и та же математическая предобработка применяется ко всем спектрам образцов градуировочного набора, после проведения процедуры нормализации преобразованные спектральные данные имеют ярко выраженные характерные черты даже при незначительных изменениях анализируемых свойств [6]. В качестве критериев для оценки качества градуировочной модели используются различные статистические характеристики, например, стандартная ошибка градуировки (SEC), стандартная ошибка валидации (SEV) и стандартная ошибка кросс- валидации (SECV) [1]. Наиболее распространенным видом математической предобработки при спектральном анализе является нахождение средневзвешенных значений спектральных данных [6], уменьшающее на 1 число степеней свободы в градуировочной модели. В данной предобработке находят усредненный по всему градуировочному набору спектр и вычитают его из каждого спектра градуировочных образцов. Аналогично находят средневзвешенные значения референтных данных. Тогда при анализе неизвестного образца перед применением построенной градуировочной модели из измеренного спектра также вычитается усредненный по градуировочному набору спектр.
Для проведения градуировки отбирается набор образцов, представительный к тем образцам, которые будут в дальнейшем анализироваться. При выборе градуировочных образцов для спектрального анализа используются следующие критерии [1]: а) образцы должны содержать все химические компоненты, которые планируется анализировать; б) диапазон изменения концентрации анализируемых компонентов в образцах градуировочного набора должен превышать диапазон изменения в анализируемых неизвестных образцах; в) величины изменения в концентрации химических компонентов от образца к образцу должны быть равномерно распределены по всему диапазону изменений; г) число образцов должно обеспечивать нахождение при помощи статистических методов математических соотношений между спектроскопическими данными и концентрацией отдельных химических компонентов. Образцы, выпадающие из градуировочного набора, определяются с помощью статистического анализа выбросов, например, путем вычисления расстояния Махаланобиса [1], которое определяется как:
D2 = r'(R -R')+r (1) где R - матрица спектральных данных полного градуировочного набора, г - вектор, соответствующий спектру одного образца. Расстояние Махаланобиса показывает, сколько степеней свободы вносит в градуировочную модель данный образец. В среднем каждый градуировочный образец должен вносить k/m, где k - число переменных в регрессии, m - число образцов в градуировочном наборе. Образцы с D2>3k/m должны быть исключены из градуировочного набора. Большое значение расстояния Махалонобиса означает, что спектр данного образца почти полностью определяет один из регрессионных коэффициентов, что делает модель неустойчивой. Такое может произойти, когда, например, нарушена равномерность распределения анализируемых свойств rрадуировочных образцов по диапазону, в котором они изменяются, т.е. когда состав образца существенно отличается от остальных образцов в градуировочном наборе. Из градуировочного набора должны быть также исключены образцы, значения анализируемых свойств которых, определенные при помощи построенной модели, значительно отличаются от значений, которые дает референтный метод. Данные образцы определяются из расхождений по Стьюденту, рассчитанных по формуле:
Figure imgf000020_0001
здесь βj - отличие полученной с помощью градуировочной модели величины концентрации химического компонента или анализируемого свойства от референтного значения для i-го градуировочного образца, SEC - стандартная ошибка градуировки [1], Dj2 - расстояние Махаланобиса для i-го градуировочного образца. Расхождения по Стьюденту должны быть равномерно распределены по нормальному закону. Величина расхождения сравнивается с коэффициентом Стьюдента, для доверительной вероятности
0,95 и числа степеней свободы m-k. В случае если величина расхождения превышает коэффициент, образец исключается из градуировочного набора.
Анализируемые свойства rрадуировочных образцов должны быть заранее известны или определяются при помощи референтных методов. Полученные численные значения считаются реальными значениями анализируемых свойств, поэтому точность референтного анализа определяет точность всей градуировки, и к этой процедуре предъявляются достаточно высокие требования. Повысить точность референтного анализа можно, проводя многократное повторение анализов с последующим усреднением результатов.
Спектры образцов градуировочного набора регистрируются на приборе, после чего может быть создана многомерная градуировочная модель. Для этого используются мультивариантные математические методы, такие как, анализ принципиальных компонентов (PCA), метод дробных наименьших квадратов (PLS) и т.п. В дальнейшем, измеряя спектр неизвестного образца на данном приборе и используя созданную градуировочную модель, можно определить анализируемые свойства, например, химический состав. Точность подобных измерений будет достаточно высокой до тех пор, пока свойства, оказывающие влияние на результаты измерений прибора, будут оставаться постоянными. Но, как уже упоминалось ранее, в процессе эксплуатации, при проведении ремонтных работ или замены отдельных типовых элементов конструкции могут изменяться, например, технические характеристики прибора, что в свою очередь может привести к снижению точности предсказаний и необходимости построения новой градуировочной модели. Следовательно, для обеспечения достаточной точности предсказания свойств неизвестных образцов на приборе при условии возможных многократных изменений влияющих на результаты измерений факторов необходимо иметь градуировочную модель, устойчивую к различиям в характеристиках прибора, изменениям условий эксплуатации и прочих свойств, влияющих на результаты измерений прибора, при этом дающую заданную точность анализа, позволяющую оценить применимость модели для анализа неизвестного образца, причем процесс создания подобной градуировки не должен быть чрезмерно длительным и трудоемким.
Заявляемый способ позволяет создавать многомерную градуировочную модель для градуируемого прибора, устойчивую к возможным изменениям одного или нескольких факторов, влияющих на результаты его измерений, не проводя измерения полного набора градуировочных образцов на приборе, когда указанный фактор или факторы изменились, причем это справедливо даже в случае многократных изменений. Для создания градуировки используются спектральные данные, измеренные на градуируемом приборе, дополненные спектральными данными, полученными в результате преобразования исходных данных к виду, эквивалентному измерениям на приборе в состоянии, когда изменено по крайней мере одно из влияющих на результаты измерений свойств, например, таким образом, чтобы охватить ожидаемый диапазон подобных изменений для приборов данного типа в процессе эксплуатации. Область применимости и устойчивость градуировочной модели анализируются на основании совокупности этих спектральных данных. Так как все данные сохраняются в компьютере градуируемого прибора, это позволяет производить оценку выпадающих данных при анализе неизвестных образцов, например, с помощью статистики Махаланобиса.
Для коррекции спектральных данных градуировочного набора, проводимой с целью учета изменений свойств, влияющих на результаты измерений прибора, используется специально подобранный набор образцов, именуемый в дальнейшем как набор для расчета корректирующих соотношений, число образцов в котором намного меньше, чем в полном градуировочном наборе, при этом их свойства могут быть неизвестны, важно лишь то, что этот набор образцов обеспечивает значительные вариации в измеряемых спектральных данных, позволяющие найти выражения преобразования. Для этого спектр каждого образца набора измеряется на приборе в состоянии, соответствующем измерению спектров градуировочных образцов, и этом же приборе, когда по крайней мере в одно из свойств, оказывающих влияние на результаты измерений, внесены некоторые изменения («измeнeннoм» приборе). Проводя корреляцию полученных спектральных данных, находят соотношения, позволяющие преобразовать спектры, измеренные на приборе до внесения изменений, к виду, как если бы измерения проводились на измененном приборе, и учитывающие различия в результатах измерений одних и тех же образцов, возникающие вследствие изменения одного или нескольких факторов, влияющих на результаты измерений прибора. С целью получения оптимальной градуировки спектральные данные могут быть подвергнуты нормализации, заключающейся в проведении одинаковых математических преобразований для всех измеренных спектров. Это обеспечивает выявление явных отличий в спектральных данных, измеренных на приборе до и после внесения изменений, что в свою очередь способствует более точному определению выражений для преобразования спектральных данных. Необходимо отметить, что найденные подобным образом для какого-либо одного прибора корректирующие соотношения, учитывающие влияние на результаты измерений прибора изменений определенного свойства или свойств, могут быть использованы на любом другом приборе данного типа при построении градуировочной модели, устойчивой к указанным изменениям. Совокупность перечисленных признаков позволяет снизить трудоемкость и длительность процесса создания градуировочной модели, устойчивой к изменениям свойств, влияющих на результаты измерений прибора. Градуировочную модель для градуируемого прибора создают по совокупности исходных и преобразованных спектральных данных градуировочных образцов, используя стандартные математические методы многомерного регрессионного анализа (MLR, PCA, PLS и т.п. [6]), причем из градуировки должны быть исключены выпадающие образцы, что гарантирует устойчивость созданной модели, после чего ее можно использовать для определения свойств неизвестного образца.
В качестве практического пояснения заявляемого изобретения рассмотрим пример создания многомерной градуировочной модели, устойчивой к изменениям технических параметров спектрометра ИнфраЛЮМ ФT-10, которые могут произойти вследствие старения, выполнения ремонта или замены отдельных типовых элементов конструкции. Принцип работы данного прибора основан на Фурье спектроскопии в ближней инфракрасной области спектра (БИК). В качестве объекта исследований выбрана продовольственная пшеница, для которой с помощью построенной градуировки на данном приборе определялись основные показатели качества. Подчеркнем при этом еще раз, что данный пример используется лишь для практической иллюстрации предлагаемого способа. Для создания градуировки отбирается представительный набор образцов, спектры которых измеряются на спектрометре. Полученные спектральные данные могут быть подвергнуты процедуре нормализации, учитывающей особенности приборов, работающих на пропускание и использующих принципы Фурье спектроскопии [6]. При построении градуировочной модели в рамках данных исследований использовались следующие предобработки: выравнивание базовой линии, нормирование спектров по среднеквадратичному отклонению, расчет средневзвешенных значений [6].
Как упоминалось ранее, вид математической обработки определяется таким образом, чтобы обеспечить минимальную погрешность определения анализируемых свойств. Если спектральные данные подвергнуты процедуре нормализации, то соотношения градуировки определяются на основе сопоставления этих данных с известными, также прошедшими нормализацию, свойствами образцов градуировочного набора, определенными референтными методами.
Для оценки качества созданной градуировочной модели используются различные статистические характеристики [1].
Одним из таких параметров является стандартная ошибка градуировки (SEC), которая служит для оценки ожидаемой согласованности между величинами, полученными, используя модель градуировки, и величинами, измеренными референтным методом.
Figure imgf000024_0001
где у - вектор значений определяемых показателей для градуировочных образцов, полученных при применении этой градуировки; у - вектор референтных данных; d=п-k - число степеней свободы градуировочной модели, п - число градуировочных образцов, к - число переменных величин в градуировочной модели, зависящее от математического метода, использующегося для построения модели.
Стандартная ошибка кросс валидации (SECV) вычисляется при перекрестной проверке градуировочной модели и позволяет оценить ее устойчивость [1].
Figure imgf000024_0002
где у - вектор, содержащий оценки перекрестной валидации.
Стандартная ошибка валидации (SEV) характеризует ошибку отклонения между референтными и предсказанными по градуировочному уравнению значениями для образцов дополнительного набора, не входивших в градуировочный.
Figure imgf000024_0003
где dv- общее число референтных значений анализируемого показателя для всех спектров дополнительного набора, у. - референтные значения анализируемого показателя для i-го спектра дополнительного набора, у.- предсказанные значения анализируемого показателя для i-го спектра дополнительного набора.
Основные статистические характеристики исходной градуировочной модели, созданной по образцам продовольственной пшеницы на исследуемом спектрометре без искусственного внесения изменений в свойства, влияющие на результаты измерений, приведены в Taб.1. Градуировки были построены на два основных показателя качества - протеин и клейковину. Для определения SEV использовался валидационный набор, который выбирался таким образом, чтобы содержание в нем соответственно протеина или клейковины было равномерно распределено по всему диапазону возможных значений концентрации этих показателей. Таблица 1. Результаты градуировки исходного прибора.
Figure imgf000025_0001
Чтобы создать на спектрометре многомерную градуировочную модель, учитывающую возможные изменения свойств, оказывающих влияние на результаты измерений прибора, например, отличия технических параметров прибора, возникающие вследствие старения, выполнения ремонта или замены отдельных типовых элементов конструкции, необходимо в соответствии с заявляемым изобретением сформировать набор образцов для расчета корректирующих соотношений между результатами, получаемыми на приборе в исходном состоянии, и в состоянии, когда внесены изменения. В данном случае из градуировочного набора было отобрано 10 образцов с максимальным и минимальным значениями параметра sсоrе для каждого из анализируемых показателей. Однако отметим, что образцы в наборе для расчета корректирующих соотношений в общем случае могут не принадлежать градуировочному набору. Спектры выбранных образцов регистрировались на градуируемом приборе. Затем в конструктивные особенности данного прибора искусственно вносились некоторые изменения, например, те, которые возможны при ремонте анализатора. В частности, в приводимом примере производилась замена светоделителя в интерферометре на другой, покрытие которого отличается от исходного настолько, что вызывает максимально возможные отклонения по распределению светового пучка среди светоделителей данного класса. После этого выполнялись измерения образцов для расчета корректирующих соотношений на измененном таким образом приборе и, путем корреляции спектральных данных, полученных на приборе до и после внесения изменений, находились выражения для преобразования результатов измерений на приборе в исходном состоянии к виду, соответствующему измерениям на приборе после внесения изменений. В простейшем виде эти соотношения могут быть определены методом линейной регрессии.
Figure imgf000026_0001
где Rf1 - значения спектральных данных, измеренные на приборе, когда в его конструкцию внесены изменения (i-ая длина волны, j-й образец из набора для расчета корректирующих соотношений), i?,'" - аналогичные спектральные данные, измеренные на приборе до внесения изменений. Спектральные данные могут подвергаться процедуре нормализации, при этом используемая математическая обработка должна быть одинаковой для спектров, измеряемых на приборе как до, так и после внесения изменений. Для определения корректирующих соотношений методом наименьших квадратов находятся коэффициенты регрессии.
Figure imgf000026_0002
где с - количество образцов в наборе для расчета корректирующих соотношений. После нахождения корректирующих соотношений, используя найденные регрессионные коэффициенты и а" , спектральные данные для каждого образца из градуировочного набора в соответствии с формулой (6) преобразуются к виду, соответствующему измерениям на приборе после внесения изменений. Далее создается новая многомерная градуировочная модель по совокупности исходных и преобразованных данных градуировочного набора, которая проверяется с помощью стандартной процедуры валидации [1], в результате чего определяются ее основные статистические параметры. Для проверки эффективности предложенного метода на исследуемом приборе устанавливались несколько различных светоделителей, и в каждом случае регистрировались спектры валидационного набора образцов. После этого на основе стандартной процедуры валидации [1] оценивались основные статистические параметры исходной и скорректированной градуировочных моделей. В таблице 2 приведены данные для исходной градуировки, а в таблице 3 для скорректированной, устойчивой к влиянию замены светоделителя.
Из приведенных данных видно, что созданная в соответствии с заявляемым изобретением многомерная градуировочная модель позволяет с высокой точностью определять свойства неизвестных образцов и при этом значительно менее чувствительна к учтенному с помощью предложенного алгоритма возможному изменению конструкции прибора, в данном случае замене светоделителя. Подобным образом могут быть учтены любые другие изменения свойств, влияющих на результаты измерений прибора, приведенный здесь пример выбран в качестве наглядной иллюстрации заявляемого изобретения вследствие того, что подобные изменения являются одними из наиболее характерных для спектрометров данного типа, и способны привести к значительному снижению точности анализа по изначально построенной градуировочной модели прибора, если не будут учтены при ее создании. Таблица 2. Проверка исходной градуировки на приборе с различными светодел ителя м и .
Figure imgf000028_0001
Figure imgf000029_0001
Таблица 3 Проверка градуировки, скорректированной с учетом замены светоделителя, на приборе с различными светоделителями.
Figure imgf000029_0002
Figure imgf000030_0001
Была проведена проверка возможности использования найденных корректирующих соотношений на другом приборе той же модели. На этом приборе были зарегистрированы спектры выбранного для него градуировочного набора образцов продовольственной пшеницы, и на основании этих данных построена градуировочная модель, которая не учитывает возможного изменения свойств, оказывающих влияние на результаты измерений прибора.
После этого спектральные данные, полученные при измерении на градуируемом приборе градуировочного набора образцов, были скорректированы по найденным ранее корректирующим соотношениям к виду, соответствующему прибору, для которого проведена замена светоделителя. В соответствии с предложенным алгоритмом, была создана градуировочная модель по совокупности исходных и скорректированных спектральных данных градуировочных образцов прибора. Для проверки эффективности предложенного метода на градуируемом приборе устанавливались несколько различных светоделителей, и в каждом случае регистрировались спектры валидационного набора образцов. После этого на основе стандартной процедуры валидации [1] оценивались основные статистические параметры исходной градуировочной модели, построенной для данного прибора, и градуировочной модели, устойчивой к влиянию замены светоделителя, на основе корректирующих соотношений, рассчитанных ранее для другого прибора. Полученные данные приведены соответственно в таблицах 4 и 5.
Таблица 4. Проверка градуировки второго прибора на приборе с различными светоделителями.
Figure imgf000032_0001
Figure imgf000033_0001
Таблица 5. Проверка градуировки, скорректированной с учетом замены светоделителя по корректирующим соотношениям, определенным для другого прибора, на градуируемом приборе с различными светоделителями.
Figure imgf000033_0002
Figure imgf000034_0001
Результаты таблиц 4 и 5 подтверждают возможность использования корректирующих соотношений, определенных для какого-либо одного прибора и учитывающих влияние изменений его параметров или иных свойств, влияющих на результаты измерений, на любого другом приборе данной серии при построении устойчивой к указанным изменениям градуировочной модели. Данная особенность позволяет после создания в соответствии с предложенным способом градуировочной модели, компенсирующей отличия определенного влияющего на результаты измерений свойства или свойств, например, технических параметров прибора, изменения которых происходят вследствие старения, выполнения ремонта или замены отдельных типовых элементов конструкции, для какого-либо одного прибора, в дальнейшем использовать найденные корректирующие соотношения для всех остальных приборов данного типа при построении градуировочных моделей, учитывающих с помощью данных соотношений коррекции соответствующие изменения влияющих свойств. В таком случае нет необходимости выбирать и проводить измерения набора образцов для расчета корректирующих соотношений для всех остальных приборов, кроме первого, поскольку спектральные данные, полученные при регистрации градуировочных образцов на остальных приборах, преобразуются математически с использованием найденных корректирующих соотношений к виду, соответствующему прибору, когда внесены некоторые изменения в одно или несколько свойств, влияющих на результаты измерений.
В заключение еще раз отметим, что область применения заявляемого способа не ограничивается спектроскопией. Заявляемое изобретение может быть использовано для различных приборов, определяющих некоторые свойства образца на основе многократного измерения других свойств.
Источники информации 1. ASTM stапdаrd, E 1655 - 00, Рrасtiсеs fоr Iпfrаrеd Мultivаriаtе
Quапtitаtivе Апаlуsis.
2. Патентная заявка США NQ 4 944 589, МПК GOU 3/18, опубликована 31. 07. 1990
3. Патентная заявка США NQ 6 615 151 , МПК G01 N 015/06, опубликована 02. 09. 2003
4. Европейская патентная заявка EP 0 225 233, МПК G01N 021/27E публикована 28.03.2002
5. Российская патентная заявка NQ 2266523 C1 , МПК 7: G01N 021/01 , опубликована 20.12.2005 6. Руководство по эксплуатации ИнфраЛЮМ ФT-10, 152.00.00.00. РЭ.

Claims

Формула изобретения.
1. Способ создания многомерных градуировочных моделей, устойчивых к изменениям свойств, влияющих на результаты измерений прибора, включающий в себя подбор градуировочного набора образцов с известными вторичными свойствами; измерение на приборе первичных свойств каждого из образцов градуировочного набора; искусственное внесение изменений, по крайней мере, в одно из свойств, оказывающих влияние на результаты измерений прибора; измерение первичных свойств, как минимум, одного образца на приборе в измененном таким образом состоянии, отличающийся тем, что до внесения изменений формируют набор образцов для расчета корректирующих соотношений, измеряют первичные свойства каждого образца из этого набора на приборе до и после внесения изменений, и, сопоставляя при помощи методов многомерного регрессионного анализа результаты измерений первичных свойств образцов набора, полученные на приборе до внесения изменений, с результатами измерений первичных свойств тех же самых образцов, полученными на приборе в состоянии, когда внесены изменения, определяют корректирующие соотношения; осуществляют преобразование результатов измерения первичных свойств образцов градуировочного набора при помощи полученных корректирующих соотношений к виду, соответствующему измерениям на приборе после внесения изменений; скорректированными таким образом под измененное состояние прибора результатами измерения первичных свойств образцов градуировочного набора дополняют результаты измерения первичных свойств градуировочных образцов, полученные на приборе до внесения изменений, рассчитывают многомерную градуировочную модель по совокупности исходных и скорректированных результатов измерения первичных свойств образцов градуировочного набора, и производят ее проверку, определяя на основании нее количественные параметры валидации градуировки.
2. Способ по п. 1 , отличающийся тем, что, осуществляя сопоставление при помощи методов многомерного регрессионного анализа результатов измерений первичных свойств образцов набора для расчета корректирующих соотношений, полученных на приборе до внесения изменений, с результатами измерения первичных свойств тех же самых образцов, полученными на приборе в состоянии, когда внесены изменения, результаты измерений указанных образцов на приборе как до, так и после внесения изменений, предварительно подвергают процедуре нормализации, выявляя различие в результатах измерений образцов набора для расчета корректирующих соотношений в зависимости от внесенных изменений, причем, осуществляя преобразование результатов измерения первичных свойств образцов градуировочного набора при помощи полученных корректирующих соотношений к виду, соответствующему измерениям на приборе после внесения изменений, точно такой же процедуре нормализации предварительно подвергают результаты измерений образцов градуировочного набора на приборе до внесения изменений.
3. Способ по п. 1 , отличающийся тем, что в набор для расчета корректирующих соотношений подбирают образцы, характеризующие существующий диапазон различий в результатах измерений первичных свойств на градуируемом приборе как до, так и после внесения изменений.
4. Способ по п. 1 , отличающийся тем, что в набор для расчета корректирующих соотношений подбирают образцы, свойства которых характеризуют существующий диапазон изменений вторичных свойств образцов.
5. Способ по п.1 , отличающийся тем, что если при проведении проверки многомерной градуировочной модели, созданной по совокупности исходных и скорректированных результатов измерения первичных свойств образцов градуировочного набора, найденные количественные параметры валидации градуировки не соответствуют заданным критериям точности, совокупность исходных и скорректированных результатов измерения первичных свойств образцов градуировочного набора анализируют на присутствие выпадающих данных при помощи статистики предсказания выбросов, и перед определением соотношений градуировки исключают из градуировочного набора выпадающие образцы.
6. Способ по п.1 , отличающийся тем, что если при проведении проверки многомерной градуировочной модели, созданной по совокупности исходных и скорректированных результатов измерения первичных свойств образцов градуировочного набора, найденные количественные параметры валидации градуировки не соответствуют заданным критериям точности, совокупность исходных и скорректированных результатов измерения первичных свойств образцов градуировочного набора и (или) соответствующих известных вторичных свойств градуировочных образцов подвергают процедуре нормализации.
7. Способ по п.1, отличающийся тем, что созданная многомерная градуировочная модель обеспечивает линейную взаимосвязь между первичными и вторичными свойствами.
8. Способ по п.1 , отличающийся тем, что созданная многомерная градуировочная модель обеспечивает нелинейную взаимосвязь между первичными и вторичными свойствами.
9. Способ по п.1 , отличающийся тем, что в качестве градуируемого прибора используют спектрометр, причем первичными свойствами являются результаты измерения спектральных характеристик образцов, описывающие поглощение, отражение или рассеяние света при разных значениях волновых чисел.
PCT/RU2007/000317 2006-06-20 2007-05-31 Procédé de création de modèles de graduation multidimensionnels Ceased WO2008002192A1 (fr)

Priority Applications (5)

Application Number Priority Date Filing Date Title
EA200900029A EA012950B1 (ru) 2006-06-20 2007-05-31 Способ создания многомерных градуировочных моделей, устойчивых к изменениям свойств, влияющих на результаты измерений прибора
CN200780023164.4A CN101473197B (zh) 2006-06-20 2007-05-31 制造多维校准模型的方法
CA2655745A CA2655745C (en) 2006-06-20 2007-05-31 Method of formation of multidimensional calibration models stable to property changes affecting the instrument's measurment results
EP20070794046 EP2040042B1 (en) 2006-06-20 2007-05-31 Method for producing multidimensional calibrating patterns
US12/305,122 US20090204352A1 (en) 2006-06-20 2007-05-31 Method for producing multidimensional calibrating patterns

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
RU2006122456/28A RU2308684C1 (ru) 2006-06-20 2006-06-20 Способ создания многомерных градуировочных моделей, устойчивых к изменениям свойств, влияющих на результаты измерений прибора
RU2006122456 2006-06-20

Publications (1)

Publication Number Publication Date
WO2008002192A1 true WO2008002192A1 (fr) 2008-01-03

Family

ID=38845851

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/RU2007/000317 Ceased WO2008002192A1 (fr) 2006-06-20 2007-05-31 Procédé de création de modèles de graduation multidimensionnels

Country Status (8)

Country Link
US (1) US20090204352A1 (ru)
EP (1) EP2040042B1 (ru)
CN (1) CN101473197B (ru)
CA (1) CA2655745C (ru)
EA (1) EA012950B1 (ru)
RU (1) RU2308684C1 (ru)
UA (1) UA96946C2 (ru)
WO (1) WO2008002192A1 (ru)

Families Citing this family (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8735820B2 (en) * 2010-04-09 2014-05-27 Tesoro Refining And Marketing Company Direct match spectrographic determination of fuel properties
RU2462701C1 (ru) * 2011-04-08 2012-09-27 Государственное образовательное учреждение высшего профессионального образования Омский государственный университет путей сообщения Способ построения устойчивой градуировочной зависимости при определении количественного состава элементов в цинковых сплавах
GB201203108D0 (en) * 2012-02-23 2012-04-04 Beers Centenary De Ag Calibration of measuring instruments
RU2541906C1 (ru) * 2013-07-18 2015-02-20 Общество с ограниченной ответственностью "ВИНТЕЛ" Способ создания многомерных градуировочных моделей аналитического прибора
CN103837484B (zh) * 2014-02-24 2016-08-17 广西科技大学 一种用于消除光谱乘性随机误差的角度化多变量分析方法
RU2625642C1 (ru) * 2016-08-09 2017-07-17 Акционерное общество "Концерн "Центральный научно-исследовательский институт "Электроприбор" Способ одновременной калибровки трех и более однотипных устройств с измерительными функциями без опоры на эталонное устройство или эталонный испытательный сигнал
JP6863341B2 (ja) * 2018-06-28 2021-04-21 横河電機株式会社 フィールド機器、フィールド機器の診断方法および診断装置
DE102018122411A1 (de) * 2018-09-13 2020-03-19 Endress+Hauser SE+Co. KG Verfahren zur Verbesserung der Messperformance von Feldgeräten der Automatisierungstechnik
CN113588572B (zh) * 2021-08-04 2024-03-19 广州市华南自然资源科学技术研究院 一种农田重金属在线检测校正模型智能管理系统
CN115452929B (zh) * 2022-09-30 2023-04-21 上海立迪生物技术股份有限公司 一种成像质谱流式的信号校准方法
CN115656910B (zh) * 2022-12-27 2023-04-21 太原山互科技有限公司 互感器校验仪器远程校准系统、方法及装备
CN117848973B (zh) * 2024-03-07 2024-05-28 铜川市人民医院 基于抗感染临床药学的药品成分智能检测方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0225233A1 (fr) 1985-10-29 1987-06-10 Plastimo Dispositif de fixation et de guidage d'une voile comportant une ralingue sur un mât comportant un tunnel de ralingue
US4944589A (en) 1987-12-08 1990-07-31 Tecator Ab Method of reducing the susceptibility to interference of a measuring instrument
WO2002025233A2 (en) * 2000-09-18 2002-03-28 Sensys Medical, Inc. Method of characterizing spectrometers and providing calibration models
US6615151B1 (en) 2000-08-28 2003-09-02 Cme Telemetrix Inc. Method for creating spectral instrument variation tolerance in calibration algorithms
RU2266523C1 (ru) 2004-07-27 2005-12-20 Общество с ограниченной ответственностью ООО "ВИНТЕЛ" Способ создания независимых многомерных градуировочных моделей

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5459677A (en) * 1990-10-09 1995-10-17 Board Of Regents Of The University Of Washington Calibration transfer for analytical instruments
CA2146342A1 (en) * 1992-10-07 1994-04-14 Steven M. Maggard Spectroscopic instrument calibration
US5610836A (en) * 1996-01-31 1997-03-11 Eastman Chemical Company Process to use multivariate signal responses to analyze a sample
DE19810917A1 (de) * 1998-03-13 1999-09-16 Buehler Ag Automatisches Kalibrationsverfahren
US6157041A (en) * 1998-10-13 2000-12-05 Rio Grande Medical Technologies, Inc. Methods and apparatus for tailoring spectroscopic calibration models
US20030135547A1 (en) * 2001-07-23 2003-07-17 Kent J. Thomas Extensible modular communication executive with active message queue and intelligent message pre-validation

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0225233A1 (fr) 1985-10-29 1987-06-10 Plastimo Dispositif de fixation et de guidage d'une voile comportant une ralingue sur un mât comportant un tunnel de ralingue
US4944589A (en) 1987-12-08 1990-07-31 Tecator Ab Method of reducing the susceptibility to interference of a measuring instrument
US6615151B1 (en) 2000-08-28 2003-09-02 Cme Telemetrix Inc. Method for creating spectral instrument variation tolerance in calibration algorithms
WO2002025233A2 (en) * 2000-09-18 2002-03-28 Sensys Medical, Inc. Method of characterizing spectrometers and providing calibration models
RU2266523C1 (ru) 2004-07-27 2005-12-20 Общество с ограниченной ответственностью ООО "ВИНТЕЛ" Способ создания независимых многомерных градуировочных моделей

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
"Practices for Infrared Multivariate Quantitative Analysis", ASTM STANDARD, E 1655 - 00
See also references of EP2040042A4 *

Also Published As

Publication number Publication date
EP2040042B1 (en) 2015-05-13
EA200900029A1 (ru) 2009-06-30
CA2655745C (en) 2015-07-28
RU2308684C1 (ru) 2007-10-20
EP2040042A4 (en) 2010-09-08
EP2040042A1 (en) 2009-03-25
US20090204352A1 (en) 2009-08-13
EA012950B1 (ru) 2010-02-26
UA96946C2 (ru) 2011-12-26
CN101473197A (zh) 2009-07-01
CN101473197B (zh) 2016-01-20
CA2655745A1 (en) 2008-01-03

Similar Documents

Publication Publication Date Title
WO2008002192A1 (fr) Procédé de création de modèles de graduation multidimensionnels
Workman A review of calibration transfer practices and instrument differences in spectroscopy
RU2266523C1 (ru) Способ создания независимых многомерных градуировочных моделей
Ye et al. Non-destructive prediction of protein content in wheat using NIRS
Cogdill et al. Process analytical technology case study, part III: calibration monitoring and transfer
EP2710353B1 (en) SPECTROSCOPIC APPARATUS AND METHOD of DETERMINING COMPONENTS PRESENT IN A SAMPLE
EP0954744B1 (en) Calibration method for spectrographic analyzing instruments
EP0552291A1 (en) Method of estimating property and/or composition data of a test sample
JPH08505221A (ja) 分光装置の較正
Westerhaus et al. Quantitative analysis
CN111563436B (zh) 一种基于ct-cdd的红外光谱测量仪器标定迁移方法
Broad et al. Guidelines for the development and validation of near-infrared spectroscopic methods in the pharmaceutical industry
Fan et al. Direct calibration transfer to principal components via canonical correlation analysis
Allegrini et al. Generalized error-dependent prediction uncertainty in multivariate calibration
Mark et al. Validation of a near-infrared transmission spectroscopic procedure, part A: validation protocols
US8548764B1 (en) Method and system for increasing optical instrument calibration and prediction accuracy within and across different optical instrument platforms
Chen et al. Investigation of sample partitioning in quantitative near-infrared analysis of soil organic carbon based on parametric LS-SVR modeling
JPH07151677A (ja) 濃度計
Momeni-Boroujeni et al. Systematic error detection in laboratory medicine
Workman Jr The essential aspects of multivariate calibration transfer
Nieuwoudt et al. Routine monitoring of instrument stability in a milk testing laboratory with ASCA: A pilot study
CN118777255B (zh) 一种基于近红外光谱的草甘膦产品及其副产物检测方法
Yun Model Evaluation
Reggente et al. Analysis of functional groups in atmospheric aerosols by infrared spectroscopy: functional group quantification in US measurement networks
Igne et al. 17 Calibration Transfer

Legal Events

Date Code Title Description
WWE Wipo information: entry into national phase

Ref document number: 200780023164.4

Country of ref document: CN

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

Ref document number: 07794046

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 12305122

Country of ref document: US

ENP Entry into the national phase

Ref document number: 2655745

Country of ref document: CA

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 200900029

Country of ref document: EA

WWE Wipo information: entry into national phase

Ref document number: 2007794046

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: RU