CN105447455B - Laser Modulation absorbs gas sensor signal preprocess method - Google Patents

Laser Modulation absorbs gas sensor signal preprocess method Download PDF

Info

Publication number
CN105447455B
CN105447455B CN201510788547.7A CN201510788547A CN105447455B CN 105447455 B CN105447455 B CN 105447455B CN 201510788547 A CN201510788547 A CN 201510788547A CN 105447455 B CN105447455 B CN 105447455B
Authority
CN
China
Prior art keywords
signal
array
concentration signal
high concentration
low concentration
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.)
Active
Application number
CN201510788547.7A
Other languages
Chinese (zh)
Other versions
CN105447455A (en
Inventor
贾渠
高翔
常洋
李彦林
宁召科
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.)
Beijing Aerospace Tianhong Intelligent Equipment Technology Co ltd
Original Assignee
Beijing Aerospace Yilian Science and Technology Development Co Ltd
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 Beijing Aerospace Yilian Science and Technology Development Co Ltd filed Critical Beijing Aerospace Yilian Science and Technology Development Co Ltd
Priority to CN201510788547.7A priority Critical patent/CN105447455B/en
Publication of CN105447455A publication Critical patent/CN105447455A/en
Application granted granted Critical
Publication of CN105447455B publication Critical patent/CN105447455B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/02Preprocessing
    • G06F2218/04Denoising
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2218/00Aspects of pattern recognition specially adapted for signal processing
    • G06F2218/12Classification; Matching

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The invention discloses a kind of Laser Modulations to absorb gas sensor signal preprocess method, and the sensor analog signals of acquisition are converted to digital signal, digital signal is preserved in a manner of one-dimensional fixed long number group;The array of preservation is carried out blocking optical path signal according to the threshold value of setting and is deleted, preliminary concentration is calculated after will be less than the number removal of threshold value, distinguishes high concentration signal and low concentration signal;Low concentration signal noise removal process is carried out when for low concentration signal;High concentration signal noise removal process is carried out when for high concentration signal;Signal after noise removal process is sent into subsequent processing steps, pretreatment is completed.By this method, for different interference types, respective preprocess method is proposed respectively, is realized the pretreatment after low concentration, high concentration signal sensing, is measured caused by signal interference is greatly lowered inaccurate.

Description

