CN113447701B - 雷电流峰值标定方法及装置 - Google Patents
雷电流峰值标定方法及装置 Download PDFInfo
- Publication number
- CN113447701B CN113447701B CN202110542392.4A CN202110542392A CN113447701B CN 113447701 B CN113447701 B CN 113447701B CN 202110542392 A CN202110542392 A CN 202110542392A CN 113447701 B CN113447701 B CN 113447701B
- Authority
- CN
- China
- Prior art keywords
- lightning current
- lightning
- current recording
- data
- magnetic stripe
- 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
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/04—Measuring peak values or amplitude or envelope of ac or of pulses
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Digital Magnetic Recording (AREA)
Abstract
本公开提供一种雷电流峰值标定方法及装置,涉及雷电测量技术领域,能够解决现有斜率法精度不高以及无法全量程覆盖的问题。具体技术方案为:通过计算雷电流记录卡磁条数据的包络数据,结合最小二乘拟合算法,利用雷电实验室冲击电流发生器产生的标准雷电冲击电流进行模拟试验来完成雷电流峰值标定的高精度标定方法。本发明用于雷电流峰值测量。
Description
技术领域
本公开涉及雷电测量技术领域,尤其涉及雷电流峰值标定方法及装置。
背景技术
雷电是伴有闪电和雷鸣的一种雄伟壮观而又有点令人生畏的放电现象,其电压可高达数1000MV~10MV,其电流可达数10kA~300kA,具有极大破坏力,因此常常造成人畜伤亡、建筑物被毁、引发***、火灾及电力***停电、电气设备遭受破坏等严重事故。雷电流幅值的测量是进行雷电研究的基础工作,对研究雷电特性、分析雷害事故、探讨防雷对策均具有十分重要的意义。
现有雷电流记录卡读卡器,检测雷电流记录卡的雷电流峰值的算法原理是斜率法,即通过采集磁条数据形成的磁场图形,取固定位置的两个点计算图形斜率,来确定雷电流峰值的大小。这种算法在雷电流低于80kA时,检测值相对准确。由于在自然雷电全范围内,雷电流峰值与磁条信息不构成线性关系,因而不适用于大电流的检测,因此现有雷电流记录卡读卡器精度低、量程小、体积大、智能化较低,无法覆盖0kA~300kA的宽量程。
发明内容
本公开实施例提供一种雷电流峰值标定方法及装置,能够解决现有斜率法精度不高以及无法全量程覆盖的问题。所述技术方案如下:
根据本公开实施例的第一方面,提供一种雷电流峰值标定方法,该方法包括:
获取P个雷电流记录卡,P个雷电流记录卡为未被雷电冲击处于满磁状态的雷电流记录卡;
将P个雷电流记录卡分别进行不同幅值的雷电冲击模拟实验,获取每个雷电流记录卡的雷电流峰值以及对应的磁条数据,并计算每个雷电流记录卡的磁条数据的包络数据;
采用M个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据,拟合得到雷电流记录卡的雷电流峰值与包络数据的初始拟合函数;
采用N个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对初始拟合函数中的参数进行标定,得到目标拟合函数,其中,P=M+N。
在一个实施例中,将P个雷电流记录卡分别进行不同幅值的雷电冲击模拟实验包括:
获取雷电冲击电流的幅值范围,并将雷电冲击电流的幅值范围按照预设间隔进行划分,得到M个子间隔;
将P个雷电流记录卡中的M个雷电流记录卡按照M个子间隔进行不同幅值的雷电冲击模拟实验,将M个子间隔的间隔端点最大值确定为每个雷电流记录卡的雷电流峰值;
将P个雷电流记录卡中的N个雷电流记录卡分为两组,将第一组中的每个雷电流记录卡分别进行8/20μs的雷电冲击模拟试验,将第二组中的每个雷电流记录卡分别进行10/350μs的雷电冲击模拟试验,获取N个雷电流记录卡中每个雷电流记录卡的雷电流峰值。
在一个实施例中,获取每个雷电流记录卡的磁条数据包括:通过雷电流记录卡读卡器获取每个雷电流记录卡的磁条数据。
在一个实施例中,包络数据包括包络面积或包络长度。
在一个实施例中,计算每个雷电流记录卡的磁条数据的包络数据包括:
在包络数据为包络面积时,采用第一公式计算每个雷电流记录卡的磁条数据的包络数据,磁条数据为曲边图形,包络数据为曲边图形的面积;
第一公式为:其中,S表示包络数据,n表示对曲边图形划分的区间个数,Δx表示划分的区间间隔的时间单位,f(δi)表示在δi处的值。
在一个实施例中,采用M个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据,拟合得到雷电流峰值与包络数据的初始拟合函数包括:
采用最小二乘法和滤波算法对M个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据进行拟合,得到雷电流峰值与包络数据的初始拟合函数。
在一个实施例中,初始拟合函数为三阶函数。
在一个实施例中,采用N个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对初始拟合函数中的参数进行标定,得到目标拟合函数包括:
采用第一组中每个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对初始拟合函数中的参数进行标定,得到第一目标拟合函数;
采用第二组中每个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对初始拟合函数中的参数进行标定,得到第二目标拟合函数,第一目标拟合函数的参数与第二目标拟合函数的参数不同。
根据本公开实施例的第二方面,提供一种雷电流峰值标定装置,包括:
第一获取模块,用于获取P个雷电流记录卡,P个雷电流记录卡为未被雷电冲击处于满磁状态的雷电流记录卡;
第二获取模块,用于将P个雷电流记录卡分别进行不同幅值的雷电冲击模拟实验,获取每个雷电流记录卡的雷电流峰值以及对应的磁条数据;
计算模块,用于计算每个雷电流记录卡的磁条数据的包络数据,包络数据包括包络面积或包络长度;
拟合模块,用于采用M个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据,拟合得到雷电流记录卡的雷电流峰值与包络数据的初始拟合函数;
标定模块,用于采用N个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对初始拟合函数中的参数进行标定,得到目标拟合函数,其中,P=M+N。
在一个实施例中,第二获取模块,用于获取雷电冲击电流的幅值范围,并将雷电冲击电流的幅值范围按照预设间隔进行划分,得到M个子间隔;将P个雷电流记录卡中的M个雷电流记录卡按照M个子间隔进行不同幅值的雷电冲击模拟实验,将M个子间隔的间隔端点最大值确定为每个雷电流记录卡的雷电流峰值;将P个雷电流记录卡中的N个雷电流记录卡分为两组,将第一组中的每个雷电流记录卡分别进行8/20μs的雷电冲击模拟试验,将第二组中的每个雷电流记录卡分别进行10/350μs的雷电冲击模拟试验,获取N个雷电流记录卡中每个雷电流记录卡的雷电流峰值。
在一个实施例中,第二获取模块,用于通过雷电流记录卡读卡器获取每个雷电流记录卡的磁条数据。
在一个实施例中,计算模块,用于在包络数据为包络面积时,采用第一公式计算每个雷电流记录卡的磁条数据的包络数据,磁条数据为曲边图形,包络数据为曲边图形的面积;
第一公式为:其中,S表示包络数据,n表示对曲边图形划分的区间个数,Δx表示划分的区间间隔的时间单位,f(δi)表示在δi处的值。
在一个实施例中,拟合模块,用于采用最小二乘法和滤波算法对M个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据进行拟合,得到雷电流峰值与包络数据的初始拟合函数。
在一个实施例中,初始拟合函数为三阶函数。
在一个实施例中,标定模块,用于采用第一组中每个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对初始拟合函数中的参数进行标定,得到第一目标拟合函数;
标定模块,用于采用第二组中每个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对初始拟合函数中的参数进行标定,得到第二目标拟合函数,第一目标拟合函数的参数与第二目标拟合函数的参数不同。
根据本公开实施例的第三方面,提供一种雷电流峰值标定装置,雷电流峰值标定装置包括处理器和存储器,存储器中存储有至少一条计算机指令,指令由处理器加载并执行以实现第一方面以及第一方面的任一实施例所描述的雷电流峰值标定方法中所执行的步骤。
根据本公开实施例的第四方面,提供一种计算机可读存储介质,存储介质中存储有至少一条计算机指令,指令由处理器加载并执行以实现第一方面以及第一方面的任一实施例所描述的雷电流峰值标定方法中所执行的步骤。
通过计算雷电流记录卡磁条数据的包络数据,结合最小二乘拟合算法,利用雷电实验室冲击电流发生器产生的8/20μs和10/350μs的标准雷电冲击电流进行模拟试验来完成雷电流峰值标定的高精度标定方法。通过计算磁条数据的包络数据,提高了雷电流记录卡的幅值检测范围,实现了0~300kA全量程检测。同时由于利用最小二乘拟合和滤波算法,及实验室标准数据标定,全量程精度提高到≤±10kA。进而提高了雷电流记录卡幅值读取的准确性和可靠性,为客户提供准确有效的雷击数据。
应当理解的是,以上的一般描述和后文的细节描述仅是示例性和解释性的,并不能限制本公开。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本公开的实施例,并与说明书一起用于解释本公开的原理。
图1是本公开实施例提供的一种雷电流峰值标定方法的流程图;
图2是本公开实施例提供的一种雷电流峰值标定方法的流程图;
图3是本公开实施例提供的一种50kA雷电流记录卡磁条数据图形;
图4是本公开实施例提供的一种雷电流峰值标定装置的结构图。
具体实施方式
这里将详细地对示例性实施例进行说明,其示例表示在附图中。下面的描述涉及附图时,除非另有表示,不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本公开相一致的所有实施方式。相反,它们仅是与如所附权利要求书中所详述的、本公开的一些方面相一致的装置和方法的例子。
本公开实施例提供一种雷电流峰值标定方法,如图1所示,该雷电流峰值标定方法包括以下步骤:
101、获取P个雷电流记录卡。
P个雷电流记录卡为未被雷电冲击处于满磁状态的雷电流记录卡。在本公开实施例中,雷电流记录卡采用磁条式雷电流记录卡,磁条式雷电流记录卡记录雷电流大小的主要原理是剩磁法。具体的,当导体流过电流时,它周围会产生强磁场,处于磁场中的磁性粒子会按磁力线方向排布,外加磁场消失后,磁性粒子排布方向会变的不一致,磁条形成自身磁性,即所谓的剩磁法。
102、将P个雷电流记录卡分别进行不同幅值的雷电冲击模拟实验,获取每个雷电流记录卡的雷电流峰值以及对应的磁条数据,并计算每个雷电流记录卡的磁条数据的包络数据。
在本公开实施例中,将P个雷电流记录卡分别进行不同幅值的雷电冲击模拟实验是指一个雷电流记录卡对应一次雷电冲击模拟实验。具体的,将P个雷电流记录卡分别进行不同幅值的雷电冲击模拟实验包括:
获取雷电冲击电流的幅值范围,并将雷电冲击电流的幅值范围按照预设间隔进行划分,得到M个子间隔;
将P个雷电流记录卡中的M个雷电流记录卡按照M个子间隔进行不同幅值的雷电冲击模拟实验;
将P个雷电流记录卡中的N个雷电流记录卡分为两组,将第一组中的每个雷电流记录卡分别进行8/20μs的雷电冲击模拟试验,将第二组中的每个雷电流记录卡分别进行10/350μs的雷电冲击模拟试验,获取N个雷电流记录卡中每个雷电流记录卡的雷电流峰值。
雷电冲击电流的幅值范围即自然界中雷电冲击电流,自然界中雷电冲击电流可达到0~300kA,将雷电冲击电流的幅值范围按照预设间隔进行划分,得到M个子间隔,示例性的,若预设间隔为10kA,则划分为30个子间隔;若预设间隔为20kA,则划分为15个子间隔;进而,将P个雷电流记录卡中的M个雷电流记录卡按照M个子间隔进行不同幅值的雷电流冲击模拟实验,即一个雷电流记录卡对应一个子间隔,同时,将P各雷电流记录卡中的N个雷电流记录卡分为两组,将第一组中的每个雷电流记录卡分别进行8/20μs的雷电冲击模拟试验,将第二组中的每个雷电流记录卡分别进行10/350μs的雷电冲击模拟试验,获取N个雷电流记录卡中每个雷电流记录卡的雷电流峰值。其中,第一组中雷电流记录卡的数量和第二组中雷电流记录卡的数量可以相同,也可以不同;8/20μs和10/350μs是两种不同的雷击波形,10/350μs为直接雷电冲击试验所用波形,8/20μs为间接雷击冲击试验所用波形。需要说明的是,雷电冲击模拟实验中雷电通过冲击电流发生器产生标准雷电冲击电流进行模拟实验。
在将雷电流记录卡进行雷电冲击模拟实验后,获取每个雷电流记录卡的雷电流峰值以及对应的磁条数据。其中,对于M个雷电流记录卡,将M各子间隔的间隔端点最大值确定为每个雷电流记录卡的雷电流峰值,以预设间隔为10kA为例,雷电冲击电流的幅值范围0~300kA分为30个子间隔,则需要30个雷电流记录卡进行雷电冲击模拟实验,这30个雷电流记录卡所记录的雷电流峰值分别为10kA、20kA、30kA、……、290kA、300kA。并且,通过雷电流记录卡读卡器读取每个雷电流记录卡的磁条数据,并计算每个雷电流记录卡的磁条数据的包络数据,包络数据包括包络面积或包络长度。
在包络数据为包络面积(即磁面积)时,对于计算每个雷电流记录卡的磁条数据的包络数据包括:采用第一公式计算每个雷电流记录卡的磁条数据的包络数据,磁条数据为曲边图形,包络数据为曲边图形的面积;其中,第一公式为:S表示包络数据,n表示对曲边图形划分的区间个数,Δx表示划分的区间间隔的时间单位,f(δi)表示在δi处的值。
当然,也可以计算雷电流记录卡磁条数据的包络长度,进而拟合包络长度与雷电流峰值的最小二乘函数,检测出雷电流峰值的大小。但是,由于磁条数据为不规则图形,与包络面积相比,包络长度不易准确算法,从而影响雷电流峰值的精度。
103、采用M个雷电流记录可的雷电流峰值以及对应的磁条数据的包络数据,拟合得到雷电流记录卡的雷电流峰值与包络数据的初始拟合函数。
在本公开实施例中,采用最小二乘法和滤波算法对M个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据进行拟合,得到雷电流峰值与包络数据的初始拟合函数,初始拟合函数为三阶函数。具体的,将包络数据S作为自变量,y为雷电流峰值,拟合函数y=f(S)。最小二乘多项式曲线拟合是假设给定数据由k次多项式函数生成,选择最有可能产生这些数据的k次多项式函数,即在k次多项式函数中选择一个对已知数据以及位置数据有很好预测能力的函数。给定函数y=f(S),在点S1、S2、S3、……、Sk处的函数值分别为y1、y2、y3、……、yk,求多项式y=a0+a1S+a2S2+……+akSk,依次选择一阶、二阶、三阶、四阶函数拟合,对拟合曲线图进行对比,最优为三阶函数,即y=a0+a1S+a2S2+a3S3,a0、a1、a2和a3为拟合函数参数。
104、采用N个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对初始拟合函数中的参数进行标定,得到目标拟合函数。
根据步骤103中所描述的,在确定后拟合函数后,需要对拟合函数中的参数进行标定,因此,结合步骤102中将第一组中的每个雷电流记录卡分别进行8/20μs的雷电冲击模拟试验,记录第一组中每个雷电流记录卡的雷电流峰值,并计算每个雷电流记录卡的磁条数据对应的包络数据;以及,将第二组中的每个雷电流记录卡分别进行10/350μs的雷电冲击模拟试验,记录第二组中每个雷电流记录卡的雷电流峰值,并计算每个雷电流记录卡的磁条数据对应的包络数据。进而,采用第一组中每个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对初始拟合函数中的参数进行标定,得到第一目标拟合函数;采用第二组中每个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对初始拟合函数中的参数进行标定,得到第二目标拟合函数,第一目标拟合函数的参数与第二目标拟合函数的参数不同。示例性的,采用第一组中的每个雷电流记录卡的数据,对初始拟合函数中的参数进行多次标定校准,反推出a0、a1、a2和a3的值,如a0=67.003、a1=-0.998、a2=0.00887、a3=-0.0000312。这样,在确定好拟合函数及拟合函数的参数后(即得到目标拟合函数),通过获取雷电流记录卡的磁条数据的包络数据,可实现对雷电流记录卡幅值的读取。
本公开实施例提供的雷电流方法,通过获取P个雷电流记录卡,将P个雷电流记录卡分别进行不同幅值的雷电冲击模拟实验,获取每个雷电流记录卡的雷电流峰值以及对应的磁条数据,并计算每个雷电流记录卡的磁条数据的包络数据;采用M个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据,拟合得到雷电流记录卡的雷电流峰值与包络数据的初始拟合函数;采用N个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对初始拟合函数中的参数进行标定,得到目标拟合函数。这样,弥补了现有斜率法中精度不高的确定,同时,可以将检测量程达到0~300kA,覆盖整个自然雷电的量程,可全方位实现雷电流记录卡的幅值读取,提高雷电流检测的准确性。
基于上述图1对应的实施例提供的雷电流峰值方法,本公开另一实施例提供一种雷电流峰值标定方法,该方法应用于雷电流记录卡读卡器,雷电流记录卡读卡器由采集电路、数据处理电路、自动读卡器、功能按键和液晶显示屏组成。参考图2所示,对雷电流记录卡读卡器进行初始化,将雷电流记录卡***该读卡器中并进行读卡按键操作,此时,采集电路进行AD采样采集磁条数据,数字化处理后转化为十六进制数据,数据处理电路将采集的磁条数据进行滤波和放大处理,计算出磁条数据的包络面积,即磁面积。图3所示为50kA雷电流记录卡磁条数据图形,横轴表示采集点数,纵轴表示磁条数据由于磁条数据为不规则图形,计算磁面积即曲边图形的面积,对于曲边图形的面积可近似地用矩形面积代替,分成间隔1个时间单位的小曲边图形,然后用矩形面积代替后求和。磁面积S的近似值:当Δx→0,/>n→∞。若取Δx=1,n=3000,求3000个小曲边图形面积得出磁面积。
定义满磁状态下,雷电流为0kA。通过冲击电流发生器产生标准雷电冲击电流进行模拟试验,将记录模拟雷电流幅值的雷电流记录卡作为标准卡,分别记录0~300kA范围内间隔10kA的各雷电流值(意味着需要30个雷电流记录卡)。通过读取存储有不同雷电流峰值的雷电流记录卡的磁面积,用最小二乘法和滤波算法拟合出磁面积与雷电流峰值对应关系的函数。依次选择一阶、二阶、三阶、四阶函数拟合。拟合曲线图对比,最优为三阶函数。最后,通过实验室冲击电流发生器产生的8/20μs和10/350μs的标准雷电冲击电流进行模拟试验,将记录模拟雷电流幅值的雷电流记录卡作为标准卡,对雷电流记录卡读卡器进行拟合函数参数标定,反推出拟合函数参数,并进行数据存储和显示。
本公开实施例提供了一种通过计算雷电流记录卡磁条数据的包络面积,结合最小二乘拟合算法,利用雷电实验室冲击电流发生器产生的8/20μs和10/350μs的标准雷电冲击电流进行模拟试验来完成雷电流峰值标定的高精度标定方法。该方法突出优势在于通过计算磁条数据的包络面积,提高了雷电流记录卡的幅值检测范围,实现了0~300kA全量程检测。同时由于利用最小二乘拟合和滤波算法,及实验室标准数据标定,全量程精度提高到≤±10kA。进而提高了雷电流记录卡幅值读取的准确性和可靠性,为客户提供准确有效的雷击数据。
基于上述图1对应的实施例中所描述的雷电流峰值标定方法,下述为本公开装置实施例,可以用于执行本公开方法实施例。
本公开实施例提供一种雷电流峰值标定装置,如图4所示,该雷电流峰值标定装置40包括:
第一获取模块401,用于获取P个雷电流记录卡,P个雷电流记录卡为未被雷电冲击处于满磁状态的雷电流记录卡;
第二获取模块402,用于将P个雷电流记录卡分别进行不同幅值的雷电冲击模拟实验,获取每个雷电流记录卡的雷电流峰值以及对应的磁条数据;
计算模块403,用于计算每个雷电流记录卡的磁条数据的包络数据;
拟合模块404,用于采用M个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据,拟合得到雷电流记录卡的雷电流峰值与包络数据的初始拟合函数;
标定模块405,用于采用N个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对初始拟合函数中的参数进行标定,得到目标拟合函数,其中,P=M+N。
在一个实施例中,第二获取模块402,用于获取雷电冲击电流的幅值范围,并将雷电冲击电流的幅值范围按照预设间隔进行划分,得到M个子间隔;将P个雷电流记录卡中的M个雷电流记录卡按照M个子间隔进行不同幅值的雷电冲击模拟实验,将M个子间隔的间隔端点最大值确定为每个雷电流记录卡的雷电流峰值;将P个雷电流记录卡中的N个雷电流记录卡分为两组,将第一组中的每个雷电流记录卡分别进行8/20μs的雷电冲击模拟试验,将第二组中的每个雷电流记录卡分别进行10/350μs的雷电冲击模拟试验,获取N个雷电流记录卡中每个雷电流记录卡的雷电流峰值。
在一个实施例中,第二获取模块402,用于通过雷电流记录卡读卡器获取每个雷电流记录卡的磁条数据。
在一个实施例中,包络数据包括包络面积或包络长度
在一个实施例中,计算模块403,用于在包络数据为包络面积时,采用第一公式计算每个雷电流记录卡的磁条数据的包络数据,磁条数据为曲边图形,包络数据为曲边图形的面积;
第一公式为:其中,S表示包络数据,n表示对曲边图形划分的区间个数,Δx表示划分的区间间隔的时间单位,f(δi)表示在δi处的值。
在一个实施例中,拟合模块404,用于采用最小二乘法和滤波算法对M个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据进行拟合,得到雷电流峰值与包络数据的初始拟合函数。
在一个实施例中,初始拟合函数为三阶函数。
在一个实施例中,标定模块,用于采用第一组中每个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对初始拟合函数中的参数进行标定,得到第一目标拟合函数;
标定模块,用于采用第二组中每个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对初始拟合函数中的参数进行标定,得到第二目标拟合函数,第一目标拟合函数的参数与第二目标拟合函数的参数不同。
本公开实施例提供的雷电流峰值标定装置,通过计算雷电流记录卡磁条数据的包络数据,结合最小二乘拟合算法,利用雷电实验室冲击电流发生器产生的8/20μs和10/350μs的标准雷电冲击电流进行模拟试验来完成雷电流峰值标定的高精度标定方法。通过计算磁条数据的包络数据,提高了雷电流记录卡的幅值检测范围,实现了0~300kA全量程检测。同时由于利用最小二乘拟合和滤波算法,及实验室标准数据标定,全量程精度提高到≤±10kA。进而提高了雷电流记录卡幅值读取的准确性和可靠性,为客户提供准确有效的雷击数据。
本公开实施例还提供了一种雷电流峰值标定装置,该雷电流峰值标定装置包括接收器、发射器、存储器和处理器,该发射器和存储器分别与处理器连接,存储器中存储有至少一条计算机指令,处理器用于加载并执行至少一条计算机指令,以实现上述图1对应的实施例中所描述的雷电流峰值标定方法。
基于上述图1对应的实施例中所描述的雷电流峰值标定方法,本公开实施例还提供一种计算机可读存储介质,例如,非临时性计算机可读存储介质可以是只读存储器(英文:Read Only Memory,ROM)、随机存取存储器(英文:Random Access Memory,RAM)、CD-ROM、磁带、软盘和光数据存储装置等。该存储介质上存储有计算机指令,用于执行上述图1对应的实施例中所描述的雷电流峰值标定方法,此处不再赘述。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
本领域技术人员在考虑说明书及实践这里公开的公开后,将容易想到本公开的其它实施方案。本申请旨在涵盖本公开的任何变型、用途或者适应性变化,这些变型、用途或者适应性变化遵循本公开的一般性原理并包括本公开未公开的本技术领域中的公知常识或惯用技术手段。说明书和实施例仅被视为示例性的,本公开的真正范围和精神由下面的权利要求指出。
Claims (8)
1.一种雷电流峰值标定方法,其特征在于,所述方法包括:
获取P个雷电流记录卡,所述P个雷电流记录卡为未被雷电冲击处于满磁状态的雷电流记录卡;
将所述P个雷电流记录卡分别进行不同幅值的雷电冲击模拟实验,获取所述每个雷电流记录卡的雷电流峰值以及对应的磁条数据,并计算所述每个雷电流记录卡的磁条数据的包络数据;
采用M个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据,拟合得到雷电流记录卡的雷电流峰值与包络数据的初始拟合函数;
采用N个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对所述初始拟合函数中的参数进行标定,得到目标拟合函数,其中,P=M+N;
所述将所述P个雷电流记录卡分别进行不同幅值的雷电冲击模拟实验,获取每个雷电流记录卡的雷电流峰值包括:
获取雷电冲击电流的幅值范围,并将所述雷电冲击电流的幅值范围按照预设间隔进行划分,得到M个子间隔;
将所述P个雷电流记录卡中的M个雷电流记录卡按照所述M个子间隔进行不同幅值的雷电冲击模拟实验,将所述M个子间隔的间隔端点最大值确定为所述每个雷电流记录卡的雷电流峰值;
将所述P个雷电流记录卡中的N个雷电流记录卡分为两组,将第一组中的每个雷电流记录卡分别进行8/20μs的雷电冲击模拟试验,将第二组中的每个雷电流记录卡分别进行10/350μs的雷电冲击模拟试验,获取N个雷电流记录卡中每个雷电流记录卡的雷电流峰值。
2.根据权利要求1所述的方法,其特征在于,所述获取所述每个雷电流记录卡的磁条数据包括:
通过雷电流记录卡读卡器获取所述每个雷电流记录卡的磁条数据。
3.根据权利要求1所述的方法,其特征在于,所述包络数据包括包络面积或包络长度。
4.根据权利要求1所述的方法,其特征在于,所述计算所述每个雷电流记录卡的磁条数据的包络数据包括:
在包络数据为包络面积时,采用第一公式计算所述每个雷电流记录卡的磁条数据的包络数据,所述磁条数据为曲边图形,所述包络数据为曲边图形的面积;
所述第一公式为:其中,S表示包络数据,n表示对所述曲边图形划分的区间个数,Δx表示划分的区间间隔的时间单位,f(δi)表示在δi处的值。
5.根据权利要求1所述的方法,其特征在于,所述采用M个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据,拟合得到雷电流峰值与包络数据的初始拟合函数包括:
采用最小二乘法和滤波算法对所述M个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据进行拟合,得到雷电流峰值与包络数据的初始拟合函数。
6.根据权利要求5所述的方法,其特征在于,所述初始拟合函数为三阶函数。
7.根据权利要求1所述的方法,其特征在于,所述采用N个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对所述初始拟合函数中的参数进行标定,得到目标拟合函数包括:
采用第一组中每个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对所述初始拟合函数中的参数进行标定,得到第一目标拟合函数;
采用第二组中每个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对所述初始拟合函数中的参数进行标定,得到第二目标拟合函数,所述第一目标拟合函数的参数与所述第二目标拟合函数的参数不同。
8.一种雷电流峰值标定装置,其特征在于,包括:
第一获取模块,用于获取P个雷电流记录卡,所述P个雷电流记录卡为未被雷电冲击处于满磁状态的雷电流记录卡;
第二获取模块,用于将所述P个雷电流记录卡分别进行不同幅值的雷电冲击模拟实验,获取所述每个雷电流记录卡的雷电流峰值以及对应的磁条数据;
计算模块,用于计算所述每个雷电流记录卡的磁条数据的包络数据,所述包络数据包括包络面积或包络长度;
拟合模块,用于采用M个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据,拟合得到雷电流记录卡的雷电流峰值与包络数据的初始拟合函数;
标定模块,用于采用N个雷电流记录卡的雷电流峰值以及对应的磁条数据的包络数据对所述初始拟合函数中的参数进行标定,得到目标拟合函数,其中,P=M+N;
所述第二获取模块,用于获取雷电冲击电流的幅值范围,并将所述雷电冲击电流的幅值范围按照预设间隔进行划分,得到M个子间隔;将所述P个雷电流记录卡中的M个雷电流记录卡按照所述M个子间隔进行不同幅值的雷电冲击模拟实验,将所述M个子间隔的间隔端点最大值确定为所述每个雷电流记录卡的雷电流峰值;将所述P个雷电流记录卡中的N个雷电流记录卡分为两组,将第一组中的每个雷电流记录卡分别进行8/20μs的雷电冲击模拟试验,将第二组中的每个雷电流记录卡分别进行10/350μs的雷电冲击模拟试验,获取N个雷电流记录卡中每个雷电流记录卡的雷电流峰值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110542392.4A CN113447701B (zh) | 2021-05-18 | 2021-05-18 | 雷电流峰值标定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110542392.4A CN113447701B (zh) | 2021-05-18 | 2021-05-18 | 雷电流峰值标定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113447701A CN113447701A (zh) | 2021-09-28 |
CN113447701B true CN113447701B (zh) | 2023-09-01 |
Family
ID=77810061
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110542392.4A Active CN113447701B (zh) | 2021-05-18 | 2021-05-18 | 雷电流峰值标定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113447701B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118114027B (zh) * | 2024-04-26 | 2024-07-19 | 西安爱邦电磁技术有限责任公司 | 雷电流峰值采集方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012026826A (ja) * | 2010-07-22 | 2012-02-09 | Sanyo Electric Industries Co Ltd | 雷サージ電流検出装置 |
JP2012037448A (ja) * | 2010-08-10 | 2012-02-23 | Tokyo Electric Power Co Inc:The | 雷インパルス耐電圧試験システムおよび基準波形算出プログラム、並びに雷インパルス耐電圧試験方法 |
CN103616553A (zh) * | 2013-12-05 | 2014-03-05 | 上海高试电气科技有限公司 | 短路大电流波形拟合方法 |
CN104239739A (zh) * | 2014-09-26 | 2014-12-24 | 华南理工大学 | 一种雷电流幅值概率分布函数的分段拟合方法与*** |
CN105866558A (zh) * | 2016-03-23 | 2016-08-17 | 清华大学深圳研究生院 | 一种小区域雷电监测定位方法及*** |
CN107290664A (zh) * | 2017-06-27 | 2017-10-24 | 中国船舶重工集团公司第七�三研究所 | 一种发电机组波形记录仪的瞬态参数自动计算方法 |
CN108169540A (zh) * | 2017-12-13 | 2018-06-15 | 内蒙古工业大学 | 一种风力发电机组电压闪变的测量方法 |
CN108693397A (zh) * | 2018-06-12 | 2018-10-23 | 云南电网有限责任公司昆明供电局 | 一种雷电流峰值估算方法及装置 |
CN110161301A (zh) * | 2019-06-24 | 2019-08-23 | 四川电安智能科技有限公司 | 一种抗干扰型雷电流波形采集算法 |
CN110380659A (zh) * | 2019-06-27 | 2019-10-25 | 国电南瑞科技股份有限公司 | 一种永磁同步电机快速标定方法、***和控制器 |
CN112114185A (zh) * | 2020-09-01 | 2020-12-22 | 华帝股份有限公司 | 基于导数算法的电网峰值采样方法 |
CN112362947A (zh) * | 2020-10-23 | 2021-02-12 | 熊秀 | 一种雷电供电峰值电流记录卡 |
-
2021
- 2021-05-18 CN CN202110542392.4A patent/CN113447701B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012026826A (ja) * | 2010-07-22 | 2012-02-09 | Sanyo Electric Industries Co Ltd | 雷サージ電流検出装置 |
JP2012037448A (ja) * | 2010-08-10 | 2012-02-23 | Tokyo Electric Power Co Inc:The | 雷インパルス耐電圧試験システムおよび基準波形算出プログラム、並びに雷インパルス耐電圧試験方法 |
CN103616553A (zh) * | 2013-12-05 | 2014-03-05 | 上海高试电气科技有限公司 | 短路大电流波形拟合方法 |
CN104239739A (zh) * | 2014-09-26 | 2014-12-24 | 华南理工大学 | 一种雷电流幅值概率分布函数的分段拟合方法与*** |
CN105866558A (zh) * | 2016-03-23 | 2016-08-17 | 清华大学深圳研究生院 | 一种小区域雷电监测定位方法及*** |
CN107290664A (zh) * | 2017-06-27 | 2017-10-24 | 中国船舶重工集团公司第七�三研究所 | 一种发电机组波形记录仪的瞬态参数自动计算方法 |
CN108169540A (zh) * | 2017-12-13 | 2018-06-15 | 内蒙古工业大学 | 一种风力发电机组电压闪变的测量方法 |
CN108693397A (zh) * | 2018-06-12 | 2018-10-23 | 云南电网有限责任公司昆明供电局 | 一种雷电流峰值估算方法及装置 |
CN110161301A (zh) * | 2019-06-24 | 2019-08-23 | 四川电安智能科技有限公司 | 一种抗干扰型雷电流波形采集算法 |
CN110380659A (zh) * | 2019-06-27 | 2019-10-25 | 国电南瑞科技股份有限公司 | 一种永磁同步电机快速标定方法、***和控制器 |
CN112114185A (zh) * | 2020-09-01 | 2020-12-22 | 华帝股份有限公司 | 基于导数算法的电网峰值采样方法 |
CN112362947A (zh) * | 2020-10-23 | 2021-02-12 | 熊秀 | 一种雷电供电峰值电流记录卡 |
Also Published As
Publication number | Publication date |
---|---|
CN113447701A (zh) | 2021-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Taylor | Summarizing multiple aspects of model performance in a single diagram | |
KR102068028B1 (ko) | 피뢰기 저항성 누설전류 검출 장치 및 검출 방법 | |
US10097000B2 (en) | Tool employing homotopy-based approaches in finding the controlling unstable equilibrium point in the electric power grid | |
CN113447701B (zh) | 雷电流峰值标定方法及装置 | |
US20150177279A1 (en) | Electromagnetic wave identification method and identification device | |
JP6873794B2 (ja) | 中性子強度モニタリングシステムおよび方法 | |
CN105806787A (zh) | 长空气间隙雷击放电电子密度的测量***及方法 | |
EP3012643A1 (en) | Method and apparatus for identifying causes for cable overcurrent | |
Elter et al. | Performance of Higher Order Campbell methods, Part II: calibration and experimental application | |
JP2007315911A (ja) | 落雷電界量算出装置及び方法、プログラム、コンピュータ読み取り可能な記録媒体 | |
Yousof et al. | Study of transformer winding deformation by frequency response analysis | |
CN110907879B (zh) | 一种匝间耐压测试仪的校准方法 | |
KR101420729B1 (ko) | 접지망 진단 장치 및 방법 | |
Frei et al. | A transient field measurement system to analyze the severity and occurrence rate of electrostatic discharge (ESD) | |
KR20190064075A (ko) | 일사량 데이터 평가 장치 및 방법 | |
JP3462257B2 (ja) | 部分放電検出方法 | |
Henry et al. | Transmission-Line Pulse ESD Testing of ICs: A New Beginning | |
CN117647694B (zh) | 一种适用于智能手表加工过程的质量检测方法 | |
TWI467212B (zh) | 建築物樓層之地震即時分析系統及其方法與儲存媒體 | |
US20030060992A1 (en) | Method and apparatus for accurately determining the crossing point within a logic transition of a differential signal | |
CN115877112A (zh) | 基于特高频电磁信号的金属化膜电容器自愈特征监测方法 | |
Clandfield | Analysing the current noise equalization in VeloPix | |
Townsend et al. | A 3D fault location program for use in high-voltage substations using multilateration, geolocation and machine learning | |
Lavatelli et al. | A probabilistic method to assess uncertainty in vision based modal analysis techniques | |
Denholm-Price | Can an ensemble give anything more than Gaussian probabilities? |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230919 Address after: 710075 building B4, civil military integration industrial park, 176 Biyuan 2nd Road, high tech Zone, Xi'an City, Shaanxi Province Patentee after: XI'AN AIRBORNE ELECTROMAGNETIC TECHNOLOGY Co.,Ltd. Address before: 710075 building B4, civil military integration industrial park, 176 Biyuan 2nd Road, high tech Zone, Xi'an City, Shaanxi Province Patentee before: Xiong Xiu |