CN109520941A - The receptance function bearing calibration of online spectrum measurement instruments - Google Patents
The receptance function bearing calibration of online spectrum measurement instruments Download PDFInfo
- Publication number
- CN109520941A CN109520941A CN201811381410.XA CN201811381410A CN109520941A CN 109520941 A CN109520941 A CN 109520941A CN 201811381410 A CN201811381410 A CN 201811381410A CN 109520941 A CN109520941 A CN 109520941A
- Authority
- CN
- China
- Prior art keywords
- signal
- spectral
- absorption
- correction
- spectrum
- 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.)
- Granted
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 83
- 238000005259 measurement Methods 0.000 title claims abstract description 45
- 238000012937 correction Methods 0.000 claims abstract description 77
- 230000003595 spectral effect Effects 0.000 claims abstract description 76
- 239000000126 substance Substances 0.000 claims abstract description 60
- 238000000034 method Methods 0.000 claims abstract description 42
- 238000005316 response function Methods 0.000 claims abstract 15
- 238000010521 absorption reaction Methods 0.000 claims description 73
- 239000011159 matrix material Substances 0.000 claims description 12
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 9
- 230000003287 optical effect Effects 0.000 claims description 9
- 238000012935 Averaging Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 4
- 229910052757 nitrogen Inorganic materials 0.000 claims description 4
- 238000009499 grossing Methods 0.000 claims description 2
- 238000001675 atomic spectrum Methods 0.000 claims 1
- 150000003384 small molecules Chemical class 0.000 abstract description 2
- 239000007789 gas Substances 0.000 description 10
- 238000000862 absorption spectrum Methods 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- WEVYAHXRMPXWCK-UHFFFAOYSA-N Acetonitrile Chemical compound CC#N WEVYAHXRMPXWCK-UHFFFAOYSA-N 0.000 description 4
- 238000001514 detection method Methods 0.000 description 4
- 238000013329 compounding Methods 0.000 description 3
- 230000007613 environmental effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- 238000004611 spectroscopical analysis Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000000295 emission spectrum Methods 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 229910001873 dinitrogen Inorganic materials 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000001285 laser absorption spectroscopy Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011946 reduction process Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/27—Colour; 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/274—Calibration, base line adjustment, drift correction
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本公开提供一种在线光谱测量仪器的响应函数校正方法。利用标准物质在数据库中标准光谱线形数据，对实验测量中获得的混合物质光谱信号中的标准物质的线形进行对比，从而实现利用小分子物质光谱线形的多维信息，对混合物质信号整体线形以及特殊形状光谱进行校正，从而实现对光谱测量仪器响应函数的校正，无需额外设备且误差小，适应性强。
The disclosure provides a response function calibration method of an online spectrum measuring instrument. Use the standard spectral line shape data of the standard substance in the database to compare the line shape of the standard substance in the spectral signal of the mixed substance obtained in the experimental measurement, so as to realize the multidimensional information of the spectral line shape of the small molecule substance, and analyze the overall line shape of the mixed substance signal and the special The shape spectrum is corrected, so as to realize the correction of the response function of the spectrum measuring instrument, no additional equipment is required, the error is small, and the adaptability is strong.
Description
Technical field
This disclosure relates to spectral measurement analysis technical field more particularly to a kind of receptance function of online spectrum measurement instruments
Bearing calibration.
Background technique
Spectrum measurement instruments be it is a kind of can analyze light to obtain the common measuring device of spectral information, be now widely used in
Various fields, such as gas detection, biotechnology, laser technology field.But due to spectrum measurement instruments or the response letter of system
The factors such as non-linear of the noises of several power swing, environment or instruments and system with input signal, detector and signal path
It is closely related, so as to cause spectral signal there are baseline drift and the fluctuation of amplitude or distortion, cause to detection accuracy and
The influence of sensitivity etc..Thus it is guaranteed that and improve the detection accuracy and sensitivity of spectrum measurement instruments or system, need pair
Its receptance function is corrected.But due to the receptance function of spectrum measurement instruments cannot directly be obtained and school only,
Utilize the correction realized to the linear correction of measure spectrum to the receptance function of spectrum measurement instruments.
Present most bearing calibration needs to be arranged some empirical parameters, and this reduces correction side to a certain extent
The generality and adaptability of method, and be unable in correction signal influence spectrum bring due to the non-linear of detector itself, it is difficult
To meet to the high-precision demand of spectrum.For example, Fast Fourier Transform (FFT) method needs artificial determining bound, and spend a large amount of
Time calculated.Under many real-time applications, it is difficult to handle original spectrum one by one in face of a large amount of data.
Common polynomial fitting assumes that baseline can be approximated to lower order polynomial expressions form, goes to intend with a suitable multinomial
Data are closed, gives up any point for being greater than a specific criteria difference in matched curve, iterates to remaining point, until not having
There are data point needs to be rejected, and this method effect in the case where low signal-to-noise ratio or low signal background ratio is poor；In addition, very
In more application processes, the spectral line of spectrum is rendered as various shapes.The method of simple fitting of a polynomial is not able to satisfy special spectrum
Baseline fitting requirement.The automatic Iterative method of moving average is in the case where multimodal or overlap peak, and effect not smooth enough to baseline correction
Fruit is not fine；ALS method and air PLS method are all based on the smooth baseline correction of Whittaker, baseline correction effect compared with
It is good, but parameter optimization is required, and adjusting optimized parameter is the process for needing experience；And traditional Rubber-
If spectrum is divided into stem portion by band method, the minimum point of each part is considered as the position of baseline, then by with linearly
Interpolation or spline interpolation combine these points to complete the baseline of estimation, and this method needs artificial adjusting segmentation width, corrects result
Quality to a certain extent depend on operator experience.
And existing patent and paper have and are much referred to different spectral correction methods.For example, Taiwan super-micro optical
Hong Jianxiang of company et al. (108318137 A of CN) establishes stray light data bank using the response of more narrow frequency spectrographic pulse, with
It is established afterwards using the data in data bank and repairs the spectrum that formula removes amendment spectroradio standard light, to obtain measure spectrum number
According to generate correction coefficient formula.Although this method effectively, needs to establish database, therefore heavy workload and need by
It is gradually perfect, it is not used to online real-time measurement and application range receives the limitation of database.Hangzhou distant place photoelectric information is public
The advanced spectroscope function that Pan Jiangen et al. (1021559971 A of CN) calculates spectral instrument using photoradiometer is taken charge of, the letter is utilized
Number can be modified spectral instrument measurement result.This method can effectively improve that precision, speed are fast, design is simple and user
Just, but it is only limitted to the calibration spectrum measurement of the free wave band of shortwave, application range limitation is big.The auspicious environmental science and technology company of Shenzhen Kapp
Ou Yangbin et al. (107607144 A of CN) obtains functional relation using aiding sensors and carries out school to the baseline drift of sensor
Just.This method is only capable of the influence of baseline drift in removal measurement result and needs additional aiding sensors during processing
It assists.Cyrillic Bo Dibu gas of university, Bordeaux et al. (107076664 A of CN) in multiple absorption spectrums using absorbing
The interval spectral coverage for the position that amount is zero determines baseline correction curve.This method uses signal self-correcting, is only capable of preferably removing width
Baseline problem caused by the apparent fluctuation of value.Therefore in the prior art, there are no can handle spectral signal simultaneously in real time online
Middle noise, baseline drift and detector or signal path cause the bearing calibration of nonlinear influence.
Disclosure
(1) technical problems to be solved
Based on above-mentioned technical problem, the disclosure provides a kind of receptance function bearing calibration of online spectrum measurement instruments, real
Now to the correction of the receptance function of spectrum measurement instruments, spectrum can not be handled simultaneously in real time online in the prior art to alleviate
The technical issues of noise, baseline drift and detector or signal path cause nonlinear influence in signal.
(2) technical solution
The disclosure provides a kind of receptance function bearing calibration of online spectrum measurement instruments, by adding in measured matter
Standard substance after screening realizes the response letter of spectrum measurement instruments using the correction linear to measure spectrum Plays substance
Several correction, comprising: step A: using photodetector acquisition absorption signal and background signal, and make in two groups of spectral signals
It include one section of signal segment without laser emitting；Wherein, in the absorption signal comprising sample to be tested spectral signal and with this
The spectral signal for the standard substance that sample to be tested does not chemically react；Step B: using in the absorption signal and the background signal
Signal segment and photodetector without laser emitting eliminate environmental background light without the output signal under laser input state
Fluctuate the influence to spectroscopic data；Step C: being respectively weighted fitting to step B treated absorption signal and background signal,
And the optical stripe in absorption signal is eliminated using the regression criterion in background signal；Step D: the standard substance acquired is utilized
Spectral signal is matched with the spectral information in database, and to step C, treated that absorption signal is corrected, and is corrected
The spectral signal of sample to be tested afterwards.
In some embodiments of the present disclosure, the step A includes: step A1: known concentration being added in sample to be tested
Standard substance；Wherein, the concentration of the standard substance, so that its feature peak amplitude is less than or equal to the amplitude of sample to be tested
It is quasi-；Step A2: it is mixed with the spectral signal of the sample to be tested of standard substance using photodetector acquisition, and utilizes identical system
Nitrogen data of the parameter acquisition without absorption of uniting are as background signal；It wherein, include laser tuning in two groups of spectral signals
Current threshold output signal below；And step A3: multicycle signal is used to the two groups of spectral signals acquired in step A2
Average method is smoothed.
In some embodiments of the present disclosure, in the step A1: the standard substance contains in tested absorption bands
Characteristic peak that is independent and not interfered mutually with measured matter, and have detailed spectrum parameter in the database.
In some embodiments of the present disclosure, the database includes: HITRAN database, PNNL database or NIST number
According to library.
In some embodiments of the present disclosure, the step B includes: step B1: determining two groups of spectrum letters on a timeline
Number start position, and be aligned two groups of spectral signals on a timeline on the basis of the start position of absorption signal；Step B2: it adopts
Collect output signal of the photodetector under no laser input state and calculates signal averaging；Step B3: background signal is calculated
With absorption signal laser tuning current threshold output signal below average value；Step B4: step B3 is resulting flat
Mean value makes the difference with the resulting average value of step B2 respectively, and using respective difference as offset, moves background signal and suction respectively
Position of the collection of letters number on y direction, makes background signal and absorption signal and photodetector under no laser input state
The aligned in position of output signal.
In some embodiments of the present disclosure, the step C includes: step C1: to background signal entire data segment into
Row weight is 1 fitting；Step C2: the fitting that weight is 0, nothing are carried out to the data segment in absorption signal containing characteristic information
The data segment of characteristic information carries out the fitting that weight is 1；Step C3: the regression criterion after step C1 and step C2 fitting is carried out
Comparison matching, so that mutually be overlapped with the regression criterion of background signal without the corresponding regression criterion of characteristic information data section in absorption signal
It closes；Step C4: the regression criterion of absorption signal and background signal is made the difference, and removes the influence of optical stripe in absorption signal.
In some embodiments of the present disclosure, in the step D: utilize following updating formula, to step C treated inhale
The collection of letters number carries out a baseline correction:
I′0(v)=WI0(v)+Q(v)
Wherein, I '0(v) indicate that the spectrum after baseline correction is linear, second order polynomial fitting are as follows: I '0(v)=c+bv+
av2, I0(v) indicate that the spectrum before baseline correction is linear, second order polynomial fitting are as follows: I0(v)=m+nv+tv2,For
The gamma correction factor, For deviation constant matrix.
In some embodiments of the present disclosure, the step D includes: step D1: setting W initial value is the initial square of 1, Q (v)
Battle array value is all 0；Step D2: the spectral signal and number of the standard substance that Beer-Lambert law acquires photodetector are utilized
Linear matching is carried out according to the spectral signal in library, until the variance yields of the two is minimum；Step D3: when acquisition variance yields minimum
The gamma correction factor W and deviation constant matrix Q (v).Step D4: minimum two is carried out to deviation constant matrix Q (v) using following formula
Multiplication iteration acquires deviation constant Q:
N=Σ (Qi-Qj)2
Wherein, when N obtains minimum value, the spectral signal for the standard substance that deviation constant Q acquires photodetector
And the linear deviation of the spectral signal in database is minimum；
Step D5: gamma correction factor W and deviation constant Q are brought into updating formula, and using updating formula to light
Spectrum signal is corrected.
In some embodiments of the present disclosure, the step D2 includes: step D2a: obtaining absorption signal Plays substance
Spectrum it is linear；Step D2b: the database shape difference linear with the spectrum of absorption signal Plays substance is calculated；Step D2c:
The variance of the shape difference obtained in step D2b is calculated, and judges whether variance yields is minimum at this time, if variance is minimum, carries out step
Otherwise D3 carries out step D2d；Step D2d: updating the value of gamma correction factor W, is got the bid using updating formula to absorption signal
The linear carry out baseline correction of the spectrum of quasi- substance, and return step D2b.
In some embodiments of the present disclosure, in which: the receptance function bearing calibration institute of the online spectrum measurement instruments
Applicable spectrum includes: absorption spectrum, emission spectrum, molecular spectrum, atom spectrum；The response of the online spectrum measurement instruments
The project that function correction method can correct includes: the non-linear, dry of baseline drift, channel existing for spectral signal and detector
It disturbs and noise.
(3) beneficial effect
It can be seen from the above technical proposal that the receptance function bearing calibration for the online spectrum measurement instruments that the disclosure provides
It one of has the advantages that or in which a part:
(1) disclosure realizes the correction to spectrum measurement instruments receptance function by the correction linear to measure spectrum,
This bearing calibration is not necessarily to extras and error is small, and linear correction can be carried out to special shape spectrum, adaptable；
(2) the receptance function bearing calibration for the online spectrum measurement instruments that the disclosure provides directly utilizes standard substance in number
Linear according to the linear standard substance in the compounding substances spectral signal that obtains in experiment measurement of library Plays spectrum carries out pair
Than amendment, to realize the school using the linear multidimensional information of the spectrum of small-molecule substance to compounding substances signal integral linear
Just, to remove the non-linear and baseline drift in signal, the total regression to the entire signal path including spectrum baseline is realized
It corrects, reaches correction while realizing multiple measurement parameters using correction course of a standard substance；
(3) the receptance function bearing calibration for the online spectrum measurement instruments that the disclosure provides is simple and convenient for operation, effective,
There is wide applicability in the fields of measurement such as molecular spectrum measurement and content of material.
Detailed description of the invention
Fig. 1 is the step flow chart of the receptance function bearing calibration for the online spectrum measurement instruments that the disclosure provides.
Fig. 2 is spectroscopic data correcting process in the receptance function bearing calibration for the online spectrum measurement instruments that the disclosure provides
Figure.
Fig. 3 is the direct absorption spectrum measuring system schematic diagram based on TLAS technology.
Fig. 4 is the initial data schematic diagram of photodetector acquisition.
Fig. 5 is the weighted fitting result and residual error comparison diagram of background signal and absorption signal.
Fig. 6 be using the disclosure provide line spectrum measurement instruments receptance function bearing calibration correction after baseline with
Baseline comparison diagram before correction.
Fig. 7 is the comparison diagram and correcting value before and after spectrum correction.
Specific embodiment
The receptance function bearing calibration for the online spectrum measurement instruments that the disclosure provides directly utilizes standard substance in data
Plays spectrum linear the linear of standard substance in the compounding substances spectral signal obtained in experiment measurement in library compares
Amendment is not necessarily to extras and error is small to realize the correction to spectrum measurement instruments receptance function, adaptable.
For the purposes, technical schemes and advantages of the disclosure are more clearly understood, below in conjunction with specific embodiment, and reference
The disclosure is further described in attached drawing.
The disclosure provides a kind of receptance function bearing calibration of online spectrum measurement instruments, by adding in measured matter
Standard substance after screening realizes the response letter of spectrum measurement instruments using the correction linear to measure spectrum Plays substance
Several corrections, as shown in Figure 1 to Figure 2, comprising:
Step A: using photodetector acquisition absorption signal and background signal, and make include in two groups of spectral signals
One section of signal segment without laser emitting；
Wherein, in the absorption signal comprising sample to be tested spectral signal and do not chemically reacted with the sample to be tested
Standard substance spectral signal；
Step B: using in the absorption signal and the background signal without laser emitting signal segment and photodetector without
Output signal under laser input state eliminates influence of the fluctuation of environmental background light to spectroscopic data；
Step C: fitting is respectively weighted to step B treated absorption signal and background signal, and is believed using background
Regression criterion in number eliminates the optical stripe in absorption signal；
Step D: it is matched using the spectral signal of the standard substance of acquisition with the spectral information in database, to step
Treated that absorption signal is corrected by C, the spectral signal of the sample to be tested after being corrected.
In some embodiments of the present disclosure, step A includes:
Step A1: the standard substance of known concentration is added in sample to be tested；
Wherein, the concentration of the standard substance, so that its feature peak amplitude is less than or equal to subject to the amplitude of sample to be tested；
Step A2: it is mixed with the spectral signal of the sample to be tested of standard substance using photodetector acquisition, and utilizes phase
Same nitrogen data of the system parameter acquisition without absorption are as background signal；
It wherein, include laser tuning current threshold output signal below in two groups of spectral signals；
Step A3: the two groups of spectral signals acquired in step A2 are smoothly located using the method for multicycle signal averaging
Reason.
In some embodiments of the present disclosure, step A1: standard substance in tested absorption bands containing it is independent and not with
The characteristic peak that measured matter mutually interferes, and have detailed spectrum parameter in the database.
In some embodiments of the present disclosure, database includes: HITRAN database, PNNL database or NIST database.
In some embodiments of the present disclosure, as shown in " signal alignment " part in Fig. 2, step B includes:
Step B1: the start position (s and r) of two groups of spectral signals is determined, and on a timeline with the starting point of absorption signal
Two groups of spectral signals are aligned on the basis of the r of position on a timeline；
Step B2: in acquisition longer period of time, output signal and meter of the photodetector under no laser input state
Signal averaging is calculated, d is denoted as；
Step B3: background signal and absorption signal being averaged in laser tuning current threshold output signal below are calculated
Value, is denoted as f and e respectively；
Step B4: the resulting average value f and e of step B3 is made the difference with the resulting average value d of step B2 respectively, and with respective
Difference be offset, move the position of background signal and absorption signal on y direction respectively, make background signal and absorption
The aligned in position of the output signal of signal and photodetector under no laser input state.
In some embodiments of the present disclosure, step C includes:
Step C1: the fitting that weight is 1 is carried out in entire data segment to background signal；
Step C2: carrying out the fitting that weight is 0 to the data segment in absorption signal containing characteristic information, no characteristic information
Data segment carries out the fitting that weight is 1；
Step C3: the regression criterion after step C1 and step C2 fitting is compared into matching, so that nothing in absorption signal
The regression criterion of the corresponding regression criterion of characteristic information data section and background signal coincides；
Step C4: the regression criterion of absorption signal and background signal is made the difference, and removes the shadow of optical stripe in absorption signal
It rings.
In some embodiments of the present disclosure, step D: following updating formula is utilized, to step C treated absorption signal
Carry out baseline correction:
I′0(v)=WI0(v)+Q(v)
Wherein, I '0(v) indicate that the spectrum after baseline correction is linear, second order polynomial fitting are as follows: I '0(v)=c+bv+
av2, It(v) indicate that the spectrum before baseline correction is linear, second order polynomial fitting are as follows: I0(v)=m+nv+tv2,For
The gamma correction factor, For deviation constant matrix.
In some embodiments of the present disclosure, step D includes:
Step D1: setting W initial value is that 1, Q (v) initial matrix value is all 0；
Step D2: the spectral signal and data of the standard substance that Beer-Lambert law acquires photodetector are utilized
Spectral signal in library carries out linear matching, until the variance yields of the two is minimum；
Step D3: the gamma correction factor W and deviation constant matrix Q (v) when variance yields minimum are obtained.
Step D4: least square method iteration is carried out to deviation constant matrix Q (v) using following formula and acquires deviation constant Q:
N=Σ (Qi-Qj)2
Wherein, when N obtains minimum value, the spectral signal for the standard substance that deviation constant Q acquires photodetector
And the linear deviation of the spectral signal in database is minimum；
Step D5: gamma correction factor W and deviation constant Q are brought into updating formula, and using updating formula to light
Spectrum signal is corrected.
In some embodiments of the present disclosure, as shown in " extension Internal standard correction methods " part in Fig. 2, step D2 includes:
Step D2a: the spectrum for obtaining absorption signal Plays substance is linear；
Step D2b: the database shape difference linear with the spectrum of absorption signal Plays substance is calculated；
Step D2c: calculating the variance of the shape difference obtained in step D2b, and judges whether variance yields is minimum at this time, if side
It is poor minimum, step D3 is carried out, step D2d is otherwise carried out；
Step D2d: the value of gamma correction factor W is updated, using updating formula to the light of absorption signal Plays substance
Spectral line shape carries out baseline correction, and return step D2b.
In some embodiments of the present disclosure, in which: the receptance function correction for the online spectrum measurement instruments that the disclosure provides
The spectrum that method is applicable in includes: absorption spectrum, emission spectrum, molecular spectrum, atom spectrum；The disclosure provide in line spectrum
The project that the receptance function bearing calibration of measuring instrument can correct includes: baseline drift, channel existing for spectral signal and inspection
Survey the non-linear of device, interference and noise.
Below with a specific embodiment, the receptance function bearing calibration for the online spectrum measurement instruments that the disclosure provides is verified
Validity:
As shown in figure 3, being based on tunable laser absorption spectrum (Tunable Laser
AbsorptionSpectroscopy, TLAS) technology constructs spectral measurement system, and SG is signal generator in Fig. 3, and LDC is sharp
CD-ROM driver, DFB-ICL are distributed Feedback interband cascade lasers, and PD is photodetector, and DAQ is data collecting card, PC
Acetonitrile (CH is selected according to the existence conditions that experiment has for computer3CN it) is used as measured matter, and is swashed in Binding experiment system
The scanning range of light device determines that detection wave band is 3037.5-3043cm-1。
Step 1: screening criteria substance: the requirement of combination technology scheme, the standard substance finally chosen are methane.
Step 2: preparing sample to be tested: for the validity for verifying bearing calibration, being mixed with the calibrating gas of known concentration
Close the proportion of gaseous sample.Experiment utilizes 10.2ppm methane (CH4) as a kind of standard substance addition 160.8ppm CH3CN sample
In product.
Step 3: acquisition data: by proportioned mixed gas sample under conditions of normal temperature and pressure (23.6 DEG C, 1atm)
It is passed through in the gas cell of 10.24455m light path, and absorbs information containing mixed gas using what capture card acquired 20 periods
Signal It.Then, the pure nitrogen gas for being passed through enough time blows mixed gas, and is believed using identical sampling parameter acquisition nitrogen
Number as background signal such as.Collected original signal is as shown in figure 4, without absorption peak be background signal I0, contain absorption peak
Be mixed gas absorption signal It。
Step 4: noise reduction process: two groups of signals will be obtained in step 3 and carry out monocycle signal extraction respectively, then to 20
The cycle data progress period is average, to realize signal smoothing.
Step 5: signal alignment: as shown in " signal alignment " part in Fig. 2, determining that two groups of signals are each on a timeline first
From initial point position (s and r), and be aligned two signals on a timeline on the basis of the start position r of absorption signal.Then,
Background signal I is obtained respectively0With absorption signal ItWhen below laser tuning current threshold (in 0.002 second), output signal
Average value be respectively f and e.Then, the photodetector in one hour is acquired under the natural working condition that no laser inputs
Output signal and calculate signal averaging be d.Later, background signal I is utilized0With absorption signal ItResulting average value f, e points
The average value d not obtained with photodetector makes the difference, and using respective difference as offset, moves background signal I respectively0And suction
Collection of letters ItPosition on y direction, so that background signal I0With absorption signal ItIt is inputted with photodetector in no laser
Natural working condition under output signal aligned in position, to eliminate the influence of background light beats.
Step 6: weighted fitting: to by step 5 treated background signal I0With absorption signal ItIt is respectively weighted quasi-
It closes.To background signal I0After entire data segment carries out the fitting that weight is 1, the available complete light in regression criterion
Striped and noise are learned, as shown in Fig. 5 (a), due to absorption signal ItIt is middle to there is absorption information, therefore using sectionally weighting fitting
Method will be set as 1 without the position weight for absorbing information, other positions 0, to curve matching be come out, such as Fig. 5 (b) institute
Show, then, by background signal I0With absorption signal ItRegression criterion compare, it can be found that background signal from Fig. 5 (c)
I0With absorption signal ItOptical stripe can be well matched with, so as to utilize background signal I0In regression criterion disappear
Except absorption signal ItMiddle optical stripe.
Step 7: baseline correction: will be by the absorption signal I of step 6 removal optical stripetCarry out baseline correction, processing stream
Journey sets 1 for gamma correction factor W first, utilizes the mixed of experiment acquisition as shown in " extension Internal standard correction methods " part in Fig. 2
Close gas absorption signal Plays substance C H4Absorption lineshape matched with the linear of CH4 in HITRAN database.Correction is calculated
Method is iterated by least square method, when the value of W is adjusted to the two variance minimum, acquires the final gamma correction factor
W and deviation constant matrix Q (v), since Q (v) is an array, it contain between two curves with all independent variable v phases
Corresponding offset.Therefore, the smallest offset of variance that can make two curves is obtained by the iteration of least square method
Constant Q, i.e., when quadratic sum N is minimum, corresponding QiThe deviation of two curves can be made minimum, Q=Qi, after corrected
For baseline as shown in fig. 6, top curve is before correcting, lower curve is after correcting.Finally, utilizing the resulting gamma correction factor
W and deviation constant Q combines the signal expression before correction, the signal after correction can be obtained, as shown in Figure 7 the mixing of correction front and back
The absorption spectrum of gas has obvious gap, and can intuitively find out the value of correcting value.
In conclusion the receptance function bearing calibration for the online spectrum measurement instruments that the disclosure provides can remove in signal
Non-linear and baseline drift, realization the total regression of the entire signal path including spectrum baseline is corrected, reach and utilize one
Correction course of a standard substance corrects while realizing multiple measurement parameters, and final realize responds letter to spectrum measurement instruments
Several corrections.
It should also be noted that, the direction term mentioned in embodiment, for example, "upper", "lower", "front", "rear", " left side ",
" right side " etc. is only the direction with reference to attached drawing, not is used to limit the protection scope of the disclosure.Through attached drawing, identical element by
Same or similar appended drawing reference indicates.When may cause understanding of this disclosure and cause to obscure, conventional structure will be omitted
Or construction.
And the shape and size of each component do not reflect actual size and ratio in figure, and only illustrate the embodiment of the present disclosure
Content.In addition, in the claims, any reference symbol between parentheses should not be configured to the limit to claim
System.
Similarly, it should be understood that in order to simplify the disclosure and help to understand one or more of each open aspect,
Above in the description of the exemplary embodiment of the disclosure, each feature of the disclosure is grouped together into single implementation sometimes
In example, figure or descriptions thereof.However, the disclosed method should not be interpreted as reflecting the following intention: i.e. required to protect
The disclosure of shield requires features more more than feature expressly recited in each claim.More precisely, such as front
Claims reflect as, open aspect is all features less than single embodiment disclosed above.Therefore,
Thus the claims for following specific embodiment are expressly incorporated in the specific embodiment, wherein each claim itself
All as the separate embodiments of the disclosure.
Particular embodiments described above has carried out further in detail the purpose of the disclosure, technical scheme and beneficial effects
Describe in detail it is bright, it is all it should be understood that be not limited to the disclosure the foregoing is merely the specific embodiment of the disclosure
Within the spirit and principle of the disclosure, any modification, equivalent substitution, improvement and etc. done should be included in the guarantor of the disclosure
Within the scope of shield.
Claims (10)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811381410.XA CN109520941B (en) | 2018-11-20 | 2018-11-20 | Response function correction method for on-line spectrometer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811381410.XA CN109520941B (en) | 2018-11-20 | 2018-11-20 | Response function correction method for on-line spectrometer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109520941A true CN109520941A (en) | 2019-03-26 |
| CN109520941B CN109520941B (en) | 2021-02-09 |
Family
ID=65776744
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811381410.XA Expired - Fee Related CN109520941B (en) | 2018-11-20 | 2018-11-20 | Response function correction method for on-line spectrometer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109520941B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110108694A (en) * | 2019-05-17 | 2019-08-09 | 广西科技大学 | Method for Improving Measurement Accuracy of Raman Spectrometer Wavenumber Migration |
| CN114384028A (en) * | 2021-12-14 | 2022-04-22 | 安徽皖仪科技股份有限公司 | Peak drift correction method for continuous flow analyzer |
| CN117270024A (en) * | 2023-11-20 | 2023-12-22 | 合肥综合性国家科学中心人工智能研究院(安徽省人工智能实验室) | Method and device for correcting energy spectrum response function, computer equipment and storage medium |
| CN118425098A (en) * | 2024-07-04 | 2024-08-02 | 鞍山天汇科技有限公司 | Distributed laser methane detection method and system |
| CN118566173A (en) * | 2024-08-05 | 2024-08-30 | 华夏天信传感科技(大连)有限公司 | Noise suppression effect inspection method and device |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0727703A (en) * | 1993-07-08 | 1995-01-31 | Hitachi Ltd | Quantitative analysis method for multi-component substances |
| CN1149334A (en) * | 1994-05-27 | 1997-05-07 | 伊斯特曼化学公司 | Raman spectrometry apparatus and method |
| CN101101338A (en) * | 2007-07-10 | 2008-01-09 | 廊坊开发区中油油田科技工贸有限责任公司 | Exploration method and system for oil and gas |
| CN102155991A (en) * | 2011-03-18 | 2011-08-17 | 杭州远方光电信息股份有限公司 | Advanced spectrum correcting method |
| CN205821460U (en) * | 2016-08-03 | 2016-12-21 | 苏州热工研究院有限公司 | A kind of portable reference electrode calibration system being applicable to run mediator depletion |
| CN106769906A (en) * | 2016-12-14 | 2017-05-31 | 重庆川仪自动化股份有限公司 | Spectrometer data wander compensation method |
-
2018
- 2018-11-20 CN CN201811381410.XA patent/CN109520941B/en not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0727703A (en) * | 1993-07-08 | 1995-01-31 | Hitachi Ltd | Quantitative analysis method for multi-component substances |
| CN1149334A (en) * | 1994-05-27 | 1997-05-07 | 伊斯特曼化学公司 | Raman spectrometry apparatus and method |
| CN101101338A (en) * | 2007-07-10 | 2008-01-09 | 廊坊开发区中油油田科技工贸有限责任公司 | Exploration method and system for oil and gas |
| CN102155991A (en) * | 2011-03-18 | 2011-08-17 | 杭州远方光电信息股份有限公司 | Advanced spectrum correcting method |
| CN205821460U (en) * | 2016-08-03 | 2016-12-21 | 苏州热工研究院有限公司 | A kind of portable reference electrode calibration system being applicable to run mediator depletion |
| CN106769906A (en) * | 2016-12-14 | 2017-05-31 | 重庆川仪自动化股份有限公司 | Spectrometer data wander compensation method |
Non-Patent Citations (2)
| Title |
|---|
| NAOKI SAITO等: "Extended internal standard method for quantitative 1H NMR assisted by chromatography (EIC) for analyte overlapping impurity on 1H NMR spectra", 《TALANTA》 * |
| 董学锋等: "基于高斯函数卷积的色散型拉曼光谱仪温度校正", 《分析化学研究报告》 * |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110108694A (en) * | 2019-05-17 | 2019-08-09 | 广西科技大学 | Method for Improving Measurement Accuracy of Raman Spectrometer Wavenumber Migration |
| CN110108694B (en) * | 2019-05-17 | 2021-11-02 | 广西科技大学 | Raman Spectrometer Wavenumber Shift Measurement Accuracy Improvement Method |
| CN114384028A (en) * | 2021-12-14 | 2022-04-22 | 安徽皖仪科技股份有限公司 | Peak drift correction method for continuous flow analyzer |
| CN114384028B (en) * | 2021-12-14 | 2023-10-24 | 安徽皖仪科技股份有限公司 | Peak drift correction method for continuous flow analyzer |
| CN117270024A (en) * | 2023-11-20 | 2023-12-22 | 合肥综合性国家科学中心人工智能研究院(安徽省人工智能实验室) | Method and device for correcting energy spectrum response function, computer equipment and storage medium |
| CN117270024B (en) * | 2023-11-20 | 2024-02-20 | 合肥综合性国家科学中心人工智能研究院(安徽省人工智能实验室) | Method and device for correcting energy spectrum response function, computer equipment and storage medium |
| CN118425098A (en) * | 2024-07-04 | 2024-08-02 | 鞍山天汇科技有限公司 | Distributed laser methane detection method and system |
| CN118425098B (en) * | 2024-07-04 | 2024-08-23 | 鞍山天汇科技有限公司 | Distributed laser methane detection method and system |
| CN118566173A (en) * | 2024-08-05 | 2024-08-30 | 华夏天信传感科技(大连)有限公司 | Noise suppression effect inspection method and device |
| CN118566173B (en) * | 2024-08-05 | 2024-10-01 | 华夏天信传感科技(大连)有限公司 | Noise suppression effect inspection method and device |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109520941B (en) | 2021-02-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109520941A (en) | The receptance function bearing calibration of online spectrum measurement instruments | |
| CN101105446B (en) | Differential Absorption Spectroscopy Air Quality Detection System | |
| US10113999B2 (en) | Method and a device for detecting a substance | |
| CN109425583A (en) | Spectral analysis device, spectroscopic analysis methods and storage medium | |
| US10557792B2 (en) | Spectral modeling for complex absorption spectrum interpretation | |
| CN103940767B (en) | Gas concentration inversion method based on multi-manifold learning | |
| CN107219190A (en) | A kind of non-calibrating gas parameter measuring method based on triangular modulation | |
| CN108918446A (en) | A kind of super low concentration sulfur dioxide ultraviolet difference feature extraction algorithm | |
| Tang et al. | On-line multi-component alkane mixture quantitative analysis using Fourier transform infrared spectrometer | |
| CN116559110A (en) | An Adaptive Near Infrared Spectral Transformation Method Based on Correlation and Gaussian Curve Fitting | |
| Takahama et al. | Analysis of functional groups in atmospheric aerosols by infrared spectroscopy: sparse methods for statistical selection of relevant absorption bands | |
| CN106126879B (en) | A kind of soil near-infrared spectrum analysis prediction technique based on rarefaction representation technology | |
| CN104316511A (en) | Correcting method for spectral line interference in atomic emission spectrometric analysis of spark source | |
| CN106018331A (en) | Stability evaluation method of multi-channel spectrum system and pretreatment optimization method | |
| Chan et al. | Dispersive infrared spectroscopy measurements of atmospheric CO2 using a Fabry–Pérot interferometer sensor | |
| CN110208251A (en) | Plasma emission spectroscopy interference correction method | |
| Griffith et al. | Intercomparison of NDSC ground-based solar FTIR measurements of atmospheric gases at Lauder, New Zealand | |
| US11646186B2 (en) | System and method for optimizing peak shapes | |
| CN115541563A (en) | Element quantitative analysis method based on laser-induced breakdown spectroscopy | |
| CN119354911A (en) | Rubber additive composition analysis and quality control method and system based on infrared spectroscopy | |
| CN114791420A (en) | Calibration method of near-infrared spectrometer | |
| CN106872397A (en) | A kind of method based on existing calibration model quick detection agricultural product chemical constituent | |
| CN111579526B (en) | A method for characterizing near-infrared instrument variance and correction | |
| CN114387260A (en) | Interpolation reconstruction method based on pepper spectrum characteristic peak | |
| CN109521002A (en) | A kind of fuel characteristic measurement method of solid fuel particle stream |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210209 |