CN111887807A - 一种纳秒脉冲电场诱导热声信号处理方法及系统 - Google Patents
一种纳秒脉冲电场诱导热声信号处理方法及系统 Download PDFInfo
- Publication number
- CN111887807A CN111887807A CN202010173737.9A CN202010173737A CN111887807A CN 111887807 A CN111887807 A CN 111887807A CN 202010173737 A CN202010173737 A CN 202010173737A CN 111887807 A CN111887807 A CN 111887807A
- Authority
- CN
- China
- Prior art keywords
- target body
- target
- electric field
- conductivity
- electric
- 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
- 230000005684 electric field Effects 0.000 title claims abstract description 68
- 238000003672 processing method Methods 0.000 title claims description 20
- 239000000523 sample Substances 0.000 claims abstract description 47
- 238000003384 imaging method Methods 0.000 claims abstract description 45
- 238000009826 distribution Methods 0.000 claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 17
- 238000001514 detection method Methods 0.000 claims abstract description 6
- 238000004088 simulation Methods 0.000 claims description 23
- 238000004458 analytical method Methods 0.000 claims description 13
- 230000006698 induction Effects 0.000 claims description 10
- 238000002604 ultrasonography Methods 0.000 claims description 10
- 230000004044 response Effects 0.000 claims description 5
- 238000005316 response function Methods 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 4
- 241000486463 Eugraphe sigma Species 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000010438 heat treatment Methods 0.000 claims description 2
- 230000008569 process Effects 0.000 claims description 2
- 230000001939 inductive effect Effects 0.000 claims 1
- 230000035515 penetration Effects 0.000 abstract description 3
- 238000005259 measurement Methods 0.000 abstract description 2
- 230000002159 abnormal effect Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 8
- 230000008859 change Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 230000005284 excitation Effects 0.000 description 3
- 230000003068 static effect Effects 0.000 description 3
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000002547 anomalous effect Effects 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- 239000011575 calcium Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0093—Detecting, measuring or recording by applying one single type of energy and measuring its conversion into another type of energy
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/44—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
- A61B8/4416—Constructional features of the ultrasonic, sonic or infrasonic diagnostic device related to combined acquisition of different diagnostic modalities, e.g. combination of ultrasound and X-ray acquisitions
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Surgery (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Public Health (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种纳秒脉冲电场诱导热声信号处理方法及系统，属于热声成像检测领域，目标体和A，B两个电极的极板间填充绝缘油，采用电极施加纳秒宽的脉冲电压，通过超声探头采集声信号；利用有限元方法求解得到目标体的电标位，再求解目标体内的电场强度和电流密度；目标体内的电场强度和电流密度与目标体的电导率分布有关；因目标体的热函数与其电导率分布有关，最终获得目标体声信号与电导率的分布关系。该方法具有高对比度、高分辨率、非接触测量、结构简单和穿透深度更大的优势。
Description
技术领域
本发明属于热声成像检测领域，尤其涉及的是一种纳秒脉冲电场诱导热声信号处理方法及系统。
背景技术
针对生物组织电特性参数的成像技术具有重要意义，其有望用于生物检测领域，例如在生物组织发生异常前，能够提前检测目标体并对其进行成像。然而如何深度采集目标体电特性参数获取中间信息，并进行后续准确处理一直是电特性参数成像普遍存在的难题。
目前对生物组织电特性参数成像的技术主要有磁声成像、微波热声成像、感应式热声成像等。其中，磁声成像是将目标体置于静磁场中，用注入电流或感应电流激励目标体，电流在静磁场作用下产生洛伦兹力振动形成超声波，探测超声信号，然而该技术需要磁体设备产生静磁场，提高了成像成本。
微波热声成像采用微波信号对生物组织进行辐照产生热效应激发超声信号，以超声作为载体，反映生物组织内部微波吸收率差异，得到组织吸收微波的信息，从而进行成像，但是信号在组织中的穿透深度与频率成反比，这导致微波难以透入组织内部进行成像。
感应式热声成像通过激励线圈对导电目标体施加MHz量级的交变磁场，在导电物体内部产生感应电场，进而产生焦耳热，激发热弹性的超声信号，检测超声信号进行成像，与微波热声成像相比，采用的激励源频率更低，穿透深度有所增加，但是随着与激励线圈距离的增加，其感应出的电场就会变弱，不能深入组织内部进行成像。
发明内容
本发明为了解决现有技术需要与目标体直接接触，无法深入内部获取信息以便后续成像的问题，本发明提供一种纳秒脉冲电场诱导热声信号处理方法(ns pulse electricfield induced thermo-acoustic imaging，简称nsPEFTAI)，该方法具有高对比度、高分辨率、非接触测量、结构简单和穿透深度更大的优势。
本发明提供的技术方案如下：
一种纳秒脉冲电场诱导热声信号处理方法，目标体和A，B两个电极的极板间填充绝缘油，采用电极施加纳秒宽的脉冲电压U(t)，通过超声探头采集声信号；利用有限元方法求解得到目标体的电标位，再求解目标体内的电场强度和电流密度；因目标体内的电场强度和电流密度与目标体的电导率分布有关，且目标体的热函数与其电导率分布有关，从而最终获得目标体声信号与电导率的分布关系。
具体处理步骤优选包括：
将目标体置于目标区域内，在目标体和A，B两个电极的极板之间填充电导率为0的绝缘油，并将超声探头和目标体浸没在绝缘油当中，采用A，B两个电极对目标区域内施加一个纳秒宽的脉冲电压U(t)，通过超声探头采集声信号；
①由于目标体由非导电介质包裹，目标体中的电场是通过电解质感应产生的，因此符合安培定律：
其中，▽是哈密尔顿算子，J是电流密度，j是虚数单位，E是电场强度，D是电位移矢量，t是时间，ω是频率，ε是介电常数，σ是电导率，由于绝缘油中电导率为0，目标体中介电常数相对于电导率大小能够忽略不计，由此在绝缘油能够忽略σ*的实部，在目标体中能够忽略σ*的虚部，即σ*＝σ；
②由于目标体内中的电场是通过电解质感应产生的，感应的磁场非常小，故采用电准静态近似，纳秒脉冲电场诱导热声成像能够描述为
其中，
是电标位，U为上极板的电压，r为区域，∈表示属于关系，Σ1为高压电极板所在平面，Σ2为无限大的平面，Σ3-6为无穷远处的平面，n是边界的法向单位矢量；
③通过上式能够求得目标区域电标位
进而目标区域内的电场强度能够描述为：
目标体内的电流密度J为：
J＝σE (4)
由于有电流经过，目标体会吸收焦耳热，热函数H为
H(r,t)＝σE2(r,t) (5)
④通过上式能够构建出热函数图像H(r,t)，目标体吸收焦耳热膨胀会产生超声信号，声压满足的波动方程为
其中，p为目标体激发的声源图像，c0为绝缘油中的声速，β为目标体的体积热膨胀系数，Cp为目标体的比热容；由公式(6)得目标体发出的超声信号与其电导率的分布有关，通过声信号反演能够得到目标体的电导率。
优选的，在纳秒脉冲热声成像中，利用有限元方法求解公式(2)得到目标体的电标位，再利用公式(3)和(4)求解目标体内的电场强度和电流密度；目标体内的电场强度和电流密度与目标体的电导率分布有关，电导率越高的区域，其电场强度和电流密度就越高，目标体内电流密度的值与电导率相对应；由公式(5、6)可得目标体的热函数与其电导率分布有关，内部电导率越高的区域，其热函数越大，激发的声源也越大，故声源的强度与目标体中电导率的分布相对应。
优选的，在目标体中产生感应电场需要的脉冲电压幅值为10-30KV，脉宽在10-500ns之间。
优选的，利用纳秒脉冲热声成像方法分别建立简单模型和复杂模型两仿真模型进行仿真模拟分析，为纳秒脉冲热声成像的实际检测提供理论指导。
更优选的，在目标区域内感应产生一个电场，目标体内部受热膨胀产生声压p(r,t)，此热声信号被放置在目标体周围的超声探头接收，由感应式热声成像中声场正问题的分析可知，超声探头接收到的热声信号(r,t)为超声换能器位置的声压p(r,t)与超声探头的脉冲响应函数h(t)的卷积，即
式(7)表明，超声探头接收到的热声信号ω(r,t)既能反映原始声场的波形，也包含了超声探头的响应特性。
本发明还提供一种纳秒脉冲电场诱导热声成像系统，具体采用上述热声信号处理方法进行处理。
本发明综合技术方案及综合效果包括：
纳秒脉冲电场诱导热声成像继承了传统热声成像无创、高对比度的优点。
1、与磁声成像相比，本发明方法不需要静磁场，省去了磁体设备，进一步降低了设备成本和复杂度。更重要的问题是在重建过程中，避开了声源奇异性问题，更易实现快速成像。
2、与微波热声成像相比，本发明的nsPEFTAI允许更低的功率进行高效的成像，并且具有便携式成像的潜力。原因在于本领域目标体多是非磁性的，并且对磁场的响应接近真空，同时，激励源的频率降低，使得磁场穿透组织更深，不受含钙区域等高阻屏蔽效应，避免了其他辐射方法那样大的反射(反射系数>70％)，也避免了高介电常量的组织内部沉积电场的影响，提高了能量转换效率。
3、与感应式热声成像相比，采用高压脉冲电流作为激励源，避免了随着与激励源距离的增加感应电场减小的情况，穿透深度更大。
附图说明
图1为本发明实施例1纳秒脉冲电场诱导热声信号处理方法中的成像系统结构示意图。
图2为本发明实施例1纳秒脉冲电场诱导热声信号处理方法中的纳秒脉冲电压示意图。
图3为本发明实施例2纳秒脉冲电场诱导热声信号处理方法的仿真分析中，目标体电导率仿真模型图，(a)简单模型，(b)复杂模型。
图4本发明实施例纳秒脉冲电场诱导热声信号处理方法的仿真分析中，目标体仿真模型电场强度示意图，(a)简单模型，(b)复杂模型。
图5本发明实施例纳秒脉冲电场诱导热声信号处理方法的仿真分析中，目标体仿真模型电流密度示意图，(a)简单模型，(b)复杂模型。
图6本发明实施例纳秒脉冲电场诱导热声信号处理方法的仿真分析中，目标体仿真模型的热函数在1μs时刻的分布图，(a)简单模型，(b)复杂模型。
图7本发明实施例纳秒脉冲电场诱导热声信号处理方法的仿真分析中，目标体仿真模型的热函数在t＝1μs时刻，y＝0.035直线上随着x变化的分布图，(a)简单模型，(b)复杂模型。
图8本发明实施例纳秒脉冲电场诱导热声信号处理方法的仿真分析中，仿真的原始声信号图，(a)简单模型，(b)复杂模型。
图9本发明实施例纳秒脉冲电场诱导热声信号处理方法的仿真分析中，仿真超声探头接收到的信号图，(a)简单模型，(b)复杂模型。
具体实施方式
以下结合附图和具体实施例，对本发明进行详细说明。
实施例1
一种纳秒脉冲电场诱导热声信号处理方法，如图1所示，目标体和A，B两个电极的极板间填充绝缘油，采用电极施加纳秒宽的脉冲电压U(t)，通过超声探头采集声信号；利用有限元方法求解得到目标体的电标位，再求解目标体内的电场强度和电流密度；因目标体内的电场强度和电流密度与目标体的电导率分布有关，且目标体的热函数与其电导率分布有关，从而最终获得目标体声信号与电导率的分布关系。
具体处理步骤包括：
将目标体置于目标区域内，在目标体和A，B两个电极的极板之间填充电导率为0的绝缘油，并将超声探头和目标体浸没在绝缘油当中，采用A，B两个电极对目标区域内施加一个纳秒宽的脉冲电压U(t)，通过超声探头采集声信号；
①由于目标体由非导电介质包裹，目标体中的电场是通过电解质感应产生的，因此符合安培定律：
其中，▽是哈密尔顿算子，J是电流密度，j是虚数单位，E是电场强度，D是电位移矢量，t是时间，ω是频率，ε是介电常数，σ是电导率，由于绝缘油中电导率为0，目标体中介电常数相对于电导率大小能够忽略不计，由此在绝缘油能够忽略σ*的实部，在目标体中能够忽略σ*的虚部，即σ*＝σ；
②由于目标体内中的电场是通过电解质感应产生的，感应的磁场非常小，故采用电准静态近似，纳秒脉冲电场诱导热声成像能够描述为
其中，
是电标位，U为上极板的电压，r为区域，∈表示属于关系，Σ1为高压电极板所在平面，Σ2为无限大的平面，Σ3-6为无穷远处的平面，n是边界的法向单位矢量；
③通过上式能够求得目标区域电标位
进而目标区域内的电场强度能够描述为：
目标体内的电流密度J为：
J＝σE (4)
由于有电流经过，目标体会吸收焦耳热，热函数H为
H(r,t)＝σE2(r,t) (5)
④通过上式能够构建出热函数图像H(r,t)，目标体吸收焦耳热膨胀会产生超声信号，声压满足的波动方程为
其中，p为目标体激发的声源图像，c0为绝缘油中的声速，β为目标体的体积热膨胀系数，Cp为目标体的比热容；由公式(6)得目标体发出的超声信号与其电导率的分布有关，通过声信号反演能够得到目标体的电导率。
在纳秒脉冲热声成像中，利用有限元方法求解公式(2)得到目标体的电标位，再利用公式(3)和(4)求解目标体内的电场强度和电流密度；目标体内的电场强度和电流密度与目标体的电导率分布有关，电导率越高的区域，其电场强度和电流密度就越高，目标体内电流密度的值与电导率相对应；由公式(5、6)可得目标体的热函数与其电导率分布有关，内部电导率越高的区域，其热函数越大，激发的声源也越大，故声源的强度与目标体中电导率的分布相对应。
在目标体中产生感应电场需要的脉冲电压幅值为10-30KV，脉宽在10-500ns之间。
实施例2
利用实施例1所述纳秒脉冲热声成像方法，分别建立简单模型和复杂模型两仿真模型进行仿真分析。
具体的，本实施例建立的仿真模型如图3所示，图3(a)和3(b)分别是简单模型和复杂模型：
图3(a)和图3(b)外围的圆形模拟正常组织，其半径为0.08m，圆心为(0.06,0.06)m，电导率为1S/m上极板施加的电压U(t)＝100g(t)V，
其中b,c为系统常数，设置为b＝5*10-7,c＝1*10-7，下极板接地；
图3(a)内部的正方形模拟异常组织，中心位置为(0.06,0.06)m，边长为0.03m，电导率为5S/m。
图3(b)内部的椭圆形和正方形模拟异常组织，椭圆的中心位置为(0.03,0.06)m,椭圆的长半轴和短半轴分别是0.04m和0.02m，其电导率为5S/m，正方形的中心为(0.1,0.06)m，边长为0.03m，其电导率为5S/m。
在求解电场强度时，利用有限元方法对式(2)求解，能够得到目标体的电标位，然后利用式(3)和式(4)求解目标体内的电场强度和电流密度，在t＝1μs时，计算得到简单模型和复杂模型电场强度如图(4)中(a)和(b)所示，此时简单模型和复杂模型中电流密度如图(5)中(a)和(b)所示。
从图(4)和图(5)中能够看出，随着目标体的内部结构变得复杂，其电场强度和电流密度分布也趋于复杂化，内部电导率高的区域其电流密度也大，并且目标体中电流密度的模值和目标体电导率的变化相对应，电导率比较大的区域中电流密度也大，电导率小的区域中电流密度也变小，电导率均匀的区域中电流密度并不是均匀的。
仿真得到的简单模型和复杂模型中热函数在1μs时刻的分布分别如图(6)中(a)和(b)所示。图(7)中(a)和(b)分别是简单模型和复杂模型的热函数在t＝1μs时刻，y＝0.035直线上随着x变化的分布图。
由图(6)和图(7)能够看出，对于目标体内部电导率复杂的模型，热函数的分布也趋于复杂，目标体中热函数与目标体电导率的变化相对应，内部电导率比较高的区域，其热函数也比较大，激发的声源也大。因此，声源的强度与目标体中电导率的分布存在着对应关系，但是电导率均匀的区域中热函数并不是均匀的。
在脉冲电场诱导热声成像中，在目标区域内感应产生一个电场，目标体内部受热膨胀产生声压p(r,t)，此热声信号被放置在目标体周围的超声探头接收，由感应式热声成像中声场正问题的分析可知，超声探头接收到的热声信号(r,t)为超声换能器位置的声压p(r,t)与超声探头的脉冲响应函数h(t)的卷积，即
式(7)表明，超声探头接收到的热声信号ω(r,t)既能反映原始声场的波形，也包含了超声探头的响应特性。
下面分别是对图(3)中的简单模型和复杂模型激发的声信号进行的模拟仿真：
将目标体和超声探头均放置于绝缘油中，其电导率为0，目标体和绝缘油声速为1404m/s。先通过纳秒脉冲电场诱导热声成像正问题求得声压，然后利用式(7)将仿真得到的声压与超声探头的脉冲响应进行卷积，能够求解得超声探头最终测得的信号。
超声换能器的中心频率为1MHz，将超声探头放置于距远点0.08m的位置，此时的声压波形如图(8)所示，其中图(a)所示为简单模型产生声压，图(b)所示为复杂模型产生的声压；超声探头接收到的声信号随着时间的变化曲线如图(9)所示，其中图(a)所示为简单模型产生的声信号，图(b)所示为复杂模型产生的声信号。
从仿真的声信号分布曲线能够看出，声信号幅值发生变化的位置与仿真模型电导率变化的位置相对应：
图8(a)和图9(a)为简单模型的声信号波形，简单模型激发的四个脉冲信号通过与超声探头卷积后得到四个波簇，这四个波簇分别与目标体电导率变化的位置相对应。第一个波簇对应于声波从模拟正常组织的外围圆形的边界(0.04,0)m传播到超声探头位置的时间，即(0.08-0.04)/1404≈28.5μs；第二个波簇对应声波从模拟异常正方形靠近超声探头边界(0.015,0)m传播到超声探头位置的时间，即(0.08-0.015)/1404≈46.3μs；第三个波簇对应声波从模拟异常正方形远离超声探头边界(-0.015,0)m传播到超声探头位置的时间，即(0.08+0.015)/1404≈67.7μs；第四个波簇对应声波从模拟正常组织圆形的第二个边界(-0.04,0)m传播到超声探头位置的时间，即(0.08+0.04)/1404≈85.5μs。
图8(b)和图9(b)为复杂模型的声信号波形，复杂模型激发的四个个脉冲信号通过与超声探头卷积后得到四个波簇，这四个波簇分别与目标体电导率变化的边界相对应。第一个波簇对应声波从模拟异常正方形靠近超声探头边界(0.045,0)m传播到超声探头位置的时间，即(0.08-0.045)/1404≈25μs。第二个波簇对应声波从模拟异常正方形的第二个边界(0.015,0)m传播到超声探头位置的时间，即(0.08-0.015)/1404≈46.3μs；第三个波簇对应声波从模拟异常椭圆形靠近超声探头边界(-0.02,0)m传播到超声探头位置的时间，即(0.08+0.02)/1404≈71.2μs；第四个波簇对应声波从模拟异常椭圆的第二个边界(0.01,0)m传播到超声探头位置的时间，即(0.08+0.04)/1404≈85.5μs。
因此，超声探头测得的声信号也能够反映模型电导率的变化情况，此仿真结果能为纳秒脉冲热声成像的实际检测提供理论指导。
虽然本发明已作了详细描述，但对本领域技术人员来说，在本发明精神和范围内的修改将是显而易见的。此外，应当理解的是，本发明记载的各方面、不同具体实施方式的各部分、和列举的各种特征可被组合或全部或部分互换。在说明书的各个具体实施方式中，那些参考另一个具体实施方式的实施方式可适当地与其它实施方式组合，这是将由本领域技术人员所能理解的。此外，本领域技术人员将会理解，前面的描述仅是示例的方式，并不旨在限制本发明。
Claims (7)
1.一种纳秒脉冲电场诱导热声信号处理方法，其特征在于，目标体和A，B两个电极的极板间填充绝缘油，采用电极施加纳秒宽的脉冲电压U(t)，通过超声探头采集声信号；利用有限元方法求解得到目标体的电标位，再求解目标体内的电场强度和电流密度；因目标体内的电场强度和电流密度与目标体的电导率分布有关，且目标体的热函数与其电导率分布有关，从而最终获得目标体声信号与电导率的分布关系。
2.根据权利要求1所述的处理方法，其特征在于，具体处理步骤包括：
将目标体置于目标区域内，在目标体和A，B两个电极的极板之间填充电导率为0的绝缘油，并将超声探头和目标体浸没在绝缘油当中，采用A，B两个电极对目标区域内施加一个纳秒宽的脉冲电压U(t)，通过超声探头采集声信号；
①由于目标体由非导电介质包裹，目标体中的电场是通过电解质感应产生的，因此符合安培定律：
其中，
是哈密尔顿算子，J是电流密度，j是虚数单位，E是电场强度，D是电位移矢量，t是时间，ω是频率，ε是介电常数，σ是电导率，由于绝缘油中电导率为0，目标体中介电常数相对于电导率大小能够忽略不计，由此在绝缘油能够忽略σ*的实部，在目标体中能够忽略σ*的虚部，即σ*＝σ；
②由于目标体内中的电场是通过电解质感应产生的，感应的磁场非常小，故采用电准静态近似，纳秒脉冲电场诱导热声成像能够描述为
其中，
是电标位，U为上极板的电压，r为区域，∈表示属于关系，Σ1为高压电极板所在平面，Σ2为无限大的平面，Σ3-6为无穷远处的平面，n是边界的法向单位矢量；
③通过上式能够求得目标区域电标位
进而目标区域内的电场强度能够描述为：
目标体内的电流密度J为：
J＝σE (4)
由于有电流经过，目标体会吸收焦耳热，热函数H为
H(r,t)＝σE2(r,t) (5)
④通过上式能够构建出热函数图像H(r,t)，目标体吸收焦耳热膨胀会产生超声信号，声压满足的波动方程为
其中，p为目标体激发的声源图像，c0为绝缘油中的声速，β为目标体的体积热膨胀系数，Cp为目标体的比热容；由公式(6)得目标体发出的超声信号与其电导率的分布有关，通过声信号反演能够得到目标体的电导率。
3.根据权利要求2所述的处理方法，其特征在于，在纳秒脉冲热声成像中，利用有限元方法求解公式(2)得到目标体的电标位，再利用公式(3)和(4)求解目标体内的电场强度和电流密度；目标体内的电场强度和电流密度与目标体的电导率分布有关，电导率越高的区域，其电场强度和电流密度就越高，目标体内电流密度的值与电导率相对应；由公式(5、6)可得目标体的热函数与其电导率分布有关，内部电导率越高的区域，其热函数越大，激发的声源也越大，故声源的强度与目标体中电导率的分布相对应。
4.根据权利要求1所述的处理方法，其特征在于，其特征在于，在目标体中产生感应电场需要的脉冲电压幅值为10-30KV，脉宽在10-500ns之间。
5.根据权利要求1-4任一项所述的处理方法，其特征在于，利用纳秒脉冲热声成像方法分别建立简单模型和复杂模型两仿真模型进行仿真模拟分析，为纳秒脉冲热声成像的实际检测提供理论指导。
6.根据权利要求5所述的处理方法，其特征在于，在目标区域内感应产生一个电场，目标体内部受热膨胀产生声压p(r,t)，此热声信号被放置在目标体周围的超声探头接收，由感应式热声成像中声场正问题的分析可知，超声探头接收到的热声信号(r,t)为超声换能器位置的声压p(r,t)与超声探头的脉冲响应函数h(t)的卷积，即
式(7)表明，超声探头接收到的热声信号ω(r,t)既能反映原始声场的波形，也包含了超声探头的响应特性。
7.一种纳秒脉冲电场诱导热声成像系统，其特征在于，采用如权利要求1-6任一项所述热声信号处理方法进行处理。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010173737.9A CN111887807B (zh) | 2020-03-13 | 2020-03-13 | 一种纳秒脉冲电场诱导热声信号处理方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010173737.9A CN111887807B (zh) | 2020-03-13 | 2020-03-13 | 一种纳秒脉冲电场诱导热声信号处理方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111887807A true CN111887807A (zh) | 2020-11-06 |
| CN111887807B CN111887807B (zh) | 2021-08-17 |
Family
ID=73169821
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010173737.9A Expired - Fee Related CN111887807B (zh) | 2020-03-13 | 2020-03-13 | 一种纳秒脉冲电场诱导热声信号处理方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111887807B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112535469A (zh) * | 2021-01-14 | 2021-03-23 | 华南师范大学 | 基于磁场调控的差分热声成像方法与装置 |
| CN112694974A (zh) * | 2020-11-26 | 2021-04-23 | 中国石油大学(华东) | 针对纳秒脉冲电场消融动态监测系统构建及监测方法 |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102223840A (zh) * | 2008-09-10 | 2011-10-19 | 安德拉有限公司 | 光声成像设备 |
| CN102269717A (zh) * | 2011-04-27 | 2011-12-07 | 华南师范大学 | 超短脉冲微波热声成像方法及其装置 |
| CN105395167A (zh) * | 2014-09-05 | 2016-03-16 | 佳能株式会社 | 被检体信息获取装置 |
| CN106580249A (zh) * | 2016-11-08 | 2017-04-26 | 中国科学院电工研究所 | 一种注入电流式热声成像方法 |
| CN106821334A (zh) * | 2017-03-28 | 2017-06-13 | 中国石油大学(华东) | 一种低温等离子体诱导热声成像方法 |
| CN106885842A (zh) * | 2017-02-20 | 2017-06-23 | 中国科学院电工研究所 | 一种注入电流式热声成像电阻率重建方法 |
| CN107049315A (zh) * | 2017-02-20 | 2017-08-18 | 中国科学院电工研究所 | 基于最优化迭代方法的注入电流式热声电阻率图像重建方法 |
| CN107064302A (zh) * | 2017-02-20 | 2017-08-18 | 中国科学院电工研究所 | 一种注入电流式热声成像电导率重建方法 |
| CN110333294A (zh) * | 2019-08-13 | 2019-10-15 | 中国石油大学(华东) | 一种注入电流式热声成像方法 |
| CN110613431A (zh) * | 2019-10-21 | 2019-12-27 | 四川华微康医疗科技有限公司 | 一种便携式微波热声、超声双模态乳腺成像装置及其方法 |
-
2020
- 2020-03-13 CN CN202010173737.9A patent/CN111887807B/zh not_active Expired - Fee Related
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102223840A (zh) * | 2008-09-10 | 2011-10-19 | 安德拉有限公司 | 光声成像设备 |
| CN102269717A (zh) * | 2011-04-27 | 2011-12-07 | 华南师范大学 | 超短脉冲微波热声成像方法及其装置 |
| CN105395167A (zh) * | 2014-09-05 | 2016-03-16 | 佳能株式会社 | 被检体信息获取装置 |
| CN106580249A (zh) * | 2016-11-08 | 2017-04-26 | 中国科学院电工研究所 | 一种注入电流式热声成像方法 |
| CN106885842A (zh) * | 2017-02-20 | 2017-06-23 | 中国科学院电工研究所 | 一种注入电流式热声成像电阻率重建方法 |
| CN107049315A (zh) * | 2017-02-20 | 2017-08-18 | 中国科学院电工研究所 | 基于最优化迭代方法的注入电流式热声电阻率图像重建方法 |
| CN107064302A (zh) * | 2017-02-20 | 2017-08-18 | 中国科学院电工研究所 | 一种注入电流式热声成像电导率重建方法 |
| CN106821334A (zh) * | 2017-03-28 | 2017-06-13 | 中国石油大学(华东) | 一种低温等离子体诱导热声成像方法 |
| CN110333294A (zh) * | 2019-08-13 | 2019-10-15 | 中国石油大学(华东) | 一种注入电流式热声成像方法 |
| CN110613431A (zh) * | 2019-10-21 | 2019-12-27 | 四川华微康医疗科技有限公司 | 一种便携式微波热声、超声双模态乳腺成像装置及其方法 |
Non-Patent Citations (2)
| Title |
|---|
| ZARAFSHANI等: "Real-time in-situ monitoring of electrotherapy process using electric pulse-induced acoustic tomography(EpAT)", 《SPIE》 * |
| 郭亮: "Non-equilibrium plasma jet inducedthermo-acoustic resistivity imaging for higher contrast and resolution", 《NATURE SCIENTIFIC REPORTS》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112694974A (zh) * | 2020-11-26 | 2021-04-23 | 中国石油大学(华东) | 针对纳秒脉冲电场消融动态监测系统构建及监测方法 |
| CN112535469A (zh) * | 2021-01-14 | 2021-03-23 | 华南师范大学 | 基于磁场调控的差分热声成像方法与装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111887807B (zh) | 2021-08-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102894974B (zh) | 一种磁声电成像系统及成像方法 | |
| CN111887807B (zh) | 一种纳秒脉冲电场诱导热声信号处理方法及系统 | |
| Wang et al. | Analysis and design of coil-based electromagnetic-induced thermoacoustic for rail internal-flaw inspection | |
| Wang et al. | Novel coil transducer induced thermoacoustic detection of rail internal defects towards intelligent processing | |
| KR20070102400A (ko) | 전계 제어장치 및 검출장치 | |
| Sun et al. | Lorentz force electrical-impedance tomography using linearly frequency-modulated ultrasound pulse | |
| CN111948291B (zh) | 一种纳秒脉冲电场诱导热声成像的电导率重建方法 | |
| CN104434101B (zh) | 一种磁热声成像方法及其成像系统 | |
| US10368750B2 (en) | Shear wave imaging method and installation for collecting information on a soft solid | |
| CN114646818A (zh) | 空间电荷非接触测试系统及方法 | |
| Xu et al. | Magneto-acousto-electrical tomography: A new imaging modality for electrical impedance | |
| CN106885842B (zh) | 一种注入电流式热声成像电阻率重建方法 | |
| Shunqi et al. | Magneto-acoustic imaging by continuous-wave excitation | |
| CN110710972A (zh) | 一种线圈检测式磁声电层析成像的信息处理方法 | |
| Desjardins et al. | Transient response of a driver coil in transient eddy current testing | |
| Hedia et al. | A full time domain methodology based on near field time reversal for equivalent source identification | |
| CN107049315B (zh) | 基于最优化迭代方法的注入电流式热声电阻率图像重建方法 | |
| Zhang et al. | Design of bulk wave EMAT using a pulsed electromagnet | |
| Yang et al. | Conductivity reconstruction for magnetically mediated thermoacoustic imaging | |
| Samosir et al. | Propagation analysis of electromagnetic wave radiated by partial discharge in single phase and three phase 150 kV Gas Insulated Switchgear | |
| Jacobs et al. | Low-power microwave induced thermoacoustic imaging: Experimental study and hybrid FEM modeling | |
| Stobbs et al. | A supersonic underwater discharge as a high-power ultrasound source | |
| Liu et al. | Analysis of short pulse impacting on microwave induced thermo-acoustic tomography | |
| Ahmed et al. | Modeling of an EMP simulator using a 3-D FDTD code | |
| Hiwaki et al. | The property of nerve excitation elicited by magnetic stimulation of peripheral nerve |
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 |
Granted publication date: 20210817 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |