The Portable Gas Analyzer Based on the Spectrum

Similar documents
Pressure Sensitive Paint (PSP) / Temperature Sensitive Paint (TSP) Part 1

A Coal Mine Multi-Point Fiber Ethylene Gas Concentration Sensor

FLUORESCENCE DETERMINATION OF OXYGEN

MATRIX-MG Series. Innovation with Integrity. Automated High-Performance Gas Analyzers FT-IR

Journal of Chemical and Pharmaceutical Research, 2014, 6(3): Research Article

Real-time Analysis of Industrial Gases with Online Near- Infrared Spectroscopy

Alphasense Application Note AAN NDIR: Gas Concentration Calculation Overview

Section 1 Types of Waves. Distinguish between mechanical waves and electromagnetic waves.

New Generation System M, leading the World in the Non-Invasive Measurement of Critical Real-Time Parameters.

Section 4.2. Travelling Waves

Laser Spectrometers for Online Moisture Measurement in Natural Gas. Ken Soleyn GE M&C

Waves. harmonic wave wave equation one dimensional wave equation principle of wave fronts plane waves law of reflection

HIGH ACCURACY MULTI-GAS MONITORING USING AUTOMATED SELF- CALIBRATION

2600T Series Pressure Transmitters Plugged Impulse Line Detection Diagnostic. Pressure Measurement Engineered solutions for all applications

LOW PRESSURE EFFUSION OF GASES revised by Igor Bolotin 03/05/12

Chapter 20 - Waves. A wave - Eg: A musician s instrument; a cell phone call & a stone thrown into a pond A wave carries from one place to another.

DEVICES FOR FIELD DETERMINATION OF WATER VAPOR IN NATURAL GAS Betsy Murphy MNM Enterprises 801 N. Riverside Drive Fort Worth, Texas 76111

Level MEASUREMENT 1/2016

Basketball field goal percentage prediction model research and application based on BP neural network

Title: Standard Operating Procedure for Measurement of Ethylene (C 2 H 4 ) in Ambient Air by Reduced Gas Detection (RGD)

Slide 2 / 28 Wave Motion. A wave travels along its medium, but the individual particles just move up and down.

Title: Standard Operating Procedure for Dasibi Model 5008 Gas Dilution Calibrator

White Paper. Chemical Sensor vs NDIR - Overview: NDIR Technology:

Parts of Longitudinal Waves A compression

Outline Chapter 7 Waves

ANALYSIS OF THE POSITIVE FORCES EXHIBITING ON THE MOORING LINE OF COMPOSITE-TYPE SEA CAGE

Wave. 1. Transverse 2. Longitudinal 3. Standing

Chapter 7. SCIIB Pressure Sensor Performance Evaluations: Experiments, Results and Discussions

LOW PRESSURE EFFUSION OF GASES adapted by Luke Hanley and Mike Trenary

Implementation of Height Measurement System Based on Pressure Sensor BMP085

Chapter # 08 Waves. [WAVES] Chapter # 08

20.1 Waves. A wave is an oscillation that travels from one place to another. Because waves can change motion, they are a travelling form on energy.

HCL MONITORING. ICAC MARAMA WORKSHOP May 18-19th, Presented by: Ty Smith. CEMTEK

Dr. Tetsuzo KITAGAWA, Tokihiko YONEDA, Shigefumi NAKAJIMA. 1. Preface

Dissolved Oxygen and measurement possibilities. Berno Lüpkes, 15 th March 2017

Hardware Triaxial and Consolidation Testing Systems Pressure Measurement and Control

Columbus Instruments

Wave Motion. interference destructive interferecne constructive interference in phase. out of phase standing wave antinodes resonant frequencies

Chapter 15 Wave Motion. Copyright 2009 Pearson Education, Inc.

Acoustic Emission Testing of The Shell 0f Electromagnetic Valve

RESPIRATORY PHYSIOLOGY, PHYSICS AND

RESEARCH OF BLOCKAGE SEGMENT DETECTION IN WATER SUPPLY PIPELINE BASED ON FLUID TRANSIENT ANALYSIS ABSTRACT

Installation and Operation Manual

ACCURATE PRESSURE MEASUREMENT FOR STEAM TURBINE PERFORMANCE TESTING