Laser Modulation absorbs gas sensor signal preprocess method
Technical field
The invention belongs to sensor signal treatment technology, more particularly to it is pre- that a kind of Laser Modulation absorbs gas sensor signal Processing method is before one kind carrying out gas concentration calculating for TDLAS test methods, to carry out Signal Pretreatment algorithm.
Background technology
TDLAS technologies are gas sensing techniques of new generation, and gas sensor is under low concentration, high concentration made of the technology Various signal interferences are easily received, influence evaluation when meter gaseous, and then influence gas sensing test concentrations.
Invention content
The purpose of the present invention is to propose to a kind of Laser Modulations to absorb gas sensor signal preprocess method, be it is a kind of for Before TDLAS test methods carry out gas concentration calculating, Signal Pretreatment algorithm is carried out, according to multiple test, is found out under low concentration Noise behavior is formed automatic threshold, is removed jump signal caused by noise according to threshold value using first order derivative;In high concentration Under, standard curve is found out out, related algorithm is carried out by standard curve and acquisition signal curve, it, will according to automatic threshold method It is removed including signal interference caused by curve mutation, the various noises such as noise is big.
To achieve the goals above, the technical scheme is that:
A kind of Laser Modulation absorption gas sensor signal preprocess method, the sensor analog signals of acquisition are converted to Digital signal, method are:
Digital signal is preserved in a manner of one-dimensional fixed long number group;
The array of preservation is carried out blocking optical path signal according to the threshold value of setting and is deleted, after the number removal that will be less than threshold value Preliminary concentration is calculated, high concentration signal and low concentration signal are distinguished;
Low concentration signal noise removal process is carried out when for low concentration signal;High concentration is carried out when for high concentration signal Signal noise removal is handled;
Signal after noise removal process is sent into subsequent processing steps, pretreatment is completed.
Scheme is further:The differentiation high concentration signal and low concentration signal are:Preliminary concentration is mono- higher than 1500ppmV Position is high concentration signal, is then low concentration signal less than this value.
Scheme is further:The long array of fixation is the array of continuous 100 points, and the analog signal is converted to number The A/D conversion digits of signal are 12, and the threshold value is 10.
Scheme is further:The method of the low concentration signal noise removal process is:
It is first determined whether there are spikes, if spike is not present, continue subsequent processing steps without processing;If in the presence of Spike, the then value that original signal is acquired to the maximum previous point in multiple positions of array sort substitute the point value;It then looks for The array maximum point position, if subsequent processing steps are continued in maximum point position in the centre position of array;If maximum point position Not in centre position, then the array exists abnormal, deletes the array;
The method of the high concentration signal noise removal process is:
High concentration signal and standard curve are subjected to the curve after related algorithm obtains related algorithm first, if related algorithm There are a minor peaks in curve afterwards, and peak value is in the centre position of curve;Or there is multiple peak value, and multiple peak value in curve Interpeak be in the centre position of curve, then normally to continue subsequent processing steps;Otherwise the array is deleted;
The multiple point is 5% to 10% always to count;
The standard curve is the curve made previously according to ideal data;
The method for judging whether spike is:If fixed long array is the array of continuous 100 points, to continuous 100 low concentration signals obtain the first order derivative array of one-dimension array, the length of array is 99 by the way of first order derivative A, which takes maximum 10 points after its numerical ordering, then calculate remaining 89 points average value and maximum point 10 A average value;Compare the two average value size, when maximum 10 points cell means are higher than remaining 8 times of 89 average value More than, then there is " spike ", spike is otherwise not present.
The beneficial effects of the invention are as follows:By this method, for different interference types, respective pre- place is proposed respectively Reason method realizes the pretreatment after low concentration, high concentration signal sensing, is measured caused by signal interference is greatly lowered inaccurate Really.
Invention is described in detail with reference to the accompanying drawings and examples.
Description of the drawings
Fig. 1 is low concentration abnormal signal curve;
Fig. 2 is that treated signal curve;
Fig. 3 is high concentration signal contrast standard curve used in connection with;
Fig. 4 is a kind of abnormal signal under high concentration;
Fig. 5 be to Fig. 4 signals after relevant treatment result;
Fig. 6 is another abnormal signal under high concentration;
Fig. 7 be to Fig. 6 signals after relevant treatment result.
Specific implementation mode
A kind of Laser Modulation absorption gas sensor signal preprocess method, the sensor analog signals of acquisition are converted to Digital signal, method are:
Digital signal is preserved in a manner of one-dimensional fixed long number group;
The array of preservation is carried out blocking optical path signal according to the threshold value of setting and is deleted, after the number removal that will be less than threshold value Preliminary concentration is calculated, high concentration signal and low concentration signal are distinguished;
Low concentration signal noise removal process is carried out when for low concentration signal;High concentration is carried out when for high concentration signal Signal noise removal is handled;
Signal after noise removal process is sent into subsequent processing steps, pretreatment is completed.
In embodiment:The differentiation high concentration signal and low concentration signal are:Preliminary concentration is higher than 1500ppmV units High concentration signal is then low concentration signal less than this value.
Scheme is further:The long array of fixation is the array of continuous 100 points, and the analog signal is converted to number The A/D conversion digits of signal are 12, and the threshold value is 10.
In embodiment:The method of the low concentration signal noise removal process is:
It is first determined whether there are spikes, if spike is not present, continue subsequent processing steps without processing;If in the presence of Spike, the then value that original signal is acquired to the maximum previous point in multiple positions of array sort substitute the point value;It then looks for The array maximum point position, if subsequent processing steps are continued in maximum point position in the centre position of array;If maximum point position Not in centre position, then the array exists abnormal, deletes the array;
The method of the high concentration signal noise removal process is:
The curve that high concentration signal is formed is compared with standard curve first, related algorithm is obtained using related algorithm Curve afterwards, if a minor peaks occurs in the curve after related algorithm, and peak value is in the centre position of curve;Or curve occurs Multiple peak value, the and repeatedly interpeak of peak value is in the centre position of curve, then normally to continue subsequent processing steps;Otherwise Delete the array;
The multiple point is 5% to 10% always to count, when the long array of fixation is the array of continuous 100 points in the present embodiment When, the multiple point is selected as 10 points;
The standard curve is the curve made previously according to ideal data;
The method for judging whether spike is:If fixed long array is the array of continuous 100 points, to continuous 100 low concentration signals obtain the first order derivative array of one-dimension array, the length of array is 99 by the way of first order derivative A, which takes maximum 10 points after its numerical ordering, then calculate remaining 89 points average value and maximum point 10 A average value;Compare the two average value size, when maximum 10 points cell means are higher than remaining 8 times of 89 average value More than, then there is " spike ", spike is otherwise not present.
Fig. 1 to Fig. 7 illustrates different experiment curvs and different handling results.
The present embodiment signal is transformed into digital signal by TDLAS technology opto-electronic conversions, and by the A/D signals of ARM, with Under be to further illustrate this method by specific data:
It is preserved in a manner of one-dimensional fixed long number group by the signal collected, length is usually 100 points.Signal magnitude Corresponding A/D acquisitions, the A/D digits of ARM are 12, then signal corresponds to 0 ~ 4096.
It carries out judging whether that TDLAS signal probe light paths are blocked by foreign matter first, numerical values recited ratio is carried out to one-dimensional data It is right, according to test experience, when blocking light path, for digit it is that 12 signal amplitudes are generally below 10, and normal signal is then not It will appear this signal amplitude, therefore after the numerical value appearance in one-dimension array less than 10, this time acquisition signal is directly deleted.
If not being less than the signal of 10 amplitudes, its concentration of primary Calculation, when gas concentration is higher than 1500ppmV units, It is then high concentration signal, is then low concentration signal less than this value.
Low concentration signal processing:
The noise that low concentration signal often occurs is under low concentration, and gas flow rates are too fast, Gas concentration distribution is uneven Caused by signal curve there is " spike ", and " spike " amplitude is far more than signal amplitude, it is therefore desirable to carry out noise to it and go It removes.
By the way of first order derivative, when array length is usually at 100, the first order derivative number of one-dimension array is obtained Group length is 99 points, which takes maximum 10 points after its numerical ordering, then calculates the average value for being left 89 points;Most 10 big points equally calculate average value;Compare the two average value size, after tested, when interference noise occurs, kurtosis is universal Higher than 10 times of signal or more, when maximum 10 points cell means are higher than remaining 8 times of 89 average value or more, then explanation is deposited At " spike ", if spike is not present in anti-regular explanation, without pretreatment, then signal continuation is follow-up calculates.For sampled data, ask Leading the formula that formula can be reduced to is:
Wherein:N >=1, x (n) indicate the discrete data array of acquisition, and n indicates number of the data in array, such as from 1 The array after derivation is indicated to 99, y (n).
If there are spike, substituted using the previous value in maximum 10 positions of array sort in original signal acquisition array The position numerical value.Then signal peaks are removed.
Then look for the array maximum position, be peak position, if peak position between 50 ~ 60, the signal without Problem continues signal and subsequently calculates;
If peak position exceeds the part, illustrate laser absorption peak not in central area, surface laser device temperature control occurs Deviate, then the data exist abnormal, then the signal is deleted.
High concentration signal processing:
The noise that high concentration signal often occurs is that absorption peak deviates center or curve is distorted, and the distortion is general It is that sensor probe light path is abnormal refraction or receives position deviation.
High concentration signal and standard curve are subjected to cross correlation algorithm first, standard curve is normally believed in extraction test Number, and by it be filtered after obtained normalized curve.Related algorithm is carried out therewith, can obtain the curve after related algorithm.
Standard curve is at 20 DEG C, in the case of 50%RH, using humidity device1 by environment set under this condition, profit Data acquisition is carried out with TDLAS sensors, after acquiring 100 times, data is carried out averagely and filtering obtains canonical function.
The formula of related algorithm:
Wherein:Xn indicates that the array of canonical function, y (m) indicate that the array of data after acquisition, n indicate the number of canonical function According to the number in array, m indicates number of the data after acquisition again in array.
If the signal curve of acquisition is normal, which only obtains peak value in central area;If curve occurs various Distortion, then after related algorithm, it may appear that a variety of situations such as multiple small leaks, central peak deviation.
After related algorithm, in related algorithm one-dimension array, peak-seeking algorithm is carried out, that is, judges data maximums, will obtain The peak of curve position judges then to illustrate that signal is normal in central area, if the not specification exception if, signal deletion;Meanwhile Using first order derivative combination second derivative mode, each peak position in cross correlation algorithm is marked out, if only there is a main peak Illustrate that signal is normal, if occurring being more than 1 small peak, illustrates abnormal signal, which deletes.
Matlab algorithms, first order derivative combination second derivative judge the position at multiple peaks
t = 1:length(n);
%t1= 0:length(n)-1;
Lmin = diff(sign(diff(n)))== 2; % logic vector for the local min value
Lmin = [false; Lmin; false];
tmin = t (Lmin); % locations of the local min elements
vmin = n (Lmin); % values of the local min elements
After Signal Pretreatment, it may be determined that whether signal is abnormal, and abnormal signal is deleted, normal signal or through place Signal after reason continues to calculate in next step.

