CN103630750A - 一种基于电阻抗成像的凝胶电导率测量方法 - Google Patents

一种基于电阻抗成像的凝胶电导率测量方法 Download PDF

Info

Publication number
CN103630750A
CN103630750A CN201310597254.1A CN201310597254A CN103630750A CN 103630750 A CN103630750 A CN 103630750A CN 201310597254 A CN201310597254 A CN 201310597254A CN 103630750 A CN103630750 A CN 103630750A
Authority
CN
China
Prior art keywords
conductivity
gel
solution
electrical impedance
sigma
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
Application number
CN201310597254.1A
Other languages
English (en)
Other versions
CN103630750B (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.)
Sofia medical technology (Qidong) Co., Ltd
Original Assignee
Shanghai Jiaotong University
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 Shanghai Jiaotong University filed Critical Shanghai Jiaotong University
Priority to CN201310597254.1A priority Critical patent/CN103630750B/zh
Publication of CN103630750A publication Critical patent/CN103630750A/zh
Application granted granted Critical
Publication of CN103630750B publication Critical patent/CN103630750B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明涉及一种基于电阻抗成像的凝胶电导率测量方法,包括以下步骤:1)将凝胶放置于盛有导电溶液的电阻抗实验水槽中;2)用电阻抗成像数据采集系统进行电流激励和电压采集,并使用灵敏度算法进行电导率图像重建;3)根据得到的成像结果,判断凝胶的电导率是否等于溶液的电导率,如为否,执行步骤4),否则,通过液体电导率仪测量溶液的电导率,即得到凝胶的电导率;4)调整溶液电导率,并返回步骤3)。与现有技术相比,本发明具有对于任意形状的凝胶,只要其电导率分布均匀,测量电极与凝胶不接触,就可以方便地测量一定频率的电流激励下凝胶的电导率值等优点。

Description

一种基于电阻抗成像的凝胶电导率测量方法
技术领域
本发明涉及一种凝胶电导率测量方法,尤其是涉及一种基于电阻抗成像的凝胶电导率测量方法。
背景技术
目前,实验室测量凝胶的一般技术是利用阻抗分析仪测出凝胶阻值R,再用长度测量工具测量得到凝胶的形状参数,如圆柱体凝胶的底面直径d、高度h等,然后利用下面的数学公式计算得到凝胶的电阻率,再求倒数得到凝胶电导率σ。以圆柱体凝胶为例:
R = ρ · L S
ρ = R · S L = R · 0.25 · π · d 2 h
σ = 1 ρ = 4 · h π · R · d 2
然而,这种方法不仅成本很高(高精度的阻抗分析仪价格昂贵),而且无法测量形状不规则的凝胶的电导率。而且,上述技术无法判断凝胶电导率是否分布均匀,只能测得整块凝胶的平均电导率。另外,在上述测量过程中,电极与凝胶表面接触,电极会对凝胶表面施加一定压力,这不足以完全消除接触阻抗,还造成凝胶电阻抗的变化,从而引入误差,接触阻抗和凝胶电阻抗变化的大小与压力大小有关。
发明内容
为了克服直接测量技术不能测量不规则形状的凝胶,不能验证凝胶电导率的均匀性以及受到接触阻抗等因素影响的不足,本发明提供了一种基于电阻抗成像的凝胶电导率测量方法,对于任意形状的凝胶,只要其电导率分布均匀,测量电极与凝胶不接触,就可以方便地测量一定频率的电流激励下凝胶的电导率值。
本发明可以通过以下技术方案来实现:
一种基于电阻抗成像的凝胶电导率测量方法,其特征在于,包括以下步骤:
1)将凝胶放置于盛有导电溶液的电阻抗实验水槽中;
2)用电阻抗成像数据采集系统进行电流激励和电压采集,并使用灵敏度算法进行电导率图像重建;
3)根据得到的成像结果,判断凝胶的电导率是否等于溶液的电导率,如为否,执行步骤4),否则,通过液体电导率仪测量溶液的电导率,即得到凝胶的电导率;
4)调整溶液电导率,并返回步骤3)。
凝胶被放置于用于电阻抗实验的水槽中,凝胶的高度至少与电极的上边缘平齐,凝胶的位置应尽量靠近而又不接触电极片;目标物越接近电极,测量灵敏度越高,结果越精确,但不能接触电极片,以免产生接触阻抗。
所述的导电溶液为食盐水,通过加盐或蒸馏水增大或减小食盐水的电导率。
所述的用电阻抗成像数据采集系统进行电流激励和电压采集具体为:
1)通过压控电流源使信号发生器产生恒定电流输出的正弦波信号;
2)再经过模拟开关选择敏感电极阵列中的一对电极施加激励电流,依次测量其余所有相邻电极对之间的电压,然后重复步骤2),直至所有的电极对都激励过;
3)此时得到的测量电压,经过缓冲装置以减少噪声,再经过可编程增益放大器差分放大和滤波解调后,最后经过模数转换后传输给计算机;
4)计算机根据灵敏度算法进行电导率图像重建。
根据灵敏度算法进行电导率图像重建具体为:
1)建立网格,将场域分解为M个小单元;
2)根据以下的公式,使用边界测量电压,不断迭代计算出场域内的电导率分布;
e j ( n ) = v j ′ ( σ ′ ) v j ( σ ) - u j ′ ( σ ^ k ( n ) ) u j ( σ ^ k ( 0 ) )
S ‾ j , k ( n ) = S j , k ( σ ^ k ( n ) ) Σ K = 1 M S j , k ( σ ^ k ( 0 ) )
[ γ k ( n ) ] = [ S ‾ j , k ( n ) ] - 1 · [ e j ( n ) ]
σ ^ k ( n + 1 ) = σ ^ k ( n ) + γ k ( n )
式中,V′j(σ′)和Vj(σ)分别代表了电导率分布为σ′的被重建场域和电导率分布为σ的参考场域的第j次边界测量电压,
Figure BDA0000420670810000025
代表第n次迭代的计算电导率分布,
Figure BDA0000420670810000026
为假设的初始电导率分布,
Figure BDA0000420670810000027
Figure BDA0000420670810000028
分别代表了场域电导率为
Figure BDA0000420670810000029
Figure BDA00004206708100000210
时的第j次边界计算电压,
Figure BDA0000420670810000031
表示第n次迭代的归一化后的电阻抗变化量,
Figure BDA0000420670810000032
是当电导率分布为
Figure BDA0000420670810000033
时,第k个单元的灵敏度系数,而
Figure BDA0000420670810000034
表示归一化后的灵敏度系数,
Figure BDA0000420670810000035
表示第n次迭代的归一化后的电导率修正值;
定义阈值当T<0.001时,迭代收敛,此时的
Figure BDA0000420670810000037
即为所求电导率分布计算值。
电导率的分布被表示为可视化图像,不同的电导率值用不同颜色表示。
所述的根据得到的成像结果,判断凝胶的电导率是否等于溶液的电导率具体为:
当二者不相等时,可视化图像中凝胶与溶液呈现不同颜色,需要根据颜色的差别增大或减小溶液的电导率,使其接近凝胶电导率,不断重复此过程,若能达到电导率图像中溶液和凝胶混为一体,即只显示一种颜色,说明凝胶电导率分布均匀,且凝胶的电导率等于溶液的电导率。
使用液体电导率仪前进行标定,插入溶液之前保持其干净、干燥,以防改变溶液电导率。
与现有技术相比,本发明具有以下优点:
1)这种方法可以测量任意形状的凝胶的电导率,即使凝胶的形状不规则;
2)可以检验凝胶电导率分布的均匀性;
3)电极与凝胶不接触,故凝胶电阻抗不会改变,且测量不受接触阻抗的影响;
4)它简单有效,操作简便,测量过程成本低,得到的结果无需后续计算,精度较高。
具体实施方式
下面结合具体实施例对本发明进行详细说明。
实施例
在很多电阻抗实验中,需要用一定形状的凝胶来模拟人体器官进行实验,常常要求准确测量不同形状凝胶的电导率。
以底面直径5cm,高度8cm的电导率分布均匀的圆柱体凝胶为例,用很薄的手术刀将凝胶底面切平,使凝胶能够平放在直径为30cm的大水槽中,大水槽的电极面高度范围为1.6-4.1cm;向水槽中加入一定电导率的自来水使其没过电极平面,将凝胶小心平放在靠近电极片的位置,圆柱体底面圆心距离电极片5cm,通过切割凝胶或在凝胶下方加垫,调整待测横截面的高度,使其与电极平面水平重合,用夹子等干净工具而不是双手放置凝胶可以防止汗液进入溶液改变溶液的电导率;当水面平静,凝胶也静止不动,用电阻抗成像数据采集系统向水槽通入20KHz频率的激励电流,同时测量边界电压(毫伏级),用灵敏度算法进行图像重建,完成可视化成像,重建过程使用系统自带的TJU-MEIT软件完成。
根据成像结果中溶液和凝胶区域颜色的不同,可以判断二者电导率大小的不同。观察成像结果,当背景区域显示红色或黄色而凝胶区域显示蓝色,说明溶液电导率小于凝胶,向溶液中加入食盐以提高电导率,重新进行图像重建;当背景区域显示蓝色而凝胶所在区域显示红色或黄色,说明溶液电导率大于凝胶,加入蒸馏水稀释溶液降低电导率,并调整水位与之前相等,重新进行图像重建;由于电阻抗成像灵敏度较高,加入蒸馏水或食盐都应微量,缓慢调整溶液的电导率。如上所述,反复进行图像重建,直到成像结果只显示一种颜色如蓝色,此时凝胶与溶液区域混为一体,说明了凝胶的电导率分布均匀,且其电导率等于溶液电导率。用液体电导率仪测出溶液的电导率为1.095mS/cm;然后进行验证,加微量蒸馏水,调整水位,重新图像重建,发现凝胶区域为蓝色而背景区域为黄色,测得此时溶液电导率为0.990mS/cm,说明凝胶实际电导率大于0.990mS/cm;再加入微量食盐,搅拌均匀,重新图像重建,此时凝胶区域为红色而背景区域为蓝色,测得此时溶液电导率为1.199mS/cm,说明凝胶实际电导率小于1.199mS/cm。得出结论:该圆柱体凝胶电导率分布均匀,在20KHz的激励电流下,其电导率值约为1.095mS/cm,测量误差小于±0.105mS/cm。

