CN101587150A - 溶液电导率的测量方法 - Google Patents
溶液电导率的测量方法 Download PDFInfo
- Publication number
- CN101587150A CN101587150A CNA2009101394391A CN200910139439A CN101587150A CN 101587150 A CN101587150 A CN 101587150A CN A2009101394391 A CNA2009101394391 A CN A2009101394391A CN 200910139439 A CN200910139439 A CN 200910139439A CN 101587150 A CN101587150 A CN 101587150A
- Authority
- CN
- China
- Prior art keywords
- electrode
- solution
- cable
- voltage
- 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.)
- Granted
Links
Images
Landscapes
- Measurement Of Resistance Or Impedance (AREA)
Abstract
本发明公开一种溶液电导率的测量方法,测量系统包括二极或四极电导电极、测量仪表以及连接电缆,此方法包括以下步骤:采用电流幅值一定、频率一定的交流方波电流通过电缆对电极进行激励;对电极响应的电压波形进行高速A/D变换,在波形上取得三个不同时刻的电压值或三个不同时间区段的三个电压平均值;根据已知的电极响应电压波形函数 联立方程,求得三个独立参数R,B,D,其中R是待测电极之间的电阻值,B=CP/(Cs+CP),,T=(CsCP/(Cs+CP))R,Cp代表电极极化的等效电容,Cs代表电缆的等效电容、测量仪表输入电路的电容以及电极的分布电容之和,TH代表所述交流方波电流的半周期;再结合电极常数由R获得待测的溶液电导率值。
Description
本申请是申请号为200610030555.6,申请日为2006年8月30日的名称为“溶液电导率的测量方法”的中国专利申请的分案申请
技术领域
本发明涉及溶液电导率或电阻率的测量方法,尤其涉及能消除电极极化以及长电缆影响的溶液电导率或电阻率的测量方法。
背景技术
溶液电导率的基本测量方法是将电极置入溶液,通过电极向溶液通电,测量电极两端的电压U以及流过的电流I,计算电极之间的电阻R=U/I,用公式G=C/R来计算溶液电导率,其中C为电极常数,电阻率为电导率的倒数。但是溶液中的电极通电后会产生极化,U/I会随着时间而改变,因此,会产生测量误差。防止极化引起测量误差的基本方法是通过交流电进行测量,当测量系统的连接电缆很短时,频率越高,精度越高。但是连接电缆总存在一定的等效电容,频率越高,该电容将引起越大的误差。因此两者是矛盾的。
中国专利申请号为CN 200410066147.7的文件公布一种测量方法,电极被二个频率的正弦信号进行激励,分别求得二个阻抗模|Za|和|Zb|,以及r=|Za|/|Zb|然后利用
求得溶液电导率值。该方法能消除极化的影响,但不能同时消除电极极化以及长连接电缆的影响。
中国专利号为ZL02111820.5的专利介绍一种用有用功的概念进行测量的方法,测量溶液的电压和电流,利用公式
G=C*∫I2dt/∫U*Idt
求得溶液电导率,它也能消除极化的影响,但是方法本身并不涉及消除长电缆影响的问题。
文件DE4233110A1公布了一种找到最优激励频率的方法,它用二个相邻的频率(频率约相差20%)进行测量,如果二次结果相差很小,则认为所用频率是合适的,否则改变频率,重新测量,直到找到最优频率,并且得到误差的估计。但是测量方法本身是传统的单频率方法。文件US6369579B1公布的一种方法,介绍了电极极化的模型,据此模型,发明了一种至少用二个频率来激励电极并求出电极之间的电阻值R和溶液电导率的方法,虽然该模型是一个较复杂的电极模型,方法中是利用复阻抗计算来求解的,是一种谐波激励法。但模型中没有特别考虑长电缆的影响。所以不太适合于有长电缆时的情况。
在溶液电导率测量中,给电极的激励可以分为交流方波电压、交流方波电流或正弦波激励,在激励电路或输入运放测量电路等不会饱和的情况下,交流方波电流激励法能用较低的频率测得精确的结果。但是频率不能太低,否则由于极化的原因,交流方波电流激励法会导致激励电路或输入运放电路的饱和,引起测量误差。由于在实际测量中,电缆等效电容和电极极化总是同时存在的,因此不论用何种频率进行测量,用通常的方法,总有误差存在。
发明内容
本发明所要解决的技术问题是提供一种能避免电极极化,电缆电容对测量的不利影响,获得溶液准确的电导率或电阻率值的测量方法。
本发明提出一种溶液电导率的测量方法,当连接电缆及电极在某一频率的交流方波电流的激励下,可以得到确定的响应电压波形函数
它包含三个独立的参数R,B,和D,这三个参数与电极极化的等效电容Cp,电缆的等效电容、测量仪表输入电路的电容和电极的分布电容之和Cs以及待测电极之间的电阻值R相关。这三个参数不同,响应电压波形就不同。在响应电压波形函数V上取三个不同时刻的值,或三个不同时间区段的三个电压平均值,可以证明,它们是相互独立的,则用三个电压值可以联立三个方程,可以求解得到三个未知数R,B,和D。其中R为待测量。这种方法不需要用双频来测量,但要用三个电压值,所以可以称为三电压法。
本发明的溶液电导率的测量方法相比已有的测量方法具有如下有益效果:能够用比较简单的电路和方法,即使在电极与测量装置之间的连接电缆很长,并且电极存在一定的极化时,也能精确测量溶液电导率。
附图说明
以下结合附图和具体实施方式进一步说明本发明的特征和优点。
图1是本发明溶液电导率的测量方法的第一实施例采用二极电极二线制时的原理示意图;
图2是本发明溶液电导率的测量方法的第一实施例采用二极电极四线制时的原理示意图;
图3是本发明溶液电导率的测量方法的第一实施例采用四极电极四线制时的原理示意图;
图4是本发明溶液电导率的测量方法的第一实施例的原理图的等效示意图;
图5是本发明溶液电导率的测量方法的第一实施例的步骤流程示意图;
图6是电流电压波形示意图;
图7是本发明溶液电导率的测量方法第二实施例的原理示意图;
图8是本发明溶液电导率的测量方法的第二实施例的步骤流程示意图。
具体实施方式
图1~图3示出本发明溶液电导率的测量方法的第一实施例的原理示意图。请先参阅图1,测试系统包括插入待测溶液1的二极电极2、测量仪表4和连接电缆3。测量仪表4进一步包括交流方波电流激励单元5、同步检波器6、运算单元7和双频控制单元8。图1中的电极采用二极电极二线制接法,交流方波电流激励单元5通过连接电缆对一个电极输出交流方波激励电流,同步检波器6从连接电缆上获得电极响应电压。然而需要指出的是,也可采用如图2中的二极电极四线制接法,即电极102通过四线电缆103与测量仪表104连接,同步检波器6通过运算放大器106从两个电极之间获得电极响应电压;由于运算放大器106的输入阻抗很高,联接至运算放大器106的二根导线上的电流可以忽略不计,又由于交流方波电流激励相当于恒流源激励,则长电缆线的导线电阻的影响可以忽略不计。或者也可采用如图3中的四极电极四线制接法,采用一个四极电极202与测量仪表104连接。这种接法,除四线制接法能消除长电缆的导线电阻的影响外,还能在一定程度上减少电极极化的影响。这些接法采用的测量方法可以是一样的,即以下所要说明的双频法,或三电压法。
图4为本发明溶液电导率的测量方法的第一实施例的等效图。图4中,R是待测电极之间的电阻值,R=C/G,G为溶液电导率,C为电极常数。电容Cp代表电极极化的等效电容;电容Cs代表电缆的等效电容、仪表输入电路的电容以及电极的分布电容之和。电流I是等幅的电流,即IH或IL;频率FH,FL代表二个激励频率。电压V代表电极的响应电压波形。Vr代表同步检波后的电压波形。
本实施例采用双频法进行测量,以下先近似说明双频法的基本原理。
在一定频率的交流方波电流激励下进行测量,由于长电缆引起的误差主要是在电流反向的过渡过程中产生的,过渡过程的时间长度与激励的半周期之比与测量误差成正比,即使激励的频率不同,过渡过程的时间长度在二次测量中几乎是一样的,但是二次测量中的激励的半周期不同,所以频率较低的误差较小。
设过渡过程时间为t1;频率FH的半周期为 RH为交流方波激励电流频率为FH时测得的视在电阻值;频率FL的半周期为 RL是交流方波激励电流频率为FH时测得的视在电阻值;TL/TH=FH/FL=n为频率比,因此有:
RH=R*(TH-t1)/TH,RL=R*(TL-t1)/TL
则:
如果n=2,则:R=2RL-RH。
以上只是近似解释,但是能很好地理解双频测量法的本质,下面将结合图4所示的等效电路模型、图5所示的测量方法的第一实施例的步骤流程图和图6所示的电流电压波形示意图,详细说明如何求得精确的电极之间的电阻值R。本实施例中的方法包括:
步骤S1:采用某一定频率的交流方波电流通过电缆3(103)对电极2(102,202)进行激励;对如图4所示等效电路进行分析,可以推导得到电极响应的电压波形函数V:
在测量的短时间内,假设电极之间的电阻值R,电容Cs,Cp不变。假设激励为幅值IH,频率为FH的交流方波电流。在电流的正半周(t=0到TH),i=IH,根据电容电流定律:
根据基尔霍夫电流定律:
将基尔霍夫电压定律用于电容Cp二端:
设 可以得到:
微分方程的通解是:
假设,电极及激励处于稳定状态,则:
V(t=TH-)=-V(t=0+)
Vp(t=TH-)=-Vp(t=0+)
设 可以得到:
所以得到电极的响应电压波形函数:
步骤S2:对电极的响应电压波形进行同步检波,求出第一平均电压值,除以电流幅值,求得第一视在电阻值RH;具体公式如下:
对电压波形的某一时间区段t=T1~T2进行平均:
如果对整个半波进行平均(T1=0,T2=TH),则第一平均电压值
通过同步检波,第二个半波(t=TH~2TH)的平均值与第一个半波(t=TH~2TH)的平均值是一样的,即与同步检波后的总电压平均值是一样的。定义第一视在电阻值:
步骤S3:采用另一频率的交流方波电流通过电缆3对电极2进行激励;设激励电流是幅值为IL,频率为FL的交流方波电流,如果它与第一交流方波电流之间的关系为FH/FL=n;TL/TH=n,同样可以推导得到电极响应的电压波形函数为:
步骤S4:对电极的响应电压波形进行同步检波,求出第二平均电压值,除以电流幅值,求得第二视在电阻值RL。
经过同步检波后的平均电压为:
定义第二视在电阻值:
步骤S5:根据第一、第二视在电阻值RH、RL以及频率FH、FL的比值n计算获得电极之间的电阻值R,再结合电极常数获得待测的溶液电导率。具体计算如下:
视在电阻值的相对偏差:
设:R=PnwRL+(1-Pnw)RH
通常在Pnw算式中,假设 已足够精确。则Pnw是n和w的函数:先确定n,Pn(w)函数关系可以事先用数值法,分段多项式拟合求得,具体的做法是,列出一系列的D值,根据各个D值分别求出一系列w和Pn(w)的值,利用这些值作分段多项式拟合。在测量时,先求二个视在电阻值RH、RL的相对偏差w=(RL-RH)/RH,再算出Pn(w),求得电极之间的电阻值R。利用公式G=C/R求出溶液电导率G,其中C是电极常数。
在本方法中,原理上只要频率比n不等于1即可。实际上为了减少随机干扰等影响,频率比n在1.2到4之间较好,或者其倒数,即从0.25到0.84之间较好。以下举二个频率比n的例子:
第一个例子,如果n=2:
当w>0.53时:Pnw=(0.038599w4-0.309326w3+0.950803w2-1.422355w+5.761831)/(3-w)
当0.11<w≤0.53时:Pnw=-1.566025w4+1.885248w3-0.042810w2-0.059616w+2.004903
当w≤0.11时,Pnw=2
特别地当n=2, 时,R=2RL-RH。
第二个例子,如果n=3/2=1.5:
当w>0.23时:Pnw=(2.472277w4-6.618677w3+7.091540w2-3.932823w+3.766115)/(1.25-w)
当0.067<w≤0.23时:Pnw=-65.905750w4+46.746107w3-6.836475w2+0.370457w+2.993219
当w≤0.067时,Pnw=3
特别地当n=1.5, 时,R=3RL-2RH。
上述公式的采用可以根据实际精度需要,适当简化。由于Pnw的计算不准确引起的理论相对误差可以分析如下:设Pnw变成了Pnw1,引起的相对误差为:
通过数值分析,可知在第一个例子(n=2)中,如果使用上述给出的公式,当w<0.53时,由于Pnw的计算不准确引起的理论相对误差Err<0.01%;在第二个例子(n=1.5)中,使用上述给出的公式,当w<0.76时,由于Pnw的计算不准确引起的理论相对误差Err<0.01%。
当第一、第二视在电阻值RH、RL的相对偏差w=(RL-RH)/RH小于一参考值时,Pnw=n/(n-1)是一个很好的公式。通过数值分析,当n=1.2到4时,可以得到以下估计:
如果w<-0.0209*n2+0.194*n-0.151(第一参考值),或w<0.15|n=2,或w<0.09|n=1.5,由于Pnw=n/(n-1)引起的理论相对误差为Err<0.02%;
如果w<-0.0238*n2+0.247*n-0.197(第二参考值),或w<0.20|n=2,或w<0.12|n=1.5,由于Pnw=n/(n-1)引起的理论相对误差为Err<0.1%。
可见参考值的取值可依计算精确度要求而定,参考值的取值越小,采用公式Pnw=n/(n-1)所引起的理论相对误差也越小。
应当指出,在测量过程中,使IH=IL=I也是完全可以的,而且可见某些公式的表述还可能简化一些。
在本实施例中,同步检波器6可以用硬件电路实现,例如多路切换开关,由于在双频法中,公式中只用到电压平均值,可以在同步检波器之后,A/D变换之前插入低通滤波器,这样对A/D变换的速度要求较低。或者同步检波器6也可以先对电极响应的电压波形进行高速A/D变换,再由软件来实现。
如果A/D变换的速度足够高,也可以只用单一频率的交流方波电流来激励,用以下方法完成测量,同样能消除电极极化以及长电缆对测量精度的影响。
图7示出本发明溶液电导率的测量方法第二实施例的原理示意图;图8是本发明溶液电导率的测量方法的第二实施例的步骤流程示意图。测试系统包括交流方波电流激励单元5、频率控制单元408、以及运算单元407,请结合图7和图8所示,本实施例的方法包括以下步骤:
步骤S11:采用电流幅值一定、频率一定的交流方波电流通过电缆对电极进行激励;设激励为幅值IH,频率FH的交流方波电流,在正半周,如前文所述电极响应电压波形函数为:
该电压波形有三个独立的变量D、B、R,可以用三个不同时刻的电压,或者三个不同时间区段的平均电压来求得。
步骤S12:对电极响应的电压波形进行高速A/D变换,在半波波形上取得三个不同时刻的电压值或三个不同时间区段的三个电压平均值;例如,将电压波形进行A/D变换,取得波形V中的三个时刻t=(1/4)*tH、t=(1/2)*tH和t=(3/4)*tH的三个电压V1、V2、V3。
步骤S13:根据已知的电极响应电压波形函数
其中R是待测电极之间的电阻值, TH是频率FH的半周期。根据联立方程求得三个独立参数R,B,D,再根据公式G=C/R求出溶液电导率G。
以下列举一种解法:
由B1可以求得D:
由B2,D不难求得B:
其中 然后,
利用公式G=C/R求出溶液电导率G,其中C是电极常数。
或者对三个不同的时间区段,例如t=[(1/8)~(3/8)]*tH、t=[(3/8)~(5/8)]*tH和t=[(5/8)~(7/8)]*tH进行平均以便减少随机干扰的影响,得到Va1,Va2,Va3。在用Va1,Va2,Va3求解时,用前述的Va函数式。
为了减少随机干扰的影响,V1、V2、V3或Va1,Va2,Va3可以用多个电压波形的同相位点的平均值来进行计算。符号反向后,负半周电压波形上的电压值也可以一起用于计算。
用V1、V2、V3或者Va1,Va2,Va3联立方程,可以求出三个未知数R、B、D,其中R为所需测量的电极之间的电阻值。这种方法可以称为三电压测量法。
应当指出,以上说明中提到的双频法,或单频三电压法。是对某一次测量而言的。即对某一次测量,用某一对双频进行激励,或某一个频率进行激励。本领域技术人员都知道,对不同电导率范围的溶液进行测量时,比较合适的激励频率是不同的。通常溶液电导率越低,宜用越低的激励频率,越小的激励电流。另外,由于溶液电导率通常与温度相关,因此在溶液电导率传感器(电极)中,通常还可以有温度传感器,所以在连接传感器与测量装置的电缆线中,另外还可以有连接温度传感器的导线,测量装置中也会有测量温度的电路。这些基本原理也适合于本发明的方法中。
以上说明中所用的术语,符号,公式和例子并不对该发明的应用构成限制,只是为了便于它们的说明。
以上的实施例说明仅为本发明的较佳实施例说明,本领域技术人员可依据本发明的上述实施例说明而作出其它种种等效的替换及修改。然而这些依据本发明实施例所作的种种等效替换及修改,属于本发明的发明精神及由权利要求所界定的专利范围内。
Claims (1)
1.溶液电导率的测量方法,测量系统包括二极或四极电导电极、测量仪表以及连接电缆,其特征在于,所述方法包括以下步骤:
a)采用电流幅值一定、频率一定的交流方波电流通过电缆对电极进行激励;
b)对电极响应的电压波形进行高速A/D变换,在波形上取得三个不同时刻的电压值或三个不同时间区段的三个电压平均值;
c)根据已知的电极响应电压波形函数
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101394391A CN101587150B (zh) | 2006-08-30 | 2006-08-30 | 溶液电导率的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101394391A CN101587150B (zh) | 2006-08-30 | 2006-08-30 | 溶液电导率的测量方法 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2006100305556A Division CN100541208C (zh) | 2006-08-30 | 2006-08-30 | 溶液电导率的测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101587150A true CN101587150A (zh) | 2009-11-25 |
CN101587150B CN101587150B (zh) | 2011-04-13 |
Family
ID=41371479
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009101394391A Active CN101587150B (zh) | 2006-08-30 | 2006-08-30 | 溶液电导率的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101587150B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101726518B (zh) * | 2009-12-24 | 2012-09-26 | 福建师范大学 | 溶液电导率的测量方法 |
CN102809699A (zh) * | 2012-08-24 | 2012-12-05 | 福建师范大学 | 溶液电导率测量涉及的电极分布电容动态测定方法 |
CN105588983A (zh) * | 2014-11-14 | 2016-05-18 | 佛山市顺德区美的电热电器制造有限公司 | 电导率的测试装置和用电设备 |
CN107543952A (zh) * | 2016-06-28 | 2018-01-05 | 深圳市水务科技有限公司 | 测量溶液中带电体和溶液组成的回路中微小电流的方法 |
CN110057874A (zh) * | 2010-12-17 | 2019-07-26 | 通用电气健康护理生物科学股份公司 | 预测液体混合物的电导率的方法 |
-
2006
- 2006-08-30 CN CN2009101394391A patent/CN101587150B/zh active Active
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101726518B (zh) * | 2009-12-24 | 2012-09-26 | 福建师范大学 | 溶液电导率的测量方法 |
CN110057874A (zh) * | 2010-12-17 | 2019-07-26 | 通用电气健康护理生物科学股份公司 | 预测液体混合物的电导率的方法 |
US11670402B2 (en) | 2010-12-17 | 2023-06-06 | Cytiva Sweden Ab | Method for predicting the conductivity of a liquid mixture |
CN102809699A (zh) * | 2012-08-24 | 2012-12-05 | 福建师范大学 | 溶液电导率测量涉及的电极分布电容动态测定方法 |
CN102809699B (zh) * | 2012-08-24 | 2014-09-03 | 福建师范大学 | 溶液电导率测量涉及的电极分布电容动态测定方法 |
CN105588983A (zh) * | 2014-11-14 | 2016-05-18 | 佛山市顺德区美的电热电器制造有限公司 | 电导率的测试装置和用电设备 |
CN107543952A (zh) * | 2016-06-28 | 2018-01-05 | 深圳市水务科技有限公司 | 测量溶液中带电体和溶液组成的回路中微小电流的方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101587150B (zh) | 2011-04-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100541208C (zh) | 溶液电导率的测量方法 | |
US5461321A (en) | Apparatus and method for measuring capacitance from the duration of a charge-discharge charge cycle | |
CN101726518B (zh) | 溶液电导率的测量方法 | |
CN101587150B (zh) | 溶液电导率的测量方法 | |
CN102970007B (zh) | 用于时间电流转换的方法和装置 | |
WO2007107108A1 (fr) | Procédé permettant de mesurer la conductivité d'une solution et dispositif correspondant | |
RU2013143163A (ru) | Определение электрической емкости при электрохимическом анализе с улучшенным откликом | |
Ayari et al. | Active power measurement comparison between analog and digital methods | |
CN103969614A (zh) | 一种数字万用表的校准方法 | |
Park et al. | Novel instrumentation in electrochemical impedance spectroscopy and a full description of an electrochemical system | |
CN102043084B (zh) | 一种检测避雷器阻性电流的检测方法 | |
JP2583833Y2 (ja) | パルス測定装置 | |
CN107643227B (zh) | 基于石英晶体微天平的液体密度和粘度的测量装置及方法 | |
CN107228886B (zh) | 水的电导率电阻率检测的测量装置和方法 | |
CN109581062A (zh) | 用于示波器校准仪探头的高精度阻抗测量系统 | |
Mejía-Aguilar et al. | Electrical impedance measurement using voltage/current pulse excitation | |
CN108037358B (zh) | 单片机频率测试系统及方法 | |
RU88460U1 (ru) | Ультразвуковой расходомер (варианты) | |
RU2695030C1 (ru) | Устройство двухзондового измерения фазовых сдвигов распределённой RC-структуры | |
RU2338212C1 (ru) | Способ определения угла сдвига фаз между двумя сигналами, представленными цифровыми отсчетами | |
RU2204839C2 (ru) | Измеритель параметров электрохимических объектов | |
RU2267791C2 (ru) | Измеритель амплитуды гармонических процессов (варианты) | |
RU2260190C1 (ru) | Релаксационный измеритель параметров cg-двухполюсников | |
WO2007010320A1 (en) | Method and device for liquid conductivity measurement | |
RU2490660C1 (ru) | Измеритель фазовых погрешностей масштабного преобразователя |
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 |