Claims (2)

1. Laser Modulation absorbs gas sensor signal preprocess method, the sensor analog signals of acquisition are converted into digital letter Number, it is characterised in that:
Digital signal is preserved in a manner of one-dimensional fixed long number group;
The array of preservation is carried out blocking optical path signal according to the threshold value of setting and is deleted, is calculated after will be less than the number removal of threshold value Preliminary concentration distinguishes high concentration signal and low concentration signal;
Low concentration signal noise removal process is carried out when for low concentration signal;High concentration signal is carried out when for high concentration signal Noise removal process;
Signal after noise removal process is sent into subsequent processing steps, pretreatment is completed;
The long array of fixation is the array of continuous 100 points, and the analog signal is converted to the A/D conversion digits of digital signal It it is 12, the threshold value is 10;
The method of the low concentration signal noise removal process is:
It is first determined whether there are spikes, if spike is not present, continue subsequent processing steps without processing;If there is point Peak, the then value that original signal is acquired to the maximum previous point in multiple positions of array sort substitute the point value;Then look for this Array maximum point position, if subsequent processing steps are continued in maximum point position in the centre position of array;If maximum point position is not In centre position, then the array exists abnormal, deletes the array;
The method of the high concentration signal noise removal process is:
High concentration signal and standard curve are subjected to the curve after related algorithm obtains related algorithm first, if after related algorithm There are a minor peaks in curve, and peak value is in the centre position of curve;Or there is multiple peak value in curve, and repeatedly in peak value Between peak value be in the centre position of curve, then normally to continue subsequent processing steps;Otherwise the array is deleted;
The multiple point is 5% to 10% always to count;
The standard curve is the curve made previously according to ideal data;
The method for judging whether spike is:If fixed long array is the array of continuous 100 points, to continuous 100 Point low concentration signal obtains the first order derivative array of one-dimension array by the way of first order derivative, and the length of array is 99 points, The array takes maximum 10 points after its numerical ordering, and 10 points of the average value and maximum point that then calculate remaining 89 points are flat Mean value;Both compare average value size, when maximum 10 points cell means are higher than remaining 8 times of 89 average value or more, Then there is " spike ", spike is otherwise not present.
2. preprocess method according to claim 1, which is characterized in that the differentiation high concentration signal and low concentration signal It is:Preliminary concentration is high concentration signal higher than 1500ppmV units, is then low concentration signal less than this value.
CN201510788547.7A 2015-11-17 2015-11-17 Laser Modulation absorbs gas sensor signal preprocess method Active CN105447455B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510788547.7A CN105447455B (en) 2015-11-17 2015-11-17 Laser Modulation absorbs gas sensor signal preprocess method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510788547.7A CN105447455B (en) 2015-11-17 2015-11-17 Laser Modulation absorbs gas sensor signal preprocess method

