JP3785064B2 - 電気導電率推定方法および装置 - Google Patents
電気導電率推定方法および装置 Download PDFInfo
- Publication number
- JP3785064B2 JP3785064B2 JP2001234095A JP2001234095A JP3785064B2 JP 3785064 B2 JP3785064 B2 JP 3785064B2 JP 2001234095 A JP2001234095 A JP 2001234095A JP 2001234095 A JP2001234095 A JP 2001234095A JP 3785064 B2 JP3785064 B2 JP 3785064B2
- Authority
- JP
- Japan
- Prior art keywords
- conductivity
- region
- measurement
- distribution
- sensitivity detector
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)
- Measuring Magnetic Variables (AREA)
- Measurement Of Resistance Or Impedance (AREA)
- Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
Description
【発明の属する技術分野】
本発明は、定量性のある非破壊検査技法及び装置、物質・材料特性評価技法及び装置に関する。
【0002】
【従来の技術】
図2に従来の技術による電気導電率推定装置の構成図を示す。
【0003】
従来の方法は、図2の構成による電気インピーダンスCT法である。電流源13によって対象物に積極的に電流を流し込み、それによって生じる電圧を電圧検出器12によって物体表面の多数の点にて測定し、データ処理手段8で有限差分法や要素法等を用いてモデル化し、感度理論に基づいて物体内部の導電率を推定するというものである。
【0004】
【発明が解決しようとする課題】
従来の方法では、外部から電流を強制的に流し込む必要が有り、高感度検出器も接触型を使用していたため、試料内部に既に電流場が存在した場合にその場を乱してしまうという問題が有った。また、試料の全体を有限差分法や要素法等でモデル化する必要が有り、未知試料の関心領域のみの測定から、その関心領域の電気導電率を推定することは困難であった。
【0005】
【課題を解決するための手段】
そこで本発明では、高感度な磁気検出器や電磁波検出器や電荷センサーなどの非接触検出器で電流密度分布又は電位分布を測定し、その一方の測定値によって記述される一階の空間偏微分方程式を解くことによって、既に電流場がある場合でもそれを乱すこと無く、関心領域内において与えられた参照導電率に対する相対的な導電率空間分布を推定可能であること、更に、正則化された代数方程式を用いることによって、測定データに含まれるエラー(ノイズ)データや参照領域(導電率の参照値の与えられる領域)が狭くて位置が悪い場合においても関心領域内の一方の測定データのみから電気導電率を推定可能であることを特徴とする。
【0006】
本発明の正則化された代数方程式は、以下に説明するように導出される。
【0007】
測定された値によって記述される電気導電率に関する一階空間編微分方程式に対して、有限差分近似や有限要素法を用いて離散化する。ここでは有限差分法を適用した場合について述べる。有限差分法を適用した場合も同様の手順による。
【0008】
関心領域の2つの独立した電流場に対する電流密度分布を(J1x(x,y),J1y(x,y))と(J2x(x,y),J2y(x,y))とすると、下記の一階偏微分方程式が成り立つ。
【0009】
【数1】
【0010】
電流場が1つである場合は、一つの偏微分方程式のみが成立する。
【0011】
尚、電位分布が測定される場合は、2つの独立した電位分布をそれぞれV1、V2とすると、電場はそれぞれE1=−∇V1、E2=−∇V2となり、下記の一階偏微分方程式が成り立つ。
【0012】
【数2】
【0013】
初期条件は、次の「発明の実施の形態」にて詳述するが、一般的には、参照導電率は
【0014】
【数3】
【0015】
という様に、ある領域内に与えられるが、二つの場が測定された場合には参照値は関心領域内の一点で与えられれば良い。
【0016】
この様な一階空間偏微分方程式に対し、離散デカルト座標系(x,y)〜(IΔx,JΔy)を使い有限差分近似を適用し、初期条件式3を代入することにより、
連立方程式
【0017】
【数4】
【0018】
を得る。但し、sは未知導電率の空間分布を表すベクトル、Dは一階偏微分の有限差分近似定数からなる行列、Jおよびjは各々電流密度分布から定まる行列およびその一階変微分値からなるベクトルである。
【0019】
これを最小二乗法を用いて解くことになるが、J、jは電流密度およびその空間微分値に低域通過型フィルタをかけたもので決まるため、JDの逆作用素はjに含まれる高周波数帯のノイズを増幅させてしまう。つまり、sは不安定な結果となってしまう。そこで、いわゆる正則化を応用して再構成の安定化を図る。具体的には、正則化パラメータα1およびα2(正値)を用いて、
【0020】
【数5】
【0021】
をsに関して最小化する。但し、DsおよびDT Dsは各々未知導電率の空間分布の勾配およびラプラシアンである。DsおよびDT Dsは正定値であるため、e(s)は必ず一つの最小値を持つことになる。
式5のe(s)の最小化により、正則化された正規行列
【0022】
【数6】
【0023】
が得られ、従って、解は
【0024】
【数7】
【0025】
である。
【0026】
【発明の実施の形態】
以下に請求項1〜3記載の発明に対応する実施の形態について説明する。
【0027】
図1に本発明の第一の実施形態による計測装置の全体構成図を示す。本装置は計測対象物4の関心領域の電流密度分布や電位分布を測定するための高感度検出器1と走査機構3と、測定対象物4と高感度検出器1の間隔を調整する間隔調整手段14と、計測対象物4に電流場を発生させる電流場発生手段5と、高感度検出器を駆動する駆動装置2と、それらをコントロールする計測制御手段6と測定データを記憶するデータ記録手段7と、測定データから導電率分布を推定するデータ処理手段8で構成される。
【0028】
本発明では、電流密度分布又は電圧分布を非接触で測定するために、磁気センサーや静電容量センサー、電磁波センサー等を用いる。高感度検出器1は駆動装置2に接続し、これによって駆動される。
【0029】
計測制御手段6とデータ記録手段7とデータ処理手段8は、一台の装置上で実現することも可能であり、或いは、複数の装置に分散させることも可能である。走査機構3と間隔調整手段14の位置を逆にすることも可能である。
【0030】
図4は、高感度検出器1の側に走査機構3と間隔調整手段14を設けた実施例である。走査機構3と間隔調整手段14の位置を逆にすることも可能である。
【0031】
図5は、高感度検出器1側に第一の走査機構10を設け、計測対象物4側に第二の走査機構11を設けた実施例である。この実施例では、間隔調整手段14を高感度検出器側1に設けているが、計測対象物4側でもかまわない。第一の走査機構10と間隔調整手段14の位置を逆にすることも可能である。
【0032】
次に図3のフローチャートに従って、導電率分布を推定する方法について説明する。
【0033】
まず、参照領域を適当に設定する。参照領域の設定は以下のように行う。
【0034】
測定時に独立な電流場が2つ以上設定できる場合には、つまり、2つ以上の独立した電流密度分布又は電位分布が測定可能である場合には、参照領域として少なくとも1つの参照点を設定する。参照点は導電率が既知である点である。
【0035】
測定時に電流場が1つしか設定できない場合には、電流密度分布の測定の場合は、電流が支配的に流れる方向に広く伸びるように参照領域を設定する。電位分布の測定の場合は、電流が支配的に流れる方向と広く交わるように参照領域を設定する。参照領域とは導電率が既知である領域である。
【0036】
次に、設定した参照領域と関心領域を含むように測定領域を設定し、電流密度分布或いは電位分布を測定する。関心領域とは導電率分布を知りたい領域である。
【0037】
測定は、2つの独立した電流密度分布を測定する場合次のように行う。
【0038】
電流場発生装置を用いて、関心領域に第一の電流場を生じせしめ、第一の電流密度分布を測定する。続いて、第一の電流場とは独立の第二の電流場を生じせしめ、第二の電流密度分布を測定する。既に電流場が存在する場合は、その1つの場のみを推定し、参照領域を設定すればよい。測定は、計測制御手段によってセンサー駆動装置とXYZステージをコントロールし、XYZステージでサンプルを走査しながら、位置情報と検出信号をデータ記録手段に入力する。
【0039】
データ処理手段8において、測定データに対してノイズ除去のためのフィルタリングを行い、空間的に平滑化し、式4のJの係数を求める。それを元に、式7から関心領域の導電率分布sを求めることが出来る。
【0040】
【発明の効果】
本発明によれば、未知試料の関心領域の電気導電率を、関心領域のみの測定で得た1つの物理量の分布から求めることができる。特に、試料内部に既に電流場が存在した場合にはその場を乱すことなく容易に関心領域の電気導電率を推定することが可能となる。
【図面の簡単な説明】
【図1】本発明の一実施形態による計測装置の構成を示すブロック図である。
【図2】従来の技術による計測装置の構成を示すブロック図である。
【図3】本発明の一実施形態による制御プログラムのフローチャートである。
【図4】本発明の一実施形態による計測装置の構成を示すブロック図である。
【図5】本発明の一実施形態による計測装置の構成を示すブロック図である。
【符号の説明】
1・・・高感度検出器
2・・・駆動装置
3・・・走査機構
4・・・計測対象物
5・・・電流場発生装置
6・・・計測制御手段
7・・・データ記録手段
8・・・データ処理手段
9・・・ハウジング
10・・・第一の走査機構
11・・・第二の走査機構
12・・・電圧検出器
13・・・電流源
14・・・間隔調整手段
15・・・試料設置台
Claims (4)
- 測定対象物の電流密度分布を測定するための非接触の高感度検出器と、
前記測定対象物と前記高感度検出器の間隔を調整する間隔調整手段と、
前記高感度検出器を駆動する駆動装置と、
前記高感度検出器により測定した電流密度の位置情報と検出信号の測定データを記憶するデータ記録手段と、
導電率分布を算出するためのデータ処理手段と、を具備し、
電流が支配的に流れる方向に広く伸びる領域であり導電率が既知である参照領域と導電率を求める関心領域とを含む測定領域の電流密度分布を一つのみ測定することで導電率分布を求めることを特徴とする電気導電率推定装置。 - 測定対象物の電位分布を測定するための非接触の高感度検出器と、
前記測定対象物と前記高感度検出器の間隔を調整する間隔調整手段と、
前記高感度検出器を駆動する駆動装置と、
前記高感度検出器により測定した電位の位置情報と検出信号の測定データを記憶するデータ記録手段と、
導電率分布を算出するためのデータ処理手段と、を具備し、
電流が支配的に流れる方向と広く交わる領域であり導電率が既知である参照領域と導電率を求める関心領域とを含む測定領域の電位分布を一つのみ測定することで導電率分布を求めることを特徴とする電気導電率推定装置。 - 電流が支配的に流れる方向に広く伸びる領域であり導電率が既知である参照領域を設定する工程と、
前記参照領域と導電率を求める関心領域とを含む測定領域を設定する工程と
前記測定領域の電流密度分布を非接触の高感度検出器で一つのみ測定する工程と、
前記測定した電流密度の位置情報と検出信号の測定データをデータ記録手段に記憶させる工程と、
前記データ記憶手段に記録された測定データからデータ処理手段により導電率分布を求める工程と、からなる電気導電率推定方法。 - 電流が支配的に流れる方向と広く交わる領域であり導電率が既知である参照領域を設定する工程と、
前記参照領域と導電率を求める関心領域とを含む測定領域を設定する工程と
前記測定領域の電位分布を非接触の高感度検出器で一つのみ測定する工程と、
前記測定した電位の位置情報と検出信号の測定データをデータ記録手段に記憶させる工程と、
前記データ記憶手段に記録された測定データからデータ処理手段により導電率分布を求める工程と、からなる電気導電率推定方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2001234095A JP3785064B2 (ja) | 2000-08-07 | 2001-08-01 | 電気導電率推定方法および装置 |
Applications Claiming Priority (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2000-239036 | 2000-08-07 | ||
| JP2000239036 | 2000-08-07 | ||
| JP2001234095A JP3785064B2 (ja) | 2000-08-07 | 2001-08-01 | 電気導電率推定方法および装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2002156362A JP2002156362A (ja) | 2002-05-31 |
| JP3785064B2 true JP3785064B2 (ja) | 2006-06-14 |
Family
ID=26597510
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2001234095A Expired - Fee Related JP3785064B2 (ja) | 2000-08-07 | 2001-08-01 | 電気導電率推定方法および装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3785064B2 (ja) |
-
2001
- 2001-08-01 JP JP2001234095A patent/JP3785064B2/ja not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2002156362A (ja) | 2002-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111615355B (zh) | 通道完整性检测和电生理信号的重构 | |
| CN101517436A (zh) | 磁感应断层成像的方法和设备 | |
| CN103926276B (zh) | 一种在线油液磨粒监测装置及检测方法 | |
| Kourunen et al. | Suitability of a PXI platform for an electrical impedance tomography system | |
| Marashdeh et al. | Electrical capacitance tomography | |
| CN110990757B (zh) | 利用无相位数据解决高度非线性电磁逆散射问题的方法 | |
| JP2002532723A (ja) | ある容量内の電磁場特性を決定するための方法および装置 | |
| Soleimani et al. | Dynamic imaging in electrical capacitance tomography and electromagnetic induction tomography using a Kalman filter | |
| CN102187252A (zh) | 用于磁感应断层成像的方法和系统 | |
| US20060239328A1 (en) | Thermal properties measurement apparatus | |
| Kluth et al. | Model uncertainty in magnetic particle imaging: Nonlinear problem formulation and model-based sparse reconstruction | |
| Cui et al. | Planar electrical capacitance tomography dynamic imaging for non-destructive test | |
| JP4286664B2 (ja) | 誘電率又は導電率推定装置 | |
| JP4460808B2 (ja) | 電流密度ベクトル推定装置および電気導電率推定装置 | |
| US11543443B2 (en) | Impedance measuring apparatus | |
| Henn et al. | Improving model-based MPI image reconstructions: Baseline recovery, receive coil sensitivity, relaxation and uncertainty estimation | |
| JP3785064B2 (ja) | 電気導電率推定方法および装置 | |
| CN114373024B (zh) | 一种基于序贯蒙特卡洛原理的emt图像重建方法 | |
| JP2003222664A (ja) | 電流密度ベクトル推定装置および電気導電率推定装置 | |
| US10775346B2 (en) | Virtual channels for eddy current array probes | |
| WO2017008172A1 (en) | A method and device for determining the wear of a carbon ceramic brake disc in a vehicle by impedance measurements | |
| JP5441292B2 (ja) | 誘電率又は導電率推定装置 | |
| JP7769405B2 (ja) | 検査装置及び検査方法 | |
| CN117664383A (zh) | 基于同步检波的温度测量方法和装置 | |
| Suppan et al. | A Kalman filter approach for the application of electrical capacitance tomography in dynamic processes using a state reduction |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| RD01 | Notification of change of attorney |
Free format text: JAPANESE INTERMEDIATE CODE: A7421 Effective date: 20040427 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20040428 |
|
| A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A712 Effective date: 20040611 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20050802 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20050930 |
|
| A02 | Decision of refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A02 Effective date: 20051206 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20060206 |
|
| A911 | Transfer of reconsideration by examiner before appeal (zenchi) |
Free format text: JAPANESE INTERMEDIATE CODE: A911 Effective date: 20060216 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20060314 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20060316 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 Ref document number: 3785064 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100324 Year of fee payment: 4 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110324 Year of fee payment: 5 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110324 Year of fee payment: 5 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120324 Year of fee payment: 6 |
|
| R360 | Written notification for declining of transfer of rights |
Free format text: JAPANESE INTERMEDIATE CODE: R360 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120324 Year of fee payment: 6 |
|
| R370 | Written measure of declining of transfer procedure |
Free format text: JAPANESE INTERMEDIATE CODE: R370 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120324 Year of fee payment: 6 |
|
| S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120324 Year of fee payment: 6 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120324 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130324 Year of fee payment: 7 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20140324 Year of fee payment: 8 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |