CN105403567B - 光度法自动滴定的自动判定与计算方法 - Google Patents
光度法自动滴定的自动判定与计算方法 Download PDFInfo
- Publication number
- CN105403567B CN105403567B CN201511010600.7A CN201511010600A CN105403567B CN 105403567 B CN105403567 B CN 105403567B CN 201511010600 A CN201511010600 A CN 201511010600A CN 105403567 B CN105403567 B CN 105403567B
- Authority
- CN
- China
- Prior art keywords
- titration
- mrow
- curve
- automatic
- signal
- 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
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/77—Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated by observing the effect on a chemical indicator
- G01N21/78—Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated by observing the effect on a chemical indicator producing a change of colour
- G01N21/79—Photometric titration
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 Or Analyzing Non-Biological Materials By The Use Of Chemical Means (AREA)
Abstract
本发明公开了光度法自动滴定的自动判定与计算方法，属于自动分析仪器领域，包括滴定终止判定方法、滴定终点判定与计算方法。本发明根据滴定中的颜色突变，选择突变最大波长作为光源波长，在初试平均值的基础上判定滴定终止，在信号微分的基础上判定滴定终点。通过本方法能够实现光度法自动滴定中自动判定滴定终止，能够自动判定滴定突变点，并在此基础上计算出试剂消耗量，甚至直接计算出待测量。本方法精密度高，有利于自动滴定装置的小型化、自动化和智能化。
Description
技术领域
本发明涉及滴定的判定与计算，具体涉及一种光度法自动滴定的自动判定与计算方法。
背景技术
滴定是一种将标准溶液逐量加入被分析溶液中，以待测组分、滴定液、反应产物在滴定过程中引起的颜色、沉淀或电导率变化等来确定反应终点，分析溶液成分的方法。按滴定终点判定的原理不同，滴定可分为光度滴定、电位滴定、库伦滴定等；按反应原理的不同，可分为酸碱、氧化还原、沉淀或络合滴定等。
光度滴定是将一定量的标准溶液滴定到待测溶液中，同时测定待测溶液体系在适当波长处的吸光度，通过吸光度对滴定剂用量作图(即光度滴定曲线)来确定反应终点的方法。一般采用突变前、突变中和突变后的数据做三条直线，三条直线交点的中点即是反应终点。但是对于水样成分复杂的水样来说，在滴定突变前的数据是多变的，并且有时候是不成直线的。因此该方法在实际应用中是有局限性的。
实验室中滴定反应对滴定过程要求高，人工操作常常滴定过量，而且消耗的试剂量很大。现在常常采用自动滴定装置来完成滴定过程，但是光度法的自动滴定装置中自动判定滴定终止和滴定终点以及自动计算滴定结果还没有全面实现，这限制了光度法自动滴定装置的有效推广和应用。
发明内容
为了解决上述问题，本发明提供了一种光度法自动滴定的自动判定与计算方法，该方法可以自动判定滴定终止和滴定终点，并自动计算滴定结果。
上述目的是通过如下技术方案实现的：
一种光度法自动滴定的自动判定与计算方法，包括滴定终止判定方法、滴定突变点判定与滴定结果计算方法：
采用光电设备测量光源透过待测液后的信号N，并记录得到时间信号曲线Q；
滴定终止判定方法为：滴定开始后，即时记录滴定过程中的信号N，并比较N值与滴定开始后A秒内信号平均值O的大小，若N≥O+P，则继续滴定，若N＜O+P，则继续滴定B秒后停止滴定，滴定终止；其中，1≤A≤10；-N≤P≤N，使得O+P在突变中点值与突变终端值之间；5≤B≤20，使时间信号曲线Q在滴定终止前的部分为水平直线；
滴定终止后，自动判定滴定突变点并计算滴定结果，滴定突变点判定与滴定结果计算方法包括如下步骤：
(1)滴定终止后，对得到的时间信号曲线Q进行滤波平滑，得到曲线R；
(2)按如下公式计算得到曲线R的微分曲线S，对微分曲线S滤波平滑后得到曲线U，
(3)找出曲线U中的最大值C以及C相邻全部正值，得到数域V；
(4)通过时间轴找出数域V对应于曲线Q中的原始数据QW；
(5)对QW进行多项式拟合，得到信号与时间的多项式函数X；
(6)按如下公式计算得到多项式函数X的微分曲线Z，找出曲线Z中的最大值D，D所对应的时间E即为滴定突变点，
(7)根据滴定时间与滴定剂消耗量之间的关系，得到滴定突变点E时刻滴定剂的消耗量F，得出滴定结果。
进一步地，所述信号N为光电压或吸光度。
进一步地，所述滤波平滑方法为均值滤波
本发明的有益效果：
(1)本发明提供的光度法自动滴定的自动判定与计算方法可以自动判定滴定终止和滴定突变点，并自动计算滴定结果。
(2)将本发明提供的光度法自动滴定的自动判定与计算方法用于自动滴定仪，能够实现光度法滴定过程的完全自动化和智能化，有利于滴定结果精密度的提高，进一步减小试剂消耗量，实现滴定过程的小型化、自动化和智能化。
附图说明
图1为光度法自动滴定的自动判定与计算方法流程图；
图2为一种自动滴定装置图；
图3为光度法COD-Cr滴定信号曲线Q(毛刺信号部分)及其滤波平滑后的曲线R(毛刺信号中间的一条平滑曲线，灰度图中不易看出)；
图4为R的微分曲线S(毛刺信号部分)及其滤波平滑后的曲线U(毛刺信号中间的一条平滑曲线，灰度图中不易看出)；
图5为QW关于时间轴的拟合多项式X；
图6为拟合多项式X的微分曲线Z和其最大值D(227.1s时最大值D为0.01782)；
其中，1、支架；2、步进电机；3、传动轴；4、搅拌桨；5、滴定池；6、滴定池外壳；7、光路通道；8、进排液管；9、进排液孔；10、支架座；11、滴定孔；12、搅拌头。
具体实施方式
下面以滴定COD为施例详细说明本发明的技术方案。
如图1流程如所示的一种光度法自动滴定的自动判定与计算方法，包括滴定终止判定方法、滴定终点判定与计算方法：
采用光电设备测量光源透过待测液后的信号N，并记录得到时间信号曲线Q；
滴定终止判定方法为：滴定开始后，即时记录滴定过程中的信号N，并比较N值与滴定开始后A秒内信号平均值O的大小，若N≥O+P，则继续滴定，若N＜O+P，则继续滴定B秒后停止滴定，滴定终止；其中，1≤A≤10；-N≤P≤N，使得O+P在突变中点值(如图3所示的突变初始值和突变终端值得平均值)与突变终端值(如图3所示)之间；5≤B≤20，使时间信号曲线Q在滴定终止前的部分为水平直线；
滴定终止后，自动判定滴定突变点并计算滴定结果，滴定突变点判定与滴定结果计算方法包括如下步骤：
(1)滴定终止后，对得到的时间信号曲线Q进行滤波平滑，得到曲线R(如图3)；
(2)按如下公式计算得到曲线R的微分曲线S，对微分曲线S滤波平滑后得到曲线U(如图4)，
(3)找出曲线U中的最大值C以及C相邻全部正值(如图4)，得到数域V；
(4)通过时间轴找出数域V对应于时间信号曲线Q中的原始数据QW；
(5)对QW进行多项式拟合，得到信号与时间的多项式函数X(如图5)；
(6)按如下公式计算得到多项式函数X的微分曲线Z，找出曲线Z中的最大值D，D所对应的时间E即为滴定突变点(如图6)，
(7)根据滴定时间与滴定剂消耗量之间的关系，得到滴定突变点E时刻滴定剂的消耗量F，得出滴定结果。
上述方法中，信号N为光电压或吸光度；所述滤波平滑方法可以为均值滤波。数据处理计算过程通过计算机处理单元实现。
该方法可以用于如图2所示的自动连续滴定系统的机械搅拌系统，包括支架1、步进电机2、传动轴3、搅拌桨4、滴定池5和不透光的滴定池外壳6。所述滴定池5位于滴定池外壳6内，滴定池5为石英材质，滴定池外壳6设有能让光线穿过滴定池下部的光路通道7，所述滴定池5底部设有进排液管8，进排液管8穿过滴定池外壳6底部。所述支架1具有不透光的支架座10，支架座10封盖在滴定池外壳6上形成不透光的腔体，支架座10上设有滴定孔11。所述步进电机2设于支架1上，步进电机2的输出轴通过传动轴3与搅拌桨4连接，搅拌桨4穿过支架座10使搅拌头12悬于滴定池中。所述搅拌头12位于光路通道7水平上方。所述不透光的滴定池外壳6和不透光的支架座10用黑色尼龙材料制成。搅拌桨4和搅拌头12均为石英材质。滴定池5底部中心凹，进排液管8设于底部中心，从而有利于将滴定池中的液体排空，利于实验结束后的清洗工作，避免污染。
微流滴定系统包括精密注射泵、三通阀、连接细管和滴定管。精密注射泵与三通阀的公共端连接，三通阀另外两个接口中的一个接口与通入滴定剂中的连接细管连通，另一个接口通过连接细管与滴定管上端连通，滴定管通过滴定孔插入滴定池中。
滴定时，滴定系统中三通阀先切换至公共端与连通滴定剂的接口连通，精密注射泵吸取滴定剂，再将三通阀切换至公共端与滴定管连通的接口连通，滴定管插入滴定池中液面下；搅拌系统开始搅拌，缓慢加速至设定速度2转/s，稳定后开始滴定；以0.0061ml/s的速度推动精密注射泵实现微流滴定。
下面以光度法滴定COD-Cr为例阐述本发明的方法，方法流程图见图1，采用的自动连续滴定装置结构示意图见图2，通过滴定孔的毛细滴定管内径为1mm，长度为4cm。
各种试剂如下：
水样(邻苯二甲酸氢钾)；氧化剂(重铬酸钾)：取5.00ml铬标准贮备液(3.4)置于500ml容量瓶，用水稀释至标线，摇匀，此溶液1mL含1.00μg六价铬；催化剂(硫酸银-硫酸)。
试亚铁灵指示液：称取0.7g七水合硫酸亚铁溶于50mL水中，加入1.5g邻菲啰啉，搅拌溶解，移入100mL容量瓶中，定容稀释至标线，贮存于棕色瓶内。
硫酸亚铁铵标准溶液：称取39g硫酸亚铁铵溶于水中，缓慢加入20mL浓硫酸，待其冷却后移入1000mL容量瓶中，定容稀释至标线，摇匀，得(NH4)2Fe(SO4)2·6H2O浓度约为0.1mol/L的硫酸亚铁铵标准溶液。
本实施例中，采用的光度法自动滴定的自动判定与计算的方法，其流程图如图1所示，其中：N为光电压信号，单位为mv；A为3秒；P为-600mv；B为5秒；Q为滴定信号N与时间t的曲线。数据采集频率为50次/秒，机械搅拌为8圈/秒，比色池检测过程40mm，以恒速8.3μl/s推入滴定剂。
图3显示的是光度法COD-Cr滴定信号线Q，以及滤波平滑后的曲线R。图4显示的是R的微分曲线S，以及S滤波平滑后的曲线U。找出最大值C为0.02059，与C相邻的全部正值数域V对应的时间段为[224.24s,230.5s]，V对应Q中的原始数据QW。图5显示的是QW关于时间轴的拟合多项式X。图6显示的是X的微分曲线Z，其中最大值D为0.01782，对应的时间E为227.1s，E就是突变时间。根据突变时间，滴定剂消耗量和测量值之间的关系，将E带入计算出COD-Cr测量值。
为了验证本发明方法的准确度，对实验室准确配制的浓度分别为30、50、100mg/L的COD浓度溶液进行测定，每个浓度分别进行三次平行测定。采用新配的0.1mol/L硫酸亚铁铵标准溶液作为滴定溶液进行滴定，试验结果见下表：
结果表明：测量值与实际值的偏差均在5％以内，准确度高，该算法具有一定的适用性。
上述实施例的作用在于说明本发明的实质性内容，但并不以此限定本发明的保护范围。本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的实质和保护范围。
Claims (3)
1.一种光度法自动滴定的判定与计算方法，包括滴定终止判定方法、滴定突变点判定与滴定结果计算方法，其特征在于：
采用光电设备测量光源透过待测液后的信号N，并记录得到时间信号曲线Q；
滴定终止判定方法为：滴定开始后，即时记录滴定过程中的信号N，并比较N值与滴定开始后A秒内信号平均值O的大小，若N≥O+P，则继续滴定，若N＜O+P，则继续滴定B秒后停止滴定，滴定终止；其中，1≤A≤10；-N≤P≤N，使得O+P在突变中点值与突变终端值之间；5≤B≤20，使时间信号曲线Q在滴定终止前的部分为水平直线；
滴定终止后，进一步判定滴定突变点并计算滴定结果，滴定突变点判定与滴定结果计算方法包括如下步骤：
(1)滴定终止后，对得到的时间信号曲线Q进行滤波平滑，得到曲线R；
(2)按如下公式计算得到曲线R的微分曲线S，对微分曲线S滤波平滑后得到曲线U，
<mrow>
<mi>s</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
(3)找出曲线U中的最大值C以及C相邻全部正值，得到数域V；
(4)通过时间轴找出数域V对应于时间信号曲线Q中的原始数据QW；
(5)对QW进行多项式拟合，得到信号与时间的多项式函数X；
(6)按如下公式计算得到多项式函数X的微分曲线Z，找出曲线Z中的最大值D，D所对应的时间E即为滴定突变点，
<mrow>
<mi>z</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
(7)根据滴定时间与滴定剂消耗量之间的关系，得到滴定突变点E时刻滴定剂的消耗量F，进一步得出滴定结果。
2.根据权利要求1所述的判定与计算方法，其特征在于：所述信号N为光电压或吸光度。
3.根据权利要求1所述的判定与计算方法，其特征在于：所述滤波平滑方法为均值滤波。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511010600.7A CN105403567B (zh) | 2015-12-29 | 2015-12-29 | 光度法自动滴定的自动判定与计算方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201511010600.7A CN105403567B (zh) | 2015-12-29 | 2015-12-29 | 光度法自动滴定的自动判定与计算方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105403567A CN105403567A (zh) | 2016-03-16 |
| CN105403567B true CN105403567B (zh) | 2017-12-08 |
Family
ID=55469187
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201511010600.7A Active CN105403567B (zh) | 2015-12-29 | 2015-12-29 | 光度法自动滴定的自动判定与计算方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105403567B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106198532B (zh) * | 2016-09-05 | 2018-09-04 | 江苏德林环保技术有限公司 | 一种快速光度滴定的双波长自动控制方法及装置 |
| CN106404775B (zh) * | 2016-09-05 | 2019-11-12 | 江苏德林环保技术有限公司 | 一种快速光度滴定的单波长自动控制方法及装置 |
| CN108120796A (zh) * | 2017-11-20 | 2018-06-05 | 太原鹏跃电子科技有限公司 | 一种检测铁路蓄电池co32-时ph值突变点的测定计算方法 |
| CN109342408B (zh) * | 2018-10-08 | 2020-11-17 | 辽宁科技大学 | 基于图像颜色信息提取的工业化学智能滴定方法及系统 |
| CN112505239A (zh) * | 2020-09-10 | 2021-03-16 | 浙江伯利恒仪器设备有限公司 | 用于粮食脂肪酸值测定分析仪器的滴定终点自动判定方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3730685A (en) * | 1971-02-24 | 1973-05-01 | Du Pont | Computer controlled method of titration |
| JPS62110144A (ja) * | 1985-11-08 | 1987-05-21 | Agency Of Ind Science & Technol | 温度の異なる電極を用いた電位差滴定法 |
| CN1363834A (zh) * | 2001-11-22 | 2002-08-14 | 华南热带农产品加工设计研究所 | 氯化天然橡胶氯含量快速测定方法 |
| CN1392409A (zh) * | 2002-06-27 | 2003-01-22 | 上海精密科学仪器有限公司 | 能自动判断终点的自动电位滴定法 |
| WO2005003754A2 (en) * | 2003-06-24 | 2005-01-13 | Exxonmobil Research And Engineering Company | Computer-controlled automated titration apparatus for optically determining incompatibility of petroleum oils |
| CN101451974A (zh) * | 2008-12-24 | 2009-06-10 | 浙江赞宇科技股份有限公司 | 脂肪酸甲酯磺酸钠中二钠盐含量的测定方法 |
| CN102590318A (zh) * | 2012-02-24 | 2012-07-18 | 深圳市特种设备安全检验研究院 | 连续分析pH值、酚酞碱度、总碱度和氯离子浓度的分析方法及装置 |
-
2015
- 2015-12-29 CN CN201511010600.7A patent/CN105403567B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3730685A (en) * | 1971-02-24 | 1973-05-01 | Du Pont | Computer controlled method of titration |
| JPS62110144A (ja) * | 1985-11-08 | 1987-05-21 | Agency Of Ind Science & Technol | 温度の異なる電極を用いた電位差滴定法 |
| CN1363834A (zh) * | 2001-11-22 | 2002-08-14 | 华南热带农产品加工设计研究所 | 氯化天然橡胶氯含量快速测定方法 |
| CN1392409A (zh) * | 2002-06-27 | 2003-01-22 | 上海精密科学仪器有限公司 | 能自动判断终点的自动电位滴定法 |
| WO2005003754A2 (en) * | 2003-06-24 | 2005-01-13 | Exxonmobil Research And Engineering Company | Computer-controlled automated titration apparatus for optically determining incompatibility of petroleum oils |
| CN101451974A (zh) * | 2008-12-24 | 2009-06-10 | 浙江赞宇科技股份有限公司 | 脂肪酸甲酯磺酸钠中二钠盐含量的测定方法 |
| CN102590318A (zh) * | 2012-02-24 | 2012-07-18 | 深圳市特种设备安全检验研究院 | 连续分析pH值、酚酞碱度、总碱度和氯离子浓度的分析方法及装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105403567A (zh) | 2016-03-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105403567B (zh) | 光度法自动滴定的自动判定与计算方法 | |
| CN102012389B (zh) | 使用传感器的定量分析方法和定量分析装置 | |
| US3551109A (en) | Method and apparatus for the titration of chloride and bicarbonate in serum | |
| CN107703202A (zh) | 一种高锰酸盐指数在线监测自动滴定判定方法及其装置 | |
| JPS6142211B2 (zh) | ||
| CN101793902A (zh) | 一种流动注射快速分析水质余氯的装置及其分析方法 | |
| CN1176372C (zh) | 一种电位滴定仪 | |
| CN111272751A (zh) | 一种滴定法水质在线自动监测系统及方法 | |
| CN107884535A (zh) | 一种水质在线分析仪自动校准装置及校准方法 | |
| CN111323531A (zh) | 一种自动滴定过程中快慢结合的滴定装置及滴定方法 | |
| CN203572811U (zh) | 快速测定非水溶液中碱度的装置 | |
| CN105510518B (zh) | 一种自动连续微流滴定系统和方法 | |
| CN206618713U (zh) | 一种水质氨氮在线监测装置 | |
| Gaião et al. | An automated flow-injection titrator for spectrophotometric determinations of total acidity in wines, using a single standard solution and gradient calibration | |
| CN207248817U (zh) | 用于测定化学镀液稳定性的自动滴定电位仪 | |
| CN105866219A (zh) | 羟乙基/羟丙基甲基纤维素中氯化钠含量快速测定方法 | |
| JPH02268271A (ja) | 溶液中遊離酸および第２鉄イオンの同時定量分析方法および装置 | |
| JP4474628B2 (ja) | 自動滴定装置 | |
| CN114924029A (zh) | 一种自动滴定系统及其实现方法 | |
| CN108508138A (zh) | 一种重量计量自动滴定装置及其使用方法 | |
| CN211785161U (zh) | 一种滴定法水质在线自动监测系统 | |
| CN114236039A (zh) | 一种与Zeta电位仪连用的自动滴定仪及其使用方法 | |
| CN205581070U (zh) | 纤维素醚中氯化钠含量快速测定装置 | |
| CN105954343A (zh) | 羟乙基/羟丙基甲基纤维素中氯化钠含量快速测定装置 | |
| JPH11108917A (ja) | 滴定制御方法 |
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 |