CN105956242B - 基于生物体触电阻抗参数计算的多端口阻抗模型构建方法 - Google Patents

基于生物体触电阻抗参数计算的多端口阻抗模型构建方法 Download PDF

Info

Publication number
CN105956242B
CN105956242B CN201610261319.9A CN201610261319A CN105956242B CN 105956242 B CN105956242 B CN 105956242B CN 201610261319 A CN201610261319 A CN 201610261319A CN 105956242 B CN105956242 B CN 105956242B
Authority
CN
China
Prior art keywords
electric shock
impedance
organism
skin
signal
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.)
Active
Application number
CN201610261319.9A
Other languages
English (en)
Other versions
CN105956242A (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.)
China Agricultural University
Original Assignee
China Agricultural 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 China Agricultural University filed Critical China Agricultural University
Priority to CN201610261319.9A priority Critical patent/CN105956242B/zh
Publication of CN105956242A publication Critical patent/CN105956242A/zh
Application granted granted Critical
Publication of CN105956242B publication Critical patent/CN105956242B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/30Circuit design
    • G06F30/36Circuit design at the analogue level
    • G06F30/367Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明提供一种基于生物体触电阻抗参数计算的多端口阻抗模型构建方法,包括:采集活体试验对象在完整皮肤和破损皮肤状况下的多组低压配电网生物体触电信号;根据所述生物体触电信号,获取生物体触电阻抗参数;根据所述生物体触电阻抗参数,利用非线性时变电阻和电容元件搭建试验生物体触电阻抗电路模型,对生物体触电的暂态过渡过程进行描述;根据所述试验生物体触电阻抗电路模型,基于外推法获取生物体触电多端口阻抗模型。该方法能计算生物体触电阻抗参数并建立响应精度较高的生物体触电多端口阻抗模型,实现在低压配电网工频状况下对生物体触电暂态过渡过程的精确描述。

Description

