CN106053438B - A kind of synthetic biological toxicity in water remote auto analyzer - Google Patents
A kind of synthetic biological toxicity in water remote auto analyzer Download PDFInfo
- Publication number
- CN106053438B CN106053438B CN201610301333.7A CN201610301333A CN106053438B CN 106053438 B CN106053438 B CN 106053438B CN 201610301333 A CN201610301333 A CN 201610301333A CN 106053438 B CN106053438 B CN 106053438B
- Authority
- CN
- China
- Prior art keywords
- data
- curve
- coefficient
- luminous intensity
- value
- 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
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 37
- 231100000419 toxicity Toxicity 0.000 title claims abstract description 32
- 230000001988 toxicity Effects 0.000 title claims abstract description 32
- 238000012545 processing Methods 0.000 claims abstract description 16
- 239000010865 sewage Substances 0.000 claims abstract description 13
- 238000012544 monitoring process Methods 0.000 claims abstract description 10
- 241000193830 Bacillus <bacterium> Species 0.000 claims abstract description 5
- 239000000463 material Substances 0.000 claims abstract description 5
- 230000001473 noxious effect Effects 0.000 claims abstract description 5
- 230000005540 biological transmission Effects 0.000 claims abstract description 4
- 238000005415 bioluminescence Methods 0.000 claims abstract description 4
- 230000029918 bioluminescence Effects 0.000 claims abstract description 4
- 230000005764 inhibitory process Effects 0.000 claims abstract description 4
- 239000012488 sample solution Substances 0.000 claims abstract description 4
- NRZZLYODXDSLEK-UHFFFAOYSA-N (6-ethoxy-6-oxohexyl) 3,5-diacetamido-2,4,6-triiodobenzoate Chemical compound CCOC(=O)CCCCCOC(=O)C1=C(I)C(NC(C)=O)=C(I)C(NC(C)=O)=C1I NRZZLYODXDSLEK-UHFFFAOYSA-N 0.000 claims description 27
- 238000000034 method Methods 0.000 claims description 27
- 230000000694 effects Effects 0.000 claims description 10
- 231100000167 toxic agent Toxicity 0.000 claims description 9
- 239000003440 toxic substance Substances 0.000 claims description 9
- 230000006872 improvement Effects 0.000 claims description 7
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 claims description 6
- 238000003860 storage Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 239000000203 mixture Substances 0.000 claims description 4
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000003190 augmentative effect Effects 0.000 claims description 3
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 abstract description 11
- 238000005516 engineering process Methods 0.000 abstract description 7
- 238000012546 transfer Methods 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 45
- 239000007788 liquid Substances 0.000 description 17
- 230000006870 function Effects 0.000 description 15
- 241000894006 Bacteria Species 0.000 description 14
- 239000000243 solution Substances 0.000 description 12
- 230000004913 activation Effects 0.000 description 8
- 238000004458 analytical method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 230000007812 deficiency Effects 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 5
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 239000000523 sample Substances 0.000 description 4
- 230000007059 acute toxicity Effects 0.000 description 3
- 231100000403 acute toxicity Toxicity 0.000 description 3
- 239000007853 buffer solution Substances 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000002347 injection Methods 0.000 description 3
- 239000007924 injection Substances 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 239000002351 wastewater Substances 0.000 description 3
- 241000607565 Photobacterium phosphoreum Species 0.000 description 2
- 230000003139 buffering effect Effects 0.000 description 2
- 238000004140 cleaning Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 239000004615 ingredient Substances 0.000 description 2
- 230000013011 mating Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000005622 photoelectricity Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 231100000331 toxic Toxicity 0.000 description 2
- 230000002588 toxic effect Effects 0.000 description 2
- 241000589614 Pseudomonas stutzeri Species 0.000 description 1
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000013019 agitation Methods 0.000 description 1
- 238000007630 basic procedure Methods 0.000 description 1
- 239000012267 brine Substances 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 230000008014 freezing Effects 0.000 description 1
- 238000007710 freezing Methods 0.000 description 1
- 238000011534 incubation Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- 230000015654 memory Effects 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 239000002574 poison Substances 0.000 description 1
- 231100000614 poison Toxicity 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 239000000843 powder Substances 0.000 description 1
- 230000005180 public health Effects 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
- 239000010703 silicon Substances 0.000 description 1
- HPALAKNZSZLMCH-UHFFFAOYSA-M sodium;chloride;hydrate Chemical compound O.[Na+].[Cl-] HPALAKNZSZLMCH-UHFFFAOYSA-M 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 239000008400 supply water Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 231100000820 toxicity test Toxicity 0.000 description 1
- 231100000027 toxicology Toxicity 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- 238000004457 water analysis Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/75—Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated
- G01N21/76—Chemiluminescence; Bioluminescence
- G01N21/763—Bioluminescence
Landscapes
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Plasma & Fusion (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
The present invention provides a kind of synthetic biological toxicity in water remote auto analyzer, including synthetic biological toxicity in water analytical equipment, data sending device and data receiver and display device, inhibition level of the synthetic biological toxicity in water analytical equipment for monitoring the noxious material in sample solution to bioluminescence bacillus luminous intensity；Data sending device is used to wirelessly for luminous intensity information to be sent to setting data receiver at the far end and display device, data receiver and display device determine the class of pollution of sewage toxicity after carrying out data processing to luminous intensity information, and result information is shown.Compared with the existing technology, the present invention is transmitted using wireless remote data, transmission circuit failure problems when avoiding previous wired data transfer circuit, it is easy to maintain, simultaneously remotely the data of acquisition can be analyzed to obtain the situation of live water quality, more severe remote auto detection is can be suitably used for, can be used as single machine, can also act as the partial node of detection network system.
Description
Technical field
The present invention relates to water quality monitoring technical fields more particularly to a kind of synthetic biological toxicity in water remote auto to analyze
Instrument.
Background technology
In order to adapt to the needs of environmental monitoring, solve cumbersome, time-consuming, costly and repeated existing for biological toxicity tests
The problem of with sensitivity etc..The country is continuously improving always the test method of bio-toxicity, keeps its quick, sensitive, economical
It is easy, from middle nineteen sixties begin one's study with photobacteria monitoring of environmental pollute, in the seventies method become comparative maturity
Measuring technology.More with scholar's research application of the country such as the U.S., Britain, the Soviet Union, Japan, Australia, some projects are
Through having applied for patent, since the eighties the Beckman in the U.S. together the research of company in this respect achieve it is breakthrough into
Exhibition.They have developed special test equipment and photobacteria experiment is mating, and the life of noxious material is measured in application photobacteria
Many successful researchs are made in terms of object toxicity.So that the popularization that this method is developed quickly.In recent years, Chinese scholar is to shining
The research of bacterium and application made extensive work, and the special test equipment of development and production establishes Luminous bacteria biology poison
Property test system (LB systems), monitor popularization and application in China for Photobacterium Phosphoreum Toxicity and lay the foundation.
But there are following some problems for the above-mentioned previous test equipment made using photobacteria：As Bake Mann gives birth to
The special test equipment of production and photobacteria experiment are mating, but belong in experiment and use, at high price and to the requirement of environment
It is excessively harsh.The on-line automatic analysis of North America instrument production is COD and BOD methods using dividing, and also only rests on laboratory use
Bottom on.Domestic biology department of East China Normal University and Nanjing Transistor Factory once produced the tester made using photobacteria
Device, but it is to belong to manual sampling, manually reads data, intermittence test mode, it is non-it is online can only be reality together with
Freezing equipment work in combination.Public Health College, Harbin Medical Univ Hygiene Toxicology teaching and research room once reports production acute toxicity
Photobacteria sensor and acute toxicity high speed tester, a kind of urgency with photobacterium phosphoreum and silicon photo diode composition of design
Property the toxic bacteria luminescence sensor and quick determination of the environment Acute Toxicity of Pollutants analyzer of flow type Microcomputer, it is fixed to realize
Property test purpose, can not achieve on-line testing function and reduced concentration analysis experiment judge the toxic grade of sewage.Shenyang seven
The CWDZF-1 type municipal sewage bio-toxicity on-line automatic analyzer (Patent No. of color science and technology engineering company development and production
2004100213555) have many advantages, such as live on-line testing and to anti-adverse environment, it is comprehensive particularly suitable for municipal sewage biology
Close toxicity detection, but it there are following deficiency and disadvantages：
1, CWDZF-1 types in-line analyzer detects required shine each time due to not activating device online
Bacterium is re-introduced into instrument after prior activation is good, because the service life of bacterium is limited, MaLS is 117 hours, so
The time longest that follow-on test does not need human intervention is exactly 5 days, shows photogen according to photogen life curve result of detection
Best fluorescent lifetime is that activation starts 20 hours to 36 hour this periods, and CWDZF-1 type analysis instrument does not solve best raw
The problem of life period and extend as far as possible continuously monitor.
2, CWDZF-1 types in-line analyzer pipeline uses medical catheter and containers for medical use bottle, real using the principle for hanging brine
The injection of existing photogen, it may appear that the problem of occurring other bacteriums inside pipeline and striving resource with photogen.To influence to examine
Survey the accuracy of result.
3, CWDZF-1 types in-line analyzer is original only there are one warm area, cannot ensure that the activity of photogen is in most very well
Good state influences to detect quality so the luminous intensity of the photogen in test process is insufficient.
4, the original test cell of CWDZF-1 types in-line analyzer is made using the click structure for moving up and down photomultiplier
It obtains test cell rotation and moves up and down common presence with photomultiplier, the testing time is long, the number that each test tube is measured
According to few.
5, after the test cell rotation of CWDZF-1 types in-line analyzer, into test tube, injection liquid is not easy accurately, especially
When the non-head of photomultiplier is not directly corresponding with test tube, since what photomultiplier was clicked moves up and down, photomultiplier transit can be damaged
Pipe.Influence the reliability of detection.
6, CWDZF-1 types in-line analyzer can not achieve long-range detection and analysis, it is necessary to data could be obtained at the scene,
Especially bad weather when, staff has to stay at acquisition and work, and acquisition is inconvenient, cannot receive long-range and
Analyze data,.
7, CWDZF-1 types in-line analyzer can not realize that multinode acquires, and can not realize the networking of data monitoring.
Therefore for drawbacks described above present in currently available technology, it is really necessary to be studied, to provide a kind of scheme,
Solve defect existing in the prior art.
Invention content
In view of this, the purpose of the present invention is to provide a kind of accuracy height, stability is good, reliability is high, can be long-range
The synthetic biological toxicity in water remote auto analyzer of the monitoring of real-time, to solve the above problems.
In order to overcome the deficiencies of existing technologies, technical scheme is as follows：
A kind of synthetic biological toxicity in water remote auto analyzer, including synthetic biological toxicity in water analytical equipment, data
Sending device and data receiver and display device, wherein synthetic biological toxicity in water analytical equipment is for monitoring in sample solution
Noxious material to the inhibition level of bioluminescence bacillus luminous intensity；The data sending device is used for luminous intensity information
Setting data receiver at the far end and display device, the data receiver and display device are wirelessly sent to shining
Strength information determines the class of pollution of sewage toxicity after carrying out data processing, and result information is shown；
The synthetic biological toxicity in water analytical equipment include Instrument shell, the cultivation unit of photogen, photogen shine
Intensity collection and storage unit, photogen luminous intensity transmission unit, data acceptance unit and data display unit.
Preferably, the data receiver and display device determine sewage toxicity after carrying out data processing to luminous intensity information
The class of pollution include the following steps：
Step 1：Data processing is carried out to luminous intensity information by the curve fitting algorithm based on numerical value phase approximately principle, i.e.,
The parametric equation of curve matching is obtained using least square solution：X=[ATA]-1ATB obtains model y1, and exact value is found out,
Middle X is least square method parameter, and A, B are fitting coefficient；
Step 2：Data processing is carried out to luminous intensity information by the curve fitting algorithm based on form principle of similarity：
(1) time-concerning impact factor is set as λ, reacts amount of the photogen activity by the time-to-live；If the similar matched curve of form
For：y2(x)=a0+a1x+a2x2+...+anxn, wherein aiFor each level number；
(2) by matched curve combination shape coefficient equationWherein x12iTo be fitted point value, S12For deviator
Value, can obtain：
(3) augmented objective function with parameter is constructed, when meeting constraints, equal sign both sides are set up；When being unsatisfactory for
When constraints range, a fully big number μ is taken>0, construct following function：
A={ a in formula1,a2,...,an,λ}；
(4) pass through the solved function extreme value of powell algorithms：
Step a：Selected initial point x(0), n linearly independent vector group, the first search direction system { p of composition0,p1,....pn -1, given accuracy ε ＞ 0 set k=0, and k is initial extreme value；
Step b：Enable y0=xk, successively along { p0,p1,....pn-1In direction carry out a pile search,
It is corresponding to obtain auxiliary iteration point y1,y2,....yn, i.e.,
β in formulaj-1For along pj-1The step-length in direction；
Such as | y(k)-x(k)|<ε is set up, then stops calculating, no to then follow the steps c；
Step c：Acceleration direction is constructed, p is enabledn=yn-y0If | | pn| |≤ε then stops iteration, exports xk+1=yn, otherwise
Go to step d；
Step d：Determine adjustment direction：Find out m, wherein m is minimum point so that
f(ym-1)-f(ym)=max | f (ym-1)-f(ym)|1≤j≤n}
If following formula is set up：
f(y0)-2f(yn)+f(2yn-y0) [f (the y of ＜ 2m-1)-f(ym)],
Step f is gone to, e is otherwise gone to step；
Step e：Enable xk+1=yn+βnpn,Meanwhile enabling { p0,p1,....pn-1}k+1=
{p0,...,pm-1,pm+1,....pn-1,pnK=k+1 goes to step b；
Step f：Enable xk+1=yn, set k=k+1 and go to step b；
(5) by powell algorithms, a is found out1,a2,......,an, it is contemplated that a0For DC component, fitting is not interfered with
The form of curve substitutes into formula y2(x)=a0+a1x+a2x2+...+anxnIn, a is asked using least square method0, obtain being based on form phase
Like the curve fitting algorithm y of criterion2；Accuracy is found out using formula of correlation coefficient：
Step 3：From the two above steps, cubic fitting model y is obtained respectively1, y2, construct improved curve matching mould
Type is：Y (x)=w1y1(x)+w2y2(x), wherein 0 ＜ w1≤ 1,0 ＜ w2≤ 1, w1、w2According to the precision value found out；
Step 4：Further optimize each term coefficient found out using Fabonacci methods, it is made to be more nearly optimal value；
It is as follows：
1) set the deviation of respective value in j-th of discrete point and matched curve as
It can be found out according to above formula, the maximum positively biased in n discrete point is not good enoughIt is not good enough with maximum negative biasThen m powers coefficient amInitial section [A1,B1] be
2) first judge aiIn initial section [A1,B1] it is unimodal function, the best a found outTIt is exactly aiIn section [A1,B1] in
Approximate minimum or maximum, i.e. aiIn [A1,aT] section answer strictly decreasing or be incremented by, in [aT,B1] on answer strictly increasing or
Successively decrease, it is optimized using Fabonacci methods；
3) according to the positive and negative inclined absolute value of the difference approximately equal of maximum of optimum fit curve, then deviation takes sum of the two
Half, i.e.,
It can thus be concluded that best constant coefficient is
Step 5：Utilize the expert system to prestore in obtained improvement matched curve y combination this system databases, march
Lines matching respectively obtains the fitting coefficient B for improving matched curve and photogen in expert system and variety classes, the toxicity of concentration
The fitting coefficient A of each reaction mechanism curve of substance seeks the distance between fitting coefficient B and each fitting coefficient A, coefficient minimum
By the toxicant to be predicted.
Compared with the existing technology, the present invention is transmitted using wireless remote data, when avoiding previous wired data transfer circuit
Transmission circuit failure problems, it is easy to maintain, while remotely the data of acquisition can be analyzed to obtain the situation of live water quality,
More severe remote auto detection is can be suitably used for, can be used as single machine, can also act as the partial node of detection network system.
Data processing advanced optimizes improved curve matching mould by being implanted into improved curve fitting algorithm with Fabonacci methods
Each term coefficient of type realizes the prediction of toxicant ingredient and concentration；Simultaneously present invention further contemplates that magnetic agitation time, culture temperature
Under the action of degree, incubation time, exposure duration and pH value range are to result measured value, each condition is analyzed to photogen bio-toxicity
The influence of experiment realizes remote online monitoring to improve while analyzer measurement accuracy, stability.
Description of the drawings
Fig. 1 is the general structure schematic diagram of the synthetic biological toxicity in water remote auto analyzer of the present invention；
Fig. 2 is the structural schematic diagram of synthetic biological toxicity in water analytical equipment in the present invention；
Fig. 3 is the structure diagram of wireless base station apparatus in the present invention；
Fig. 4 is the structure diagram of data receiver and display device in the present invention；
Comparisons of the Fig. 5 between follow-on cubic fitting model.
Specific implementation mode
Following is a specific embodiment of the present invention in conjunction with the accompanying drawings, technical scheme of the present invention will be further described,
However, the present invention is not limited to these examples.
Referring to Fig. 1, it is shown the general structure schematic diagram of the synthetic biological toxicity in water remote auto analyzer of the present invention,
Including synthetic biological toxicity in water analytical equipment, data sending device and data receiver and display device, wherein water quality integrates
The inhibition level that bio-toxicity analytical equipment is used to monitor the noxious material in sample solution to bioluminescence bacillus luminous intensity；Institute
Data sending device is stated for luminous intensity information to be wirelessly sent to setting data receiver at the far end and display
Device, the data receiver and display device carry out luminous intensity information the pollution etc. of determining sewage toxicity after data processing
Grade, and result information is shown.
Referring to Fig. 2, it is shown the synthetic biological toxicity in water analytical equipment of the present invention, test philosophy is as follows:
System changes the position of test tube by controlling the rotation of motor 13.When test tube disk 57 is in the initial position time
Electric Hall switch 7 exports a low level signal, and it is high level signal otherwise to export.Before test starts every time and test
When, system can all detect the output signal of photoelectricity Hall switch 7.It is calculated using rational positioning if signal is high level
Method rotary electric machine 13 to position for test tube disk, is returned to initial position.Energy is used in combination in photoelectricity Hall switch 7 and stepper motor
Positioning test tube disk 57 is to realize the positioning of test tube well.Photomultiplier 8 is located at the side of No. 4 test tubes 9, it can be acquired
The optical signal that photogen sends out, and faint electric signal is converted optical signals to for control system acquisition process.In test tube disk
57 when being in initial position, and the alignment of No. 1 test tube 61 adds 56, No. 2 test tubes 59 of liquid pipe to be measured and is directed at 55, No. 3 test tubes 10 of water pipe
Alignment buffering 11, No. 4 test tubes 9 of liquid pipe are directed at bacterium solution pipe 12, while No. 4 test tubes are also targeted by photomultiplier 8.Add bacterium solution pipe
12 are connected with photogen storage tank 17 after activation, and system adds bacterium solution flow control electromagnetic valve to be added into test tube by opening and closing
Quantitative bacterium solution.Buffering liquid pipe 11 is connected with buffer solution storage tank 36, is added into test tube by control flow solenoid valve quantitative
Buffer solution.Water pipe 55 is connected with water pump pipe or running water pipe to supply water at any time, and system controls solenoid valve by switch traffic
To add clear water to test tube, for rinsing test tube.Liquid pipe 56 to be measured is connected with fluid cylinder 51 to be measured, and system can pass through flow control
Solenoid valve is opened and closed to add quantitative prepare liquid into test tube.Fluid cylinder to be measured is connected with water pump, and water pump can be always to prepare liquid
It is newest that prepare liquid is sent in cylinder to keep the prepare liquid in fluid cylinder to be measured.There are overflow port and overflow pipe in the upside of fluid cylinder 51 to be measured
53 are connected, and prepare liquid can flow to other places to prevent the prepare liquid of prepare liquid cylinder excessive by overflow pipe 53.Delivery pipe 2,
3,64,65 have wastewater trough 1 below, test the liquid of the waste liquid in test tube and cleaning test tube and will all be discharged in wastewater trough 1.
Analyzer internal upper part is provided with the containers such as freeze-dried powder storage container, buffer solution kettle, nutrient tank, for storing activation
Required raw materials and reagents when photogen.Bacterium solution stores kettle and sewage cylinder after being provided with activation bacterium kettle, activation in the middle part of instrument, respectively
For bacterium liquid activation, the storage of test bacterium solution and testing liquid.Instrument lower part is provided with oxicity analysis test device, main to wrap
Include the two phase mixing stepper motor for precision positioning, sample container disk, photoelectric detection system and test container cleaning device
Deng.It is connected by various pipelines between entire instrument, the addition of liquid and the control of additive amount are completed by flow solenoid valve to control.
It is constant at 20 DEG C to start sample to be tested temperature in thermoelectric refrigerator and temperature detection sensor adjustment measurement darkroom
Left and right.It using the time-sharing automatic acquisition water quality of sewage pump multi-point sampling and is uniformly mixed, passes through water conservancy diversion solenoid valve control sewage sample
Into darkroom is measured, by the luminous bacillus of auto-control ration drug injector automatic ration filling activation, to measuring darkroom
Interior sewage sample and the mixing of luminous bacillus nitrogen injection start photomultiplier progress 15min and continuously measure, and measurement data is interim
It preserves and data receiver and display device, data receiver and display device pair is wirelessly transmitted to by data sending device
Data carry out analyzing processing and non-volatile memories, with the figure and analysis result after the display processing of lattice lcd module scene,
Analyzer progress pipeline cleans waiting and measures next time automatically to be started.
Referring to Fig. 3, it show as the structure diagram of wireless base station apparatus in the present invention, synthetic biological toxicity in water analysis
The data of collected reflection water quality situation are sent to data and send control module 46 by device, and control module 46 sends out sensor
The data sent are handled, and control water quality signal lamp 45 is bright, bright more of lamp, indicate that water quality is more bad.By treated
Data are sent to WiFi receiving modules 47 by WiFi sending modules 43, signal lamp 42 indicate current data send whether work
Make normally, if working properly, signal lamp 42 is bright, and signal lamp 44 shows whether sensing data receives normal, if passing
Sensor data receiver is normal, then signal lamp 44 is bright.
Referring to Fig. 4, it is shown the structure diagram of data receiver and display device in the present invention, WiFi receiving modules 47 receive
The data that WiFi sending modules 43 are sent, whether ultraviolet lamp 48 receives normally for display data, purple if data receiver is normal
Outer lamp 48 is bright.Power switch 49 controls entire data and receiving module switch.It control section 41 will be at the data that received
The data of each node are carried out analyzing processing by reason if it is multinode, and analysis and processing result is shown on LCD display 40
It shows to come, and display system working condition.
Wherein, the step of carrying out data processing to luminous intensity information in data receiver and display device is as follows：
1, it is based on the curve fitting algorithm of numerical value close (least square method) combination similar to form, utilizes two kinds of algorithms
Precision value introduces time-concerning impact factor and powell algorithms as the weight in innovatory algorithm model.Solves numerical value phase
The closely fitting speed of (least square method) principle, time factor are avoided to photogen activity influence and using weighted value advantage
The preceding n direction of search of powell algorithms must linear independence the problems such as.The basic procedure of the innovatory algorithm is as follows：
Curve fitting algorithm based on numerical value phase approximately principle：
The parametric equation of curve matching is obtained using least square solution：X=[ATA]-1ATB obtains model y1, and find out
Exact value.
Curve fitting algorithm based on form principle of similarity：
If time-concerning impact factor is λ, amount of the reaction photogen activity by the time-to-live.If the similar matched curve of form is：
y2(x)=a0+a1x+a2x2+...+anxn。
By matched curve combination shape coefficient equationIt can obtain：
The solution of the problem is converted by solution Unconstrained Optimization Problems using the outer point method of penalty function.It constructs with parameter
Augmented objective function, when meeting constraints, equal sign both sides are set up, and when being unsatisfactory for constraints range, equation value is very
Greatly.Take a fully big number μ>0, construct following function：
A={ a in formula1,a2,...,an,λ}
The step of being the solved function extreme value of powell algorithms below：
Step a：Selected initial point x(0), n linearly independent vector group, the first search direction system { p of composition0,p1,....pn -1, given accuracy ε ＞ 0 set k=0, and k is initial extreme value；
Step b：Enable y0=xk, successively along { p0,p1,....pn-1In direction carry out a pile search,
It is corresponding to obtain auxiliary iteration point y1,y2,....yn, i.e.,
β in formulaj-1For along pj-1The step-length in direction；
Such as | y(k)-x(k)|<ε is set up, then stops calculating, no to then follow the steps c；
Step c：Acceleration direction is constructed, p is enabledn=yn-y0If | | pn‖≤ε then stops iteration, exports xk+1=yn, otherwise turn
Step d；
Step d：Determine adjustment direction：Find out m, wherein m is minimum point so that
f(ym-1)-f(ym)=max | f (ym-1)-f(ym)|1≤j≤n}
If following formula is set up：
f(y0)-2f(yn)+f(2yn-y0) [f (the y of ＜ 2m-1)-f(ym)],
Step f is gone to, e is otherwise gone to step；
Step e：Enable xk+1=yn+βnpn,Meanwhile enabling { p0,p1,....pn-1}k+1=
{p0,...,pm-1,pm+1,....pn-1,pnK=k+1 goes to step b；
Step f：Enable xk+1=yn, set k=k+1 and go to step b；
According to powell algorithms, a is found out1,a2,......,an, it is contemplated that a0For DC component, matched curve is not interfered with
Form, substitute into formula y2(x)=a0+a1x+a2x2+...+anxnIn, a is asked using least square method0.Therefore it obtains being based on form phase
Like the curve fitting algorithm y of criterion2.Accuracy is found out using formula of correlation coefficient：
From the two above steps, cubic fitting model y is obtained respectively1, y2, constructing improved cubic fitting model is：y
(x)=w1y1(x)+w2y2(x), wherein 0 ＜ w1≤ 1,0 ＜ w2≤ 1, w1w2According to the precision value found out.
2, optimization of the Fabonacci methods to improvement curve fitting algorithm
With improved curve fitting algorithm fit come effect be preferable, but be not necessarily best.Therefore, may be used
On the basis of improved curve fitting algorithm, further optimizes each term coefficient found out using Fabonacci methods, make
It is more nearly optimal value, with y=a0+a1x+a2x2+...+anxnFor, it first has to find out anInterval (original area
Between), then its interval is optimized using Fabonacci methods, determines optimum value.Under the conditions of this optimum value, ask
Go out an-1 optimum value.The optimum value of other coefficients is according to said method found out one by one.
Algorithm steps are as follows：
amInitial section [A1,B1] algorithm
If the deviation of respective value is in j-th of discrete point and matched curve
It can be found out according to above formula, the maximum positively biased in n discrete point is not good enoughIt is not good enough with maximum negative biasThen m powers coefficient amInitial section [A1,B1] be
Optimum coefficient
Fabonacci methods are adapted to unimodal function, it is therefore necessary to first judge aiIn initial section [A1,B1] it is unimodal function.
Known according to the definition of unimodal function, the best a found outTIt is exactly aiIn section [A1,B1] in minimum (big) value of approximation, i.e. ai
[A1,aT] section answers strictly decreasing (increasing), in [aT,B1] on answer strictly increasing (subtracting).Judge aiAfter unimodal function, use
Fabonacci methods optimize it.
Constant coefficient a0Algorithm
According to the positive and negative inclined absolute value of the difference approximately equal of maximum of optimum fit curve, then deviation takes the one of sum of the two
Half, i.e.,
It can thus be concluded that best constant coefficient is
3, the prediction of toxicant ingredient and concentration
To realize the prediction of toxicant type and concentration, the special of obtained improvement matched curve y combination this system is utilized
Family's system, carries out Curve Matching.Respectively obtain improve fitting coefficient B and the photogen in expert system of matched curve with it is not of the same race
The fitting coefficient A of each reaction mechanism curve of toxicant of class, concentration, seeks between fitting coefficient B and each fitting coefficient collection A
Distance, coefficient it is minimum by the toxicant to be predicted.
The above results information, user can be inquired by client.Client/server is with the side such as chart, curve
Formula is presented to the user, and foundation is provided for administrative staff's monitoring, data analysis, decision.In the client, user can only be to currently examining
What survey, toxicant library and history detected checks.After login, information that can be in automatic download server business to client
End, to realize checking for user.Application and development based on android system belongs to user front end function.Its function is main
Have：Realize detection, popular science knowledge, individual center, system configuration, user interaction, search, the binding etc. of analyzer.
Photogen different time points in waste water to be measured are collected by water quality biological toxicity in-line analyzer data
Corresponding photogen luminous intensity numerical value, and improved curve fitting algorithm is combined, graphing, as shown in Figure 5.
Here also due to x=1:1:50, span is larger, therefore military order t=(x-51)/51, the value model of independent variable
It encloses for [- 1,1], in Figure 5, the value range of abscissa is [- 1,0], therefore obtained prediction model is：
Y (t)=6641.463t10+29039.504t9+54058.737t8+55957.231t7+35297.68t6
+13988.965t5+3466.473t4+518.092t3+43.192t2+1.614t+1.243
=0.41586y1+0.5724y2
Wherein w1=0.41586, w2=0.5724 matched curve made is closest to real curve.
The curve values of fitting are with the error of corresponding actual value：
As seen from Figure 5, the matched curve of least square is flatter, and the matched curve of form similarity criterion is anti-well
Variation and the jitter conditions of real curve have been answered, but upper deficiency is fitted in the precision of numerical value, and improved matched curve is fitted
Effect is best.Innovatory algorithm is the combination of two kinds of curve matching thought, and algorithm idea is to maintain to match with real curve form
While numerical value it is also close, this algorithm idea has more advantage in more complicated curve, and embodies the life of photogen
Activity is influenced by time-concerning impact factor.It can be seen that according to figure, the matched curve of least square is flatter, form similarity criterion
Variation and the jitter conditions of real curve have been reacted in matched curve well, but are fitted upper deficiency in the precision of numerical value, and are changed
Into matched curve fitting effect it is best.Innovatory algorithm is the combination of two kinds of curve matching thought, algorithm idea be to maintain with very
Numerical value is also close while solid-line curve form matches, and this algorithm idea has more advantage, and body in more complicated curve
Reveal the vital activity of photogen is influenced by time-concerning impact factor.
Due to powell algorithms in iteration before n direction of search must linear independence, otherwise asking not optimal solution
Topic, therefore powell algorithms fail, and can effectively avoid the deficiency of powell algorithms using weighted value in innovatory algorithm.As y (x)
=w1y1(x)+w2y2(x) w in2When very little, improved curve fitting algorithm will be least square method y (x)=w at this time1y1(x)+k
(wherein k is the numerical value of very little), avoids deficiency of the powell algorithms without optimal solution in innovatory algorithm.
In improvement curve algorithm, the weight of two formula is determined using accuracy, utilizes the calculating of powell algorithms
Fireballing feature improves the fitting speed of innovatory algorithm.
The explanation of above example is only intended to facilitate the understanding of the method and its core concept of the invention.It should be pointed out that pair
For those skilled in the art, without departing from the principle of the present invention, the present invention can also be carried out
Some improvements and modifications, these improvement and modification are also fallen within the protection scope of the claims of the present invention.To these embodiments
A variety of modifications are it will be apparent that General Principle defined herein can be for those skilled in the art
It is realized in other embodiments in the case of not departing from the spirit or scope of the present invention.Therefore, the present invention is not intended to be limited to
These embodiments shown in the application, and be to fit to consistent with principle disclosed in the present application and features of novelty widest
Range.
Claims (1)
1. a kind of synthetic biological toxicity in water remote auto analyzer, which is characterized in that analyzed including synthetic biological toxicity in water
Device, data sending device and data receiver and display device, wherein synthetic biological toxicity in water analytical equipment is for monitoring
The inhibition level of noxious material in sample solution to bioluminescence bacillus luminous intensity；The data sending device will be for that will shine
Strength information is wirelessly sent to setting data receiver at the far end and display device, and the data receiver and display fill
It sets and determines the class of pollution of sewage toxicity after carrying out data processing to luminous intensity information, and result information is shown；
The synthetic biological toxicity in water analytical equipment includes Instrument shell, the cultivation unit of photogen, photogen luminous intensity
Acquisition and storage unit, photogen luminous intensity transmission unit, data acceptance unit and data display unit；
The data receiver and display device carry out luminous intensity information the class of pollution of determining sewage toxicity after data processing
Include the following steps：
Step 1：Data processing is carried out to luminous intensity information by the curve fitting algorithm based on numerical value phase approximately principle, that is, is utilized
Least square solution obtains the parametric equation of curve matching：X=[ATA]-1ATB obtains model y1, and exact value is found out, wherein X
For least square method parameter, A, B are fitting coefficient；
Step 2：Data processing is carried out to luminous intensity information by the curve fitting algorithm based on form principle of similarity：
(1) time-concerning impact factor is set as λ, reacts amount of the photogen activity by the time-to-live；If the similar matched curve of form is：y2
(x)=a0+a1x+a2x2+...+anxn, wherein aiFor each level number；
(2) by matched curve combination shape coefficient equationWherein x12iTo be fitted point value, S12It, can for offset value
：
(3) augmented objective function with parameter is constructed, when meeting constraints, equal sign both sides are set up；When being unsatisfactory for constraining
When condition and range, a fully big number μ is taken>0, construct following function：
A={ a in formula1,a2,...,an,λ}；
(4) pass through the solved function extreme value of powell algorithms：
Step a：Selected initial point x(0), n linearly independent vector group, the first search direction system { p of composition0,p1,....pn-1, it gives
Determine precision ε ＞ 0, set k=0, k is initial extreme value；
Step b：Enable y0=xk, successively along { p0,p1,....pn-1In direction carry out a pile search,
It is corresponding to obtain auxiliary iteration point y1,y2,....yn, i.e.,
β in formulaj-1For along pj-1The step-length in direction；
Such as | y(k)-x(k)|<ε is set up, then stops calculating, no to then follow the steps c；
Step c：Acceleration direction is constructed, p is enabledn=yn-y0If | | pn| |≤ε then stops iteration, exports xk+1=yn, otherwise turn to walk
Rapid d；
Step d：Determine adjustment direction：Find out m, wherein m is minimum point so that
f(ym-1)-f(ym)=max | f (ym-1)-f(ym)|1≤j≤n}
If following formula is set up：
f(y0)-2f(yn)+f(2yn-y0) [f (the y of ＜ 2m-1)-f(ym)],
Step f is gone to, e is otherwise gone to step；
Step e：Enable xk+1=yn+βnpn,Meanwhile it enabling
{p0,p1,....pn-1}k+1={ p0,...,pm-1,pm+1,....pn-1,pnK=k+1 goes to step b；
Step f：Enable xk+1=yn, set k=k+1 and go to step b；
(5) by powell algorithms, a is found out1,a2,......,an, it is contemplated that a0For DC component, matched curve is not interfered with
Form substitutes into formula y2(x)=a0+a1x+a2x2+...+anxnIn, a is asked using least square method0, obtain being based on form similarity criterion
Curve fitting algorithm y2；Accuracy is found out using formula of correlation coefficient：
Step 3：From the two above steps, cubic fitting model y is obtained respectively1, y2, constructing improved cubic fitting model is：
Y (x)=w1y1(x)+w2y2(x), wherein 0 ＜ w1≤ 1,0 ＜ w2≤ 1, w1、w2According to the precision value found out；
Step 4：Further optimize each term coefficient found out using Fabonacci methods, it is made to be more nearly optimal value；
It is as follows：
1) set the deviation of respective value in j-th of discrete point and matched curve as
It can be found out according to above formula, the maximum positively biased in n discrete point is not good enoughIt is not good enough with maximum negative biasThen m powers coefficient amInitial section [A1,B1] be
2) first judge aiIn initial section [A1,B1] it is unimodal function, the best a found outTIt is exactly aiIn section [A1,B1] in it is close
Like minimum or maximum, i.e. aiIn [A1,aT] section answer strictly decreasing or be incremented by, in [aT,B1] on answer strictly increasing or successively decrease,
It is optimized using Fabonacci methods；
3) according to the positive and negative inclined absolute value of the difference approximately equal of maximum of optimum fit curve, then deviation takes the half of sum of the two,
I.e.
It can thus be concluded that best constant coefficient is
Step 5：Using the expert system to prestore in obtained improvement matched curve y combination this system databases, curve is carried out
Match, respectively obtains the fitting coefficient B for improving matched curve and photogen in expert system and variety classes, the toxicant of concentration
The fitting coefficient A of each reaction mechanism curve, seeks the distance between fitting coefficient B and each fitting coefficient A, and coefficient minimum will be
The toxicant predicted.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610301333.7A CN106053438B (en) | 2016-05-07 | 2016-05-07 | A kind of synthetic biological toxicity in water remote auto analyzer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610301333.7A CN106053438B (en) | 2016-05-07 | 2016-05-07 | A kind of synthetic biological toxicity in water remote auto analyzer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106053438A CN106053438A (en) | 2016-10-26 |
| CN106053438B true CN106053438B (en) | 2018-09-11 |
Family
ID=57176201
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610301333.7A Active CN106053438B (en) | 2016-05-07 | 2016-05-07 | A kind of synthetic biological toxicity in water remote auto analyzer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106053438B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108375570A (en) * | 2017-12-27 | 2018-08-07 | 深圳市有为环境科技有限公司 | A kind of water quality biological comprehensive toxicity Quality Control instrument and its method |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101477131A (en) * | 2009-01-13 | 2009-07-08 | 浙江理工大学 | On-line automatic analyzer for synthetic biological toxicity in water and its use method |
| CN201344934Y (en) * | 2009-01-13 | 2009-11-11 | 浙江理工大学 | Water quality comprehensive biological toxicity on-line automatic analyzer |
| WO2010075121A1 (en) * | 2008-12-26 | 2010-07-01 | General Electric Company | Control system for monitoring localized corrosion in an industrial water system |
| CN102234609A (en) * | 2010-04-27 | 2011-11-09 | 杭州绿洁水务科技有限公司 | System and method of comprehensive toxicity detection of water body |
| CN102252990A (en) * | 2011-05-05 | 2011-11-23 | 中国科学院合肥物质科学研究院 | Automatic calibration method of water quality monitoring system with application of ultraviolet absorption method |
-
2016
- 2016-05-07 CN CN201610301333.7A patent/CN106053438B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010075121A1 (en) * | 2008-12-26 | 2010-07-01 | General Electric Company | Control system for monitoring localized corrosion in an industrial water system |
| CN101477131A (en) * | 2009-01-13 | 2009-07-08 | 浙江理工大学 | On-line automatic analyzer for synthetic biological toxicity in water and its use method |
| CN201344934Y (en) * | 2009-01-13 | 2009-11-11 | 浙江理工大学 | Water quality comprehensive biological toxicity on-line automatic analyzer |
| CN102234609A (en) * | 2010-04-27 | 2011-11-09 | 杭州绿洁水务科技有限公司 | System and method of comprehensive toxicity detection of water body |
| CN102252990A (en) * | 2011-05-05 | 2011-11-23 | 中国科学院合肥物质科学研究院 | Automatic calibration method of water quality monitoring system with application of ultraviolet absorption method |
Non-Patent Citations (1)
| Title |
|---|
| A simple piecewise cubic spline method for approximation of highly nonlinear data;Mehdi Zamani;《Natural Science》;20111231;第79-83页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106053438A (en) | 2016-10-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110352344A (en) | Can field deployment Multi-channel type sampling and monitoring device and germ contamination measurement method | |
| CN105973877B (en) | A kind of water quality remote on-line monitoring method based on curve matching and toxicity test algorithm | |
| CN201344934Y (en) | Water quality comprehensive biological toxicity on-line automatic analyzer | |
| CN101477131B (en) | On-line automatic analyzer for synthetic biological toxicity in water and its use method | |
| CN104977263A (en) | Multi-parameter water quality monitor and monitoring method | |
| CN102131914A (en) | Apparatus and method for bacteriological testing on plasma | |
| Petropoulos et al. | Re-modeling ELISA kits embedded in an automated system suitable for on-line detection of algal toxins in seawater | |
| CN105277535B (en) | Ammonia nitrogen field fast detection method in a kind of water for eliminating reagent blank influence | |
| CN102175661A (en) | Online analyzer of Escherichia coli | |
| CN201707297U (en) | Escherichia coli fluorescence fast detection system | |
| CN106053438B (en) | A kind of synthetic biological toxicity in water remote auto analyzer | |
| CN115290615B (en) | High-flux method for detecting microplastic content in marine organism | |
| CN103336001B (en) | Method for rapidly evaluating water quality biological stability of drinking water | |
| CN208902601U (en) | Seawater total nitrogen concentration test device | |
| KR101810810B1 (en) | Method and Apparatus for Determining Microorganisms in a Water Sample | |
| Amrita et al. | Analysing the water quality parameters from traditional to modern methods in aquaculture | |
| EP0304406A2 (en) | Apparatus for automatically counting the microorganisms possibly present in liquids, particularly in waters for human use | |
| CN211122516U (en) | Wide-spectrum five-day biochemical oxygen demand online monitoring device | |
| CN207958369U (en) | A kind of microbial detection device | |
| CN209296582U (en) | The rapid quantitative detection device of viable bacteria in a kind of sample | |
| CN101625317A (en) | Full-automatic online chemical oxygen demand (BOD) and biological oxygen demand (BOD) monitor and use method thereof | |
| CN105277542B (en) | A kind of water nitrite field fast detection method for eliminating reagent blank influence | |
| CN108593867A (en) | A kind of water quality automatic detection analysis instrument | |
| CN1243232C (en) | Photochemical and biochemical microorganic membrane dynamic responding sensor for COD with optical fiber | |
| CN202421070U (en) | Detecting unit for medical sewage |
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 |