CN113050175A - 直升机航空电磁收录装置及发射源参数识别方法 - Google Patents
直升机航空电磁收录装置及发射源参数识别方法 Download PDFInfo
- Publication number
- CN113050175A CN113050175A CN202110250033.1A CN202110250033A CN113050175A CN 113050175 A CN113050175 A CN 113050175A CN 202110250033 A CN202110250033 A CN 202110250033A CN 113050175 A CN113050175 A CN 113050175A
- Authority
- CN
- China
- Prior art keywords
- data
- receiver
- current
- value
- sequence
- 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
- 238000000034 method Methods 0.000 title claims abstract description 41
- 230000005284 excitation Effects 0.000 claims abstract description 15
- 238000001514 detection method Methods 0.000 claims abstract description 14
- 230000001360 synchronised effect Effects 0.000 claims abstract description 13
- 238000005070 sampling Methods 0.000 claims description 17
- 230000003321 amplification Effects 0.000 claims description 12
- 238000005311 autocorrelation function Methods 0.000 claims description 12
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 12
- 230000005540 biological transmission Effects 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 10
- 238000013519 translation Methods 0.000 claims description 7
- 238000005314 correlation function Methods 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 4
- 239000004033 plastic Substances 0.000 claims description 4
- 229920003023 plastic Polymers 0.000 claims description 4
- 230000007246 mechanism Effects 0.000 claims description 3
- 230000001681 protective effect Effects 0.000 claims description 3
- 230000004044 response Effects 0.000 claims description 3
- 125000004122 cyclic group Chemical group 0.000 claims description 2
- 238000010606 normalization Methods 0.000 claims description 2
- 230000006698 induction Effects 0.000 abstract description 9
- 230000008859 change Effects 0.000 abstract description 4
- 239000000306 component Substances 0.000 description 29
- 239000007787 solid Substances 0.000 description 11
- 230000008569 process Effects 0.000 description 7
- 230000000875 corresponding effect Effects 0.000 description 6
- 238000002955 isolation Methods 0.000 description 6
- 238000004804 winding Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 4
- 230000003750 conditioning effect Effects 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 239000010410 layer Substances 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000001052 transient effect Effects 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 238000012850 discrimination method Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 230000005355 Hall effect Effects 0.000 description 1
- 239000004677 Nylon Substances 0.000 description 1
- 239000004698 Polyethylene Substances 0.000 description 1
- 239000004809 Teflon Substances 0.000 description 1
- 229920006362 Teflon® Polymers 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 230000005674 electromagnetic induction Effects 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005290 field theory Methods 0.000 description 1
- 230000017525 heat dissipation Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000000696 magnetic material Substances 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229920001778 nylon Polymers 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- -1 polyethylene Polymers 0.000 description 1
- 229920000573 polyethylene Polymers 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 230000008054 signal transmission Effects 0.000 description 1
- 238000001228 spectrum Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
- G01V3/083—Controlled source electromagnetic [CSEM] surveying
- G01V2003/085—Receivers
Abstract
本发明属于航空电磁探测领域，涉及一种直升机航空电磁收录装置及发射源参数识别方法，通过直升机平台搭载，包括磁场传感器和接收机，所述磁场传感器固定于直升机下方吊舱，通过长距离信号线与机舱内接收机连接，所述接收机的控制器内运行高频m序列编码发射源参数识别软件，接收机通过控制器启动采集卡采集数据；采集卡板载时钟通过PXIe背板路由至计数卡，经计数卡分频后作为同步信号同步发射机；所述发射机内部使用电流传感器收录激发电流数据，通过线缆传输至接收机，所述接收机运行高频m序列编码发射源参数识别软件通过采集的发射机的电流数据，对发射机进行参数识别。能够准确感知其二次感应场的变化，精确快速地收录海量数据。
Description
技术领域
本发明属于航空电磁探测领域，具体而言涉及适用于直升机平台的针对高频m序列发射源的宽频带三分量磁场传感器和高采样率、高速电磁法接收机的直升机航空电磁收录装置及发射源参数识别方法。
背景技术
直升机航空电磁法是航空电磁法(AEM)的一种。直升机电磁法勘探系统以直升机为载体，在飞行过程中，通过吊舱上发射线圈产生的大发射磁矩激发地下介质，地下介质因涡流效应产生二次感应场，利用接收装置对二次场信号进行接收，根据电磁场理论实现对地下介质电阻率结构的解释。直升机航空电磁法收录系统包括接收机和磁场传感器，是直升机航空电磁法勘探系统的核心组成部分。
传统的航空时间域电磁法一般采用基频为25Hz、75Hz、125Hz等的双极性多边形发射电流，由于发射源频率较低对地下浅层异常体丢失了探测精度，由于接收机采样率限制、一次场干扰等影响勘探结果，存在浅层探测盲区。
m序列，全称最长线性移位寄存器序列，作为目前应用最广泛的一种伪随机序列已被应用于通信领域，如扩频通信、卫星通信的码分多址、数字数据中的加密等。在地球物理领域，伪随机编码作为发射波形已被应用于地面和海洋的多道瞬变电磁系统(MTEM)，实现了矿产资源的高分辨探测，说明其拥有更强的分辨能力，展现出了浅层探测的潜力。近年来首次提出将高频伪随机码源应用于直升机航空平台，旨在将高频伪随机序列码源低噪声浅层探测能力与直升机航空电磁探测平台低成本、高灵活性、更适合复杂地形调查的优势相结合。该系统采用200kHz码元频率的伪随机编码发射电流，对收录系统从传感器带宽、接收机采样率到实时数据存储速率的要求均远高于以往。
直升机航空电磁探测系统单架次飞行3小时，其间工作人员无法干涉仪器设备运行。发射机与接收机不同源导致的两系统时钟漂移、负载对发射电流的影响和发射机故障等可能性必须考虑，若设备工作异常，将造成较大的损失。因此，飞行前，在地面对系统工作状态做出准确判断，及时做出修复和飞行中及时发现工作异常，及时建议返航挽回损失都极为重要。传统的直升机航空电磁探测系统采用的双极性多边形发射电流(如双极性三角波、双极性梯形波等)有明显的时域特征，观测方便；利用感应信号抽道实时显示配合短暂闭合安装于吊舱的人工异常线圈的方法可了解接收系统状态。而伪随机编码作为航空电磁探测发射源，难于直接在时域观测，且不适用常规抽道方法处理，因此无法直接沿用以往判断系统工作状态的方法。
发明内容
本发明所要解决的技术问题在于提供一种直升机航空电磁收录装置及发射源参数识别方法，适用于高频m序列激发源的，能够准确感知其二次感应场的变化，精确快速地收录海量数据。
本发明是这样实现的，
种高频m序列编码发射源参数识别方法，包括如下步骤：
A)循环录入长度不断增加的m波形电流数据进行线性自相关运算，直至出现相邻两个明显尖峰，其间隔为m波形电流的周期，对单周期m波形电流数据做循环自相关运算，取1/5T～4/5T内数据取均值得到循环自相关函数最小值min的估计值；
B)采用最小值min的估计值，利用公式
求得阶数估计值N，对阶数估计值取整，得到阶数；
C)根据公式
并取整，求得码元频率估计值；
D)根据已知的采样率，以及求取的阶数和码元频率，使用当前阶数下所有可能的反馈系数，生成可能的m序列，与待测m波形进行对比，判断待测序列最可能采用的反馈系数；
E)计算待测m波形与参考m序列重合所需的平移量，定义为其相位差，以参考m序列为参考，右移为平移正参考方向，将参考m序列和待测m波形做互相关运算，尖峰位置的互相关函数偏移量，为所求相位；
F)上述步骤所求得包括周期、阶数、码元频率、反馈系数和相位的发射参数，对照发射机理论参数，综合评价其合理性，对航空电磁探测发射系统状态的基本判断。
进一步地，所述步骤A具体包括：
A1)获取N点发射机电流波形发射时间内)双极性电流数据，并取平均值，发射电流数据直接减去该均值以校正发射电流偏置，获得无偏置双极性发射电流；
A2)将无偏置双极性电流数据翻转至正向，并取均值，电流幅值在均值以下的数据点视为低幅值噪声，将其全部置零，将去除低幅值噪声的数据求得各极值点，并取均值，该均值数据被视为双极性发射电流幅值I的估计值；
A3)截取适当点数步骤A1中的无偏置双极性发射电流数据；
A4)调整阈值，使用滞回比较方法对电流数据进行二值化；
A5)二值化后数据线性自相关并归一化；
A6)小噪声数据置零，求极值，存入极值数组；
A7)判断极值数组长度是否大于1，若大于则进行下一步；
A8)取两个尖峰，其横坐标间隔为m波形周期T；
A9)按照步骤A8求取的周期，截取一个周期步骤A1中的无偏置双极性发射电流数据；
A10)对步骤A9截取的数据做循环自相关；
A11)对T/5～4T/5内数据取均值，作为循环自相关函数最小值min的估计值.
进一步地，步骤D包括:
D1)读取反馈系数表；
D2)由已知采样频率Fs，求得的码元频率，阶数n，读取n对应的有限m个反馈系数，生成m序列；
D3)生成的m序列与步骤A9中截取的一个周期无偏置双极性待测发射电流做等长化处理；
D4)生成的m序列与步骤D5中截取的待测电流做互相关处理；
D5)互相关函数取绝对值后求极值，并减去互相关函数绝对值的均值，差值为Dmm；
D6)判断阶数n对应的反馈系数是否遍历完毕；若是则进行下一步；
D7)找到差值Dmm的最大值对应的反馈系数；
D8)得出待测电流对应的反馈系数。
一种基于权利要求1所述方法的直升机航空电磁收录装置，通过直升机平台搭载，包括磁场传感器和接收机，所述磁场传感器固定于直升机下方吊舱，通过长距离信号线与机舱内接收机连接，所述接收机的控制器内运行高频m序列编码发射源参数识别软件，接收机通过控制器启动采集卡采集数据；采集卡板载时钟通过PXIe背板路由至计数卡，经计数卡分频后作为同步信号同步发射机；所述发射机内部使用电流传感器收录激发电流数据，通过线缆传输至接收机，所述接收机运行高频m序列编码发射源参数识别软件通过采集的发射机的电流数据，对发射机进行参数识别。
进一步地，所述磁场传感器包括三分量空心接收线圈、前置放大电路、保护固定结构以及传输线缆，三分量接收线圈由相互垂直的三个空心线圈组成，其中Z分量的空心线圈为圆形，X分量的空心线圈与Y分量的空心线圈组成一个正方体，其中X分量的空心线圈包括一组平行的正方形差分空心线圈，另一组正方形差分空心线圈组成Y分量的空心线圈，X分量的空心线圈与Y分量的空心线圈平行于吊舱平面，使用塑料螺丝将线圈安装于保护结构内，整体固定于吊舱中心处，使接收线圈与发射线圈为同心圆；或固定于吊舱偏心处，使接收线圈与发射线圈不构成同心圆；信号线跟随吊舱绳索向上延伸，经过固定于绳索上的前置放大电路后，连接至直升机机舱内的接收机。
进一步地，发射机内使用电流传感器感知激发电流信号与接收机接收的电磁响应信号一同经采集卡转换，各通道数据通过总线传递至控制器。
本发明与现有技术相比，有益效果在于：本发明方法采用高频m序列发射时，针对m序列编码的发射波形采用参数识别方法判断系统工作状态，并且参数可供进一步数据分析。本装置传感器收录三个分量磁场信息，并大幅度提升有效带宽，避免高频信息丢失，适应高频m序列激发场下的磁场信息收录；接收机针对需求下的1.25MSa/s/ch(每秒每通道一兆个采样点)的高频6通道采样产生的的海量数据。
附图说明
图1是直升机航空电磁法勘察系统组成图；
图2是直升机航空电磁法三分量磁场传感器结构图；
图3是直升机航空电磁法接收机内部结构框图；
图4是直升机航空电磁法接收机收录软件界面；
图5是系统发射参数识别方法流程图；
图6是伪随机发射电流数据；
图7是参数识别结果。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
参见图1结合图2所示，高频m序列编码发射源的直升机航空电磁法收录系统由磁场传感器和接收机两部分组成。其中，三分量磁场传感器由三分量空心接收线圈4和前置放大电路2两个基本部分组成，并包含必需的传感器第一保护固定结构43、第二保护固定结构44和传输线缆3，其中。三分量接收线圈由相互垂直的三个空心线圈组成，其中Z分量(垂直于吊舱平面方向)的空心线圈41为圆形，X分量空心线圈绕制于由第一线槽421、第二线槽422、第三线槽423和第四线槽424组成的线框以及对称的四个线槽组成的骨架结构上，构成差分线圈结构；同理，Y分量空心线圈绕制于由第五线槽425、第六线槽426、第七线槽427和第八线槽428组成的线框以及对称的四个线槽组成的骨架结构上，构成差分线圈结构，X、Y分量平行于吊舱平面。各线圈的谐振频率参数应通过绕线匝数等方式调整一致。使用塑料螺丝将线圈安装于保护结构内，整体固定于吊舱中心处，使接收线圈4与发射线圈5为同心圆，也可以固定于吊舱偏心处，使接收线圈4与发射线圈5不构成同心圆。信号线3跟随吊舱绳索向上延伸，经过固定于绳索上的前置放大电路2后，连接至直升机机舱内的接收机1，该信号线缆总长度在30至40米之间。
三分量磁场传感器空心线圈结构采用无磁材料制作。Z分量接收线圈41为直径0.5米的圆形空心线圈，采用中心抽头绕线方式的铁氟龙镀银铜线绕制而成，Z分量接收线圈尼龙骨架凿有12个等间隔圆形通孔，用于固定。X和Y轴接收线圈42为边长0.2米的正方形空心线圈，二者相互垂直，组成一个镂空的正方体形状。底部四个棱凿有4个等间距的圆形通孔，用于固定。三个垂直分量空心线圈利用圆形通孔，使用航空塑料螺丝可靠固定于圆形的高分子聚乙烯材料底座44，并使用半球形亚克力透明保护罩43与底座44紧密固定，将接收线圈完全包裹，避免飞行过程中气流对线圈的扰动。为抑制长距离信号传输过程中产生的共模噪声，空心线圈均采取的中心抽头绕线方式形成了对称的差分结构。为提升传感器带宽，在吊舱体积、重量等限制下，通过控制空心线圈匝数、绕线层数、绕线方式和增加隔层等方式，将三个分量的空心线圈谐振频率均调整至150kHz以上，减小信号高频成分畸变。
三分量空心接收线圈三路信号单独放大，采用两级三运放结构的仪用前置放大电路2。前级选用低电压电流噪声的通用放大器，后级采用全差分放大器，使用精密电阻保证电桥平衡。在匹配差分空心线圈三端输出结构的同时，能够提供极高共模和差模输入阻抗，很低的输出阻抗，稳定的、较高的放大倍数和极高的共模抑制比。前置放大电路的设计考虑空心线圈高带宽特性，在保证较高增益前提下使空心线圈带宽内具有平稳的幅频特性，避免放大环节产生高频信号畸变。前置放大电路前端设计有匹配电路，使用精密电阻与空心线圈输出端并联，用以调整线圈的阻尼比，避免信号时域震荡。
参见图3，高采样率直升机航空电磁法接收机1由直升机搭载，固定于机舱内标准机架上相应位置，由直升机发电机提供28V直流电源。在接收机外部，前面板内嵌高亮工业显示器用于保证野外强光环境下的人机交互，内嵌按键矩阵用于控制接收机总电源、前置放大电路电源、同步信号、触发信号、辅助功能等开关；后面板内嵌连接各信号线、数据线与电源线的航空连接器；接收机机箱内部，固定着PXIe机箱、电源转接电路板和信号转接电路板。PXIe机箱内包含3U尺寸的控制器13、采集卡14、计数卡15，及直升机航空条件下必需的28V直流转ATX电源模块，各板卡通过PXIe背板总线12实现信息高速传递；电源管理19转接电路板分别供给PXIe机箱、工业显示器、数字隔离与信号调理模块以及磁场传感器前置放大电路所需的28伏、12伏、5伏和正负15伏直流电源；信号转接板包含了必要的通信接口转换、开关控制、光耦隔离、信号单端化调理功能和必要的保护电路。
高采样率直升机航空电磁法收录系统实现对二次场感应信息的完整收录，包括以下步骤：
参见图3，高采样率直升机航空电磁法接收机与操作人员通过VGA接口连接的内嵌工业显示器161和通过USB接口连接的鼠标163键盘162实现人机交互功能。仪器启动后，采集卡15开始采集，且产生的与设置的采样率关联的时钟信号，经PXIe背板触发总线12路由至计数卡，分频至50Hz后作为同步信号，通过光耦隔离，经同步线缆同步发射机。发射系统正常工作后，包含一次激发场与二次感应场的磁场信号被差分结构的三分量空心线圈获取，感应电压由传输线缆接入信号调理模块18，分别经过滤波、放大和单端化处理后，由电磁接收机内高采样率采集卡14以1.25MSa/s/ch(每秒每通道一兆个采样点)的速率完成高速模数转换。同时，发射机内使用电流传感器感知激发电流信号，与电磁响应信号和辅助装置17信号一同经采集卡14转换，各通道数据通过背板总线12传递至控制器13内存，由软件控制实时存储至内部固态硬盘及外置固态硬盘。
直升机提供的28伏直流电源经短接保护、过流保护、无源滤波和DC-DC隔离处理后，使用合适功率的多块DC-DC模块，分别给PXIe机箱、工业显示器、数字隔离与信号调理和接收传感器前置放大电路所需的28伏、12伏、5伏和正负15伏直流电源。接收机机箱内所用到的DC-DC隔离电源模块均设计了对应的散热装置，保证运行于适宜的工作温度环境。
开始采集前，需要经过以下操作：
参见图4，接收机1经十分钟以上预热后，可准备开始工作。通过接收机前面板内嵌工业显示器161和外接键盘162、鼠标163，对图4所示的收录软件进行操作。根据实际勘测任务需要、发射系统参数等，在“工作设置”选项卡中，设置发射信号基频、采样频率，电流增益、主存储路径(内置固态硬盘)和从存储路径(方便航次间数据转移的外置固态硬盘)及其文件格式、存储目录名、工作描述、航次间隔、测线间隔、背景线间隔，并可输入TEM装置结构简要说明文字。
2.待发射系统就绪后，通过接收前面板上开关给传感器前置放大电路供电。点击收录软件上“启动采集”按钮，按钮由红变绿，此时，采集卡与计数卡开始工作，若此时查看“核对波形”选项卡，可见各通道外部噪声信号波形。在采集卡工作的同时，由采集卡板载时钟产生的时钟滴答信号经背板路由至计数卡，经分频后产生50Hz方波信号，经过信号转接电路光耦隔离后，输出幅值12V的50Hz方波，作为同步信号同步发射机。发射机收到同步信号后，方能开始发射瞬变电流。
3.点击“启动采集”按钮后，接收机采集卡和计数卡开始工作，但数据尚未开始存储。在开启发射机前，应点击收录软件上“开始飞行”按钮(或F2快捷键)，按钮由红变绿，接收机开始按1秒钟一个文件的方式存储采集的数据至固态硬盘，内置固态硬盘与外置固态硬盘数据内容完全一致，外置固态硬盘可选择使用数据量更小的格式。
4.在发射机正常发射瞬变电流后，一次激发场和二次感应场被三分量磁场传感器以三个相互垂直的方向感知，根据法拉第电磁感应定律，产生与磁场变化率相关的感应电动势，三个通路的微弱感应电压信号经前置放大电路放大，传输至直升机机舱内接收机；与此同时，发射机内使用基于霍尔效应的电流传感器获取发射电流波形，也作为一个信号通道，汇总至接收机，由采集卡以24位、1.25MSa/s/ch的参数进行高速高精度模数转换，并由收录软件控制，以可靠方式存储至固态硬盘中。
5.在实际勘察工作中，直升机拉升至900米高空时，应点击“开始背景”按钮，存储背景场数据，便于地面后期数据处理中去除背景场干扰。在正常飞行高度时，点击“开始测线”按钮，开始单条测线的数据收录。
6.直升机进入每条测线都应重新点击“开始测线”按钮，在内存中按文件夹组织存储测线数据。当本次飞行勘察任务结束，飞机降落在机场时，应将外置固态硬盘弹出系统，连接至地面工作站处理，并为接收机替换备用固态硬盘，勘察系统可及时开展下一航次飞行任务。
针对m序列发射波形的系统工作状态判别方法，从接收系统角度出发，在仅已知采样频率的条件下对发射机发射的电流波形数据进行实时的参数识别，实测参数与发射电流理论参数进行对照，给出发射机至接收机整个系统工作状态的提示。特别地，接收机内处理器资源将大部分倾斜于高采样率的数据全波收录过程，因此，参数识别方法应当避免过于复杂的运算，权衡执行效率与效果。
进一步的，系统工作状态判别方法基于m序列具有良好的自相关特性，其自相关函数类似白噪声的自相关函数(δ函数)。该方法对电流采样数据进行分析，实时识别发射参数包括周期、阶数、码元频率、反馈系数和相位差，具体而言参见图5：
A、周期，即m序列生成时各寄存器的初始值。参与自相关运算的m波形长度大于一个周期时，其自相关函数将会呈现一个以上的尖峰，据此，截取1至2个周期适当长度的m波形数据进行自相关运算，相邻两个明显尖峰的间隔即为其周期。实际电流波形因为噪声的存在，这一过程将直接导致其自相关函数各尖峰附近出现抖动，难以确定各尖峰横坐标。考虑到循环截取序列长度过程中，自相关函数尖峰附近的噪声频率也处于动态的变化之中，常规的滤波方法不易成为简洁高效的消噪选择。参见图6，本发明提出，首先以滞回比较的方法对电流波形进行二值化处理，对二值化后的波形做自相关，回避对自相关函数曲线消噪的环节。合适的阈值使滞回比较后产生的二值波形不再存在边沿处的震荡，从根源上去除了自相关函数尖峰附近抖动对求极值产生的干扰，达到简单快速求取各尖峰位置的目的。在二值化的过程中，滞回比较的阈值的选择应在合理范围内，阈值过小，将会产生误跳变，边沿处的二值化抖动的效应同样会体现在其自相关函数的尖峰周围；阈值过大，跳变条件过于严苛，将会带来二值化后相位的偏移，码元宽度失真，将带来二值化误差。
B、阶数，即m序列生成时移位寄存器个数。对实测电流波形做循环自相关运算，考虑到实际电流的噪声带来理论最小值min的波动，在求min时，宜用截取合适长度取均值的方法处理。结合对电流幅值的估计值I，利用公式
求得阶数估计值N。对阶数估计值取整，得到阶数。
C、码元频率，即移位寄存器移位时钟频率，其倒数称为码元宽度，简称码宽。根据公式
并取整，求得码元频率估计值。(取整数有利于生成参考m序列，但码元频率并不一定是整数)。
D、反馈系数是确定的移位寄存器个数(阶数)下的反馈抽头分布。对于确定的阶数n，构成最长线性反馈移位寄存器序列(m序列)的反馈系数是有限的，其必须满足本原多项式。因此，在采样率已知，阶数、码元频率已经求得的情况下，可以使用当前阶数下所有可能的反馈系数，生成可能的m序列，与待测m波形进行对比，判断待测序列最可能采用的反馈系数。生成的对比序列(0，1序列)在转换成对称幅值序列时，有两种逻辑，在无法知晓待测电流波形是何种逻辑时，互相关的对比方式可以很好解决这个问题。根据互相关的性质，两个序列越相似，越正相关的位置将会出现越大的幅值，越负相关的位置，将会出现越大的负向幅值。因此，将互相关结果取绝对值，无论待测m波形与对比m序列采用何种转换逻辑，都可以认为，出现最明显峰值的即是最可能的反馈系数。
E、从发射角度而言，生成序列时必须赋予各移位寄存器以初始值，即初始状态。而对于接收者，在随机截取一段观测数据的条件下，无法推测出初始状态.为此，本实施例中引入参考m序列的概念，即接收者生成这样一个m序列：其阶数、码元频率、采样率、反馈系数均与待测m波形相同，只有初始状态由接收者自定义，以此为参考，计算待测m波形与参考m序列重合所需的平移量，定义为其相位差。以参考m序列为参考，右移为平移正参考方向。基于以上定义，将参考m序列和待测m波形做互相关运算，尖峰位置的互相关偏移量，即为所求相位。
参见图7，为一次参数识别的显示结果。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (6)
1.一种高频m序列编码发射源参数识别方法，其特征在于，包括如下步骤：
A)循环录入长度不断增加的m波形电流数据进行线性自相关运算，直至出现相邻两个明显尖峰，其间隔为m波形电流的周期，对单周期m波形电流数据做循环自相关运算，取1/5T～4/5T内数据取均值得到循环自相关函数最小值min的估计值；
B)采用最小值min的估计值，利用公式
求得阶数估计值N，对阶数估计值取整，得到阶数；
C)根据公式
并取整，求得码元频率估计值；
D)根据已知的采样率，以及求取的阶数和码元频率，使用当前阶数下所有可能的反馈系数，生成可能的m序列，与待测m波形进行对比，判断待测序列最可能采用的反馈系数；
E)计算待测m波形与参考m序列重合所需的平移量，定义为其相位差，以参考m序列为参考，右移为平移正参考方向，将参考m序列和待测m波形做互相关运算，尖峰位置的互相关函数偏移量，为所求相位；
F)上述步骤所求得包括周期、阶数、码元频率、反馈系数和相位等发射参数，对照发射机理论参数，综合评价其合理性，得出对航空电磁探测发射系统状态的基本判断。
2.按照权利要求1所述的方法，其特征在于，所述步骤A具体包括：
A1)获取N点发射机电流波形发射时间内)双极性电流数据，并取平均值，发射电流数据直接减去该均值以校正发射电流偏置，获得无偏置双极性发射电流；
A2)将无偏置双极性电流数据翻转至正向，并取均值，电流幅值在均值以下的数据点视为低幅值噪声，将其全部置零，将去除低幅值噪声的数据求得各极值点，并取均值，该均值数据被视为双极性发射电流幅值I的估计值；
A3)截取适当点数步骤A1中的无偏置双极性发射电流数据；
A4)调整阈值，使用滞回比较方法对电流数据进行二值化；
A5)二值化后数据线性自相关并归一化；
A6)小噪声数据置零，求极值，存入极值数组；
A7)判断极值数组长度是否大于1，若大于则进行下一步；
A8)取两个尖峰，其横坐标间隔为m波形周期T；
A9)按照步骤A8求取的周期，截取一个周期步骤A1中的无偏置双极性发射电流数据；
A10)对步骤A9截取的数据做循环自相关；
A11)对T/5～4T/5内数据取均值，作为循环自相关函数最小值min的估计值。
3.按照权利要求1所述的方法，其特征在于，步骤D包括:
D1)读取反馈系数表；
D2)由已知采样频率Fs，求得的码元频率，阶数n，读取n对应的有限m个反馈系数，生成m序列；
D3)生成的m序列与步骤A9中截取的一个周期无偏置双极性待测发射电流做等长化处理；
D4)生成的m序列与步骤D5中截取的待测电流做互相关处理；
D5)互相关函数取绝对值后求极值，并减去互相关函数绝对值的均值，差值为Dmm；
D6)判断阶数n对应的反馈系数是否遍历完毕；若是则进行下一步；
D7)找到差值Dmm的最大值对应的反馈系数；
D8)得出待测电流对应的反馈系数。
4.一种基于权利要求1所述方法的直升机航空电磁收录装置，其特征在于，通过直升机平台搭载，包括磁场传感器和接收机，所述磁场传感器固定于直升机下方吊舱，通过长距离信号线与机舱内接收机连接，所述接收机的控制器内运行高频m序列编码发射源参数识别软件，接收机通过控制器启动采集卡采集数据；采集卡板载时钟通过PXIe背板路由至计数卡，经计数卡分频后作为同步信号同步发射机；所述发射机内部使用电流传感器收录激发电流数据，通过线缆传输至接收机，所述接收机运行高频m序列编码发射源参数识别软件通过采集的发射机的电流数据，对发射机进行参数识别。
5.按照权利要求4所述的装置，其特征在于，所述磁场传感器包括三分量空心接收线圈、前置放大电路、保护固定结构以及传输线缆，三分量接收线圈由相互垂直的三个空心线圈组成，其中Z分量的空心线圈为圆形，平行于吊舱平面固定，X分量的空心线圈与Y分量的空心线圈组成一个正方体，其中X分量的空心线圈包括一组平行的正方形差分空心线圈，另一组正方形差分空心线圈组成Y分量的空心线圈，X分量的空心线圈与Y分量的空心线圈垂直于于吊舱平面固定，使用塑料螺丝将线圈安装于保护结构内，整体固定于吊舱中心处，使接收线圈与发射线圈构成同心圆；或固定于吊舱偏心处，使接收线圈与发射线圈不构成同心圆；信号线跟随吊舱绳索向上延伸，经过固定于绳索上的前置放大电路后，连接至直升机机舱内的接收机。
6.按照权利要求4所述的装置，其特征在于，发射机内使用电流传感器感知激发电流信号与接收机接收的电磁响应信号一同经采集卡转换，各通道数据通过总线传递至控制器。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110250033.1A CN113050175B (zh) | 2021-03-08 | 2021-03-08 | 直升机航空电磁收录装置及发射源参数识别方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110250033.1A CN113050175B (zh) | 2021-03-08 | 2021-03-08 | 直升机航空电磁收录装置及发射源参数识别方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113050175A true CN113050175A (zh) | 2021-06-29 |
| CN113050175B CN113050175B (zh) | 2022-05-13 |
Family
ID=76510618
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110250033.1A Expired - Fee Related CN113050175B (zh) | 2021-03-08 | 2021-03-08 | 直升机航空电磁收录装置及发射源参数识别方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113050175B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113917549A (zh) * | 2021-09-08 | 2022-01-11 | 成都理工大学 | 基于光纤传感技术的航空电磁数据采集系统及采集方法 |
Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4629990A (en) * | 1982-12-17 | 1986-12-16 | Shell Internationale Research Maatschappij, B.V. | Method and apparatus for correcting the relative motion of a transmitter and a receiver in airborne electromagnetic prospecting |
| US20100237870A1 (en) * | 2009-03-17 | 2010-09-23 | Geo Equipment Manufacturing Limited | Geophysical Prospecting Using Electric And Magnetic Components Of Natural Electromagnetic Fields |
| US20110037639A1 (en) * | 2008-01-04 | 2011-02-17 | Mario Manuel Duran Toro | System and method for detecting, locating and identifying objects located above the ground and below the ground in a pre-referenced area of interest |
| CN102721982A (zh) * | 2012-06-26 | 2012-10-10 | 中国科学院电工研究所 | 基于扩频编码技术的地面电磁勘探方法及其探测系统 |
| CN103630941A (zh) * | 2013-01-30 | 2014-03-12 | 中国科学院电子学研究所 | 长线源伪随机编码发射和阵列式接收瞬变电磁系统及方法 |
| CN104865608A (zh) * | 2015-05-22 | 2015-08-26 | 吉林大学 | 时间域航空电磁法运动噪声检测装置及抑制方法 |
| CN105607131A (zh) * | 2016-01-11 | 2016-05-25 | 甘肃省有色地质调查院 | 编码源电磁测深法获取大地全区视电阻率谱的方法及装置 |
| US20170068014A1 (en) * | 2014-02-28 | 2017-03-09 | Cgg Services Sa | Systems and methods for a composite magnetic field sensor for airborne geophysical surveys |
| CN108919366A (zh) * | 2018-05-15 | 2018-11-30 | 吉林大学 | 一种直升机磁共振与瞬变电磁联合探测装置及探测方法 |
| CN109597134A (zh) * | 2019-01-28 | 2019-04-09 | 吉林大学 | 基于绝热脉冲激发源的核磁共振地下水探测装置及其方法 |
| CN109946744A (zh) * | 2019-04-01 | 2019-06-28 | 吉林大学 | 一种基于可控源补偿的瞬变电磁勘探系统及方法 |
| CN110471117A (zh) * | 2019-09-26 | 2019-11-19 | 国科（重庆）仪器有限公司 | 一种航空电磁探测系统及方法 |
| CN110515129A (zh) * | 2019-08-08 | 2019-11-29 | 吉林大学 | 城市地下空间多线圈拖曳式偶极电磁探测装置及方法 |
| CN110703343A (zh) * | 2019-10-28 | 2020-01-17 | 吉林大学 | 基于pwm调控技术的宽配谐磁共振探测装置及探测方法 |
| CN111965717A (zh) * | 2020-08-12 | 2020-11-20 | 彩虹无人机科技有限公司 | 一种固定翼无人机时间域航空电磁勘查系统 |
-
2021
- 2021-03-08 CN CN202110250033.1A patent/CN113050175B/zh not_active Expired - Fee Related
Patent Citations (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4629990A (en) * | 1982-12-17 | 1986-12-16 | Shell Internationale Research Maatschappij, B.V. | Method and apparatus for correcting the relative motion of a transmitter and a receiver in airborne electromagnetic prospecting |
| US20110037639A1 (en) * | 2008-01-04 | 2011-02-17 | Mario Manuel Duran Toro | System and method for detecting, locating and identifying objects located above the ground and below the ground in a pre-referenced area of interest |
| US20100237870A1 (en) * | 2009-03-17 | 2010-09-23 | Geo Equipment Manufacturing Limited | Geophysical Prospecting Using Electric And Magnetic Components Of Natural Electromagnetic Fields |
| CN102721982A (zh) * | 2012-06-26 | 2012-10-10 | 中国科学院电工研究所 | 基于扩频编码技术的地面电磁勘探方法及其探测系统 |
| CN103630941A (zh) * | 2013-01-30 | 2014-03-12 | 中国科学院电子学研究所 | 长线源伪随机编码发射和阵列式接收瞬变电磁系统及方法 |
| US20170068014A1 (en) * | 2014-02-28 | 2017-03-09 | Cgg Services Sa | Systems and methods for a composite magnetic field sensor for airborne geophysical surveys |
| CN104865608A (zh) * | 2015-05-22 | 2015-08-26 | 吉林大学 | 时间域航空电磁法运动噪声检测装置及抑制方法 |
| CN105607131A (zh) * | 2016-01-11 | 2016-05-25 | 甘肃省有色地质调查院 | 编码源电磁测深法获取大地全区视电阻率谱的方法及装置 |
| CN108919366A (zh) * | 2018-05-15 | 2018-11-30 | 吉林大学 | 一种直升机磁共振与瞬变电磁联合探测装置及探测方法 |
| CN109597134A (zh) * | 2019-01-28 | 2019-04-09 | 吉林大学 | 基于绝热脉冲激发源的核磁共振地下水探测装置及其方法 |
| CN109946744A (zh) * | 2019-04-01 | 2019-06-28 | 吉林大学 | 一种基于可控源补偿的瞬变电磁勘探系统及方法 |
| CN110515129A (zh) * | 2019-08-08 | 2019-11-29 | 吉林大学 | 城市地下空间多线圈拖曳式偶极电磁探测装置及方法 |
| CN110471117A (zh) * | 2019-09-26 | 2019-11-19 | 国科（重庆）仪器有限公司 | 一种航空电磁探测系统及方法 |
| CN110703343A (zh) * | 2019-10-28 | 2020-01-17 | 吉林大学 | 基于pwm调控技术的宽配谐磁共振探测装置及探测方法 |
| CN111965717A (zh) * | 2020-08-12 | 2020-11-20 | 彩虹无人机科技有限公司 | 一种固定翼无人机时间域航空电磁勘查系统 |
Non-Patent Citations (2)
| Title |
|---|
| 嵇艳鞠 等: ""基于异常线圈的时间域AEM系统测试和标定方法研究"", 《地球物理学报》 * |
| 王世隆 等: ""同心补偿式直升机时间域航空电磁法吊舱校准装置研究"", 《地球物理学报》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113917549A (zh) * | 2021-09-08 | 2022-01-11 | 成都理工大学 | 基于光纤传感技术的航空电磁数据采集系统及采集方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113050175B (zh) | 2022-05-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10746894B2 (en) | Magnetotelluric measurement system | |
| CN102053280B (zh) | 带有参考线圈的核磁共振地下水探测系统及探测方法 | |
| CN102565862B (zh) | 一种瞬变电磁响应信号梯度测量方法及观测装置 | |
| CN111257951B (zh) | 地空短偏移电磁探测系统、多源发射信号分离方法 | |
| CN103823244B (zh) | 磁共振三分量消噪装置及消噪方法 | |
| US7977940B2 (en) | Low temperature squid transient electromagnetic receiver system | |
| CN203502602U (zh) | 消除工频谐波干扰的核磁共振探测装置 | |
| CN113050175B (zh) | 直升机航空电磁收录装置及发射源参数识别方法 | |
| CN109001997B (zh) | 一种防干扰的环境信息采集系统 | |
| CN110501665B (zh) | 一种磁场辐射敏感度测试设备的比对装置 | |
| Adami et al. | HPM detection system for mobile and stationary use | |
| CN207908433U (zh) | 瞬变电磁探伤仪噪声抑制装置 | |
| CN109521477A (zh) | 一种基于squid的岸基水下磁异常探测装置及其探测方法 | |
| CN113009290A (zh) | 一种特高压gis设备局部放电种类识别和定位方法及系统 | |
| CN202870310U (zh) | 一种适用于矿用本安型瞬变电磁仪的信号采集电路 | |
| CN110068871B (zh) | 一种车载时域电磁感应-极化效应的微型测量系统及方法 | |
| CN105629315B (zh) | 主动场补偿式数字超低频电磁传感器 | |
| CN105572482A (zh) | 一种密闭空间电场测量设备 | |
| CN206757026U (zh) | 一种电子式电能表电快速瞬变脉冲群抗扰度测试系统 | |
| CN104730589A (zh) | 一种航磁测量数据采集系统和装置 | |
| Macnae et al. | Joint sensing of B and dB/dt responses | |
| CN110441716A (zh) | 低频磁场测量装置 | |
| CN109857099A (zh) | 一种磁循迹传感器以及应用该传感器的传感方法 | |
| CN112532245B (zh) | 一种基于宽带模数转换器件的多通道磁信号采集模块及方法 | |
| CN107024724A (zh) | 反卷积快速优化直升机航空电磁探测早期数据的方法 |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220513 |