基于生物体触电阻抗参数计算的多端口阻抗模型构建方法
技术领域
本发明涉及电气工程技术领域,尤其涉及一种基于生物体触电阻抗参数计算的多端口阻抗模型构建方法。
背景技术
触电事故是较为常见的恶性事故。目前,触电现象的仿真缺少对生物体触电阻抗模型的研究和描述。在低压配电网中,触电事故中考虑的直接物理因素为流过生物体的工频触电电流,而触电电流幅值直接受到生物体触电阻抗的影响。发生触电事故时,生物体可以等效为电路组合,通常利用两种阻抗来描述触电路径阻抗:内部阻抗和接触皮肤阻抗,两种阻抗的矢量和为触电路径总阻抗。当前的生物体触电电路模型通常由定阻抗元件进行构建,其适用范围具有一定的局限性,同时也无法满足生物体触电暂态过渡过程描述的精度要求。
鉴于此,需提供一种生物体触电阻抗参数计算方法,并建立响应精度较高的生物体触电多端口阻抗模型,以实现在低压配电网工频状况下,对生物体触电暂态过渡过程的精确描述。
发明内容
为解决上述的技术问题,本发明提供一种基于生物体触电阻抗参数计算的多端口阻抗模型构建方法,能够计算生物体触电阻抗参数并建立响应精度较高的生物体触电多端口阻抗模型,实现在低压配电网工频状况下,对生物体触电暂态过渡过程的精确描述。
第一方面,本发明提供一种基于生物体触电阻抗参数计算的多端口阻抗模型构建方法,包括:
采集活体试验对象在完整皮肤和破损皮肤状况下的多组低压配电网生物体触电信号;
根据所述生物体触电信号,获取生物体触电阻抗参数;
根据所述生物体触电阻抗参数,利用非线性时变电阻和电容元件搭建试验生物体触电阻抗电路模型,对生物体触电的暂态过渡过程进行描述;
根据所述试验生物体触电阻抗电路模型,基于外推法获取生物体触电多端口阻抗模型。
可选地,所述采集活体试验对象在完整皮肤和破损皮肤状况下的多组低压配电网生物体触电信号,包括:
通过低压配电网,利用电压激励、测量电流的方式,采集活体试验对象在完整皮肤和破损皮肤状况下的多组生物体触电信号;
其中,所述生物体触电信号,包括:触电电压信号序列U(m)和触电电流信号序列I(m),其中,m=0,1,…,M-1;采样参数包括:采样频率fS,采样周期TS,采样点数M,总采样时间TP
可选地,所述根据所述生物体触电信号,获取生物体触电阻抗参数,包括:
采用短时傅里叶变换STFT对所述生物体触电信号进行分析,得到所述生物体触电信号的时频特性;
从所述时频特性中提取频率为50Hz的生物体触电工频信号的幅值和相角,所述生物体触电工频信号包括:工频电压信号分量U1(m)和工频电流信号分量I1(m);
将提取的生物体触电工频信号进行分段线性化,分段区间为[m,m+1],在分段区间内利用欧姆定律进行计算,获得工频状况下生物体的触电总阻抗序列Z(m)=R(m)+j X(m),m=0,1,…,M-2,R(m)是电阻分量,X(m)是电抗分量;
根据所述触电总阻抗序列,获取生物体触电时的皮肤电阻RS(m)、皮肤电容CS(m)和内电阻R0(m);
令t=0为触电起始时刻,t=tS为触电信号的采样结束时刻,通过曲线拟合方法获取对离散数学序列RS(m)、CS(m)和R0(m)的函数表达式。
可选地,所述根据所述触电总阻抗序列,获取生物体触电时的皮肤电阻RS(m)、皮肤电容CS(m)和内电阻R0(m),包括:
根据所述触电总阻抗序列Z(m),获取皮肤阻抗ZS(m)与内阻抗Z0(m),m=0,1,…,M-2;
以弗莱贝尔格电路作为参数计算的拓扑电路,将内阻抗视为恒定电阻,取内阻抗Z0(m)的模的期望值作为内电阻,通过第一公式,求取生物体触电时的皮肤电阻RS(m)、皮肤电容CS(m)和内电阻R0(m);
其中,所述第一公式为:
其中,m=0,1,…,M-2,XC(m)=1/jωCs(m),ω为电流角频率,XC(m)为皮肤容抗。
可选地,内电阻R0(m)为常函数,函数表达式为:
R0(t)=R0
皮肤电阻RS(m)和皮肤电容CS(m)的函数表达式为:
其中,0≤t≤tS,R1、R2、C1、C2为曲线拟合解得的常数,τ1、τ2为拟合所得的衰减时间常数。
可选地,所述根据所述生物体触电阻抗参数,利用非线性时变电阻和电容元件搭建试验生物体触电阻抗电路模型,包括:
根据RS(m)、CS(m)和R0(m)的函数表达式,利用非线性时变电阻和电容元件搭建试验生物体触电阻抗电路模型。
可选地,所述根据所述试验生物体触电阻抗电路模型,基于外推法获取生物体触电多端口阻抗模型,包括:
根据生物体触电原理,求取生物体触电阻抗幅值与体重的映射关系曲线;
根据所述映射关系曲线,对所述试验生物体触电阻抗电路模型补充增益环节K,外推得到生物体触电多端口阻抗模型。
可选地,在所述采集活体试验对象在完整皮肤和破损皮肤状况下的多组低压配电网生物体触电信号之后,在所述根据所述生物体触电信号,获取生物体触电阻抗参数之前,还包括:
基于改进的FIR数字滤波器对所述生物体触电信号进行预处理;
相应地,所述根据所述生物体触电信号,获取生物体触电阻抗参数,具体为:
根据预处理后的生物体触电信号,获取生物体触电阻抗参数。
可选地,所述基于改进的FIR数字滤波器对所述生物体触电信号进行预处理,包括:
采用加窗FIR低通数字滤波器对所述生物体触电信号U(m)和I(m)进行滤波;
采用曲线匹配拟合方法消除滤波后的生物体触电信号产生的相位延迟,得到预处理后的电压信号序列U'(m)和电流信号序列I'(m)。
由上述技术方案可知,本发明的基于生物体触电阻抗参数计算的多端口阻抗模型构建方法,能够计算生物体触电阻抗参数并建立响应精度较高的生物体触电多端口阻抗模型,实现在低压配电网工频状况下,对生物体触电暂态过渡过程的精确描述。
附图说明
图1为本发明一实施例提供的基于生物体触电阻抗参数计算的多端口阻抗模型构建方法的流程示意图;
图2为本发明另一实施例提供的基于生物体触电阻抗参数计算的多端口阻抗模型构建方法的流程示意图;
图3为本发明实施例提供的弗莱贝尔格电路的结构示意图;
图4为本发明实施例提供的试验生物体触电阻抗电路模型的结构示意图;
图5为本发明实施例提供的对所述试验生物体触电阻抗电路模型和定值阻抗电路模型进行仿真分析,得到仿真触电电流信号的绝对误差对比图;
图6为本发明实施例提供的生物体触电多端口阻抗模型的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他的实施例,都属于本发明保护的范围。
图1示出了本发明一实施例提供的基于生物体触电阻抗参数计算的多端口阻抗模型构建方法的流程示意图,如图1所示,本实施例的基于生物体触电阻抗参数计算的多端口阻抗模型构建方法,包括:
101、采集活体试验对象在完整皮肤和破损皮肤状况下的多组低压配电网生物体触电信号。
在具体应用中,所述步骤101,可包括:
通过低压配电网,利用电压激励、测量电流的方式,采集活体试验对象在完整皮肤和破损皮肤状况下的多组生物体触电信号;
其中,所述生物体触电信号,包括:触电电压信号序列U(m)和触电电流信号序列I(m),其中,m=0,1,…,M-1;采样参数包括:采样频率fS,采样周期TS,采样点数M,总采样时间TP
102、根据所述生物体触电信号,获取生物体触电阻抗参数。
在具体应用中,所述步骤102,可包括图中未示出的步骤P1-P5:
P1、采用短时傅里叶变换STFT对所述生物体触电信号进行分析,得到所述生物体触电信号的时频特性。
P2、从所述时频特性中提取频率为50Hz的生物体触电工频信号的幅值和相角,所述生物体触电工频信号包括:工频电压信号分量U1(m)和工频电流信号分量I1(m)。
P3、将提取的生物体触电工频信号进行分段线性化,分段区间为[m,m+1],在分段区间内利用欧姆定律进行计算,获得工频状况下生物体的触电总阻抗序列Z(m)=R(m)+j X(m),m=0,1,…,M-2,R(m)是电阻分量,X(m)是电抗分量。
可理解的是,当采样频率足够大时,可认为生物体触电工频信号分量序列在[m,m+1]区间内是线性的,此时将提取的生物体触电工频信号进行分段线性化,可获得工频状况下生物体的触电总阻抗序列。
应说明的是,触电路径总阻抗Z(m)由皮肤阻抗ZS(m)和内阻抗Z0(m)组成。皮肤完整时,计算得到的触电总阻抗Z1(m)为皮肤阻抗ZS(m)与内阻抗Z0(m)之和,即Z1(m)=ZS(m)+Z0(m);皮肤破损时,认为皮肤阻抗为零,即ZS(m)≈0,此时,总阻抗Z2(m)≈Z0(m),联立求解上述两个方程式Z1(m)=ZS(m)+Z0(m)=ZS(m)+Z2(m),可获得数据序列ZS(m)、Z0(m)的参数值,其中,m=0,1,…,M-2。
P4、根据所述触电总阻抗序列,获取生物体触电时的皮肤电阻RS(m)、皮肤电容CS(m)和内电阻R0(m)。
在具体应用中,所述步骤P4,可以包括:
根据所述触电总阻抗序列Z(m),获取皮肤阻抗ZS(m)与内阻抗Z0(m),m=0,1,…,M-2;
以弗莱贝尔格电路(如图3所示)作为参数计算的拓扑电路,将内阻抗视为恒定电阻,取内阻抗Z0(m)的模的期望值作为内电阻,通过第一公式,求取生物体触电时的皮肤电阻RS(m)、皮肤电容CS(m)和内电阻R0(m);
其中,所述第一公式为:
其中,m=0,1,…,M-2,XC(m)=1/jωCs(m),ω为电流角频率,XC(m)为皮肤容抗。
P5、令t=0为触电起始时刻,t=tS为触电信号的采样结束时刻,通过曲线拟合方法获取对离散数学序列RS(m)、CS(m)和R0(m)的函数表达式。
其中,内电阻R0(m)为常函数,函数表达式为:
R0(t)=R0
皮肤电阻RS(m)和皮肤电容CS(m)的函数表达式为:
其中,0≤t≤tS,R1、R2、C1、C2为曲线拟合解得的常数,τ1、τ2为拟合所得的衰减时间常数。
103、根据所述生物体触电阻抗参数,利用非线性时变电阻和电容元件搭建试验生物体触电阻抗电路模型,对生物体触电的暂态过渡过程进行描述。
在具体应用中,所述步骤103中的“根据所述生物体触电阻抗参数,利用非线性时变电阻和电容元件搭建试验生物体触电阻抗电路模型”可具体为:
根据RS(m)、CS(m)和R0(m)的函数表达式,利用非线性时变电阻和电容元件搭建试验生物体触电阻抗电路模型(如图4所示)。
104、根据所述试验生物体触电阻抗电路模型,基于外推法获取生物体触电多端口阻抗模型。
在具体应用中,所述步骤104,可包括图中未示出的步骤A1和A2:
A1、根据生物体触电原理,求取生物体触电阻抗幅值与体重的映射关系曲线;
A2、根据所述映射关系曲线,对所述试验生物体触电阻抗电路模型补充增益环节K,外推得到生物体触电多端口阻抗模型,如图6所示。
应说明的是,根据生物体触电原理,在生物体触电过程中,生物体的接触触电电流与生物体的体重近似成正比。
本实施例的基于生物体触电阻抗参数计算的多端口阻抗模型构建方法,能够计算生物体触电阻抗参数并建立响应精度较高的生物体触电多端口阻抗模型,实现在低压配电网工频状况下,对生物体触电暂态过渡过程的精确描述,为工频条件下生物体触电暂态过程的数学建模和仿真分析提供借鉴,同时,对实际工况中生物体触电短路电流的计算和检测具有一定的参考价值。
在具体应用中,在所述步骤101与102之间,本实施例所述方法还可以包括:
基于改进的FIR数字滤波器对所述生物体触电信号进行预处理;
相应地,所述步骤102则具体为:
根据预处理后的生物体触电信号,获取生物体触电阻抗参数。
在具体应用中,所述基于改进的FIR数字滤波器对所述生物体触电信号进行预处理,可包括:
采用加窗FIR低通数字滤波器对所述生物体触电信号U(m)和I(m)进行滤波;
采用曲线匹配拟合方法消除滤波后的生物体触电信号产生的相位延迟(延迟数据点数为N),得到预处理后的电压信号序列U'(m)和电流信号序列I'(m)。
可理解的是,根据预处理后的生物体触电信号获取的生物体触电阻抗参数更加精确,相应本实施例最后获得的生物体触电多端口阻抗模型也更加精确。
图2示出了本发明另一实施例提供的基于生物体触电阻抗参数计算的多端口阻抗模型构建方法的流程示意图,如图2所示,本实施例的基于生物体触电阻抗参数计算的多端口阻抗模型构建方法,包括:
S1、采集活体试验对象在完整皮肤和破损皮肤状况下的多组低压配电网生物体触电信号。
所述信号通过生物体触电物理试验平台进行采集。以220V/50Hz低压配电网为电源,采用电压激励,测量电流的方式进行触电物理试验,试验电路中设有限流电阻。设定的采样参数包括:采样频率fS=10kHz,采样周期TS=1/fS=0.1ms,采样点数M=3000,总采样时间TP=300ms。考虑最危险的触电方式,设定触电路径为左前肢-右后肢。
利用高精度录波器采集活体试验对象在完整皮肤和破损皮肤状况下的多组触电电压信号序列U(m)和触电电流信号序列I(m),其中,m=0,1,…,M-1。
S2、基于改进的FIR数字滤波器对所述生物体触电信号进行预处理。
在具体应用中,所述步骤S2,可包括:
采用加窗FIR低通数字滤波器对所述生物体触电信号U(m)和I(m)进行滤波;采用曲线匹配拟合方法消除滤波后的生物体触电信号产生的相位延迟(延迟数据点数为N),得到预处理后的电压信号序列U'(m)和电流信号序列I'(m)。
在步骤S2中,可选取窗口函数为Hamming窗,窗函数长度为65,截止频率为0.04π,其FIR数字滤波器的差分方程为:
式中,bn为常系数,N为滤波器阶数,y为输出信号,x为输入信号。通过大量的触电信号验证,这样既能有效地保留触电信号中的低频分量,又能有效地滤除噪声和高频谐波。
考虑到所述FIR数字滤波器滤波后,会使信号产生相位延迟,延迟数据点数N'=N/2,N为滤波器阶数。基于最小二乘法在信号序列中选取拟合程度最佳、数据长度为N'的子序列,对原序列进行延拓。
式中,σ越小,拟合程度越好,y为原序列,y'为子序列。在滤波后通过曲线匹配拟合方法消除相位延迟。经过改进的FIR数字滤波预处理,得到电压信号序列U'(m)和电流信号序列I'(m)。
S3、根据预处理后的生物体触电信号,获取生物体触电阻抗参数。
在具体应用中,所述步骤S3,可包括图中未示出的步骤Q1-Q5:
Q1、采用短时傅里叶变换STFT对预处理后的生物体触电信号进行分析,得到预处理后的生物体触电信号的时频特性:
式中,WM=e-j2π/M,y为信号序列,γ为窗口序列。
Q2、从所述时频特性中提取频率为50Hz的生物体触电工频信号的幅值和相角,所述生物体触电工频信号包括:工频电压信号分量U1(m)和工频电流信号分量I1(m)。
在步骤Q2中,可选用窗口长度为200个采样点的矩形窗进行计算,得到各整次频率分量随采样时间变化的幅值和相角参数;然后提取频率为50Hz的生物体触电工频信号的幅值和相角,所述生物体触电工频信号包括:工频电压信号分量U1(m)和工频电流信号分量I1(m)。
Q3、将提取的生物体触电工频信号进行分段线性化,分段区间为[m,m+1],在分段区间内利用欧姆定律进行计算,获得工频状况下生物体的触电总阻抗序列Z(m)=R(m)+j X(m),m=0,1,…,M-2,R(m)是电阻分量,X(m)是电抗分量。
可理解的是,当采样频率足够大时,可认为生物体触电工频信号分量序列在[m,m+1]区间内是线性的,此时将提取的生物体触电工频信号进行分段线性化,可获得工频状况下生物体的触电总阻抗序列。
应说明的是,触电路径总阻抗Z(m)由皮肤阻抗ZS(m)和内阻抗Z0(m)组成。皮肤完整时,计算得到的触电总阻抗Z1(m)为皮肤阻抗ZS(m)与内阻抗Z0(m)之和,即Z1(m)=ZS(m)+Z0(m);皮肤破损时,认为皮肤阻抗为零,即ZS(m)≈0,此时,总阻抗Z2(m)≈Z0(m),联立求解上述两个方程式Z1(m)=ZS(m)+Z0(m)=ZS(m)+Z2(m),可获得数据序列ZS(m)、Z0(m)的参数值,其中,m=0,1,…,M-2。
Q4、根据所述触电总阻抗序列,获取生物体触电时的皮肤电阻RS(m)、皮肤电容CS(m)和内电阻R0(m)。
在具体应用中,所述步骤P4,可以包括:
根据所述触电总阻抗序列Z(m),获取皮肤阻抗ZS(m)与内阻抗Z0(m),m=0,1,…,M-2;
以弗莱贝尔格电路(如图3所示)作为参数计算的拓扑电路,将内阻抗视为恒定电阻,取内阻抗Z0(m)的模的期望值作为内电阻,通过第一公式,求取生物体触电时的皮肤电阻RS(m)、皮肤电容CS(m)和内电阻R0(m);
其中,所述第一公式为:
其中,m=0,1,…,M-2,XC(m)=1/jωCs(m),ω为电流角频率,XC(m)为皮肤容抗。
Q5、令t=0为触电起始时刻,t=tS为触电信号的采样结束时刻,通过曲线拟合方法获取对离散数学序列RS(m)、CS(m)和R0(m)的函数表达式。
其中,内电阻R0(m)为常函数,函数表达式为:
R0(t)=R0
皮肤电阻RS(m)和皮肤电容CS(m)的函数表达式为:
其中,0≤t≤tS,R1、R2、C1、C2为曲线拟合解得的常数,τ1、τ2为拟合所得的衰减时间常数。
S4、根据所述生物体触电阻抗参数,利用非线性时变电阻和电容元件搭建试验生物体触电阻抗电路模型,对生物体触电的暂态过渡过程进行描述。
在具体应用中,所述步骤S4中的“根据所述生物体触电阻抗参数,利用非线性时变电阻和电容元件搭建试验生物体触电阻抗电路模型”可具体为:
根据RS(m)、CS(m)和R0(m)的函数表达式,利用非线性时变电阻和电容元件搭建试验生物体触电阻抗电路模型(如图4所示)。
可理解的是,通过搭建低压配电网仿真电路对所述试验生物体触电阻抗电路模型进行仿真分析,得到仿真触电电流信号波形,仿真信号的绝对误差如图5中实线所示,与试验动物触电信号相比,触电电流仿真信号的相似系数为0.9969,均方根误差为1.2127。对定值阻抗电路模型进行仿真,其仿真结果的绝对误差如图5中虚线所示。对比发现,在触电发生后(第1000个采样点时刻发生触电),本实施例的改进的试验生物体触电阻抗电路模型较为显著的提高了仿真结果和触电暂态过渡过程描述的精度。
S5、根据所述试验生物体触电阻抗电路模型,基于外推法获取生物体触电多端口阻抗模型。
在具体应用中,所述步骤S5,可包括图中未示出的步骤A1和A2:
A1、根据生物体触电原理,求取生物体触电阻抗幅值与体重的映射关系曲线。
应说明的是,根据生物体触电原理,在生物体触电过程中,生物体的接触触电电流与生物体的体重近似成正比,所以根据生物体触电原理,求取的生物体触电阻抗幅值与体重的映射关系曲线为Zm=km。
以所述活体试验对象作为参考对象,则映射关系表达式为Zm/Zk=m/mk,其中Zm为待测生物体的触电路径总阻抗(单位:Ω),m为待测生物体的体重(单位:kg)。Zk为活体试验对象的触电路径总阻抗(单位:Ω),mk为活体试验对象的体重(单位:kg)。
A2、根据所述映射关系曲线,对所述试验生物体触电阻抗电路模型补充增益环节K,外推得到生物体触电多端口阻抗模型,如图6所示。
其中,增益系数K=m/mk,根据所述生物体触电多端口阻抗模型,可得:
在所述生物体触电多端口阻抗模型中,在端口1输入激励电压信号U(t),在端口2输入生物体体重m,可在端口3得到生物体触电电流暂态信号I(t),实现对其触电暂态过渡过程较为精确地描述。
本实施例的基于生物体触电阻抗参数计算的多端口阻抗模型构建方法,能够计算生物体触电阻抗参数并建立响应精度较高的生物体触电多端口阻抗模型,实现在低压配电网工频状况下,对生物体触电暂态过渡过程的精确描述,为工频条件下生物体触电暂态过程的数学建模和仿真分析提供借鉴,同时,对实际工况中生物体触电短路电流的计算和检测具有一定的参考价值。
本领域普通技术人员可以理解:实现上述各方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成。前述的程序可以存储于一计算机可读取存储介质中。该程序在执行时,执行包括上述各方法实施例的步骤;而前述的存储介质包括:ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。

Claims (8)

1.一种基于生物体触电阻抗参数计算的多端口阻抗模型构建方法,其特征在于,包括:
采集活体试验对象在完整皮肤和破损皮肤状况下的多组低压配电网生物体触电信号;
根据所述生物体触电信号,获取生物体触电阻抗参数;
根据所述生物体触电阻抗参数,利用非线性时变电阻和电容元件搭建试验生物体触电阻抗电路模型,对生物体触电的暂态过渡过程进行描述;
根据所述试验生物体触电阻抗电路模型,基于外推法获取生物体触电多端口阻抗模型,包括:根据生物体触电原理,求取生物体触电阻抗幅值与体重的映射关系曲线;根据所述映射关系曲线,对所述试验生物体触电阻抗电路模型补充增益环节K,外推得到生物体触电多端口阻抗模型。
2.根据权利要求1所述的方法,其特征在于,所述采集活体试验对象在完整皮肤和破损皮肤状况下的多组低压配电网生物体触电信号,包括:
通过低压配电网,利用电压激励、测量电流的方式,采集活体试验对象在完整皮肤和破损皮肤状况下的多组生物体触电信号;
其中,所述生物体触电信号,包括:触电电压信号序列U(m)和触电电流信号序列I(m),其中,m=0,1,…,M-1;采样参数包括:采样频率fS,采样周期TS,采样点数M,总采样时间TP
3.根据权利要求2所述的方法,其特征在于,所述根据所述生物体触电信号,获取生物体触电阻抗参数,包括:
采用短时傅里叶变换STFT对所述生物体触电信号进行分析,得到所述生物体触电信号的时频特性;
从所述时频特性中提取频率为50Hz的生物体触电工频信号的幅值和相角,所述生物体触电工频信号包括:工频电压信号分量U1(m)和工频电流信号分量I1(m);
将提取的生物体触电工频信号进行分段线性化,分段区间为[m,m+1],在分段区间内利用欧姆定律进行计算,获得工频状况下生物体的触电总阻抗序列Z(m)=R(m)+j X(m),m=0,1,…,M-2,R(m)是电阻分量,X(m)是电抗分量;
根据所述触电总阻抗序列,获取生物体触电时的皮肤电阻RS(m)、皮肤电容CS(m)和内电阻R0(m);
令t=0为触电起始时刻,t=tS为触电信号的采样结束时刻,通过曲线拟合方法获取对离散数学序列RS(m)、CS(m)和R0(m)的函数表达式。
4.根据权利要求3所述的方法,其特征在于,所述根据所述触电总阻抗序列,获取生物体触电时的皮肤电阻RS(m)、皮肤电容CS(m)和内电阻R0(m),包括:
根据所述触电总阻抗序列Z(m),获取皮肤阻抗ZS(m)与内阻抗Z0(m),m=0,1,…,M-2;
以弗莱贝尔格电路作为参数计算的拓扑电路,将内阻抗视为恒定电阻,取内阻抗Z0(m)的模的期望值作为内电阻,通过第一公式,求取生物体触电时的皮肤电阻RS(m)、皮肤电容CS(m)和内电阻R0(m);
其中,所述第一公式为:
其中,m=0,1,…,M-2,XC(m)=1/jωCs(m),ω为电流角频率,XC(m)为皮肤容抗。
5.根据权利要求3所述的方法,其特征在于,内电阻R0(m)为常函数,函数表达式为:
R0(t)=R0
皮肤电阻RS(m)和皮肤电容CS(m)的函数表达式为:
其中,0≤t≤tS,R1、R2、C1、C2为曲线拟合解得的常数,τ1、τ2为拟合所得的衰减时间常数。
6.根据权利要求3所述的方法,其特征在于,所述根据所述生物体触电阻抗参数,利用非线性时变电阻和电容元件搭建试验生物体触电阻抗电路模型,包括:
根据RS(m)、CS(m)和R0(m)的函数表达式,利用非线性时变电阻和电容元件搭建试验生物体触电阻抗电路模型。
7.根据权利要求2所述的方法,其特征在于,在所述采集活体试验对象在完整皮肤和破损皮肤状况下的多组低压配电网生物体触电信号之后,在所述根据所述生物体触电信号,获取生物体触电阻抗参数之前,还包括:
基于改进的FIR数字滤波器对所述生物体触电信号进行预处理;
相应地,所述根据所述生物体触电信号,获取生物体触电阻抗参数,具体为:
根据预处理后的生物体触电信号,获取生物体触电阻抗参数。
8.根据权利要求7所述的方法,其特征在于,所述基于改进的FIR数字滤波器对所述生物体触电信号进行预处理,包括:
采用加窗FIR低通数字滤波器对所述生物体触电信号U(m)和I(m)进行滤波;
采用曲线匹配拟合方法消除滤波后的生物体触电信号产生的相位延迟,得到预处理后的电压信号序列U'(m)和电流信号序列I'(m)。
CN201610261319.9A 2016-04-25 2016-04-25 基于生物体触电阻抗参数计算的多端口阻抗模型构建方法 Active CN105956242B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610261319.9A CN105956242B (zh) 2016-04-25 2016-04-25 基于生物体触电阻抗参数计算的多端口阻抗模型构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610261319.9A CN105956242B (zh) 2016-04-25 2016-04-25 基于生物体触电阻抗参数计算的多端口阻抗模型构建方法