Claims (8)

1.一种基于电阻抗成像的凝胶电导率测量方法,其特征在于,包括以下步骤:
1)将凝胶放置于盛有导电溶液的电阻抗实验水槽中;
2)用电阻抗成像数据采集系统进行电流激励和电压采集,并使用灵敏度算法进行电导率图像重建;
3)根据得到的成像结果,判断凝胶的电导率是否等于溶液的电导率,如为否,执行步骤4),否则,通过液体电导率仪测量溶液的电导率,即得到凝胶的电导率;
4)调整溶液电导率,并返回步骤3)。
2.根据权利要求1所述的基于电阻抗成像的凝胶电导率测量方法,其特征在于,凝胶被放置于用于电阻抗实验的水槽中,凝胶的高度至少与电极的上边缘平齐,凝胶的位置应尽量靠近而又不接触电极片。
3.根据权利要求1所述的基于电阻抗成像的凝胶电导率测量方法,其特征在于,所述的导电溶液为食盐水,通过加盐或蒸馏水增大或减小食盐水的电导率。
4.根据权利要求1所述的基于电阻抗成像的凝胶电导率测量方法,其特征在于,所述的用电阻抗成像数据采集系统进行电流激励和电压采集具体为:
1)通过压控电流源使信号发生器产生恒定电流输出的正弦波信号;
2)再经过模拟开关选择敏感电极阵列中的一对电极施加激励电流,依次测量其余所有相邻电极对之间的电压,然后重复步骤2),直至所有的电极对都激励过;
3)此时得到的测量电压,经过缓冲装置以减少噪声,再经过可编程增益放大器差分放大和滤波解调后,最后经过模数转换后传输给计算机;
4)计算机根据灵敏度算法进行电导率图像重建。
5.根据权利要求1所述的基于电阻抗成像的凝胶电导率测量方法,其特征在于,根据灵敏度算法进行电导率图像重建具体为:
1)建立网格,将场域分解为M个小单元;
2)根据以下的公式,使用边界测量电压,不断迭代计算出场域内的电导率分布;
e j ( n ) = v j &prime; ( &sigma; &prime; ) v j ( &sigma; ) - u j &prime; ( &sigma; ^ k ( n ) ) u j ( &sigma; ^ k ( 0 ) )
S &OverBar; j , k ( n ) = S j , k ( &sigma; ^ k ( n ) ) &Sigma; K = 1 M S j , k ( &sigma; ^ k ( 0 ) )
[ &gamma; k ( n ) ] = [ S &OverBar; j , k ( n ) ] - 1 &CenterDot; [ e j ( n ) ]
&sigma; ^ k ( n + 1 ) = &sigma; ^ k ( n ) + &gamma; k ( n )
式中,V′j(σ′)和Vj(σ)分别代表了电导率分布为σ′的被重建场域和电导率分布为σ的参考场域的第j次边界测量电压,
Figure FDA0000420670800000024
代表第n次迭代的计算电导率分布,
Figure FDA0000420670800000025
为假设的初始电导率分布,
Figure FDA0000420670800000026
分别代表了场域电导率为
Figure FDA0000420670800000028
Figure FDA0000420670800000029
时的第j次边界计算电压,
Figure FDA00004206708000000210
表示第n次迭代的归一化后的电阻抗变化量,是当电导率分布为
Figure FDA00004206708000000212
时,第k个单元的灵敏度系数,而表示归一化后的灵敏度系数,表示第n次迭代的归一化后的电导率修正值;
定义阈值
Figure FDA00004206708000000215
当T<0.001时,迭代收敛,此时的
Figure FDA00004206708000000216
即为所求电导率分布计算值。
6.根据权利要求5所述的基于电阻抗成像的凝胶电导率测量方法,其特征在于,电导率的分布被表示为可视化图像,不同的电导率值用不同颜色表示。
7.根据权利要求6所述的基于电阻抗成像的凝胶电导率测量方法,其特征在于,所述的根据得到的成像结果,判断凝胶的电导率是否等于溶液的电导率具体为:
当二者不相等时,可视化图像中凝胶与溶液呈现不同颜色,需要根据颜色的差别增大或减小溶液的电导率,使其接近凝胶电导率,不断重复此过程,若能达到电导率图像中溶液和凝胶混为一体,即只显示一种颜色,说明凝胶电导率分布均匀,且凝胶的电导率等于溶液的电导率。
8.根据权利要求1所述的基于电阻抗成像的凝胶电导率测量方法,其特征在于,使用液体电导率仪前进行标定,插入溶液之前保持其干净、干燥,以防改变溶液电导率。
CN201310597254.1A 2013-11-22 2013-11-22 一种基于电阻抗成像的凝胶电导率测量方法 Active CN103630750B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310597254.1A CN103630750B (zh) 2013-11-22 2013-11-22 一种基于电阻抗成像的凝胶电导率测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310597254.1A CN103630750B (zh) 2013-11-22 2013-11-22 一种基于电阻抗成像的凝胶电导率测量方法

