CN106248901B - 一种利用鱼类活动电位功率变化监测水中有毒物质的方法 - Google Patents

一种利用鱼类活动电位功率变化监测水中有毒物质的方法 Download PDF

Info

Publication number
CN106248901B
CN106248901B CN201610692197.9A CN201610692197A CN106248901B CN 106248901 B CN106248901 B CN 106248901B CN 201610692197 A CN201610692197 A CN 201610692197A CN 106248901 B CN106248901 B CN 106248901B
Authority
CN
China
Prior art keywords
power
fish
mrow
water
frequency
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
Application number
CN201610692197.9A
Other languages
English (en)
Other versions
CN106248901A (zh
Inventor
杨凌升
纪明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Qingcheng Electronic Technology Co ltd
Original Assignee
Nanjing University of Information Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Nanjing University of Information Science and Technology filed Critical Nanjing University of Information Science and Technology
Priority to CN201610692197.9A priority Critical patent/CN106248901B/zh
Publication of CN106248901A publication Critical patent/CN106248901A/zh
Application granted granted Critical
Publication of CN106248901B publication Critical patent/CN106248901B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/18Water
    • G01N33/186Water using one or more living organisms, e.g. a fish

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Biological Materials (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明公开了一种利用鱼类生物电信号频谱变化来监测水质的方法,通过对采集的鱼类两种活动电信号进行分析,以鱼类呼吸运动和游泳运动所占的功率在总运动功率中的比例变化作为判定依据,识别有毒物质对水体的危害程度。利用本发明,可以在提高系统对毒物敏感性的同时,有效消除了由鱼的个体差异等所带来的测量误差,并且无需外加设备,进一步提高了系统的精度。

Description

一种利用鱼类活动电位功率变化监测水中有毒物质的方法
技术领域
本发明属于环境保护中水质监测的技术领域,特别涉及一种利用鱼类活动电位功率变化监测水中有毒物质的方法。
背景技术
鱼类作为水生食物链的顶层生物,它的生命行为是水质毒性监测的重要指标。鱼类对水环境的变化十分敏感,当水体中的有毒物质达到一定浓度时就会引起一系列中毒反应,这些反应可以监测水体是否受到有毒物质污染的判定依据。
在[1]中,鱼的活动电位(电压值)变化是监测系统判定水体中有有毒物质的判定依据。由于系统对生物电信号提取的精度等问题,该系统只能采用鲫鱼或锦鲤等体长特定的鱼类。体长小于3cm的鱼类由于所得电位小,体长大于15cm的鱼由于电位大,不能使用在该系统中。斑马鱼,青鳉鱼等国际通用的毒性实验标准用鱼具有与人类极为相似的毒性反应,这种鱼的体长均小于3cm,可见在发明[1]中,由于系统的精度以及判定依据的问题,发明的实用性值得商榷。
在[2]中,科学工作者利用德国的DFB多物种净水监测仪对有毒物质进行了监测,DFB系统也提出了将鱼类的运动按照频率不同加以区分,在判定水体是否受到有毒物质污染时,也采用的是监测生物行为强度的变化作为判定依据。在该方法中为了提高灵敏度,对观测电极要加外加电场,类似于AM调制,对监测结果有一定的影响。同时生物个体差异,数量差异以及与电极位置关系等带来的对监测数据的影响并不能消除。
在[3]中,杭州电子科技大学的工作者通过对鱼类的运动图像进行监测,并通过对运动情况的变化进行水质监测。鱼的运动要产生显著的变化(翻转,撞击水槽,规避毒物)必然需要毒物浓度提高,这对系统的精度是个限制。
[1]邹一平,王磊,洪鼎宙,“利用鱼的活动电位变化监测水中有毒物质的方法”,200510019149.5.
[2]赵晓艳,刘丽君,聂湘平,何秀婷,程章,“利用多物种净水监测仪在线监测水体抗生素药物及有机磷农药”,环境科学学报,第30卷第1期,180-185,2010.
[3]胡江龙,方景龙,王大全,“多目标跟踪算法在水质监测中的应用”,机电工程,第29卷第5期,2012.
发明内容
本发明的目的是提出一种利用鱼类生物电信号频谱变化来监测水质的方法,通过对采集的鱼类活动电信号进行分析能够帮助识别有毒物质对水体的危害程度。以鱼类呼吸运动和游泳运动所占的功率在总运动功率中的比例变化作为判定依据,在提高系统对毒物敏感性的同时,有效消除了由鱼的个体差异等所带来的误差,进一步提高了系统的精度。
为实现上述目的,本发明采用的技术方案为:提出一种利用鱼类活动电位功率变化监测水中有毒物质的方法,包括以下几个监测步骤:
步骤(A),对水中鱼类活动电信号f(t)进行测量,鱼类的活动按照频率分为呼吸运动以及游泳运动,设置游泳运动信号频率f1范围为[x0,x1],呼吸运动信号频率f2范围为[x2,x3],游泳运动效率恒定比为a,呼吸运动效率恒定比为b,毒素报警范围为大于θ,x0,x1,x2,x3为常数;
其中,鱼类活动电信号f(t)的信号周期为T;从f(t)中截取的一段,得到一个截短函数fT(t),对这个时间信号进行傅里叶变换可以得到FT(jω);当T→∞时,定义为f(t)的功率谱密度;
步骤(B),对功率谱密度在不同运动频率段进行积分计算得到对应运动的生物电信号功率P:在f1∈[x0,x1]的情况下,对所测得的功率谱密度进行积分可以求得对应游泳运动的生物电信号功率P1,在f2∈[x2,x3]的情况下,对所测得的功率谱密度进行积分可以求得对应呼吸运动的生物电信号功率P2;
其中,在频域下进行计算,为f(t)的功率谱密度,对功率谱密度进行积分得到生物电信号功率,如式(1)所示:
式中,角频率ω0和ω1可由频率计算得到,计算公式为ω=2πf;游泳运动时角频率ω0=2πx0、ω1=2πx1,呼吸运动时角频率ω0=2πx2、ω1=2πx3
步骤(C),对所得的游泳运动和呼吸运动的生物电信号功率求和得到总功率PT,PT=P1+P2;
步骤(D),对所得实时效率P1/PT和P2/PT进行对比测量,当时,报警器发出信号警示,否则,跳转到步骤(A)进行下一轮测量。
与现有技术相比,本发明具有以下优点:
1、将鱼类活动电位按照频率分为呼吸运动和游泳运动,从而能在两种数值的变化中监测水质,提高了系统对水中毒物敏感度;
2、以呼吸运动和游泳运动所占的功率在总运动功率中的比例变化作为判定依据,在提高系统对毒物敏感性的同时,有效消除了由于鱼类的个体差异问题等带来的监测误差,进一步提高了系统的精度;
3、由于测量值为电位信号频率,无需限定鱼类大小或种群,可采用斑马鱼、青鳉鱼等国际通用的毒性实验标准用鱼,其具有与人类极为相似的毒性反应,有效监测饮用水源的安全;
4、无需加外加电场或AM调制,避免因为生物个体、数量差异以及与电极位置关系等带来的对监测数据的影响,提高监测的准确性。
附图说明
图1是本发明的监测方法流程图;
图2是本发明具体实施例中,所述毒物进入水体前两种运动的功率比例变化示意图;
图3是本发明具体实施例中,所述毒物进入水体后两种运动的功率比例变化示意图。
具体实施方式
下面结合实施例对本发明作更进一步的说明。
如图1所示,一种利用鱼类活动电位功率变化监测水中有毒物质的方法,包括以下几个监测步骤:
步骤(A),对水中鱼类活动电信号f(t)进行测量,鱼类的活动按照频率分为呼吸运动以及游泳运动,设置游泳运动信号频率f1范围为[x0,x1],呼吸运动信号频率f2范围为[x2,x3],游泳运动效率恒定比为a,呼吸运动效率恒定比为b,毒素报警范围为大于θ,x0,x1,x2,x3为常数;
其中,鱼类活动电信号f(t)的信号周期为T;从f(t)中截取的一段,得到一个截短函数fT(t),对这个时间信号进行傅里叶变换可以得到FT(jω);当T→∞时,定义为f(t)的功率谱密度;
步骤(B),对功率谱密度在不同运动频率段进行积分计算得到对应运动的生物电信号功率P:在f1∈[x0,x1]的情况下,对所测得的功率谱密度进行积分可以求得对应游泳运动的生物电信号功率P1,在f2∈[x2,x3]的情况下,对所测得的功率谱密度进行积分可以求得对应呼吸运动的生物电信号功率P2;
其中,在频域下进行计算,为f(t)的功率谱密度,对功率谱密度进行积分得到生物电信号功率,如式(1)所示:
式中,角频率ω0和ω1可由频率计算得到,计算公式为ω=2πf;游泳运动时角频率ω0=2πx0、ω1=2πx1,呼吸运动时角频率ω0=2πx2、ω1=2πx3
步骤(C),对所得的游泳运动和呼吸运动的生物电信号功率求和得到总功率PT,PT=P1+P2;
步骤(D),对所得实时效率P1/PT和P2/PT进行对比测量,当时,报警器发出信号警示,否则,跳转到步骤(A)进行下一轮测量。
前述游泳运动效率恒定比a,呼吸运动效率恒定比b的计算步骤为:
步骤(E),在国际标准的安全水质情况下,对所测得的水中鱼类的信号频率f进行判断:在f1∈[x0,x1]的情况下,对所测得的功率谱密度进行积分可以求得对应游泳运动的生物电信号功率P1a,在f2∈[x2,x3]的情况下,对所测得的功率对应的频谱密度进行积分可以求得对应呼吸运动的生物电信号功率P2b;
步骤(F),对所得的游泳运动和呼吸运动的生物电信号功率求和得到总功率PT',PT'=P1a+P2b;
步骤(G),计算a=P1a/PT',b=P2b/PT'。
前述鱼类为斑马鱼或青鳉鱼。
前述步骤(A)中所述游泳运动信号频率f1范围为[0,2Hz]。
前步骤(A)中所述呼吸运动信号频率f2范围为[4Hz,10Hz]。
前步骤(A)中所述毒素报警值为20%。
如下给出该方法的一个具体实施例。
f(t)就是测得的鱼类活动电信号所代表的函数,这是个未知方程,可以通过实测数据进行拟合。在本发明中,直接对时域信号进行了傅里叶变换,后续数据处理都是在频域进行的。信号可以采用金属电极(如铜片,不锈钢铁片等制成),电极接上处理电路。时域信号做傅里叶变换即变为频域信号,然后根据频率区分生物电信号到底是哪部分运动产生的。测量值是频谱密度,但是首先要明确测量信号的频率,这才能求解频谱密度积分的频率范围。
如图2所示,在未投入有毒物质时,所测游泳运动效率恒定比a为P1`/PT=54.3%,呼吸运动效率恒定比b为P2`/PT=45.7%。如图所示各种运动在总功率中所占比例基本恒定,从而消除了鱼个体、数目、疲劳程度等带来的对监测水质的影响。
如图3所示,在水中投入有毒物质后,所测P1/PT=34.3%,P2/PT=74.3%分别与前测a=54.3%和b=45.7%进行比较,均大于20%的毒素报警值,则系统发出警报。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (6)

1.一种利用鱼类活动电位功率变化监测水中有毒物质的方法,其特征在于:包括以下几个监测步骤:
步骤(A),对水中鱼类活动电信号f(t)进行测量,鱼类的活动按照频率分为呼吸运动以及游泳运动,设置游泳运动信号频率f1范围为[x0,x1],呼吸运动信号频率f2范围为[x2,x3],游泳运动效率恒定比为a,呼吸运动效率恒定比为b,毒素报警范围为大于θ,x0,x1,x2,x3为常数,t为时间系数;
其中,鱼类活动电信号f(t)的信号周期为T;从f(t)中截取的一段,得到一个截短函数fT(t),对这个时间信号进行傅里叶变换可以得到FT(jω);当T→∞时,定义为f(t)的功率谱密度;
步骤(B),对功率谱密度在不同运动频率段进行积分计算得到对应运动的的生物电信号功率P:在f1∈[x0,x1]的情况下,对所测得的功率谱密度进行积分可以求得对应游泳运动的生物电信号功率P1,在f2∈[x2,x3]的情况下,对所测得的功率谱密度进行积分可以求得对应呼吸运动的生物电信号功率P2;
其中,在频域下进行计算,为f(t)的功率谱密度,对功率谱密度进行积分得到生物电信号功率,如式(1)所示:
<mrow> <mi>P</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <mi>&amp;pi;</mi> </mrow> </mfrac> <msubsup> <mo>&amp;Integral;</mo> <msub> <mi>&amp;omega;</mi> <mn>0</mn> </msub> <msub> <mi>&amp;omega;</mi> <mn>1</mn> </msub> </msubsup> <munder> <mi>lim</mi> <mrow> <mi>T</mi> <mo>&amp;RightArrow;</mo> <mi>&amp;infin;</mi> </mrow> </munder> <mfrac> <msup> <mrow> <mo>|</mo> <msub> <mi>F</mi> <mi>T</mi> </msub> <mrow> <mo>(</mo> <mi>j</mi> <mi>&amp;omega;</mi> <mo>)</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mi>T</mi> </mfrac> <mi>d</mi> <mi>&amp;omega;</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow>
式中,角频率ω0和ω1可由频率计算得到,计算公式为ω=2πf;游泳运动时角频率ω0=2πx0、ω1=2πx1,呼吸运动时角频率ω0=2πx2、ω1=2πx3
步骤(C),对所得的游泳运动和呼吸运动的生物电信号功率求和得到总功率PT,PT=P1+P2;
步骤(D),对所得实时效率P1/PT和P2/PT进行对比测量,当时,报警器发出信号警示,否则,跳转到步骤(A)进行下一轮测量。
2.根据权利要求1所述的一种利用鱼类活动电位功率变化监测水中有毒物质的方法,其特征在于:所述游泳运动效率恒定比a,呼吸运动效率恒定比b的计算步骤为:
步骤(E),在国际标准的安全水质情况下,对所测得的水中鱼类的信号频率f进行判断:在f1∈[x0,x1]的情况下,对所测得的功率谱密度进行积分可以求得对应游泳运动的生物电信号功率P1a,在f2∈[x2,x3]的情况下,对所测得的功率对应的频谱密度进行积分可以求得对应呼吸运动的生物电信号功率P2b;
步骤(F),对所得的游泳运动和呼吸运动的生物电信号功率求和得到总功率PT',PT'=P1a+P2b;
步骤(G),计算:a=P1a/PT',b=P2b/PT'。
3.根据权利要求1所述的一种利用鱼类活动电位功率变化监测水中有毒物质的方法,其特征在于:所述鱼类为斑马鱼或青鳉鱼。
4.根据权利要求1所述的一种利用鱼类活动电位功率变化监测水中有毒物质的方法,其特征在于:步骤(A)中所述游泳运动信号频率f1范围为[0,2Hz]。
5.根据权利要求1所述的一种利用鱼类活动电位功率变化监测水中有毒物质的方法,其特征在于:步骤(A)中所述呼吸运动信号频率f2范围为[4Hz,10Hz]。
6.根据权利要求1所述的一种利用鱼类活动电位功率变化监测水中有毒物质的方法,其特征在于:步骤(A)中所述的θ为20%。
CN201610692197.9A 2016-08-19 2016-08-19 一种利用鱼类活动电位功率变化监测水中有毒物质的方法 Active CN106248901B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610692197.9A CN106248901B (zh) 2016-08-19 2016-08-19 一种利用鱼类活动电位功率变化监测水中有毒物质的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610692197.9A CN106248901B (zh) 2016-08-19 2016-08-19 一种利用鱼类活动电位功率变化监测水中有毒物质的方法

Publications (2)

Publication Number Publication Date
CN106248901A CN106248901A (zh) 2016-12-21
CN106248901B true CN106248901B (zh) 2018-03-27

Family

ID=57593279

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610692197.9A Active CN106248901B (zh) 2016-08-19 2016-08-19 一种利用鱼类活动电位功率变化监测水中有毒物质的方法

Country Status (1)

Country Link
CN (1) CN106248901B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106932551A (zh) * 2017-02-20 2017-07-07 山东师范大学 环境压力下的鱼类响应模型生成方法及水体质量评估方法
CN107192805B (zh) * 2017-06-29 2023-07-21 南京信息工程大学 一种在线饮用水源水质监测系统
CN109682941A (zh) * 2018-12-21 2019-04-26 山东师范大学 利用鱼类qrs间期评价水体重金属的方法、装置及系统
CN109633112A (zh) * 2018-12-05 2019-04-16 环境保护部华南环境科学研究所 一种利用青鳉鱼行为强度的变化监测水体中多环芳烃的方法
CN112848935B (zh) * 2021-04-01 2023-04-25 南京信息工程大学 一种电动汽车大功率高效率无线充电自适应实现方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1866000A (zh) * 2005-07-21 2006-11-22 邹一平 利用鱼的活动电位变化监测水中有毒物质的方法
WO2008048141A1 (fr) * 2006-06-20 2008-04-24 Kholodkevich Sergei Viktorovic Procédé de monitorage biologique de l'environnement
CN101191792A (zh) * 2006-11-24 2008-06-04 中国科学院生态环境研究中心 基于水生生物回避行为的水质在线安全预警系统和方法
WO2011155943A1 (en) * 2010-06-11 2011-12-15 Empire Technology Development Llc Detection and decomposition of bisphenol-a
CN203732481U (zh) * 2014-01-07 2014-07-23 南京信息工程大学 基于鱼类生物电频谱的水质监测系统
EP2878948A1 (de) * 2013-12-02 2015-06-03 Airbus DS GmbH Verfahren zur Wasserqualitätsüberwachung

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1866000A (zh) * 2005-07-21 2006-11-22 邹一平 利用鱼的活动电位变化监测水中有毒物质的方法
WO2008048141A1 (fr) * 2006-06-20 2008-04-24 Kholodkevich Sergei Viktorovic Procédé de monitorage biologique de l'environnement
CN101191792A (zh) * 2006-11-24 2008-06-04 中国科学院生态环境研究中心 基于水生生物回避行为的水质在线安全预警系统和方法
WO2011155943A1 (en) * 2010-06-11 2011-12-15 Empire Technology Development Llc Detection and decomposition of bisphenol-a
EP2878948A1 (de) * 2013-12-02 2015-06-03 Airbus DS GmbH Verfahren zur Wasserqualitätsüberwachung
CN203732481U (zh) * 2014-01-07 2014-07-23 南京信息工程大学 基于鱼类生物电频谱的水质监测系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
利用鱼的活动电位监测原水有毒物质;王磊;《福建分析测试》;20051231;第14卷(第4期);第2309-2311页,转第2315页 *

Also Published As

Publication number Publication date
CN106248901A (zh) 2016-12-21

Similar Documents

Publication Publication Date Title
CN106248901B (zh) 一种利用鱼类活动电位功率变化监测水中有毒物质的方法
US10330617B2 (en) Wearable sensor badge for toxic industrial chemicals
Li et al. Electronic-tongue colorimetric-sensor array for discrimination and quantitation of metal ions based on gold-nanoparticle aggregation
CN105615881B (zh) 一种监测身体状况的便携装置、系统及方法
CN104407040A (zh) 一种土壤重金属离子浓度自动检测及修复系统
CN108108889A (zh) 一种水质监测数据在线处理方法及装置
CN205568946U (zh) 一种监测身体状况的便携装置及系统
CN106940363A (zh) 一种基于海洋生物行为反应的海洋污染预警方法
Liu et al. Insertable, scabbarded, and nanoetched silver needle sensor for hazardous element depth profiling by laser-induced breakdown spectroscopy
Jou et al. Synthesis and measurement of valve activities by an improved online clam-based behavioral monitoring system
Priya et al. IoT based automation of real time in-pipe contamination detection system in drinking water
Meghana et al. Design and development of real-time water quality monitoring system
Yunus et al. Application of independent component analysis for estimating nitrate contamination in natural water sources using planar electromagnetic sensor
CN211374726U (zh) 水质监测系统
CN104458870A (zh) 一种重金属离子浓度自动检测系统
CN104407041A (zh) 一种土壤重金属离子浓度自动检测及修复方法
Liu et al. Multiplex analysis of fundamental surface water quality indicators by multichannel electrochemical device combined with MFCSC model
CN202092984U (zh) 恶臭在线监测系统
CN203732481U (zh) 基于鱼类生物电频谱的水质监测系统
CN115980286A (zh) 一种利用电子鼻对污水处理厂不同阶段废水检测的方法
Sankhla et al. Screen-Printed Bimetallic Cobalt-Manganese Metal-Organic Framework Electrodes for Electrochemical Detection of Dichlorvos
Yunus et al. Planar electromagnetic sensor for the detection of nitrate and contamination in natural water sources using electrochemical impedance spectroscopy approach
CN205333744U (zh) 一种基于巨磁电阻的空间平面电磁场强测试仪
Liu et al. Developing on-site paper colorimetric monitoring technique for quick evaluating copper ion concentration in mineral wastewater
Bhardwaj et al. Real time assessment of potable water quality in distribution network based on low cost multi-sensor array

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240201

Address after: 211599 No.21 Chenlu Road, Xiongzhou street, Liuhe District, Nanjing City, Jiangsu Province

Patentee after: Nanjing Qingcheng Electronic Technology Co.,Ltd.

Country or region after: China

Address before: 210044, No. 219, Ning six road, Pukou District, Jiangsu, Nanjing

Patentee before: Nanjing University of Information Science and Technology

Country or region before: China