Publications (2)

Publication Number Publication Date
CN105447455A CN105447455A (en) 2016-03-30
CN105447455B true CN105447455B (en) 2018-11-06

Family

ID=55557615

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510788547.7A Active CN105447455B (en) 2015-11-17 2015-11-17 Laser Modulation absorbs gas sensor signal preprocess method

Country Status (1)

Country Link
CN (1) CN105447455B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011242222A (en) * 2010-05-17 2011-12-01 Shimadzu Corp Gas concentration measuring apparatus
CN102706834A (en) * 2012-02-17 2012-10-03 重庆同博测控仪器有限公司 Gas concentration value over-limit detection method for optical interference gas detection device
CN103090886A (en) * 2011-10-31 2013-05-08 深圳迈瑞生物医疗电子股份有限公司 Denoising method and denoising device for optical signal

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011242222A (en) * 2010-05-17 2011-12-01 Shimadzu Corp Gas concentration measuring apparatus
CN103090886A (en) * 2011-10-31 2013-05-08 深圳迈瑞生物医疗电子股份有限公司 Denoising method and denoising device for optical signal
CN102706834A (en) * 2012-02-17 2012-10-03 重庆同博测控仪器有限公司 Gas concentration value over-limit detection method for optical interference gas detection device

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于TDLAS技术的低浓度CO检测系统设计;王金明 等;《仪表技术与传感器》;20150930(第9期);第73-76页 *

Also Published As

Publication number Publication date
CN105447455A (en) 2016-03-30

Similar Documents