Publications (2)

Publication Number Publication Date
CN103630750A true CN103630750A (zh) 2014-03-12
CN103630750B CN103630750B (zh) 2016-07-06

Family

ID=50212006

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310597254.1A Active CN103630750B (zh) 2013-11-22 2013-11-22 一种基于电阻抗成像的凝胶电导率测量方法

Country Status (1)

Country Link
CN (1) CN103630750B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108333224A (zh) * 2018-01-31 2018-07-27 中国海洋石油集团有限公司 弱凝胶体系在岩心内成胶时间的测定装置及测定方法
CN110186583A (zh) * 2019-06-13 2019-08-30 南京航空航天大学 基于电阻抗成像的陶瓷基复合材料高温部件温度测量方法
CN110720915A (zh) * 2019-09-09 2020-01-24 浙江工业大学 一种基于gan的脑部电阻抗断层成像方法
CN112051309A (zh) * 2020-08-06 2020-12-08 南方电网科学研究院有限责任公司 一种基于电阻抗成像的直流电缆附件无损检测方法及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002228695A (ja) * 2001-02-02 2002-08-14 Mitsubishi Electric Corp 抵抗測定装置
CN101865876A (zh) * 2009-04-15 2010-10-20 上海培清科技有限公司 一种全自动分子电泳凝胶成像分析仪
CN102138790A (zh) * 2011-03-28 2011-08-03 中国医学科学院生物医学工程研究所 一种三维电阻抗断层成像系统性能测试方法与装置
CN102499678A (zh) * 2011-09-23 2012-06-20 中国人民解放军第四军医大学 一种便携式电阻抗成像系统的电阻抗测量装置及测量方法
CN103196981A (zh) * 2013-03-28 2013-07-10 上海交通大学 基于电导率检测的等电聚焦电泳测试结果预评估的方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002228695A (ja) * 2001-02-02 2002-08-14 Mitsubishi Electric Corp 抵抗測定装置
CN101865876A (zh) * 2009-04-15 2010-10-20 上海培清科技有限公司 一种全自动分子电泳凝胶成像分析仪
CN102138790A (zh) * 2011-03-28 2011-08-03 中国医学科学院生物医学工程研究所 一种三维电阻抗断层成像系统性能测试方法与装置
CN102499678A (zh) * 2011-09-23 2012-06-20 中国人民解放军第四军医大学 一种便携式电阻抗成像系统的电阻抗测量装置及测量方法
CN103196981A (zh) * 2013-03-28 2013-07-10 上海交通大学 基于电导率检测的等电聚焦电泳测试结果预评估的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
宫莲 等: "一种有效的电阻抗成像的图像重建算法", 《中国生物医学工程学报》, vol. 10, no. 3, 30 September 1991 (1991-09-30), pages 149 - 159 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108333224A (zh) * 2018-01-31 2018-07-27 中国海洋石油集团有限公司 弱凝胶体系在岩心内成胶时间的测定装置及测定方法
CN108333224B (zh) * 2018-01-31 2020-08-18 中国海洋石油集团有限公司 弱凝胶体系在岩心内成胶时间的测定装置及测定方法
CN110186583A (zh) * 2019-06-13 2019-08-30 南京航空航天大学 基于电阻抗成像的陶瓷基复合材料高温部件温度测量方法
WO2020248854A1 (zh) * 2019-06-13 2020-12-17 南京航空航天大学 基于电阻抗成像的陶瓷基复合材料高温部件温度测量方法
CN110720915A (zh) * 2019-09-09 2020-01-24 浙江工业大学 一种基于gan的脑部电阻抗断层成像方法
CN112051309A (zh) * 2020-08-06 2020-12-08 南方电网科学研究院有限责任公司 一种基于电阻抗成像的直流电缆附件无损检测方法及装置
CN112051309B (zh) * 2020-08-06 2022-07-22 南方电网科学研究院有限责任公司 一种基于电阻抗成像的直流电缆附件无损检测方法及装置

