CN105447455A - Laser modulation absorption gas sensor signal pretreatment method - Google Patents

Laser modulation absorption gas sensor signal pretreatment method Download PDF

Info

Publication number
CN105447455A
CN105447455A CN201510788547.7A CN201510788547A CN105447455A CN 105447455 A CN105447455 A CN 105447455A CN 201510788547 A CN201510788547 A CN 201510788547A CN 105447455 A CN105447455 A CN 105447455A
Authority
CN
China
Prior art keywords
signal
array
high concentration
curve
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.)
Granted
Application number
CN201510788547.7A
Other languages
Chinese (zh)
Other versions
CN105447455B (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 laser modulation absorption gas sensor signal pretreatment method. Collected sensor simulation signals are converted into digital signals. The digital signals are stored in a one-dimensional fixed long array mode. For the stored array, according to a set threshold value, shielding optical path signal deletion is performed. A figure which is lower than the threshold value is removed and then initial concentration is calculated. High concentration signals and low concentration signals are distinguished. When the signals are the low concentration signals, low concentration signal noise removing processing is performed. When the signals are the high concentration signals, high concentration signal noise removing processing is performed. The signals after noise removing processing are sent into a subsequent processing step and then pretreatment is completed. By using the method, for different interference types, the respective pretreatment method is provided; pretreatment after low concentration and high concentration signal sensing is realized; and inaccurate measurements caused by signal interferences are greatly reduced.

Description

Laser Modulation absorbs gas sensor signal preprocess method
Technical field
The invention belongs to sensor signal treatment technology, particularly a kind of Laser Modulation absorbs gas sensor signal preprocess method, is, before one carries out gas concentration calculating for TDLAS method of testing, to carry out Signal Pretreatment algorithm.
Background technology
TDLAS technology is gas sensing techniques of new generation, and the gas sensor that this technology is made easily receives various signal disturbing under low concentration, high concentration, affects evaluation during meter gaseous, and then affects gas sensing test concentrations.
Summary of the invention
The object of the invention is to propose a kind of Laser Modulation and absorb gas sensor signal preprocess method, before one carries out gas concentration calculating for TDLAS method of testing, carry out Signal Pretreatment algorithm, according to repeatedly testing, find out noise behavior under low concentration, adopt first order derivative, form automatic threshold, remove according to the jump signal that noise causes by threshold value; In higher concentrations, find out out typical curve, carry out related algorithm by typical curve and collection signal curve, according to automatic threshold method, the signal disturbing caused by the various noise such as comprise curve sudden change, noise is large is removed.
To achieve these goals, technical scheme of the present invention is:
A kind of Laser Modulation absorbs gas sensor signal preprocess method, and the sensor analog signals of collection is converted to digital signal, and its method is:
Digital signal is preserved in one dimension fixed length array mode;
Block optical path signal to the array of preserving according to the threshold value of setting to delete, calculate preliminary concentration after the numeral lower than threshold value being removed, distinguish high concentration signal and low concentration signal;
When for carrying out low concentration signal noise removal process during low concentration signal; When for carrying out high concentration signal noise removal process during high concentration signal;
Signal after noise removal process is sent into subsequent processing steps, and pre-service completes.
Scheme is further: described differentiation high concentration signal and low concentration signal are: preliminary concentration is high concentration signal higher than 1500ppmV unit, is then low concentration signal lower than this value.
Scheme is further: described fixed length array is the array of continuous 100 points, and it is 12 that the A/D that described simulating signal is converted to digital signal changes figure place, and described threshold value is 10.
Scheme is further: the method for described low concentration signal noise removal process is:
First judge whether to there is spike, if there is not spike, then do not carry out process and continue subsequent processing steps; If there is spike, then value original signal being gathered the previous point in multiple maximum positions of array sort substitutes this point value; Then find this array maximum point position, if maximum point position is in the centre position of array, then continue subsequent processing steps; If maximum point position is not or not centre position, then this array exists abnormal, deletes this array;
The method of described high concentration signal noise removal process is:
First high concentration signal and typical curve are carried out related algorithm and draw the curve after related algorithm, if a minor peaks appears in the curve after related algorithm, and peak value is in the centre position of curve; Or repeatedly peak value appears in curve, and repeatedly the interpeak of peak value is in the centre position of curve, then for normally to continue subsequent processing steps; Otherwise delete this array;
Described multiple point is always count 5% to 10%;
Described typical curve is in advance according to the curve that ideal data makes;
Describedly judge whether that the method that there is spike is: if fixed length array is the array of continuous 100 points, continuous 100 some low concentration signals are adopted to the mode of first order derivative, obtain the first order derivative array of one-dimension array, the length of array is 99 points, 10 points that this array is maximum after getting its numerical ordering, then calculate the remaining mean value of 89 points and 10 some mean values of maximum point; , when maximum 10 cell means of counting are higher than remaining 89 some mean values more than 8 times, then there is " spike ", otherwise there is not spike in mean value size relatively.
The invention has the beneficial effects as follows: by this method, for different interference types, propose respective preprocess method respectively, achieve the pre-service after low concentration, high concentration signal sensing, significantly reduce the measurement that causes of signal disturbing inaccurate.
Below in conjunction with drawings and Examples, invention is described in detail.
Accompanying drawing explanation
Fig. 1 is low concentration abnormal signal curve;
Fig. 2 is the signal curve after process;
Fig. 3 be high concentration signal be correlated with use contrast standard curve;
Fig. 4 is a kind of abnormal signal under high concentration;
Fig. 5 is to Fig. 4 signal result after relevant treatment;
Fig. 6 is another kind of abnormal signal under high concentration;
Fig. 7 is to Fig. 6 signal result after relevant treatment.
Embodiment
A kind of Laser Modulation absorbs gas sensor signal preprocess method, and the sensor analog signals of collection is converted to digital signal, and its method is:
Digital signal is preserved in one dimension fixed length array mode;
Block optical path signal to the array of preserving according to the threshold value of setting to delete, calculate preliminary concentration after the numeral lower than threshold value being removed, distinguish high concentration signal and low concentration signal;
When for carrying out low concentration signal noise removal process during low concentration signal; When for carrying out high concentration signal noise removal process during high concentration signal;
Signal after noise removal process is sent into subsequent processing steps, and pre-service completes.
In embodiment: described differentiation high concentration signal and low concentration signal are: preliminary concentration is high concentration signal higher than 1500ppmV unit, is then low concentration signal lower than this value.
Scheme is further: described fixed length array is the array of continuous 100 points, and it is 12 that the A/D that described simulating signal is converted to digital signal changes figure place, and described threshold value is 10.
In embodiment: the method for described low concentration signal noise removal process is:
First judge whether to there is spike, if there is not spike, then do not carry out process and continue subsequent processing steps; If there is spike, then value original signal being gathered the previous point in multiple maximum positions of array sort substitutes this point value; Then find this array maximum point position, if maximum point position is in the centre position of array, then continue subsequent processing steps; If maximum point position is not or not centre position, then this array exists abnormal, deletes this array;
The method of described high concentration signal noise removal process is:
First the curve formed by high concentration signal and typical curve compare, the curve after adopting related algorithm to draw related algorithm, if a minor peaks appears in the curve after related algorithm, and peak value is in the centre position of curve; Or repeatedly peak value appears in curve, and repeatedly the interpeak of peak value is in the centre position of curve, then for normally to continue subsequent processing steps; Otherwise delete this array;
Described multiple point is always count 5% to 10%, and in the present embodiment when fixed length array is the array of continuous 100 points, described multiple clicking is selected as 10 points;
Described typical curve is in advance according to the curve that ideal data makes;
Describedly judge whether that the method that there is spike is: if fixed length array is the array of continuous 100 points, continuous 100 some low concentration signals are adopted to the mode of first order derivative, obtain the first order derivative array of one-dimension array, the length of array is 99 points, 10 points that this array is maximum after getting its numerical ordering, then calculate the remaining mean value of 89 points and 10 some mean values of maximum point; , when maximum 10 cell means of counting are higher than remaining 89 some mean values more than 8 times, then there is " spike ", otherwise there is not spike in mean value size relatively.
Fig. 1 to Fig. 7 illustrates different experiment curv and different results.
The present embodiment signal by TDLAS technology opto-electronic conversion, and is transformed into digital signal by the A/D signal of ARM, is below further illustrate this method by concrete data:
Signal through collecting is preserved in one dimension fixed length array mode, and length is generally 100 points.Signal magnitude corresponding A/D gathers, and the A/D figure place of ARM is 12, then signal correspondence 0 ~ 4096.
First carry out judging whether that TDLAS signal probe light path is blocked by foreign matter, numerical values recited comparison is carried out to one-dimensional data, according to test experience, when blocking light path, that 12 signal amplitudes are generally lower than 10 for figure place, normal signal then there will not be this signal amplitude, and after being therefore less than the numerical value appearance of 10 in one-dimension array, this time collection signal is directly deleted.
If be not less than the signal of 10 amplitudes, then its concentration of primary Calculation, when gas concentration is higher than 1500ppmV unit, being then high concentration signal, is then low concentration signal lower than this value.
Low concentration signal process:
The noise that low concentration signal often occurs is under low concentration, and gas flow rates is too fast, " spike " appears in the uneven signal curve caused of Gas concentration distribution, and " spike " amplitude far exceedes signal amplitude, therefore needs to carry out noise remove to it.
Adopt the mode of first order derivative, when array length is generally 100 points, the first order derivative array length obtaining one-dimension array is 99 points, 10 points that this array is maximum after getting its numerical ordering, then calculates the mean value of remaining 89 points; 10 maximum same calculating mean values of point; Mean value size relatively, after tested, when there is interference noise, kurtosis is generally higher than signal more than 10 times, when maximum 10 cell means of counting are higher than remaining 89 some mean values more than 8 times, then existence " spike " is described, if anti regular illustrates to there is not spike, without the need to pre-service, then signal continues subsequent calculations.For sampled data, the formula that differentiate formula can be reduced to is:
Wherein: n >=1, x (n) represents the discrete data array gathered, and n represents the number of data in array, and such as, from 1 to 99, y (n) represents the array after differentiate.
If there is spike, gather array at original signal and then adopt array sort maximum 10 previous values in some position to substitute this position numerical value.Then signal peaks is removed.
Then find this array maximum position, it is peak, if peak is between 50 ~ 60, then this signal no problem, continues signal subsequent calculations;
If peak exceeds this part, illustrate that laser absorption peak is not in central area, departing from appears in the temperature control of surface laser device, then these data exist abnormal, then this signal is deleted.
High concentration signal process:
The noise that high concentration signal often occurs is that absorption peak departs from center, or curve distorts, and this distortion is generally that the generation anomalous refraction of sensor probe light path or receiving position depart from.
First high concentration signal and typical curve are carried out cross correlation algorithm, typical curve extracts normal signal in test, and the normalized curve will obtained after it filtering process.Carry out related algorithm with it, the curve after related algorithm can be drawn.
Typical curve is at 20 DEG C, in 50%RH situation, uses humidity device1 by environment set under this condition, utilizes TDLAS sensor to carry out data acquisition, after gathering 100 times, carry out data on average and filtering obtain canonical function.
The formula of related algorithm:
Wherein: Xn represents the array of canonical function, y (m) represents the data array after gathering, and n represents the number of the data of canonical function in array, and m represents the number of the data after collection again in array.
If the signal curve gathered is normal, then this related algorithm only obtains peak value in central area; If various distortion occurs curve, then after related algorithm, the multiple situation such as there will be multiple small leak, central peak departs from.
After related algorithm, in related algorithm one-dimension array, carry out peak-seeking algorithm, namely judge data maximums, this peak of curve position will be obtained, judge then to illustrate that in central area signal is normal, if do not exist, specification exception, this signal is deleted; Meanwhile, adopt first order derivative in conjunction with second derivative mode, mark out each peak in cross correlation algorithm, if only there is a main peak, illustrate that signal is normal, be greater than 1 small peak, then abnormal signal be described if occur, this signal is deleted.
Matlab algorithm, first order derivative judges the position at multiple peak in conjunction with second derivative
t=1:length(n);
%t1=0:length(n)-1;
Lmin=diff(sign(diff(n)))==2;%logicvectorforthelocalminvalue
Lmin=[false;Lmin;false];
tmin=t(Lmin);%locationsofthelocalminelements
vmin=n(Lmin);%valuesofthelocalminelements
After Signal Pretreatment, can determine that whether signal is abnormal, abnormal signal is deleted, and normal signal or signal after treatment continue next step and calculate.

Claims (4)

1. Laser Modulation absorbs gas sensor signal preprocess method, the sensor analog signals of collection is converted to digital signal, it is characterized in that:
Digital signal is preserved in one dimension fixed length array mode;
Block optical path signal to the array of preserving according to the threshold value of setting to delete, calculate preliminary concentration after the numeral lower than threshold value being removed, distinguish high concentration signal and low concentration signal;
When for carrying out low concentration signal noise removal process during low concentration signal; When for carrying out high concentration signal noise removal process during high concentration signal;
Signal after noise removal process is sent into subsequent processing steps, and pre-service completes.
2. preprocess method according to claim 1, is characterized in that, described differentiation high concentration signal and low concentration signal are: preliminary concentration is high concentration signal higher than 1500ppmV unit, is then low concentration signal lower than this value.
3. preprocess method according to claim 1 and 2, is characterized in that, described fixed length array is the array of continuous 100 points, and the A/D conversion figure place that described simulating signal is converted to digital signal is 12, and described threshold value is 10.
4. preprocess method according to claim 3, is characterized in that,
The method of described low concentration signal noise removal process is:
First judge whether to there is spike, if there is not spike, then do not carry out process and continue subsequent processing steps; If there is spike, then value original signal being gathered the previous point in multiple maximum positions of array sort substitutes this point value; Then find this array maximum point position, if maximum point position is in the centre position of array, then continue subsequent processing steps; If maximum point position is not or not centre position, then this array exists abnormal, deletes this array;
The method of described high concentration signal noise removal process is:
First high concentration signal and typical curve are carried out related algorithm and draw the curve after related algorithm, if a minor peaks appears in the curve after related algorithm, and peak value is in the centre position of curve; Or repeatedly peak value appears in curve, and repeatedly the interpeak of peak value is in the centre position of curve, then for normally to continue subsequent processing steps; Otherwise delete this array;
Described multiple point is always count 5% to 10%;
Described typical curve is in advance according to the curve that ideal data makes;
Describedly judge whether that the method that there is spike is: if fixed length array is the array of continuous 100 points, continuous 100 some low concentration signals are adopted to the mode of first order derivative, obtain the first order derivative array of one-dimension array, the length of array is 99 points, 10 points that this array is maximum after getting its numerical ordering, then calculate the remaining mean value of 89 points and 10 some mean values of maximum point; , when maximum 10 cell means of counting are higher than remaining 89 some mean values more than 8 times, then there is " spike ", otherwise there is not spike in mean value size relatively.
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 true CN105447455A (en) 2016-03-30
CN105447455B 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检测系统设计", 《仪表技术与传感器》 *

Also Published As

Publication number Publication date
CN105447455B (en) 2018-11-06

Similar Documents

Publication Publication Date Title
CN102654405B (en) Gait counting method and device based on acceleration sensor
RU2011151543A (en) METHOD AND DEVICE FOR DETECTING LEAKS IN THE FUEL SUPPLY SYSTEM
CN108630330B (en) Pressurized water reactor nuclear power station instrument system detector test processing method, device and system
CN106500754A (en) The detection method of sensor and the detection means of sensor
CN108535354B (en) Damage judgment and positioning method for magnetic flux leakage detection and magnetic emission detection of steel wire rope
CN102928225B (en) Self-adaptive resonance demodulation method for fault diagnosis of rolling bearing
US20110060554A1 (en) Interference Detector and Methods
CN104794839B (en) One kind is based on POTDR fiber optic intrusion recognizer
CN105508276A (en) Method for evaluating noise of main transformer fan of urban transformer substation
CN104457911B (en) A kind of liquid level detection system and liquid-level detecting method
CN103558647A (en) Method for detecting sampling probe blocking abnormity based on amplitude-frequency characteristic
CN105259086A (en) Detection method and system of dust concentration
WO2019075913A1 (en) Signal processing system applied to remove otdr noise
CN117314020B (en) Wetland carbon sink data monitoring system of plankton
CN105447455A (en) Laser modulation absorption gas sensor signal pretreatment method
US20200183028A1 (en) Seismic device and safety device employing same
CN105388441B (en) The detection method and system of a kind of magnetic sensor
CN102621170B (en) Method for automatically determining measurement time in detection of energy spectrometer
CN202649982U (en) Capacitive touch screen detecting device
CN106908169B (en) Data processing method of distributed optical fiber temperature sensor
CN105203439A (en) Air purification control method and device
CN210621418U (en) Pavement structure depth detection system
CN104535785A (en) Measurement datanoisereducing method
CN110007342A (en) A kind of direct first break picking method and system of time-frequency domain for low signal-to-noise ratio seismic signal
CN106908172A (en) The signal processing method and system of wireless temperature measurement system

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