Publications (2)

Publication Number Publication Date
CN105956242A CN105956242A (zh) 2016-09-21
CN105956242B true CN105956242B (zh) 2019-10-08

Family

ID=56915280

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610261319.9A Active CN105956242B (zh) 2016-04-25 2016-04-25 基于生物体触电阻抗参数计算的多端口阻抗模型构建方法

Country Status (1)

Country Link
CN (1) CN105956242B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10819340B2 (en) * 2016-11-03 2020-10-27 Vicwood Prosperity Technology Limited Living body detection method and apparatus
CN108256132B (zh) * 2016-12-29 2023-08-18 中国电力科学研究院 一种基于体重的生物体触电多端口阻抗建模方法
CN107678543B (zh) * 2017-09-25 2020-06-16 南昌大学 一种基于电触觉设备的人手皮肤-电极生物阻抗模型参数估计方法
CN108152601B (zh) * 2017-12-22 2020-07-03 重庆大学 电动汽车高压线束共模噪声引起的低频辐射发射预测方法
CN112332368B (zh) * 2020-10-23 2022-12-09 国网天津市电力公司 一种基于高阶sogi的漏电流观测识别方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102684141A (zh) * 2012-04-28 2012-09-19 中国农业大学 触电电流识别方法及剩余电流保护装置
CN102868139A (zh) * 2012-09-10 2013-01-09 中国农业大学 一种触电信号暂态量识别方法及剩余电流保护装置
CN102920454A (zh) * 2012-06-26 2013-02-13 北京四海华辰科技有限公司 一种人体阻抗测量方法、装置及设备
CN104638612A (zh) * 2015-02-13 2015-05-20 西安交通大学 一种人体触电事故识别与保护方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100274147A1 (en) * 2009-04-22 2010-10-28 Abhilash Patangay Detecting ischemia with nonlinear heart rate variability measures

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102684141A (zh) * 2012-04-28 2012-09-19 中国农业大学 触电电流识别方法及剩余电流保护装置
CN102920454A (zh) * 2012-06-26 2013-02-13 北京四海华辰科技有限公司 一种人体阻抗测量方法、装置及设备
CN102868139A (zh) * 2012-09-10 2013-01-09 中国农业大学 一种触电信号暂态量识别方法及剩余电流保护装置
CN104638612A (zh) * 2015-02-13 2015-05-20 西安交通大学 一种人体触电事故识别与保护方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于参数优化的最小二乘支持向量机触电电流检测方法;韩晓慧;《农业工程学报》;20141231;全文 *
基于有限脉冲反应和径向基神经网络的触电信号识别;关海鸥;《农业工程学报》;20130430;全文 *