APPLICATION OF ACOUSTIC EMISSION DETECTION TECHNOLOGY ON MULTILAYER BINDING HIGH PRESSURE VESSEL

Fail Operational Controls for an Independent Metering Valve

Thermo Scientific Model 146i Principle of Operation

Numerical simulation of three-dimensional flow field in three-line rollers and four-line rollers compact spinning systems using finite element method

Characterizers for control loops

General Information. Department of Physics, Kansas State University, 116 Cardwell Hall Manhattan, KS 66506, USA. Education

EXPERIMENTAL RESULTS OF GUIDED WAVE TRAVEL TIME TOMOGRAPHY

Dissolved Oxygen Guide

u = Open Access Reliability Analysis and Optimization of the Ship Ballast Water System Tang Ming 1, Zhu Fa-xin 2,* and Li Yu-le 2 " ) x # m," > 0

Applications Note: Use of "pentane equivalent" calibration gas mixtures

sensors ISSN

A new method for noninvasive venous blood oxygen detection

Gas Analysis Isn t Exactly Your Only Priority? Luckily For You, It s Ours. SIPROCESS GA700 effortless gas analytics from Siemens. siemens.

Constant pressure air charging cascade control system based on fuzzy PID controller. Jun Zhao & Zheng Zhang

Technical Data Sheet MF010-O-LC

THE PRESSURE SIGNAL CALIBRATION TECHNOLOGY OF THE COMPREHENSIVE TEST

Standard Operating Procedure for H 2 S and SO 2 Converter

AMS 6915 Board mount pressure sensor with digital output (I²C)

Chs. 16 and 17 Mechanical Waves

METHOD 25A - DETERMINATION OF TOTAL GASEOUS ORGANIC CONCENTRATION USING A FLAME IONIZATION ANALYZER

EXPERIMENTAL RESEARCH ON COEFFICIENT OF WAVE TRANSMISSION THROUGH IMMERSED VERTICAL BARRIER OF OPEN-TYPE BREAKWATER

Protective Atmospheres, Measurement Technologies and Troubleshooting Tools

Advanced Test Equipment Rentals ATEC (2832) OMS 600

Relative Dosimetry. Photons

Experimental Research on Oxygen-enriched Air Supply System for Gasoline Engine

A NOVEL SENSOR USING REMOTE PLASMA EMISSION SPECTROSCOPY FOR MONITORING AND CONTROL OF VACUUM WEB COATING PROCESSES

11/22/ (4) Harmonization: <846> SPECIFIC SURFACE AREA

DEPARTMENT OF THE NAVY DIVISION NEWPORT OFFICE OF COUNSEL PHONE: FAX: DSN:

Vibration Analysis and Test of Backup Roll in Temper Mill

Currents measurements in the coast of Montevideo, Uruguay

The M-Series Eletta Flow Meter High accuracy DP Flow Meter with multiple functions

4.4 WAVE CHARACTERISTICS 4.5 WAVE PROPERTIES Student Notes

Title: Standard Operating Procedure for R&R Environmental Devices Model MFC201 Gas Dilution Calibrator

Generating Calibration Gas Standards

Columbus Instruments

IDeA Competition Report. Electronic Swimming Coach (ESC) for. Athletes who are Visually Impaired

Study on Mixing Property of Circulation/Shear Coupling Mechanism-Double Impeller Configuration

Analysis of wind resistance of high-rise building structures based on computational fluid dynamics simulation technology

Title: Standard Operating Procedure for Elemental and Organic Carbon (EC and OC) using Non-Dispersive Infrared Detection (NDIR)

LFE OEM TCD - Thermal Conductivity Detector

Analysis and Research of Mooring System. Jiahui Fan*

Correction of Pressure Drop in Steam and Water System in Performance Test of Boiler

Transverse waves cause particles to vibrate perpendicularly to the direction of the wave's motion (e.g. waves on a string, ripples on a pond).

Preview. Vibrations and Waves Section 1. Section 1 Simple Harmonic Motion. Section 2 Measuring Simple Harmonic Motion. Section 3 Properties of Waves

Cover Sheet-Block 6 Wave Properties

INTRODUCTION TO THE SPECTROPHOTOMETER AND PIPETTING SKILLS

AUK Environmental Sensor. Determining the refractive index of the air. AUK Environmental Sensor

Setting a New Standard in Flash Chromatography Performance

LUMAZOTE T ECHNICAL BRO CHU RE

Pre AP Physics: Unit 7 Vibrations, Waves, and Sound. Clear Creek High School

TEK-THERMAL 1700B. Thermal Mass Flowmeter. FLOW. Technology Solutions

A Fault Diagnosis Monitoring System of Reciprocating Pump

GasSense NDIR User Manual

I.CHEM.E. SYMPOSIUM SERIES NO. 97 BUOYANCY-DRIVEN NATURAL VENTILATION OP ENCLOSED SPACES

LABORATORY EXPERIMENTS ON WAVE OVERTOPPING OVER SMOOTH AND STEPPED GENTLE SLOPE SEAWALLS

Vertical in situ profiles of nitrate and oxygen in the northern Japan Sea

Transcription:

018 International Conference on Sensor Network and Computer Engineering (ICSNCE 018) The Portable Gas Analyzer Based on the Spectrum Xu Shuping Xi an 71001, China e-mail: 563937848@qq.com; xusp686@163.com Dong Qiyu Xi an 71001, China Xu pei Xi an 71001, China Su Xiaohui Xi an 71001, China Abstract It is a problem of spectra analysis of flue gas that how to separate and calculate the concentration of different kinds of gas from continuous mixed gas absorption spectrum signal. So based on experimental data, a new iteration of the circular algorithm is put forward on the basis of Lambert-Beer's law. This algorithm, using the superposition of the absorbance, makes data nonlinear fitting. It takes the advantages of the wavelength optimum and the circular iteration to distinguishing the gas mixture of different composition. Experimental results show that: this method can at once calculate a variety of harmful gas concentration and accuracy up to ± 3%. it has strong anti-jamming capability and is suitable for practical application of engineering. Keywords-Circular Iteration; Ultraviolet Spectrum; Micro-spectrometer; Embedded Systems I. INTRODUCTION With the rapid development of economic and people living standard unceasing improvement the popularity of the development of heavy industry and transportation, with a large number of harmful substances into the environment atmosphere, changed the composition of normal air, worsen air quality[1]. Health problems, which caused by the atmospheric pollution has caused governments and people of great importance to it. Therefore, monitoring SO and NO x become one of the important subject of environmental protection[]. For the accurate monitoring of the real-time quality, ecological environment and pollution to the environment, to provide accurate basis for the environmental protection departments at various levels supervision, management and environment decision, urgent need a large number of modern environmental monitoring instruments. Through to the present situation of flue gas analyzer in the domestic and foreign research and analysis, portable flue gas on the domestic market at present the main monitoring method is electrochemical analysis instrument, infrared gas analyzer and differential absorption spectrum. Electrochemical analysis instrument each sensor can only be measured a chemical composition in flue gas, if need to measure six ingredients you need six sensors. The chemical battery principle led to zero drift in the development of the system and cross interference, sensor lifetime is also difficult to solve the problem[3]. There is no absorption peak to NO in the infrared spectrum analyzer, cannot monitor NO, but can only detect the NO content, through the assumptions to NO concentrations. At the same time as infrared spectrum analyzer adopts wheel filter method, each time only one gas measurement, can not measure a variety of gases at the same time[4]. Differential absorption measured gas to uv and visible light waves, on the basis of the differential absorption spectrum by the strength of the differential 13

