Abstract:
An oxygen analyzer which measures the level of dissolved oxygen in a liquid stream. The measuring process embodies chronoamperoemetry of pulsed potential voltammetry. The signal from the electrode is real time averaged. The analyzer does not require temperature, salinity or pressure compensation circuitry and does not have to be calibrated repeatedly.

Description:
BACKGROUND AND BRIEF SUMMARY OF THE INVENTION 
     The present invention relates to an apparatus and method for measuring the level of oxygen in a fluid stream. 
     Presently available devices for measuring dissolved oxygen in situ are subject to certain disabilities, the major disability being the requirement of calibrating the instrument or a daily basis. 
     It has been recognized that an oxygen sensor can be used with a measuring process known as chronoamperoemetry of pulsed potential voltammetry, first described by Lilley et. al., J. Electronanal. Chem., 23 (1969) 426-429; see also Smart et. al., In Situ Voltammetric Membrane Ozone Electrode, Anal. Chem. Vol. 51, No. 14, December &#39;79, p. 2315-2319, and Non-Steady State Measurement, &#34;The Measurement of Dissolved Oxygen&#34;, 1979. 
     In Lilley, the polarizing potential is applied as a pulse rather than a steady level. The pulse perturbs the sensor from equilibrium generating a large current transient which can be related to the oxygen partial pressure of the sample. Lilley recognized that this process had the important advantage over steady state measurements of increasing sensitivity and eliminating the requirement for stirring the sample solution as long as the current was read within four seconds of pulse application. 
     In spite of this teaching by Lilley in 1969, to date it is not believed that commercially available oxygen sensors embody the teachings of Lilley, but rather, that they use steady state measurements which measurements are subject to drift and the solution in which the electrode is immersed must be continually stirred to refresh the solution at the electrode surface. 
     The present invention provides a method and apparatus which electrochemically measures the quantity of a species in a fluid stream. More particularly, my invention provides a device wherein the acquisition, storage and generation of information corresponding to the oxygen level in a fluid stream is handled in a control module designed especially for that purpose. 
     Broadly, the invention comprises a sensor immersed in a fluid stream, the sensor specific for the species in the fluid stream to be measured, a switch to energize the sensor for a duration (d) and to prevent energization of the sensor for an interval (i); a circuit to real time average the signal from the sensor; an analog acquisition board and a computer adapted to control the activation of the switch, data acquisition and storage and all timing functions. In a preferred embodiment, temperature measurements are also derived from the sensor. 
     My invention achieves major economies in hardware using a microcomputer to control all important process parameters such as pulse duration, delay and repetition rate. Real time signal averaging is used to reduce the noise level and thereby improve precision. My device, when used for environmental sampling, does not require temperature salinity or pressure compensation circuitry. 
     The invention is particularly suitable for applications requiring long term measurements in environments where it is impracticable to service and recalibrate the sensor frequently for example, industrial processes, such as fermentors, secondary sewage treatment plants and oceanographic research. The major advantages of my invention are increased precision (changes as as small as 0.3 umolar can be detected, one time calibration) and the elimination of a stirring or flow rate requirement. The sensor is interfaced to the microprocessor, a complete record of time, temperature and oxygen concentration may be stored and made available for printout, plotting and statistical analysis. 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS 
     FIG. 1 is a block diagram of a system embodying the invention; 
     FIG. 2 is a schematic of the analog acquisition board and sensor circuit; and 
     FIG. 3 is a timing diagram. 
    
    
     DESCRIPTION OF THE PREFERRED EMBODIMENT 
     Hardware 
     Referring to FIG. 1, a block diagram of an apparatus embodying the invention is shown and comprises a voltage source 10, which consists of a regulated power supply for a polarizing voltage; an analog switch 12 under computer control to generate a pulse and an electrode 14-15 of the Clark membrane type which transmits its output to a multiplexor 16. The electrode used is a YSI 5720 with a built-in thermistor. The output from the multiplexor 16 flows to a sample and hold circuit 18 and then to an A/D converter 20. A programmable interval timer 24 (as will be described) controls the switch 12, the sample and hold circuit and the A/D converter. A computer 22, which in the preferred embodiment is a TECHNICO 9990, interfaces with the A/D converter 20, the timer 24, and a keyboard printer 26 such as a TI ASR 33. 
     The Computer 
     The control of the computer is accomplished through instructions. The instructions are written to produce the desired sequence of operations and events. The computer thus has stored in its memory the programs or routines corresponding to each mode or operation desired of the computer. It is well known to those skilled in the art, that the computer comprises suitable controls, storage and computational units for performing the various arithmetic and logical functions on data which it processes in digital form. Any standard computer language consistent with the capability of the computer can be used for the instructions. All subroutines are not described in detail, since they can be written in any desired notations, for mats or, depending upon the particular computer being utilized, computer language, (BASIC, FORTRAN) etc. Programs and instructions described are put in terms of structural flow. Where necessary and applicable for purposes of the invention, individual programs are described. For the specific computer of the preferred embodiment, the manufacturer&#39;s handbook sets forth the necessary preparations and steps for the loading and running of a program. 
     As is well known, programs are loaded into the microcomputer and for illustration purposes, the programs used herein are identified as 02INT, 02ELEC and 02 READ. After the programs have been loaded, the measurement process as controlled by 02ELEC is completely automated. The oxygen concentration and temperature are sampled under Interrupt control in background; during the data acquisition the computer can be put to any other use which the user may desire. 
     At any time during or after data acquisition, the user can examine the data which has already been collected by running the program &#34;02READ&#34;. This program takes the raw sensor readings and converts them to temperature in degrees Centigrade and oxygen concentration in micromoles per liter using the current calibration function stored on the diskette. When this program is run, even while data acquisition is underway, the time, temperature, oxygen concentration, and current rate of change of oxygen in inches per liter per hour will be printed out in real time as the measurements are made. The following is an example of the data generated. 
     
                                           TABLE I__________________________________________________________________________PULSE DURATION 1.49727 02READ                        sample outputPULSE INTERVAL 300SALINITY 32.8START PRINTOUT AT RECORD NO. 999    CALIBRATION DATA    = 9.371 - 2469.240PROBE 1  TEMP = 146.611 - 22.249THERMISTOR    TIME TEMP             02μM                  RATE μM/hr.sup.-1                           MEAN__________________________________________________________________________27       12.033         15.2             227.9                   0.00    63928       12.117         15.2             227.9                   0.00    63929       12.200         15.2             227.9                   0.00    63930       12.283         15.2             227.9                   0.00    63931       12.367         15.2             227.5                  -0.94    63932       12.450         15.2             227.5                  -1.41    63933       12.533         15.2             227.5                  -1.41    63934       12.617         15.2             227.5                  -0.94    63935       12.700         15.2             227.5                  -0.01    63936       12.783         15.2             227.5                  -0.01    63937       12.867         15.2             227.5                  -0.01    63938       12.950         15.2             227.5                  -0.01    63939       13.033         15.2             227.5                  -0.01    63940       13.117         15.2             227.5                  -0.01    639REINITIALIZE DATA PRINTER__________________________________________________________________________ 
    
     Programs and instructions described below are put in terms of structural flow: 
     02INT Structural Flow 
     1. Start 
     2. Input time (t)--time of day 
     3. Input duration (d)--pulse duration 
     4. Input interval (i)--interval between pulses 
     5. Timer 2 (24) interrupt in i seconds 
     6. R9=the first position of entry in the data table established in the memory 
     7. Stop 
     A graphic illustration of the data table is set forth below: 
     
         ______________________________________DATA TABLE             O.sub.2                    Temp.______________________________________R9            0         0        1R9 + 2        2         2        3R9 + 2        4         7        5R9 + 2        n         n        n + 1______________________________________ 
    
     02ELEC Structural Flow 
     1. Interrupt request from step 4 of 02INT 
     2. Activate timer 2 of Programmable timer (24) to output a low going pulse on output line OUT2 in d seconds 
     3. Output CRU 9  high to close switch 12 thereby energizing the electrode circuit 
     4. Has OUT2 gone low yet? If no then hold. If yes: 
     5. Prepare to sample from electrode amplifier 18 by selecting channel 2 of multiplexer 16. 
     6. Execute BURST subprogram. This will sample the input repetitively at a specified rate until a specified no. of points have been acquired. 
     (a) Specify no. of pts=111 
     (b) Specify rate as 3.3 KHz (0.3 milliseconds per point) therefore total sampling duration is 33.3 milliseconds which corresponds to 2 cycles of the major interfering noise frequency 60 Hz. 
     7. Average the 111 pts and store the result at the location pointed to by R9. (Real time averaging of the signal over 2 cycles of the 60 Hz noise very effectively filters out this source of noise. 
     8. Increment R9 by 1 to point to the next free location. 
     9. Output CRU 9  low to open switch 12 thereby de-energizing the electrode circuit. 
     10. Execute SCAN subprogram. This routine takes a single reading of the signal at a specified channel of multiplexor 16. 
     (a) Specify channel 3 which monitors the temperature sensor 15. 
     11. Store the raw temperature sensor reading at the address pointed to by R9. 
     12. Is the data table full? If yes, Stop. If no: 
     13. Increment R9 by 1. 
     14. Activate timer 2 of programmable timer 24 to generate another interrupt request in i seconds. 
     15. Acknowledge interrupt request serviced by outputting a low going pulse on CRU 6  which will clear the interrupt request IRQ6 of FF-1. 
     16. Return from interrupt. 
     02READ Structural Flow 
     1. Start 
     2. Input the answers to the following questions: 
     (a) Salinity? 
     (b) Start printout at record no. ? (1-240) if 999 then then start printing the most current temperature and oxygen measurements. 
     3. Starting at the specified record no., calculate O 2  concentration and temperature from the raw sensor readings stored in the data table. These calculations are: 
     Temp ° C=146.611-22.249 log e  (raw temp sensor reading) 
     From this in-situ temperature various quantities are calculated for the purpose of calculating temperature and salinity correction factors for converting the raw oxygen sensor readings to the desired concentration units. These are: 
     
         ______________________________________KELVIN = TEMP + 273.15RT = 100./KELVINT100 = KELVIN/100.*P1 = -173.249     *These constants areP2 = 249.634       taken from a paper byP3 = 143.348       Murray &amp; Riley 1969P4 = -21.8492      Deep Sea Res. 16,311.P5 = -3.3096 × 10.sup.-2P6 = 1.4259 × 10.sup.-2P7 = -1.7 × 10.sup.-3A1 = P1 + P2*RTA2 = P3*log.sub.e (T100) + P4*T100S1 = P5 + P6*T100 + P7*T100*T100F = EXP(A1 + A2 + Salinity &amp; S1)/EXP(A1 + A2)SLOPE = (1./EXP(9.371-2469.240*(1./KELVIN)*FB1 = -15*SLOPEO.sub.2 concentration μM = B1 + SLOPE*(raw O.sub.2 sensor reading) (Note: *equals times) 
    
     4. Print out data and also store in a diskette data file. 
     
         ______________________________________Time     temp °C.               O.sub.2 μM                         Rate (μM hr.sup.-1)______________________________________ 
    
      continue until done or the BREAK is typed. 
     5. Enter menu of possible operations 
     1. 02INT restarts data acquisition 
     2. Plot of O 2  concentration vs. time 
     3. Compute simple statistics, means, standard deviations, 95% confidence intervals, slope, intercept and correlation coefficient 
     4. Update calibration data 
     5. Print data 
     OPERATION 
     In the operation of the invention, an electrode 14-15 such as a Clark type membrane with thermistor is used in the preferred embodiment. The electrode 14-15 is inserted into the body of water, the oxygen level of which is to be measured. Whether inserted into a tank, bottle or the like is immaterial. For the data shown in Table I above, the water to be sampled was siphoned into a BOD bottle. The sensor, which had a taper the same as the ground glass neck of the bottle, was inserted carefully into the bottle so as not to trap any bubbles. 
     The following describes the operation of the invention with reference to the drawings. 
     The program 02INT is loaded into the computer 20 and the current time and pulse duration (d) and pulse interval (i) desired is input into the computer 22. In the specific example, the pulse duration was 1.5 seconds and the pulse interval was 5 minutes. The pulse duration could be any range depending upon the species in the fluid being measured and the characteristics of the membrane, including the recovery rate. For oxygen, a pulse duration from 0.5-3 seconds is suitable and a pulse interval of from 3-5 minutes or greater than five minutes is suitable. 
     Referring to FIGS. 2 and 3, the computer data acquisition and process control is accomplished using two INTEL 8255 programmable peripheral devices. As shown in FIG. 2, these devices are organized as three READ/WRITE registers and a WRITE ONLY control register. To the computer 22, these registers appear as bytes of internal memory. Port B is programmed for strobed input and Port C is programmed for simple output. Ten Port B input lines 8 bits of 8255 lo and 2 bits of 8255 hi are connected to the output lines of the AD 571 analog to digital converter 20. A data ready pulse originating from the AD 571 latches the data from the AD 571 into the Port B input after being delayed 600 nanoseconds by monostable 25 to allow time for the data lines to settle. The falling edge of the data ready pulse also sets a bit in Port C of 8255 which doubles as a status register signifying that the input buffer is full. Software testing to see whether this bit is set can tell if new data is available in the input buffer. Reading the input, Port B, automatically clears the input buffer full flag. 
     The input of the analog to digital converter 20 is driven by the sample and hold amplifier 18. The sample and hold amplifier 18 tracks the input waveform from the multiplexer (16) until a convert command is issued to the sample and hold circuit via the OUT 0 line of programmable timer 24. At this time, the input voltage is held constant for 25u seconds, the time it takes the converter 20 to make the analog to digital conversion. 
     The source of the input voltage to be sampled is selected by the multiplexer 16. The particular channel selected is determined by the data put out on the multiplexer address lines by bytes PC 5-7 of output Port C. 
     Channel 2 of the multiplexer 16 is connected to the oxygen sensor 14 and channel 3 is connected to the thermistor 15. The six other channels are available to monitor any other pertinent sensors as desired. 
     All time dependent data acquisition events are controlled by the INTEL 8253 programmable counter timer 24. The timer contains three completely independent timers, each operating on a different time base; two of which are used in this embodiment. The input of Timer 0 is connected to a 1.5 Mhz time base and used to control the sampling rate of the analog to digital converter 20. The input of Timer 2 is connected to a 91.5 Hz time base and Timer 2 is used to generate interrupts after a 0.011-360 second delay. A pulse on the OUT2 line is generated when Timer 2 counts down to zero which toggles the output of a flip-flop FF-1 to low sending an interrupt request via IRQ 6  to the computer 22. If interrupts have been enabled, the 02ELEC program pointed to by the interrupt vector stored at location 18 will be executed, upon completion of the 02ELEC program the interrupt request is cleared by sending pulse out on the interrupt acknowledge line CRU 6 . 
     Generation of a signal proportional to oxygen concentration begins with a well regulated power supply 10 which generates the polarizing voltage for the oxygen electrode 14. The analog switch 12 is in series with the oxygen electrode 14. The switch is opened and closed by computer control via CRU 9  so as to produce a pulse of precisely controlled duration (d). The height of the pulse is controlled by the 10 K variable resistor and is set to put the electrode in the plateau region of its I-V curve. The preferred pulse height has been found to be 1.2 volts for a duration of 1.5 seconds, which produces optimal linerality and long term stability. The current transient produced when the electrode 14 is energized is converted to a proportional voltage by the operational amplifier 17. A second stage of amplification at 18 inverts the signal so that the voltage transient 0-10 V is positive. The eleven volt Zener diode clips signals which exceed 11 volts protecting the subsequent stages. The signal is input into the multiplexer channel 2. 
     A voltage divider circuit 15 converts the resistance changes in the thermistor 15 into a proportional voltage signal which is fed to channel 3 of the analog multiplexer 16. See flow chart. 
     FIG. 3 illustrates the timing sequence used in the preferred embodiment of my invention. 
     Although my invention has been described with reference to the measurement of oxygen, any species in a fluid stream, gaseous or ionic, may be measured.