CN101636655A - 用于对流出物成分的变化进行定性的方法 - Google Patents
用于对流出物成分的变化进行定性的方法 Download PDFInfo
- Publication number
- CN101636655A CN101636655A CN200880008759A CN200880008759A CN101636655A CN 101636655 A CN101636655 A CN 101636655A CN 200880008759 A CN200880008759 A CN 200880008759A CN 200880008759 A CN200880008759 A CN 200880008759A CN 101636655 A CN101636655 A CN 101636655A
- Authority
- CN
- China
- Prior art keywords
- logic
- domain logic
- derivative
- parameter
- domain
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000005259 measurement Methods 0.000 claims abstract description 12
- 230000002547 anomalous effect Effects 0.000 claims description 6
- 239000000126 substance Substances 0.000 claims description 6
- 230000007246 mechanism Effects 0.000 claims description 4
- 230000005055 memory storage Effects 0.000 claims description 3
- 238000012544 monitoring process Methods 0.000 claims description 2
- 230000008859 change Effects 0.000 description 11
- 230000002349 favourable effect Effects 0.000 description 10
- 238000001228 spectrum Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 238000002211 ultraviolet spectrum Methods 0.000 description 6
- 238000011161 development Methods 0.000 description 5
- 238000012512 characterization method Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000002798 spectrophotometry method Methods 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 150000001875 compounds Chemical class 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 235000012206 bottled water Nutrition 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 239000003651 drinking water Substances 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- 230000000007 visual effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 239000010842 industrial wastewater Substances 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 239000010841 municipal wastewater Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 239000002351 wastewater Substances 0.000 description 1
Images
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/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/33—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using ultraviolet light
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F1/00—Treatment of water, waste water, or sewage
-
- 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/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/18—Water
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F1/00—Treatment of water, waste water, or sewage
- C02F1/008—Control or steering systems not provided for elsewhere in subclass C02F
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F2209/00—Controlling or monitoring parameters in water treatment
- C02F2209/001—Upstream control, i.e. monitoring for predictive control
-
- C—CHEMISTRY; METALLURGY
- C02—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F—TREATMENT OF WATER, WASTE WATER, SEWAGE, OR SLUDGE
- C02F2209/00—Controlling or monitoring parameters in water treatment
- C02F2209/003—Downstream control, i.e. outlet monitoring, e.g. to check the treating agents, such as halogens or ozone, leaving the process
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Engineering & Computer Science (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Environmental & Geological Engineering (AREA)
- Hydrology & Water Resources (AREA)
- Water Supply & Treatment (AREA)
- Organic Chemistry (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Treatment Of Sludge (AREA)
- Investigating Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)
- Physical Or Chemical Processes And Apparatus (AREA)
- Preparation Of Compounds By Using Micro-Organisms (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
- Activated Sludge Processes (AREA)
- Data Exchanges In Wide-Area Networks (AREA)
Abstract
本申请涉及一种用于对流出物成分的变化进行定性的方法,该方法包括随时间对流出物的至少第一参数和第二参数执行多个测量。本申请的特点在于,在每一时刻,该方法包括:确定参数的第一导数和第二导数;定义第一逻辑域和第二逻辑域;分配第一导数属于第一逻辑域的概率;分配第二导数属于第二逻辑域的概率;定义全局逻辑域;分配全局属性概率;以及根据所述全局属性概率对流出物成分的变化进行定性。
Description
技术领域
本发明涉及液体或气体流出物成分的分析和在线控制(on-line control)领域,液体或气体流出物例如为饮用水、城市废水或工业废水。
本发明尤其涉及一种检测这种液体或气体流出物成分的变化的方法。
背景技术
在本领域公知的检测方法中,在管线中(in-line)实现测量,例如在废水处理站上游侧的管道中实现测量。
这种方法检测流出物成分中存在的任何异常,特别是任何物理-化学或生物学的异常,这些异常易于在实现测量的位置的处理工艺下游造成事故。
获取管线中的这些测量结果并对它们进行数学处理,以警告操作者有成分与通常送到处理站的流出物成分不同的流出物到达。如果需要,操作者可采取适当措施,例如通过使“不同”的流出物转移到另一处理单元和/或对流出物进行取样以稍后进行分析。
检测流出物的变化的一些在线方法使用可见紫外线(VIS-UV)分光光度技术(spectrophotometry),以基于一组信号表示流出物的定性变化。
对于定性方面,VIS-UV光谱用作流出物的“指纹(fingerprint)”。
换言之,给定的流出物具有其自身的VIS-UV光谱,包括特定点(称作等吸收点,isobestic point),该光谱在理论上是不随时间变化的。
显然,只要检测到管线中测量的UV光谱中有变化,则该变化反映了流出物成分的变化。
所获得的信息向操作员示出流出物的整体质量是否保持不变,但并不提供对流出物所表现出的危害的评估。最后,并不一定所有的流出物都存在这些点,这通常使该方法无效。
如文件FR 2783322和FR 2787883中所描述的,在假设流出物的每个UV光谱可以被看作是少量已知UV光谱(称为参考光谱)的线性组合的基础上,处理UV光谱的另一公知方法使用UV光谱的去卷积,该方法具体针对被研究的流出物中的特殊化合物。
因此,该方法假定流出物的成分或为易于在流出物中找到的化合物为公知常识,这要求事先进行实验室分析。而且,必须事先建立化合物光谱库或预定成分光谱库。
文件EP 1070954涉及使用可见UV光谱中的分光光度技术,以通过对光谱进行去卷积来确定全局参数。该文件提出了随时间监测流出物,以确定可能的污染开始显现的精确时刻。因此该方法要求在实现测量的每个点,通过传统方法预先将流出物特征化,并预先校准、识别被分析的流出物特有的参考光谱。
因此,公知的定性方法完全依赖使用该方法的环境。尤其是,在VIS-UV分光光度技术中,假定流出物在波长所关注的范围内吸收。
现有技术的方法的另一缺陷是它们不能进行其它类型的测量,如pH,传导率或者特别是温度测量。
发明内容
本发明的一个目的是提出一种对流出物成分的变化进行定性(qualify)的方法,其中实现对流出物的至少一个第一物理-化学参数和一个第二物理-化学参数随时间的连续测量,该方法与使用该方法的环境无关,并且该方法能够如处理非参数测量值一样处理参数测量值。
本发明是通过在每个时间步长(time step)内执行下述动作来实现其目的的:
·确定第一参数和第二参数中每一参数的第一导数(first derivative);
·确定第一参数和第二参数之间的至少一个第二导数(secondderivative);
·为每一参数的第一导数定义第一逻辑域,该第一逻辑域包括至少一个正常第一逻辑域和一个异常第一逻辑域,所述正常第一逻辑域与所述第一导数的正常变化相对应,所述异常第一逻辑域与所述第一导数的异常变化相对应;
·对每个所述参数的所述第一导数属于所述第一逻辑域中每一逻辑域的第一概率进行分配;
·定义第二逻辑域,该第二逻辑域包括至少一个正常第二逻辑域和一个异常第二逻辑域,所述正常第二逻辑域与所述第二导数的正常变化相对应,所述异常第二逻辑域与所述第二导数的异常变化相对应;
·对所述第二导数属于所述第二逻辑域中每一逻辑域的第二概率进行分配;
·根据所述第一逻辑域和所述第二逻辑域定义全局逻辑域,该全局逻辑域包括至少一个正常全局逻辑域和一个异常全局逻辑域;
·分配所述参数的第一导数和至少一个第二导数的组合属于所述全局逻辑域中每一逻辑域的全局概率;以及
·根据所述全局概率对所述流出物成分的变化进行定性。
因此,通过本发明,事先不必知道流出物的“通常”成分就能检测其成分变化的异常。
而且,通过本发明,事先也不必为了初始化该方法而进行分析。
相反,本发明的方法只要求在时间上连续测量参数,以适于该方法使用的环境,该环境例如可以是流出物处理网络,或者是饮用水网络或者工业场所。
这意味着本发明的方法简单并可快速使用,而不需要使用公知方法时必须进行的那些预备校验步骤或实验室分析。
而且,有利地,本发明的方法基于分析光谱测量值的变化可与基于分析幅度测量值的变化差不多,例如但不仅仅是pH、氧化还原电势、传导率或者温度。
本发明的主要目的之一是立即、实时告诉操作员可能有含危险成分的流出物到达。
为此,分析所获得各参数中每一参数的变化发展过程,并且还分析第一参数和第二参数之间的相对变化的发展过程。
在本发明的含义中,第一导数是相对于相同的量(优选为时间),而第二导数表示一个参数相对于另一参数的变化,优选是随时间的变化。
如果参数取决于非时间变量,则第二导数例如可以是针对给定时间步长的参数相对于所述变量的导数。
换言之,每个参数发展变化的方式(优选是随时间发展变化的方式)是使用第一导数分析的。这会得到参数的各体方差(individual variance),第一导数是时间导数较为有利。
相反,参数之间相对于彼此的发展变化方式是使用第二导数分析的。这在本文被称作所述参数的互方差(cross variance)。
互方差特别用来确定参数是否是以相似或者至少是相关联的方式发展变化的。
参数的数目优选为明显大于2。
优选第二导数的数目与参数组合的数目一样多。
本发明的根本思想是将个体方差和互方差进行比较,从而建立对流出物成分的全局方差(global variance)的显示。
实施本发明方法的数学处理是以其它领域众所周知的模糊逻辑理论为基础的,因此上文提到的比较是通过在其它领域众所周知的推理运算来实现的。
全局方差显示的作用是可以将流出物成分的变化定性为正常或异常,或者通过任何其它适当的限定方式(qualifier)进行定性。
更准确地讲,分配每个参数的第一导数属于第一逻辑域中每一逻辑域的第一概率表征参数的个体方差,而分配第二导数属于第二逻辑域中每一逻辑域的第二概率表征参数的互方差。
最后,显然分配所述全局概率的步骤表征参数的全局方差。
本发明应用的数学处理优选采用由先前时间步长中的参数值组成的至少一个可变长度移动数据库,从而该方法自身就能适应流出物的自然发展变化过程。
移动数据库是一个包含预定数目的值的数据库,在数据库中放置新数据项使得必须删除最早存储的值。
这暗示操作员不必设置阈值来确定流出物是否变化或者不必设置阈值以使任何这类阈值适应流出物的自然发展变化,这是特别有益的。
相反,通过本发明,所述方法自动确定这些阈值作为存储在数据库中的数据的函数。
如上所述,本发明的方法不需要事先或以其它方式在实验室中对流出物进行任何测量。而且,不需要对所使用的数学工具进行预先校验。
这省去了针对流出物的某些特征进行测量的必要,因此本发明的方法提供了一种与流出物成分的变化关联的危险性风险的总体处理方法。
最后,可以选择简单的定性标准,例如“红灯”表示流出物变化很大,这暗示处理存在故障的风险高,相关联地,例如“绿灯”表示流出物变化较小或者根本没有变化。
对于每个参数的第一导数,根据标准差的倍数定义第一逻辑域是有利的,该标准差是通过在先前时间步长中确定的所述参数的第一导数值组成的集合估计出来的。
对于每个参数的第一导数,还优选根据由先前时间步长中确定的所述参数的第一导数值组成的集合的平均值或中间值定义第一逻辑域。
根据标准差的倍数定义第二逻辑域是有利的,该标准差是根据由先前时间步长中确定的所述至少一个第二导数值组成的集合估计出来的。
根据标准差的倍数定义全局逻辑域是有利的,该标准差是根据由先前时间步长中确定的参数的第一导数值和所述至少一个第二导数值组成的集合估计出来的。
在至少一个移动数据库中(如上文提到的移动数据库)存储在上述先前时间步长中确定的值是有利的。
为了改进对流出物成分的变化的定性,对于每个参数的第一导数,定义与所述第一导数的亚正常(subnormal)变化相对应的第一中间逻辑域。
对于所述至少一个第二导数,定义与所述第二导数的亚正常变化相对应的第二中间逻辑域。
此外,还可以定义与流出物成分的非常异常的变化相对应的一个或多个全局中间逻辑域。
实践中,除了以上提到的红灯和绿灯之外,可以使用“琥珀色灯”,所述琥珀色灯反映流出物在可接受的范围内的变化,风险被估计为适中,但需要提高警惕性(有朝红灯发展的潜在危险)。
在第一实施例中,根据所述第一概率和第二概率加权和计算所述全局概率。
在第二实施例中:
·定义包括至少一个正常第三逻辑域和一个异常第三逻辑域的第三逻辑域,该第三逻辑域是根据第一逻辑域定义的;
·定义包括至少一个正常第四逻辑域和一个异常第四逻辑域的第四逻辑域,该第四逻辑域是根据第二逻辑域定义的;
·对由参数的第一导数组成的集合属于第三逻辑域中每一逻辑域的第三概率进行分配;
·对至少由所述至少一个第二导数组成的集合属于第四逻辑域中每一逻辑域的第四概率进行分配;
·根据第三逻辑域和第四逻辑域定义全局逻辑域。
在第二实施例中,插入另外的第三逻辑域和第四逻辑域,以改进对流出物成分的变化的定性,并防止误报警。
优选地,正常第三逻辑域与正常第一逻辑域的并集对应,异常第三逻辑域与异常第一逻辑域的并集对应,正常第四逻辑域与正常第二逻辑域的并集对应,异常第四逻辑域与异常第二逻辑域的并集对应。
定义亚正常第三逻辑域和亚正常第四逻辑域是有利的,亚正常第三逻辑域优选与亚正常第一逻辑域的并集对应,亚正常第四逻辑域优选与亚正常第二逻辑域的并集对应。
优选定义五个全局逻辑域。
第一全局逻辑域优选与正常第三逻辑域和正常第四逻辑域的并集对应。
第二全局逻辑域优选与正常第三逻辑域和亚正常第四逻辑域的并集以及亚正常第三逻辑域和正常第四逻辑域的并集之间的并集对应。
第三全局逻辑域优选与正常第三逻辑域和异常第四逻辑域的并集与亚正常第三逻辑域和亚正常第四逻辑域的并集以及与异常第三逻辑域和正常第四逻辑域的并集之间的并集对应。
第四全局逻辑域优选与亚正常第三逻辑域和异常第四逻辑域的并集与异常第三逻辑域和亚正常第四逻辑域的并集之间的并集对应。
第五全局逻辑域优选与异常第三逻辑域和异常第四逻辑域的并集对应。
最后,在由参数的第一导数和第二导数组成的集合属于全局逻辑域的全部概率被计算之后,再进行分配。
根据所述第三概率和第四概率加权和计算所述全局概率是有利的。
所述加权和优选包括至少一个动态加权系数,例如由在先前时间步长中确定的第一导数的值组成的集合评估的标准差。
在第一变型中,对流出物成分的变化进行定性的全局逻辑域是具有最大全局概率的逻辑域。
优选地,只用第五全局逻辑域来确定流出物成分的变化。
显然,该第五全局逻辑域是流出物成分的变化特别异常的概率为最大的逻辑域。
在第二变型中,对流出物成分的变化进行定性的全局逻辑域是所述全局概率为非零的最异常逻辑域。
在第三变型中,对流出物成分的变化进行定性的全局逻辑域是所述全局概率大于预定阈值的最异常逻辑域。
在本发明的优选实施例中,第一参数是在第一波长测量的流出物的吸光率。
每个参数优选是给定波长的流出物的吸光率。
优选在每个时间步长中测量多个波长的流出物吸光率,优选波长范围在200纳米(nm)到700纳米(nm)。
第一导数为给定波长下吸光率相对于时间的导数是有利的,第二导数为在给定时间步长中计算的吸光率相对于波长的导数是有利的。
本发明还提供了一种使用本发明的方法跟踪流出物成分的变化的设备。
在阅读下文通过非限制性例子给出的本发明的实施例的描述后,本发明的其它特征和优点是更加显而易见的。
附图说明
说明书参照附图,其中:
图1是表示使用本发明的方法用来监测流出物成分的变化的设备的图;
图2是表示本发明的优选实施例的步骤的图;
图3是表示用来确定属于第一逻辑域的第一导数的概率的三角函数的曲线图;以及
图4是表示用来确定属于第二逻辑域的第二导数的概率的三角函数的曲线图。
具体实施方式
本发明方法的优选实施例的目的是评估化工流出物的成分随时间的变化,借助图1-图4在下文对该优选实施例进行详细描述。
图1示出了使用本发明方法的本发明的设备10。
通过该图可以看出,设备10包括用来测量在管道14中流动的流出物的吸光率的测量装置12,该测量装置12优选但不必然位于流出物处理站15的上游侧。
由测量装置12测量的值被存储在存储装置16中,微处理器18实现计算以对流出物成分的变化进行定性。
计算结果优选由可视界面20给出,该可视界面20包括表征流出物偏离正常水平的变化的三个有色灯,例如,绿灯、琥珀色灯和红灯。
借助图2在下文详细解释本发明的方法。
在该实施例中,实现随时间的连续测量,产生多个参数,即在多个波长中的每个波长下流出物的吸光率。
换言之,在每个时间步长中,在步骤S101期间,在波长为200nm到700nm的范围内(优选在波长230nm到500nm的范围内)测量流出物的吸光率。
优选通过以1nm的步幅改变波长来进行测量,时间步长t优选持续两小时。
图2表示在每个时间步长中执行的步骤。
因此,在该实施例中,在波长为230nm到500nm的范围内测量吸光率,每两小时同时测量流出物中的271个参数。因此,显然本发明的含义是,第一参数是波长为230nm的吸光率,第二参数是波长为231nm的吸光率,直到第271个参数,其为波长为500nm的吸光率。
显然,参数的数目不局限于271个,在不偏离本发明范围的情况下可以选择不同数目的参数。优选地但不必须,存储装置包括存储用于多个时间步长t的参数测量值的数据库。
在第二步骤S102中,通过计算相同波长下两个连续吸光率的值之间的差,再用该结果除以两个测量值之间的时间差Δt,来确定271个参数中的每个参数的第一导数。
因此,如果考虑第一参数,即流出物在230nm的吸光率Abs,则第一参数的第一导数d1,1(t)的表达式如下:
同样地,对于第二参数的第一导数d1,2(t),d1,2(t),其表达式如下:
如此直到得到第271个参数的第一导数。
在第三步骤S103中,确定参数间的至少一个第二导数。
为了更准确起见,在此例子中,定义271个参数的270个第二导数。第一个第二导数d2,1(t)根据波长为231nm的吸光率的值和波长为230nm的吸光率的值之间的差计算,即:
因此,显然第二导数d2,1(t)是参数Abs231(t)和Abs230(t)之间的相对差,从而函数t->d2,1(t)给出这两个参数之间的相对差随时间的变化。
同样地,对于第二个第二导数d2,2(t),其表达式如下:
如此直到得到第270个第二导数。
在第四步骤S104中,对每一第j个(j=1到271)第一导数(即对每个参数的第一导数)定义三个第一逻辑域(DL1 i(j),i=1到3,j=1到271),三个第一逻辑域DL1 i(j)通过以下计算确定:
和:
b)对于每一第j个(j=1到271)第一导数,其与第一导数的平均值的绝对差ε1,j(t)通过下式计算:
ε1,j(t)=|d1,j(t)-d1,j(t)|
c)然后在下一时间步长t,按如下所述定义第一逻辑域DL1 i(j,t)(i=1到3,j=1到271):
·每一第j个第一导数的第一个第一逻辑域DL1 1(j,t)“第一导数的正常变化”的特点在于:ε1,j(t)<σ1,j(t);
·每一第j个第一导数的第二个第一逻辑域DL1 2(j,t)“第一导数的亚正常变化”的特点在于:σ1,j(t)<ε1,j(t)<2σ1,j(t);
·每一第j个第一导数的第三个第一逻辑域DL1 3(j,t)“第一导数的异常变化”的特征在于:ε1,j(t)>2σ1,j(t)。
在第五步骤S105中,通过实施以下子步骤,来分配在时间t时271个第一导数属于第一逻辑域中每一逻辑域的概率P(DL1 i(j,t))(i=1到3,j=1到271):
a)通过下列对应关系,将得到的每一第j个第一导数的绝对偏差标准化到-1和1之间:
如果ε1,j(t)=0,则‖ε1,j(t)‖=-1;
如果ε1,j(t)≥3σ1,j(t).d1,j(t),则‖ε1,j‖=+1;
否则
b)接着进行模糊化,即用图3中表示的三角函数(triangle function)对每一第j个第一导数隶属于每一标准化值‖ε1,j(t)‖的第一逻辑域DL1 i(j,t)进行分配。
图3表示这样一种情况,其中对于第50个第一导数,在时间t0具有值‖ε1,50(t0)‖=0.3,得到 它们对应于第50个第一导数隶属于多个第一逻辑域的概率。
在第六步骤S106中,定义第二逻辑域,即对于每一第k个第二导数,通过以下计算确定三个第二逻辑域(DL2 i(k,t),i=1到3,k=1到270):
a)计算每个第二导数的后50个值的平均值和标准差:
和:
b)计算每一第k个第二导数与第二导数的平均值的绝对偏差ε2,k(t):
ε2,k(t)=|d2,k(t)-d2,k(t)|
c)然后以如下方式定义第二逻辑域DL2 i(k,t):
·第一个第二逻辑域DxL2 1(k,t)(“第k个第二导数的正常变化”)的特征在于ε2,k(t)<σ2,k(t);
·第二个第二逻辑域DL2 2(k,t)(“第k个第二导数的亚正常变化”)的特征在于σ2,k(t)<ε2,k(t)<2σ2,k(t);
·第三个第二逻辑域DL2 3(k,t)(“第k个第二导数的异常变化”)的特征在于ε2,k(t)>2σ2,k(t)。
在第七步骤S107中,通过实施以下子步骤,分配第k个第二导数属于第二逻辑域的概率P(DL2 i(k,t))(i=1到3,k=1到270):
a)通过下列对应关系,将得到的绝对差标准化到-1和1之间:
·如果ε2,k(t)=0则‖ε2,k(t)‖=-1;
·如果ε2,k(t)≥3σ2,k(t).d2,k(t)则‖ε2,k(t)‖=+1;
·否则
b)接着进行模糊化处理,即用图4中表示的三角函数分配每一第k个第二导数属于每个标准化值‖ε2,k(t)的第二逻辑域DL2 i(k,t)的概率。
图4表示这样一种情况,即对于第50个第二导数,在时刻t0具有值‖ε2,50(t0)‖=0.4,得到 它们对应于第50个第二导数属于第二逻辑域的概率。
在第八步骤S108中,优选但不是必要地,定义三个逻辑域,即,为第一导数定义三个第三逻辑域(DL3 i(t),i=1到3),这三个逻辑域对应于:
a)对于第一个第三逻辑域DL3 1(t),对应于正常的第一个第一本地域的并集;
b)对于第二个第三逻辑域DL3 2(t),对应于亚正常的第二个第一逻辑域的并集;
c)对于第三个第三本地域DL3 3(t),对应于异常的第三个第一本地域的并集。
在第九步骤S109中,优选但不是必要地,通过对每一第j个第一导数属于第一逻辑域的概率由第一导数的标准差加权之后进行求和,来计算所有271个第一导数属于第三逻辑域的概率,i=1到3;
在第十步骤S110中,优选但不是必要地定义第四逻辑域,即为第二导数定义三个第四逻辑域(DL4 i(t),i=1至3),这三个第四逻辑域对应于:
a)对于第一个第四逻辑域DL4 1(t),对应于正常的第一个第二逻辑域的并集;
b)对于第二个第四逻辑域DL4 2(t),对应于亚正常的第二个第二逻辑域的并集;
c)对于第三个第四逻辑域DL4 3(t),对应于异常的第三个第二逻辑域的并集。
在第十一步骤S111中,也是优选但不是必要地,通过对每一第k个第二导数属于第二逻辑域的概率由第二导数的标准差加权之后进行求和,来计算270个第二导数中每一第二导数属于三个第四逻辑域的概率,其中i=1到3:
在第十二步骤S112中,根据第三逻辑域第四逻辑域定义全局逻辑域,并在此基础上定义流出物的全局变化,优选定义5个全局逻辑域(DL5 l(t),l=1到5)。它们对应于:
a)对于第一全局逻辑域DL5 l(t),其对应于正常的第一个第三逻辑域和正常的第一个第四逻辑域的并集;
b)对于第二全局逻辑域DL5 2(t),其对应于正常的第一个第三逻辑域和亚正常的第二个第四逻辑域的并集以及亚正常的第二个第三逻辑域和正常的第一个第四逻辑域的并集;
c)对于第三全局逻辑域DL5 3(t),其对应于正常的第一个第三逻辑域和异常的第三个第四逻辑域的并集、亚正常的第二个第三逻辑域和亚正常的第二个第四逻辑域的并集以及异常的第三个第三逻辑域和正常的第一个第四逻辑域的并集;
d)对于第四全局逻辑域DL5 4(t),其对应于亚正常的第二个第三逻辑域和异常的第三个第四逻辑域的并集以及异常的第三个逻辑域和亚正常的第二个第四逻辑域的并集;
e)对于第五全局逻辑域DL5 5(t),其对应于异常的第三个第三逻辑域和异常的第三个第四逻辑域的并集。
在第十三步骤S113中,分配由参数的第一导数和第二导数组成的集合属于每个全局逻辑域(DL5 l(t),l=1到5)的全局概率(P(DL5 l(t),l=1到5),所述概率优选按以下公式进行计算:
在第十四步骤S114中,通过用属于第五全局逻辑域的每一概率除以这些概率的和来对所获得的结果进行标准化。
然后,在第十五步骤S115中,根据全局概率对流出物成分的变化进行定性,优选只使用第五全局逻辑域。例如,进行以下测试:
·如果属于第五全局逻辑域的概率大于第一阈值(例如66%),则流出物的成分变化很大。界面20显示红灯。
·如果属于第五全局逻辑域的概率例如在33%到66%之间,则流出物成分已经变化。界面20显示琥珀色灯。
·如果属于第五全局逻辑域的概率小于第二阈值(例如33%),则流出物成分没有变化。界面20显示绿灯。
该非限制性实施例最小化了误报警的情况,是有利的。对于大多数第一导数和第二导数而言,具有属于第五全局逻辑域的异常逻辑域的高概率是不可避免的,从而有很大的概率触发警报。
如上所述,步骤S107到S111对于本发明的实现不是必需的。
本发明的一个变型不使用第三逻辑域、第四逻辑域,而直接通过第一逻辑域和第二逻辑域定义全局逻辑域,分配由参数的第一导数和第二导数组成的集合属于全局逻辑域中每一逻辑域的全局概率,然后以与步骤S112和S115类似的方式根据全局概率对流出物成分的变化进行定性。
Claims (20)
1.一种对流出物成分的变化进行定性的方法,其中随着时间的进行,连续测量所述流出物的至少一个第一物理-化学参数和一个第二物理-化学参数,所述方法的特征在于,在每个时间步长进行以下动作:
·确定第一参数和第二参数中每一参数的第一导数;
·确定第一参数和第二参数之间的至少一个第二导数;
·为每个所述参数的第一导数定义第一逻辑域,该第一逻辑域包括至少一个正常第一逻辑域和一个异常第一逻辑域,所述正常第一逻辑域与所述第一导数的正常变化相对应,所述异常第一逻辑域与所述第一导数的异常变化相对应;
·对每个所述参数的第一导数属于所述第一逻辑域中每一逻辑域的第一概率进行分配;
·定义第二逻辑域,该第二逻辑域包括至少一个正常第二逻辑域和一个异常第二逻辑域,所述正常第二逻辑域与所述第二导数的正常变化相对应,所述异常第二逻辑域与所述第二导数的异常变化相对应;
·对所述第二导数属于所述第二逻辑域中每一逻辑域的第二概率进行分配;
·根据所述第一逻辑域和所述第二逻辑域定义全局逻辑域,该全局逻辑域包括至少一个正常全局逻辑域和一个异常全局逻辑域;
·对所述参数的第一导数和至少一个第二导数的组合属于所述全局逻辑域中每一逻辑域的全局概率进行分配;以及
·根据所述全局概率对所述流出物成分的变化进行定性。
2.根据权利要求1所述的定性方法,其中对于每个所述参数的第一导数,根据通过由先前时间步长中确定的所述参数的所述第一导数的值组成的集合而被估计出来的标准差的倍数,定义所述第一逻辑域。
3.根据权利要求1或权利要求2所述的定性方法,其中对于每个所述参数的第一导数,还根据由先前时间步长中确定的所述参数的所述第一导数的值组成的集合的平均值或中间值定义所述第一逻辑域。
4.根据权利要求1-3中任一项权利要求所述的定性方法,其中根据基于由先前时间步长中确定的所述至少一个第二导数的值组成的集合而被估计出来的标准差的倍数,定义所述第二逻辑域。
5.根据权利要求1-4中任一项权利要求所述的定性方法,其中根据基于由先前时间步长中确定的所述参数的所述第一导数的值和所述至少一个第二导数的值组成的集合而估计出来的标准差的倍数,定义所述全局逻辑域。
6.根据权利要求1-5中任一项权利要求所述的定性方法,其中对于每个所述参数的所述第一导数,定义与所述第一导数的亚正常变化相对应的第一中间逻辑域。
7.根据权利要求1-6中任一项权利要求所述的定性方法,其中对于所述至少一个第二导数,定义与所述第二导数的亚正常变化相对应的第二中间逻辑域。
8.根据权利要求1-7中任一项权利要求所述的定性方法,其中定义中间全局逻辑域。
9.根据权利要求1-8中任一项权利要求所述的定性方法,其中根据所述第一概率和所述第二概率的加权和计算所述全局概率。
10.根据权利要求1-8中任一项权利要求所述的定性方法,其中:
·定义包括至少一个正常第三逻辑域和一个异常第三逻辑域的第三逻辑域,该第三逻辑域是根据所述第一逻辑域定义的;
·定义包括至少一个正常第四逻辑域和一个异常第四逻辑域的第四逻辑域,该第四逻辑域是根据所述第二逻辑域定义的;
·对由所述参数的第一导数组成的集合属于所述第三逻辑域中每一逻辑域的第三概率进行分配;
·对至少由所述至少一个第二导数组成的集合属于所述第四逻辑域中每一逻辑域的第四概率进行分配;
·根据所述第三逻辑域和所述第四逻辑域定义所述全局逻辑域。
11.根据权利要求10所述的定性方法,其中根据所述第三概率和所述第四概率的加权和计算所述全局概率。
12.根据权利要求11所述的定性方法,其中所述加权和包括至少一个动态加权系数。
13.根据权利要求1-12中任一项权利要求所述的定性方法,其中对所述流出物成分的变化进行定性的全局逻辑域是具有最大全局概率的逻辑域。
14.根据权利要求1-12中任一项权利要求所述的定性方法,其中对所述流出物成分的变化进行定性的全局逻辑域是所述全局概率为非零的最异常的逻辑域。
15.根据权利要求1-12中任一项权利要求所述的定性方法,其中对所述流出物成分的变化进行定性的全局逻辑域是所述全局概率为大于预定阈值的最异常的逻辑域。
16.根据权利要求1-15中任一项权利要求所述的定性方法,其中所述第一导数是时间导数。
17.根据权利要求1-16中任一项权利要求所述的定性方法,其中所述第一参数是在第一波长测量到的所述流出物的吸光率。
18.根据权利要求1-17中任一项权利要求所述的定性方法,其中每个所述参数均是给定波长下所述流出物的吸光率。
19.根据权利要求18所述的定性方法,其中所述第一导数是给定波长的所述吸光率相对于时间的导数,所述第二导数是在给定时间步长中计算出的所述吸光率相对于波长的导数。
20.一种使用根据权利要求1-19中任何一项权利要求所述的方法监测流出物成分的变化的设备(10),该设备(10)包括测量装置(12)、存储装置(16)和微处理器(18)。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR0752882A FR2911960B1 (fr) | 2007-01-25 | 2007-01-25 | Procede de qualification de la variabilite de la composition d'un effluent. |
FR0752882 | 2007-01-25 | ||
PCT/FR2008/050110 WO2008104647A2 (fr) | 2007-01-25 | 2008-01-24 | Procede de qualification de la variabilite de la composition d'un effluent |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101636655A true CN101636655A (zh) | 2010-01-27 |
CN101636655B CN101636655B (zh) | 2013-07-03 |
Family
ID=38470134
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008800087597A Expired - Fee Related CN101636655B (zh) | 2007-01-25 | 2008-01-24 | 用于对流出物成分的变化进行定性的方法 |
Country Status (10)
Country | Link |
---|---|
US (1) | US8249818B2 (zh) |
EP (1) | EP2118648B1 (zh) |
KR (1) | KR20090121277A (zh) |
CN (1) | CN101636655B (zh) |
AT (1) | ATE471511T1 (zh) |
CA (1) | CA2676041C (zh) |
DE (1) | DE602008001558D1 (zh) |
DK (1) | DK2118648T3 (zh) |
FR (1) | FR2911960B1 (zh) |
WO (1) | WO2008104647A2 (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105026926A (zh) * | 2013-03-04 | 2015-11-04 | 株式会社岛津制作所 | 色谱数据处理装置以及处理方法 |
CN105518438A (zh) * | 2013-07-08 | 2016-04-20 | 弗劳恩霍弗应用技术研究院 | 用于识别和量化系统中的放射粒子的方法 |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2943696B1 (fr) | 2009-03-24 | 2016-08-26 | Veolia Eau - Cie Generale Des Eaux | Installation et procede de controle de la qualite de l'eau dans un reseau d'eau potable |
FR3068051B1 (fr) | 2017-06-22 | 2022-06-24 | Aquassay | Systeme de gestion des eaux usees |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0707247B1 (de) * | 1994-10-11 | 2007-02-07 | Endress + Hauser Conducta Gesellschaft für Mess- und Regeltechnik mbH + Co.KG. | Analysegerät, insbesondere für Abwasser |
DE29607093U1 (de) * | 1996-04-19 | 1996-07-18 | Sibet Gmbh Sican Forschungs Un | Adaptiv, wassergüteabhängig gesteuertes Abwasserbauwerk |
CA2194831C (en) * | 1997-01-10 | 2000-11-07 | Jean Millo | Effluent monitoring system |
US6023065A (en) * | 1997-03-10 | 2000-02-08 | Alberta Research Council | Method and apparatus for monitoring and controlling characteristics of process effluents |
FR2783322B1 (fr) * | 1998-09-11 | 2001-03-09 | Naphtachimie Sa | Procede et dispositif de controle de qualite d'effluents |
FR2787883B1 (fr) | 1998-11-30 | 2001-03-16 | Naphtachimie Sa | Procede et dispositif de controle de qualite d'effluents par spectrophotometrie |
DE19857014A1 (de) * | 1998-12-09 | 2000-06-15 | Wtw Weilheim | Verfahren zur Messung von Wasser- und Abwasserparametern |
US8958917B2 (en) * | 1998-12-17 | 2015-02-17 | Hach Company | Method and system for remote monitoring of fluid quality and treatment |
FR2796722B1 (fr) | 1999-07-20 | 2001-09-14 | Secomam Sa | Procede et dispositif pour la numerisation, l'archivage, la caracterisation et l'analyse d'un echantillon liquide |
US6845336B2 (en) * | 2002-06-25 | 2005-01-18 | Prasad S. Kodukula | Water treatment monitoring system |
-
2007
- 2007-01-25 FR FR0752882A patent/FR2911960B1/fr not_active Expired - Fee Related
-
2008
- 2008-01-24 DK DK08761976.3T patent/DK2118648T3/da active
- 2008-01-24 WO PCT/FR2008/050110 patent/WO2008104647A2/fr active Application Filing
- 2008-01-24 EP EP08761976A patent/EP2118648B1/fr not_active Not-in-force
- 2008-01-24 CA CA2676041A patent/CA2676041C/en not_active Expired - Fee Related
- 2008-01-24 AT AT08761976T patent/ATE471511T1/de not_active IP Right Cessation
- 2008-01-24 KR KR1020097015757A patent/KR20090121277A/ko not_active Application Discontinuation
- 2008-01-24 DE DE602008001558T patent/DE602008001558D1/de active Active
- 2008-01-24 US US12/523,890 patent/US8249818B2/en not_active Expired - Fee Related
- 2008-01-24 CN CN2008800087597A patent/CN101636655B/zh not_active Expired - Fee Related
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105026926A (zh) * | 2013-03-04 | 2015-11-04 | 株式会社岛津制作所 | 色谱数据处理装置以及处理方法 |
CN110068639A (zh) * | 2013-03-04 | 2019-07-30 | 株式会社岛津制作所 | 色谱数据处理装置以及处理方法 |
CN105518438A (zh) * | 2013-07-08 | 2016-04-20 | 弗劳恩霍弗应用技术研究院 | 用于识别和量化系统中的放射粒子的方法 |
CN105518438B (zh) * | 2013-07-08 | 2019-01-15 | 弗劳恩霍弗应用技术研究院 | 用于识别和量化系统中的放射粒子的方法 |
Also Published As
Publication number | Publication date |
---|---|
ATE471511T1 (de) | 2010-07-15 |
EP2118648A2 (fr) | 2009-11-18 |
EP2118648B1 (fr) | 2010-06-16 |
CN101636655B (zh) | 2013-07-03 |
KR20090121277A (ko) | 2009-11-25 |
CA2676041C (en) | 2014-10-07 |
WO2008104647A3 (fr) | 2008-11-06 |
WO2008104647A2 (fr) | 2008-09-04 |
US8249818B2 (en) | 2012-08-21 |
CA2676041A1 (en) | 2008-09-04 |
US20100125422A1 (en) | 2010-05-20 |
FR2911960A1 (fr) | 2008-08-01 |
DK2118648T3 (da) | 2010-10-11 |
FR2911960B1 (fr) | 2009-04-24 |
DE602008001558D1 (de) | 2010-07-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Sharafi et al. | Human health risk assessment for some toxic metals in widely consumed rice brands (domestic and imported) in Tehran, Iran: uncertainty and sensitivity analysis | |
Miguntanna et al. | Determination of a set of surrogate parameters to assess urban stormwater quality | |
Brighenti et al. | Two calibration methods for modeling streamflow and suspended sediment with the swat model | |
Sun et al. | Long-term stormwater quantity and quality analysis using continuous measurements in a French urban catchment | |
Tian et al. | Investigation and simulations of changes in the relationship of precipitation-runoff in drought years | |
Bertrand-Krajewski et al. | Accounting for sensor calibration, data validation, measurement and sampling uncertainties in monitoring urban drainage systems | |
CN104034684A (zh) | 一种基于紫外-可见吸收光谱的水质多指标检测方法 | |
CN101636655B (zh) | 用于对流出物成分的变化进行定性的方法 | |
Tsoumanis et al. | Monitoring and classification of wastewater quality using supervised pattern recognition techniques and deterministic resolution of molecular absorption spectra based on multiwavelength UV spectra deconvolution | |
Baoyuan et al. | An empirical soil loss equation | |
Harisuseno | Comparative study of meteorological and hydrological drought characteristics in the Pekalen River basin, East Java, Indonesia | |
Das et al. | Groundwater quality monitoring by correlation, regression and hierarchical clustering analyses using WQI and PAST tools | |
Feng et al. | Dynamic regression modeling of daily nitrate-nitrogen concentrations in a large agricultural watershed | |
Li et al. | Effects of sampling strategies and estimation algorithms on total nitrogen load determination in a small agricultural headwater watershed | |
CN110672144A (zh) | 污染源检测方法和装置 | |
Kulasova et al. | A method of computing uncertain nitrogen and phosphorus loads in a small stream from an agricultural catchment using continuous monitoring data | |
Brauer et al. | Temporal variability in water quality of agricultural tailwaters: Implications for water quality monitoring | |
Durdu | Stochastic approaches for time series forecasting of boron: a case study of Western Turkey | |
Santos et al. | Degradation prediction model for friction in highways | |
McNeil et al. | Defining the climatic signal in stream salinity trends using the Interdecadal Pacific Oscillation and its rate of change | |
Flores-Lopez et al. | A multivariate analysis of covariance to determine the effects of near-stream best management practices on nitrogen and phosphorus concentrations on a dairy farm in the New York Conservation Effects Assessment Project watershed | |
CN105320993A (zh) | 基于证据理论的水资源污染检测方法 | |
Rahmat et al. | Analysis of spatio-temporal trends using Standardised Precipitation Index (SPI) | |
Smeti et al. | Control charts for the toxicity of finished water—modeling the structure of toxicity | |
Bertrand-Krajewski et al. | Accounting for sensor calibration, concentration heterogeneity, measurement and sampling uncertainties in monitoring urban drainage systems |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130703 Termination date: 20150124 |
|
EXPY | Termination of patent right or utility model |