Also Published As

Publication number Publication date
CN103630750B (zh) 2016-07-06

Similar Documents

Publication Publication Date Title
Weigand et al. Multi-frequency electrical impedance tomography as a non-invasive tool to characterize and monitor crop root systems
Schlebusch et al. Bladder volume estimation from electrical impedance tomography
EP2343539B1 (en) Electrical network representation of a distributed system
CN102961136B (zh) 一种用于电阻抗断层成像系统的定标装置
Maundy et al. Extracting the parameters of the single-dispersion Cole bioimpedance model using a magnitude-only method
CN103630750B (zh) 一种基于电阻抗成像的凝胶电导率测量方法
Sanchez et al. Propagation of measurement errors through body composition equations for body impedance analysis
Bera et al. A multifrequency Electrical Impedance Tomography (EIT) system for biomedical imaging
JP2016510120A5 (zh)
Kim et al. Image reconstruction using voltage–current system in electrical impedance tomography
CN109758149B (zh) 一种结合时差成像和频差成像的电阻抗成像方法
JP2018524114A (ja) 流体レベルの決定
Taghizadeh et al. Bayesian inversion for electrical-impedance tomography in medical imaging using the nonlinear Poisson–Boltzmann equation
CN103065286A (zh) 一种准静态电阻抗成像中的图像重构方法
Gnecchi et al. Soil water infiltration measurements using electrical impedance tomography
CN105701824B (zh) 基于圆弧形网格剖分的电学层析图像重建方法
Garcia-Breijo et al. Development of a puncture electronic device for electrical conductivity measurements throughout meat salting
CN104067113A (zh) 低电导率的接触式电导率测试系统
Calvano et al. Fast methods for shape reconstruction in electrical resistance tomography
CN103163404B (zh) 一种基于相邻激励测量模式的电流-电压映射构造方法
Stacey Electrical impedance tomography
CN104083169A (zh) 基于均匀阻抗模型的bis系统验证方法
Georgi et al. Conductivity reconstructions using real data from a new planar electrical impedance tomography device
CN203953659U (zh) 一种用于生物阻抗频谱系统的复阻抗电路板
Bohuslávek The Measurement Method of Meat Conductivity.

Legal Events

Date Code Title Description
PB01 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
TR01 Transfer of patent right

Effective date of registration: 20200122

Address after: 226299 No.500 Linyang Road, Qidong City, Nantong City, Jiangsu Province

Patentee after: Sofia medical technology (Qidong) Co., Ltd

Address before: 200240 Dongchuan Road, Shanghai, No. 800, No.

Patentee before: SHANGHAI JIAO TONG University

TR01 Transfer of patent right