Publication Publication Date Title
CN105651325B (en) The deterioration sensor detection realized in transmitter
CN105527650B (en) Microseismic signals and p ripple first arrival automatic identification algorithms under a kind of engineering yardstick
CN109379379A (en) Based on the network inbreak detection method for improving convolutional neural networks
CN108535354B (en) Damage judgment and positioning method for magnetic flux leakage detection and magnetic emission detection of steel wire rope
CN116881745B (en) Pressure transmitter abnormality monitoring method based on big data
CN111521915B (en) High-voltage direct-current line corona onset field strength determination method and system
KR20140141051A (en) Yield management system and method for root cause analysis using manufacturing sensor data
KR101168968B1 (en) Detecting and filtering method of random noise signal
CN104794839B (en) One kind is based on POTDR fiber optic intrusion recognizer
CN106933127A (en) A kind of control method of gas sensor
CN116738353B (en) Pharmaceutical workshop air filter element performance detection method based on data analysis
CN107276779B (en) Monitoring method, system and equipment
CN111353131B (en) Code carrier deviation degree threshold value calculation method
CN117191305A (en) State evaluation method and system for highway bridge
CN115165079A (en) Equipment alarm method based on multiple sound characteristic analysis
CN105447455B (en) Laser Modulation absorbs gas sensor signal preprocess method
CN114186596B (en) Multi-window identification method and device for spectrogram peaks and electronic equipment
CN117761263B (en) Formaldehyde gas detection system and method
CN117314020B (en) Wetland carbon sink data monitoring system of plankton
CN103558647A (en) Method for detecting sampling probe blocking abnormity based on amplitude-frequency characteristic
CN109657329A (en) Based on the sensor fault and Structural Damage Identification for improving adaptive particle filter
CN105133104B (en) A kind of on-line checking and the method and its device for removing Compact-spun Yarn grid ring failure yarn defect
CN114722885B (en) Intelligent detection method and system for abnormal operation of carrier roller carrying trolley
US20050061057A1 (en) Interference detection systems and methods
CN106444689B (en) Intelligent housing central apparatus based on human body infrared inductor

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
PE01 Entry into force of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Laser modulation absorption gas sensor signal pretreatment method

Effective date of registration: 20191111

Granted publication date: 20181106

Pledgee: Beijing Yizhuang International Financing Guarantee Co.,Ltd.

Pledgor: BEIJING AEROSPACE YILIAN SCIENCE AND TECHNOLOGY DEVELOPMENT Co.,Ltd.

Registration number: Y2019990000503

PC01 Cancellation of the registration of the contract for pledge of patent right

Date of cancellation: 20211012

Granted publication date: 20181106

Pledgee: Beijing Yizhuang International Financing Guarantee Co.,Ltd.

Pledgor: BEIJING AEROSPACE YILIAN SCIENCE AND TECHNOLOGY DEVELOPMENT Co.,Ltd.

Registration number: Y2019990000503

PC01 Cancellation of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Signal preprocessing method for laser modulated absorption gas sensors

Effective date of registration: 20230519

Granted publication date: 20181106

Pledgee: CHINA ACADEMY OF AEROSPACE AERODYNAMICS

Pledgor: BEIJING AEROSPACE YILIAN SCIENCE AND TECHNOLOGY DEVELOPMENT Co.,Ltd.

Registration number: Y2023980041196

CP03 Change of name, title or address
CP03 Change of name, title or address

Address after: Building 24, 18th Courtyard, Kechuang 13th Street, Beijing Economic and Technological Development Zone, Daxing District, Beijing, 100176

Patentee after: Beijing Aerospace Tianhong Intelligent Equipment Technology Co.,Ltd.

Address before: Building 24, yard 18, Kechuang 13th Street, Daxing Economic and Technological Development Zone, Beijing 100176

Patentee before: BEIJING AEROSPACE YILIAN SCIENCE AND TECHNOLOGY DEVELOPMENT Co.,Ltd.

PC01 Cancellation of the registration of the contract for pledge of patent right
PC01 Cancellation of the registration of the contract for pledge of patent right

Granted publication date: 20181106

Pledgee: CHINA ACADEMY OF AEROSPACE AERODYNAMICS

Pledgor: BEIJING AEROSPACE YILIAN SCIENCE AND TECHNOLOGY DEVELOPMENT Co.,Ltd.

Registration number: Y2023980041196

PE01 Entry into force of the registration of the contract for pledge of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Signal preprocessing method for laser modulation absorption gas sensors

Granted publication date: 20181106

Pledgee: CHINA ACADEMY OF AEROSPACE AERODYNAMICS

Pledgor: Beijing Aerospace Tianhong Intelligent Equipment Technology Co.,Ltd.

Registration number: Y2024990000150