018 International Conference on Sensor Network and Computer Engineering (ICSNCE 018) absorption spectrum inversion gas concentrations, in the traditional algorithm of differential absorption, the application of least square method of gas concentration for calculating optimal estimation equation is complex, and easy in pathological conditions[5]. Accordingly based on Lambert Beer law and absorbance of dual stack, a mutually recursive iteration is presented more gas concentration decoding algorithm, the process of the algorithm is presented, and verify the effectiveness of the algorithm. This system using the Ocean Optic micro spectrometer combining with embedded technology company has developed a set of portable flue gas analyzer[6]. This instrument is to overcome the disadvantages of the electrochemical sensor, had been achieved and sensor life problems, and at the same time for accurate measurement of gases and no cross interference. II. WORKING PRINCIPLE OF THE SYSTEM Portable flue gas analyzer system structure consists of three parts, part data collection, data processing part and the human-computer interaction part. Data acquisition part by sampling probe, flue gas sampling pump, an ultraviolet light and miniature spectrum analyzer, collecting gas mainly accomplished by ultraviolet light and miniature spectrum analyzer spectrum before and after the data is converted into digital signal transmitted to the data processing part. Data processing part of ARM microcontroller CPU core module and peripheral circuit, main is to electrical signals were collected by using corresponding algorithm to calculate the gas concentration of all kinds of gas in the flue gas composition. The human-computer interaction part consists of LCD display and keyboard, mainly to complete the real-time display of gas concentration, and through the keyboard to control the storage of data and related parameters Settings. III. THE BASIC PRINCIPLE OF LOOP ITERATION METHOD A. Lambert-Beerlaw The mathematical model of Lambert-Beer's law [7-8] said = In the formula: A is absorbance; T is the transmittance, project the light intensity is than the incident light intensity;c is light-absorbing substance concentration; B is the absorption layer thickness. Its physical meaning is: when a beam of parallel light through a vertical uniform a material suction light scattering, the absorbance and the concentration and absorption material suction light is directly proportional to the thickness. Absorbance [16] A binary additive think: if the binary and multicomponent mixture components have absorbed A wave number, the total absorbance at the wave number is equal to the arithmetic of absorbance and at all levels. B. The Mathematical model The logarithmic of Lambert-Beer's law said: R Alg S D D Kc In the formula: A is absorbance; R is photon number (zero point); S is the number of through the photon; λ is selected for a wavelength; K is constant (associated with the length of wavelength and absorption tube); C is a single gas concentration; D is dark spectrum photon number (related to the integration time can take a fixed constant). Two kinds of mixed gas Lambert-Beer absorbance expression can be represented as: 1 1 lg R D A S1 D lg R D A S D According to the law of superposition of binary: A A1A R D R 1 D R D lg lg lg S D S1 D S D i.e. 14

018 International Conference on Sensor Network and Computer Engineering (ICSNCE 018) R D R1 D R D S D S1 D S D By formula (7) can be know in a particular incident light wavelength lambda, the mixed gas of photon number and subtract its dark photons through the photon number, equal to the selected wavelength under different gas on the number of incident photons and through the photon number minus the dark photon number of the product. In the actual gas absorption model of gas absorption and scattering all play a role in the formation of spectrum, so the experiment of gas absorption spectra contain two parts. Part of gas molecules on the absorption of photons, ultraviolet band is mainly caused by the electron transition of the gas molecules absorbed. Another part of the scattering is gas or smoke without selective absorption caused by light attenuation. Because this system adopts the removable detection method, gas extraction, you first after a pretreatment device for flue gas filtration, drying, cooling and other processing. Filtration process can be 0. Lum or above level particle filter out completely, and flue gas molecules of the diameter of the order of magnitude are below 1 nm. Instrument of the selected work between the wavelength of 190 nm - 90nm, greater than the molecules or residual particle diameter, so that the gases to be detected is mainly gas molecules Rayleigh scattering. The actual operation process we with nitrogen gas as zero, the spectral curve is called the zero gas, nitrogen gas molecules, although about 190 nm - 90 nm uv absorption, but also produces Rayleigh scattering. That is zero gas line is light after zero gas scattering spectral lines. By spectrogram, found that in the absence of absorption bands, such as nitric oxide line outside the three absorption peak is coincidence with zero gas lines, namely in the absence of gas absorption bands of two gases scattering light is the same. So the Rayleigh scattering of gas absorption interference can be ruled out by the method of using the zero gas line as a reference. According to the Lambert Beer-law, exclude mie scattering and Rayleigh scattering in gas detection, then the device applicable gas absorption model becomes: I( ) ln( ) L( i( ) Ci) I0( ) Which represents the light source through I ( ) 0 zero gas scattering by spectrometer in wavelength after the received light intensity, i ( ) represents the ith kind of gases at wavelength absorption cross section, L is the length of the pool by absorption, and C i is the ith the concentration of the gas. Which gas is equal to the total absorbance and of each gas absorption degree. IV. CIRCULAR ITERATION STEPS Circular iteration method using various gases in 190 ~ 90 nm band characteristic absorption peak, combined with the absorbance of dual stack, a feature in some gas absorption point assuming that other gas no absorption, launched the initial concentration of the gas, then switch to another characteristic absorption point, the gas absorbed photon number from the measured total absorbed photon number subtracting, get another gas initial concentration, and so on to get the initial concentration of each gas. Then return to the first, the characteristics of the gas absorption point from the total number of photons absorbed minus every other gas absorption of the photon number, again to get the first gas concentration of an iteration, and so on to get other gas concentrations of an iteration, the repeated iteration until the concentration difference between two times smaller than a certain value, the concentration of each gas to the gas concentration accurately. Algorithm steps are as follows. First step: to solve the first gas in the gas mixture of initial concentration C 10, and selects a characteristic wavelength 1, under the specific wavelength of a gas in the gas mixture has obvious characteristic absorption peak, and other gas absorption peak is small at this point, is to read the value of the absorbed photon number S 1and solve R S D D is obtained by absorbance look-up table and the concentration of the gas as the initial concentration of mixed gas in the gas. Step : solving the second gas in the gas mixture of initial concentration: select the second characteristic wavelength in the wavelength of second kind of gas has obvious absorption peak, and other gas absorption peak is weak, then 15

018 International Conference on Sensor Network and Computer Engineering (ICSNCE 018) read the absorbed photon number S considered under the band is only two kinds of mixed gas, according to formula (7) to calculate the second gas absorbance, look-up table against solving calculation, the second gas concentrations, as the second kind of initial concentration C 0 of gases. Step 3: calculating of other gases in the mixed gases, the initial concentration: the mth characteristic wavelengths selected m,the wavelength of the first m gas has obvious absorption peak, read the wavelength of absorption of the photon number S m by formula (7) and look-up table type inversion calculation of this kind of gas concentration as the initial concentration Cm0of the gas. Step 4: iterative inversion to calculate the first gas concentration of first order recursive: will the desires of all kinds of gas concentration in formula (7), again read 1 of S 1, reverse the first gas concentration C11 of first order recursive. Step 5: repeat the second and third step to calculate interfered with each other. At wavelength of 1 (73.33nnl) So and No are absorbed. Also at wavelength of (31.33nm) Soand No are absorbed. Only at wavelength of 3 (5.88nm) NO, Soand No are absorbed. At wavelength of 4 (08.3nm) NH 3 No and So are absorbed. TheSoand No have absorbed at each wavelength point, For NH 3 and NO, and, as long as the concentration of So and No first came out, and then the corresponding concentration on wavelength 3 or 4 of the absorbance minus the can get their absorbance. So should first of all, the concentration ofsoand No, then the concentration of NH 3 and NO. gas m first order recursive concentration C 1m. Step 6: calculation of gas concentration of the adjacent two iterative error; Calculating first-order differential iteration of each gas is shown in the following formula, C C 1 om m1 Select 1maximum of first-order differential iteration as the iteration error, i.e 1 max m1, m 1,,..., M Step 7: repeat the fourth, fifth and sixth step, until the iteration error is less than the given value namely. Figure 1. Mixed gas absorption spectrum A. Calculate the concentration of So and No For the concentration of NH3 and NO are the premise of know the concentration of SoandNo. However, Soand No throughout the working wave band are interfered with each other, need an algorithm to eliminate their interference. Here, on the basis of the look-up table using loop iteration method and calculation method to determine the concentration of Soand No. n G Step 8: termination of the algorithm, it will be the last time the concentration as a final concentration of gases. V. THE EXPERIMENTAL RESULTS AND ANALYSIS The obtained by experiments, So, No, the NO and NH 3 four gas absorption lines as shown in figure 1, in the range of 190 nm to 190 nm, four types of gas are So Figure. The absorption spectrum of No and 16

018 International Conference on Sensor Network and Computer Engineering (ICSNCE 018) As shown in figure, At the wavelength of 1 and at the So and No were absorbed, they differ greatly in the absorption intensity at two wavelengths, the absorbance of the absorbance of No at 73.33nm is far less than the same concentration of No in 31.33nm, and the So absorbance at 31.33nm is far less than the same concentration of So in the absorbance at 73.33nm. Then we can be absorbed in 31.33nm No called the main absorption of No, No in the absorption of 73.33nm is called So interference, So in the absorption at 73.33nm is the main absorption, in 31.33 absorption interference on No main absorption. Iterative method is through circulation calculation to gradually eliminate the interference of No on So, the true concentration approaching So and No. Cyclic iteration method to overcome the conflict there understand equations caused no solution of the problem, and this method can be easily implemented by software programming. Statistical chart concentrations of So and No calculated with iteration number as shown in the figure 3-4. The abscissa is the number of iterations, the longitudinal concentration coordinates So and No (ppm). Figure 3. Figure 4. 60 50 40 30 0 10 0 48 47 46 45 44 43 4 41 40 39 1 3 4 5 The results of calculation curves of iteration number and So concentration 1 3 4 5 The results of calculation curves of iteration number No and concentration According to the running results and statistics can be seen in figure, the process of iteration is a successive approximation process, the calculation results of So with the increase in the number of iterations decreases gradually, and tends to be stable, but No is exactly the opposite, results with increasing iteration times increase. Cyclic iterative method in two iterations later effect is not obvious, which is to say as long as two iterations. This can be calculated throughnoand So concentrations. B. Calculate the concentration of NOand NH 3 Because the NO and NH 3 absorption in the selected wavelengths of 3 and a wavelength of 4 does not interfere with each other, so it can be separated by the concentration of total absorbance and the So and NO: calculation. Hypothesis has been obtained in the mixed gas of No and So for C1 concentration and C concentration, the NO and NH 3 for C3 and C4, then at the wavelength of 3 according to the superposition of absorbance of the available, A A A A lg( ) A A I 0 3 1 I 1 A 1 is the C1 concentration of No absorbance at 5.8nm, A is the C concentration of So absorbance at 5.8nm, A is a mixture of gases in the 5.88nm total absorbance at 5.88nm, A 3 is NO of total absorbance. I0 is the spectral intensity of 5.88nm through zero gas, I is the intensity of 5.88nm through transmission into the mixed gas after, can be obtained directly by the spectrometer, A and A 3 concentrations of A and Noobtained by. Absorbance of A3 so you can get the NO in 5.88nm, then according to the concentration and absorbance at 5.88nm correspondence between the NO table to calculate the concentration of NO. The absorbance of A4 NH3 can be obtained with the same method in 08.3nm, then according to the corresponding relationship between concentration and absorbance at 08.3nm NH 3 for NH 3 concentration. VI. CONCLUSION Main harmful components for atmospheric environmental pollution monitoring requirements, using uv wavelength grating type continuous frequency measuring method, 17

018 International Conference on Sensor Network and Computer Engineering (ICSNCE 018) precision is put forward to solve the various harmful ingredient concentration of recursive iteration fast inversion algorithm, and validates the effectiveness of the algorithm. Portable flue gas analyzer based on this algorithm is based on embedded technology, sensor based on micro spectrometer data collection, using uv light source through the spectrum analysis method analysis of flue gas concentration, achieved through the use of a miniature spectrum analyzer to a variety of gas composition at the same time for the purpose of accurate measurement. The product has compact structure, high measuring accuracy, strong anti-interference, high sensitivity etc, has a broad application prospect and popularization value. ACKNOWLEDGMENT The authors wish to thank the cooperators. This research is partially funded by the Project funds in shaanxi province department of education (15JF019) and the Project funds in shaanxi province department of science industrial projects(015gy067). REFERENCE. [1] Zhu Fa-hua, Li Hui, Qiu Shu-guang. Development and Application of Continuous Em ission Monitoring Technology[J]. Administration and Technique of Environmental Monitoring 010() [] Lei Tian-xue. The present situation of the portable flue gas analyzer[j]. Administration and Technique of Environmental Monitoring. 1998(10) [3] Ju Hui, Wu Yi-hui, Development situation of micro spectrometers[j]. Optics and Precision Engineering. 004.9(4). pp37-376 [4] Shi Baosong, Sun Shouhong, Zhang Wei.Application of CCD in the portable spectrometer[j].electronic Measurement Technology. 010(11) [5] Limited ARM Development Guide 000-001. ARM DOI. 000.06. [6] Dong, Qing. Ocean wave spectrum from SAR image using D-ARMA model.ieeeinternational Geosciences and Remote Sensing Symposium, 005, 991-994 [7] Chen Zhi-gang. Discussion on Application of Lambert-Beer law in test [J]. Acta Metrologica Sinica.1985(1) [8] Pop, Paul.Embeddedsystems design:optimization challenges.lecture Notes in ComputerScience, 005, 35(4):16-0 18