Also Published As

Publication number Publication date
CN105956242A (zh) 2016-09-21

Similar Documents

Publication Publication Date Title
CN105956242B (zh) 基于生物体触电阻抗参数计算的多端口阻抗模型构建方法
WO2018161486A1 (zh) 一种基于动态参数的动力电池soc估算方法及系统
CN105264709B (zh) 用于感测电化学装置的内部温度的方法和设备
CN105004998A (zh) 电荷状态估计器及制造和使用电荷状态估计器的方法
DE102011009351B4 (de) System und Verfahren zum Überwachen einer Isolierung unter Verwendung einer variablen emulierten Induktivität
US20090157338A1 (en) Measuring method and apparatus for potentiometric measuring probes
CN103412242B (zh) 一种基于快速独立分量分析和互信息的谐波源定位方法
CN106248569B (zh) 一种基于内源式eit的接地网腐蚀诊断方法
CN107192956B (zh) 一种电池短路漏液在线监测方法和装置
CN111164436A (zh) 用于确定导电器件的阻抗的方法、设备和计算机程序
US9417290B1 (en) Methods and apparatus for dynamic characterization of electrochemical systems
WO2017008118A1 (en) Fluid level determination
Das et al. Cross spectrum aided surface condition assessment of metal oxide surge arrester employing convolutional neural network
Eisenberg Impedance measurement of the electrical structure of skeletal muscle
DE102012016279A1 (de) Verfahren zur Parametrierung eines impedanzbasierten Modells eines elektrochemischen Systems
CN106443549A (zh) 一种用于校准电池内阻测试仪的模拟交流电阻装置
CN115183884A (zh) 电热协同系统的红外测温补偿方法及装置
DE102014011397B4 (de) Verfahren zum Kalibrieren einer elektrochemischen Impedanzspektroskopie-Messvorrichtung und Impedanznormal
D'Aguanno et al. A system model for investigating passive electrical properties of neurons
CN106645965A (zh) 一种测量rlc参数的方法
CN111723454A (zh) 一种二阶pngv模型参数的在线生成方法、装置和计算机设备
CN108256132B (zh) 一种基于体重的生物体触电多端口阻抗建模方法
CN108508280A (zh) 线路零序参数计算方法、装置及电子设备
KR20050047295A (ko) 고주파 대역의 피부 임피던스 응답을 표현하는 피부임피던스 모델
Durand et al. Parameter estimation algorithms for the shunt cable model

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant