CN110057990A - 一种多参数水质剖面仪的pH校正方法 - Google Patents
一种多参数水质剖面仪的pH校正方法 Download PDFInfo
- Publication number
- CN110057990A CN110057990A CN201910177098.0A CN201910177098A CN110057990A CN 110057990 A CN110057990 A CN 110057990A CN 201910177098 A CN201910177098 A CN 201910177098A CN 110057990 A CN110057990 A CN 110057990A
- Authority
- CN
- China
- Prior art keywords
- data
- water quality
- section
- ratio
- parameter water
- 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
Classifications
-
- 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
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)
- Testing Or Calibration Of Command Recording Devices (AREA)
Abstract
本发明涉及一种多参数水质剖面仪的pH校正方法,包括下列步骤:数据获取:利用多参数水质剖面仪获取观测海域各站点的pH剖面S;数据质量控制:通过剔除仪器测定过程中的感温数据、逆压校正数据、计算局部标准差去除大振幅奇异值等步骤,对S进行质量控制,质量控制结果用Q表示;数据平均:对Q中各站位下行数据进行滑动平均,处理结果用A表示;数据校正:以经过pH标准缓冲溶液校准过的台式pH计测定的海水pH值为基准,对A进行校正处理。
Description
技术领域
本发明涉及海洋科学领域,具体涉及一种多参数水质剖面仪(简称RBR)的pH校正方法。
背景技术
pH表示水中酸性物质氢离子(H+)浓度的负对数(pH=log10[H+]),其作为水环境健康的重要指标,对海洋生物的生存和繁衍具有重要意义。关于海水pH值的测定,常规化学法多利用采水器采集水样,并使用台式pH计测定海水的pH值。由于台式pH计的探头价格较低、易于更换,故测量结果更加精准,但该方法耗时耗力且只能获取单点pH数据,具有较大的局限性。多参数水质剖面仪的问世为获取高频且连续的剖面数据提供了可能,目前已成为国际上应用最广泛的海洋调查仪器设备之一。该仪器包含温度、电导(盐度)、压力(深度)、溶解氧、有效光合作用辐照(PAR)、叶绿素、浊度、pH等多个探头,最多可以同时测定13个水下参数。由于其体积小巧,非常适合于近海水体快速调查、大面积调查、水质剖面观测、走航观测、拖体安装观测、浮标长期观测、海底长期监测等,但仪器长时间工作或位于温盐跃层处测定时,可能存在数据漂移或响应滞后现象。因此,现场调查过程中可以将两种方法相结合,以台式pH计测定的海水pH值为基准,对多参数水质剖面仪观测的pH 剖面进行校正,两者优势互补能够获取更高频且精确的pH剖面。
关于多参数水质剖面仪的pH校正,前人多采用“线性回归法”,而考虑到多参数水质剖面仪pH探头的实际测量误差,高pH更高,低pH更低,即测定过程中斜率发生变化,故该方法的相对误差较大,因此本发明采用“比值法”对多参数水质剖面仪观测的pH数据进行校正,可获取精度更高、质量更好的pH剖面,为研究海水理化性质提供良好的数据支持。
发明内容
有鉴于此,本发明提供了一种多参数水质剖面仪的pH校正方法,改善因仪器响应随时间漂移产生的测定误差,获取更高精度的pH剖面,为研究海水的理化性质提供更高频且精确的数据基础。本发明的技术方案如下:
一种多参数水质剖面仪的pH校正方法,包括下列步骤:
1)数据获取:利用多参数水质剖面仪获取观测海域各站点的pH剖面S;
2)数据质量控制:通过剔除仪器测定过程中的感温数据、逆压校正数据、计算局部标准差去除大振幅奇异值等步骤,对S进行质量控制,质量控制结果用Q表示;
3)数据平均:对Q中各站位下行数据进行滑动平均,处理结果用A表示;
4)数据校正:以经过pH标准缓冲溶液校准过的台式pH计测定的海水pH值为基准,对A进行校正处理,方法如下:
第1步:去除跃层数据,利用多参数水质剖面仪同步观测的温度、盐度、声速、水深数据,计算各站位的浮力频率N,并根据logN2>-3确定海水中温度跃层的位置,去除A中跃层处的pH数据,处理结果用T表示,浮力频率的计算公式为:
式中,c0为声速,ρ为海水密度,z为水深,g为重力加速度,g2/co 2一项量级甚小,忽略不计;
第2步:计算台式pH计测定的海水pH值与对应T中pH的比值,并将该比值定义为Ratio,即有:
式中,pHb为台式pH计测定的海水pH值,pHr指对应多参数水质剖面仪测定的海水pH值;
第3步:假定Ratio近似服从正态分布,计算Ratio的概率密度函数,并对离散较大的 Ratio进行1~3倍标准偏差剔除,处理结果用R表示;
第4步:对R与对应T中pH观测数据进行线性拟合,计算观测海域各站位pH剖面的校正系数,方法如下:令y=ax+b,参数a和b用来调整校正曲线的斜率和截距,x表示多参数水质剖面仪观测的pH剖面A,即包含跃层在内的全部pH数据,y表示不同pH对应的校正比例;
第5步:将A和y一一对应相乘,获得完整的pH校正剖面。
本发明的多参数水质剖面仪的pH校正方法,包括数据获取、数据质量控制、滑动平均处理和数据校正等步骤。首先用多参数水质剖面仪获取观测海域各站点的pH剖面;接下来进行数据质量控制,包括去除感温数据、逆压校正数据和去尖峰处理三个步骤;然后使用“滑动平均法”对质量控制后的pH剖面进行平均输出;最后以台式pH计测定海水的pH值为基准,对多参数水质剖面仪观测的pH平均输出剖面进行校正。本发明方法的特征有:1)多参数水质剖面仪的采样频率为6Hz,可以获取较为丰富的pH剖面;2)通过pH校正步骤,可以改善因仪器响应随时间漂移产生的测定误差,有效地保证数据质量,为研究海水理化性质提供更高频且精确的数据支持。
附图说明
图1是pH校正方法流程图
图2是部分实验结果,其中:
(a)B09站位的跃层剖面图
(b)Ratio的概率密度函数及频率分布直方图
(c)Ratio与pHr的线性相关曲线
(d)pHb与pHr的线性相关曲线
图3是2017年黄渤海夏季共享航次校正前后底层pH的大面分布图,其中:
(a)校正前 (b)校正后
具体实施方式
本发明主要包括:数据获取、数据质量控制、数据平均和数据校正等步骤。图1给出了所提方法的框图。包括:
1、数据获取:利用多参数水质剖面仪获取观测海域各站点的pH剖面,并用S表示。仪器在布放前已使用pH=4.01,7.00和9.21的NIST(National Institute of Standardsand Technology,USA)标准缓冲溶液在室温下对pH探头进行了校准;布放过程中保证布放船侧迎风,均匀下放,以便仪器稳定测量。多参数水质剖面仪的采样频率为6Hz,pH探头测定范围1~13,精度为0.1,可获取较为丰富的pH剖面;
2、数据质量控制:为了排除观测过程中由于船体运动、绞车钢缆摆动以及仪器探头不稳定性带来的外界干扰,需要对多参数水质剖面仪测定的原始pH剖面S进行质量控制,质量控制结果用Q表示。具体过程如下:
第1步:去除感温数据。多参数水质剖面仪入水后需在表层(~1m)停留约1min,完成仪器感温过程,此过程多为虚假数据,应予以去除;
第2步:逆压校正数据。当Pi-Pi-1≤0,剔除Pi水深处的pH数据,式中,i为某一站点的数据序列,Pi为i处的压力,Pi-1为i-1处的压力;
第3步:去尖峰处理。计算pH剖面在对应深度处的一阶导数及导数的局部标准差,将两者比值超过某一设定阈值的信号作为尖峰去除,并用尖峰两端信号的平均值替代。原理如下:
式中,i为某一站点的数据序列,zi为对应位置的水深,xi为多参数水质剖面仪测得的对应位置的pH数据,Di为对应水深的一阶导数,s为对应水深一阶导数的标准差。使得其中lim为变量,将两者比值超过lim这一阈值的信号作为尖峰去除,并用尖峰两端的平均信号代替。
3、数据平均:为了精确地表示测量结果,消除偶然变动引起的随机误差,常对动态测试数据作平滑和滤波处理。本专利采用“滑动平均法”对Q进行处理(处理结果用A表示),滤除0.2m以下的物理现象,以便下一步研究分析。滑动平均法的基本原理:
yj=fj+ej j=1,2,...,N
式中,yj为动态测试数据,fi为较平滑的测量结果,ei为随机误差。
4、数据校正:为改善因仪器响应随时间漂移产生的测定误差,获取更高精度的pH剖面,所提方法以台式pH计测定的海水pH值为基准,对多参数水质剖面仪观测的pH剖面A进行校正。具体过程如下:
第1步:去除跃层数据。温盐跃层附近由于水温、盐度突然而急剧变化,多参数水质剖面仪可能存在响应滞后现象,因此,首先需要对各站位跃层附近的pH数据进行去除。
利用多参数水质剖面仪同步观测的温度、盐度、声速、水深等数据,计算各站位的浮力频率N,并根据logN2>-3确定海水中温度跃层的位置,去除A中跃层处的pH数据,处理结果用T表示。浮力频率的计算公式为:
式中,c0为声速,ρ为海水密度,z为水深,g为重力加速度。g2/co 2一项量级甚小,可忽略不计。
第2步:计算台式pH计测定的海水pH值与对应T中pH的比值,并将该比值定义为Ratio,即有:
式中,pHb为台式pH计测定的海水pH值,pHr指对应多参数水质剖面仪测定的海水pH值。
第3步:假定Ratio近似服从正态分布,计算Ratio的概率密度函数,并对离散较大的 Ratio进行1~3倍标准偏差剔除,处理结果用R表示。基于上述假设,将有68%数值分布在距离平均值有1个标准差之内的范围,约95%数值分布在距离平均值有2个标准差之内的范围,约99.7%数值分布在距离平均值有3个标准差之内的范围。概率密度函数的计算公式为:
记作X~N(μ,σ2),
其中,x为第2步中定义的比值Ratio,μ为Ratio的期望,σ为Ratio的标准差,σ2为Ratio的方差。
第4步:对R与对应T中pH观测数据进行线性拟合,计算观测海域各站位pH剖面的校正系数。令y=ax+b,参数a和b用来调整校正曲线的斜率和截距,x表示多参数水质剖面仪观测的pH剖面A,y表示不同pH对应的校正比例。本专利以2017年黄渤海夏季共享航次为例,对离散较大的Ratio作2倍标准偏差剔除处理,拟合所得一元线性曲线为:
y=-0.0991x+1.8037 R2=0.8225
由R=0.9069可以看出,回归曲线中Ratio与pHr的相关性较好,故采用线性拟合法适用于多参数水质剖面仪的pH剖面校正。
第5步:将A和y一一对应相乘,获得完整的pH校正剖面。即:
Y=A×y
式中,A为校正前多参数水质剖面仪观测的pH剖面,Y为校正后多参数水质剖面仪观测的pH剖面。
对校正前后的两套数据进行线性拟合,以检验pH校正前后是否改变调查海域pH的空间分布,拟合曲线如图2d所示;图3显示了2017年渤黄海夏共享航次底层pH的大面分布,由此可以看出,校正前后pH的空间分布并没有发生改变,但数据质量明显提高,更符合正常海水pH约为7.9~8.1的科学规律。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种多参数水质剖面仪的pH校正方法,包括下列步骤:
1)数据获取:利用多参数水质剖面仪获取观测海域各站点的pH剖面S;
2)数据质量控制:通过剔除仪器测定过程中的感温数据、逆压校正数据、计算局部标准差去除大振幅奇异值等步骤,对S进行质量控制,质量控制结果用Q表示;
3)数据平均:对Q中各站位下行数据进行滑动平均,处理结果用A表示;
4)数据校正:以经过pH标准缓冲溶液校准过的台式pH计测定的海水pH值为基准,对A进行校正处理,方法如下:
第1步:去除跃层数据,利用多参数水质剖面仪同步观测的温度、盐度、声速、水深数据,计算各站位的浮力频率N,并根据logN2>-3确定海水中温度跃层的位置,去除A中跃层处的pH数据,处理结果用T表示,浮力频率的计算公式为:
式中,c0为声速,ρ为海水密度,z为水深,g为重力加速度,g2/co 2一项量级甚小,忽略不计;
第2步:计算台式pH计测定的海水pH值与对应T中pH的比值,并将该比值定义为Ratio,即有:
式中,pHb为台式pH计测定的海水pH值,pHr指对应多参数水质剖面仪测定的海水pH值;
第3步:假定Ratio近似服从正态分布,计算Ratio的概率密度函数,并对离散较大的Ratio进行1~3倍标准偏差剔除,处理结果用R表示;
第4步:对R与对应T中pH观测数据进行线性拟合,计算观测海域各站位pH剖面的校正系数,方法如下:令y=ax+b,参数a和b用来调整校正曲线的斜率和截距,x表示多参数水质剖面仪观测的pH剖面A,即包含跃层在内的全部pH数据,y表示不同pH对应的校正比例;
第5步:将A和y一一对应相乘,获得完整的pH校正剖面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910177098.0A CN110057990B (zh) | 2019-03-08 | 2019-03-08 | 一种多参数水质剖面仪的pH校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910177098.0A CN110057990B (zh) | 2019-03-08 | 2019-03-08 | 一种多参数水质剖面仪的pH校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110057990A true CN110057990A (zh) | 2019-07-26 |
CN110057990B CN110057990B (zh) | 2021-09-14 |
Family
ID=67316743
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910177098.0A Expired - Fee Related CN110057990B (zh) | 2019-03-08 | 2019-03-08 | 一种多参数水质剖面仪的pH校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110057990B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111680054A (zh) * | 2020-06-08 | 2020-09-18 | 中国科学院地理科学与资源研究所 | 基于多参数水质监测仪的时序数据存储方法及装置 |
CN117405975A (zh) * | 2023-12-14 | 2024-01-16 | 深圳鹏城新能科技有限公司 | 一种pv面板绝缘电阻检测方法、系统及介质 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7348783B1 (en) * | 2006-08-18 | 2008-03-25 | Chung Yuan Christian University | Multi-functional pH meter and fabrication thereof |
JP2011220980A (ja) * | 2010-04-14 | 2011-11-04 | Central Res Inst Of Electric Power Ind | 試料液の蛍光染料濃度の測定方法 |
CN103542875A (zh) * | 2012-07-09 | 2014-01-29 | 恩德莱斯和豪瑟尔测量及调节技术分析仪表两合公司 | 用于并行校准至少两个传感器的装置和方法 |
CN107607605A (zh) * | 2017-08-08 | 2018-01-19 | 安吉元融仪器仪表检测有限公司 | 精确的pH计的校准方法 |
CN108368471A (zh) * | 2015-10-30 | 2018-08-03 | 豪夫迈·罗氏有限公司 | Ph测量装置的校准偏差的鉴别 |
CN109060611A (zh) * | 2018-07-09 | 2018-12-21 | 天津大学 | 一种基于激光粒度仪的悬沙质量浓度数据处理方法 |
CN109141380A (zh) * | 2018-09-19 | 2019-01-04 | 天津大学 | 一种基于ctd仪的水文数据处理方法 |
-
2019
- 2019-03-08 CN CN201910177098.0A patent/CN110057990B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7348783B1 (en) * | 2006-08-18 | 2008-03-25 | Chung Yuan Christian University | Multi-functional pH meter and fabrication thereof |
JP2011220980A (ja) * | 2010-04-14 | 2011-11-04 | Central Res Inst Of Electric Power Ind | 試料液の蛍光染料濃度の測定方法 |
CN103542875A (zh) * | 2012-07-09 | 2014-01-29 | 恩德莱斯和豪瑟尔测量及调节技术分析仪表两合公司 | 用于并行校准至少两个传感器的装置和方法 |
CN108368471A (zh) * | 2015-10-30 | 2018-08-03 | 豪夫迈·罗氏有限公司 | Ph测量装置的校准偏差的鉴别 |
CN107607605A (zh) * | 2017-08-08 | 2018-01-19 | 安吉元融仪器仪表检测有限公司 | 精确的pH计的校准方法 |
CN109060611A (zh) * | 2018-07-09 | 2018-12-21 | 天津大学 | 一种基于激光粒度仪的悬沙质量浓度数据处理方法 |
CN109141380A (zh) * | 2018-09-19 | 2019-01-04 | 天津大学 | 一种基于ctd仪的水文数据处理方法 |
Non-Patent Citations (1)
Title |
---|
张开远 等: "基于最小二乘法的pH值温度补偿系统设计", 《传感器与微系统》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111680054A (zh) * | 2020-06-08 | 2020-09-18 | 中国科学院地理科学与资源研究所 | 基于多参数水质监测仪的时序数据存储方法及装置 |
CN111680054B (zh) * | 2020-06-08 | 2023-08-04 | 中国科学院地理科学与资源研究所 | 基于多参数水质监测仪的时序数据存储方法及装置 |
CN117405975A (zh) * | 2023-12-14 | 2024-01-16 | 深圳鹏城新能科技有限公司 | 一种pv面板绝缘电阻检测方法、系统及介质 |
CN117405975B (zh) * | 2023-12-14 | 2024-03-22 | 深圳鹏城新能科技有限公司 | 一种pv面板绝缘电阻检测方法、系统及介质 |
Also Published As
Publication number | Publication date |
---|---|
CN110057990B (zh) | 2021-09-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110146895B (zh) | 基于倒置式多波束回声仪的声速剖面反演方法 | |
Mackas et al. | Least squares multiple tracer analysis of water mass composition | |
CN105335336B (zh) | 一种传感器阵列的稳健自适应波束形成方法 | |
RiSeR et al. | Salinity in Argo: A modern view of a changing ocean | |
CN110673108B (zh) | 一种基于迭代Klett的机载海洋激光雷达信号处理方法 | |
CN110057990A (zh) | 一种多参数水质剖面仪的pH校正方法 | |
CN115342814B (zh) | 一种基于多传感器数据融合的无人船定位方法 | |
CN115327518A (zh) | 一种基于解析模型的海洋激光雷达多次散射衰减校正方法 | |
Shao et al. | Verification of echosounder measurements of thickness and spatial distribution of kelp forests | |
CN109060611B (zh) | 一种基于激光粒度仪的悬沙质量浓度数据处理方法 | |
CN113723182B (zh) | 一种训练样本受限条件下的sar图像舰船检测方法 | |
CN116864026A (zh) | 一种基于Argo温盐剖面的海洋溶解氧浓度重构模型构建方法 | |
笹野大輔 et al. | Testing a new quick response oxygen sensor," RINKO" | |
RU2769550C1 (ru) | Способ достроения измеренной от поверхности моря части профиля вертикального распределения скорости звука до дна | |
CN114818232B (zh) | 一种海水空间温度剖面的重构方法及系统 | |
CN109471113A (zh) | 基于相位法的多波束声呐海底地形测量质量实时评估方法 | |
Beaudoin et al. | Streamlining sound speed profile pre-processing: Case studies and field trials | |
RU2747739C1 (ru) | Методика верификации измерений солености морской воды автоматическими средствами измерений | |
CN112114367A (zh) | 一种利用主动电场探测海域中电学分界面的方法 | |
Fleagle | Note on the effect of air‐sea temperature difference on wave generation | |
CN104777526B (zh) | 一种ascat反演风速的订正方法 | |
CN110568417A (zh) | 基于多物理场的机载激光测深回波信号的处理方法 | |
CN118114030B (zh) | 一种基于单点单深度海流计数据的内孤立波参数提取方法 | |
Ubaidillah et al. | Architecture tools measure the levels of salt and ph of seawater using a fuzzy logic-based android | |
Fujiki et al. | Improvement in estimation stability and accuracy of wave directional spectrum by a Bayesian method for swell observation |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210914 |