一种基于参数估计的铅锌矿吸水井淤泥厚度测量方法技术

本发明专利技术公开一种基于参数估计的铅锌矿吸水井淤泥厚度测量方法，包括以下步骤：首先使用信号收集装置收集在吸水井中的回波信号，对收集到的信号降噪，提取特征，建立分类器进行分类，然后使用不同的数学模型对相应的信号进行参数估计，估计出超声波传播时间参数，使用该时间参数结合超声波在水中的传播速度计算出井底淤泥厚度。本发明专利技术不同于传统超声波测量中使用的阈值法，可以有效的避免回波中的干扰项，使测量结果更加精确。

A method based on parameter estimation for measuring thickness of silt in suction well of lead zinc mine

The invention discloses a method for measuring the silt thickness of a lead-zinc mine suction well based on parameter estimation, which comprises the following steps: firstly, the echo signals collected in the suction well are collected by the signal collection device, then the collected signals are de-noised, the features are extracted, and the classifier is established to classify, and then the corresponding mathematical models are used to classify the signals. The signal parameters are estimated to estimate the ultrasonic propagation time parameters, which are used to calculate the silt thickness at the bottom of the well. The invention is different from the threshold method used in traditional ultrasonic measurement, and can effectively avoid the interference term in the echo, so as to make the measurement result more accurate.

【技术实现步骤摘要】
一种基于参数估计的铅锌矿吸水井淤泥厚度测量方法
本专利技术属于超声波检测领域，具体涉及一种铅锌矿井下吸水井淤泥厚度测量方法。
技术介绍
矿井排水安全是影响矿区安全生产的重要因素，井下排水对井下工作人员和设备安全非常重要，是井下作业的最基本的条件。由于井下作业开采等原因，井下的水都比较浑浊，水在吸水井聚集后，容易在吸水井里产生淤泥。目前吸水井水位的测量装置是超声波液位计，通过测量探头到液位表面间的距离来给出液位高度，这种测量方式忽略了井底沉淀的淤泥厚度，随着水聚集时间越长，误差越来越大。吸水井中的水颜色较深，无法用眼估计出沉淀淤泥的泥位高低，当淤泥高度达到水泵吸水位时，导致水泵抽不到水或者水泵抽到淤泥，造成水泵的损坏。随着计算机自动化技术的进一步发展，水泵房等处无人值守是一个发展趋势，如果水位测量不准，或者是不能准确知道淤泥的厚度，可能引起严重的后果。因此获得淤泥的准确厚度，是实现井下水泵房等处无人值守可靠保障。测量吸水井下淤泥厚度通常有人工插杆法，光电法，高度可调的浓度计测量装置等。人工插杆法不能在线的实时获取淤泥厚度信息，测量繁琐，限制条件较多，在矿井下极为不便。光电法测量属于接触式测量，准确性较高，但安装复杂，设备维护麻烦。高度可调的浓度计利用浓度计的高度变化，测量出吸水井中各个高度的浓度，当在某一高度处浓度变化大时说明该处为淤泥界面，该高度便为淤泥厚度，该装置成本较高且在测量时需要不断移动浓度计，在矿下较易损坏。超声波测距基于超声波信号遇到被测物时产生反射信号，可以根据回波信号估计出超声波探头与被测物之间的距离，当前大部分超声波测距使用阈值法估计超声波传播时间，在铅锌矿下吸水井中，存在重金属粒子和其它矿下杂质干扰，使用传统阈值法难以得出准确信息。因此设计一种基于高斯回波信号模型的估计方法来估计超声波传播时间可以极大的提高测量精度，而不同的回波信号对应的回波数学模型不一，将不同模型下的回波信号分类出来是进行准确测量的前提。取得回波信号的包络可以简化回波数学模型而不影响测量精度，对于不同的回波信号采用不同的数学模型，在超声波测距中最关键的参数是超声波传播时间。
技术实现思路
本专利技术的目的提供一种使用超声波测量矿下吸水井淤泥厚度的方法，首先对信号进行去噪，分类，然后使用不同的数学模型对相应的回波信号进行参数估计，估计出超声波传播时间，根据超声波在水中的传输速度计算出相应的淤泥厚度。本专利技术的技术方案具体步骤如下：S1：将超声波换能器探头安装在距离吸水井底部固定高度的位置，探头距离池底高度为h，使用超声波换能器和信号采集系统，获取超声波经吸水井底部泥水界面反射后的回波信号，在不同工况下收集p组回波信号，p满足50≤p≤500；S2：对于超声波回波信号进行9层小波分解，获得各层小波系数，使用Hesusure阈值选取方法进行降噪处理，使用软阈值函数进行系数处理，使用处理后的小波系数进行信号重构，提高信号信噪比；S3：对于去噪后的信号进行特征提取，然后建立随机森林分类器进行分类，将信号分为三种类型；S4：对于不同类别的回波信号，采用不同的高斯回波数学模型，分别为单回波模型、双回波模型和三回波模型；S5：对信号做希尔伯特变换，取绝对值得到信号包络；S6：对于信号包络，使用简化后的高斯回波模型进行参数估计，高斯单回波模型为其中β为幅度，a为带宽因子，τ为传播时间，p(t)为噪声，t为时间变量，双回波模型为其中β21、β22为两个回波分别的幅度，a21、a22为两个回波分别的带宽因子，τ21、τ22为两个回波分别的传播时间，三回波模型为其中β31、β32、β33为三个回波分别的幅度，a31、a32、a33为三个回波分别的带宽因子，τ31、τ32、τ33为三个回波分别的传播时间；S7：使用人工蜂群算法对包络信号进行参数估计，分别估计出单回波飞行时间τ，双回波飞行时间τ21，τ22，三回波飞行时间τ31，τ32，τ33；S8：对于单回波信号，计算超声波传输距离其中d1为单回波超声波传输距离，v为超声波在水中的传输速度，设hx为淤泥厚度，则hx＝h-d1，对于双回波信号，超声波传输距离为或其中d21、d22分别为第一回波和第二回波的传输距离，淤泥厚度为hx＝h-d21或hx＝h-d22，对于三回波信号，舍弃第三回波对应的参数，超声波传输距离为或淤泥厚度为hx＝h-d31或hx＝h-d32，其中d31、d32分别为第一回波、第二回波的传输距离。所述S1中将超声波换能器探头安装在水下固定高度位置，所选探头发射角为6°，探头距离池底高度为h，历史最高泥位为hm，h满足h-hm≥30cm，以探头为中心点，探头面为水平面，向下至池底，半径为h*tan3°+h′的圆柱体空间内为测量空间，测量空间内应排除可能产生干扰回波的固体障碍物，h′为冗余量，满足所述S2中使用小波函数对原始信号进行9层小波分解，获得各层小波系数，使用Hesusure阈值选取规则，结合固定阈值规则和无偏似然估计自适应阈值规则，其具体步骤如下：S21：首先计算出固定阈值规则规定的阈值N为信号长度，σ为噪声的标准差；S22：然后计算无偏似然估计自适应阈值Rigrsure阈值，小波系数平方向量X＝[x1，x2，…，xn]，其中x1≤x2≤…≤xn，n为该层小波系数的个数，设一风险向量为：S＝[s1，s2，…，si，…，sn]，风险向量中各元素为：从风险向量中找出最小值smin作为风险值，找出最小风险值下标系数对应的xmin，则自适应阈值为：S23：结合上述两种阈值选取规则，设K为小波分解后某层小波系数平方和，定义参数参数混合型阈值λh取值规则为当确定阈值后，选用软阈值函数对小波系数进行处理，软阈值函数为：其中xn为经软阈值函数处理后的新小波系数；S24：对小波系数进行阈值处理后，进行信号重构，获得降噪后的信号。所述S3将铅锌矿吸水井泥水界面回波信号分为三种类别，分别对应三种不同的泥水界面情形，分别为单回波信号、双回波信号及三回波信号，能够更好的估计泥水界面位置和沉淀物组成成分。所述S5对信号做希尔伯特变换其中s(t)为原始信号，然后对经希尔伯特变换后的信号取绝对值得到g(t)为信号上半部分包络，在超声波测距中，最关键的参数为超声波在水中传播的时间，即信号幅度最高点所对应的时间参数，对原始信号求包络后，信号幅度最高点对应的时间参数不变，而相应的数学模型极大的简化了，更有利于后续的参数估计。所述S6使用简化的回波模型进行参数估计，单回波模型为双回波模型为三回波模型为使用简化后的回波模型可以减少待估计的参数个数，未简化单回波模型中需要估计五个参数，双回波参数为十个，三回波参数为十五个，简化后单回波参数为三个，双回波为六个，三回波为九个，极大的减少了估计参数时的计算量。所述S7使用人工蜂群算法对三种模型的信号分别进行参数估计，以单回波模型为例，可分为以下步骤：S71：单回波模型参数维度为3，假设li＝[li1，li2，li3]表示第i个蜜源，蜜源生成公式为：lij＝lbj+rand(0，1)(luj-lbj)，其中lbj为第j个参数的下界，luj为第j各参数的上界，luj为新生成的蜜源，rand(0，1)为0到1之间的随机数，利用蜜源生成公式生成初始蜜源，即对应的可行解；S本文档来自技高网...

【技术保护点】
1.一种基于参数估计的铅锌矿吸水井淤泥厚度测量方法，其特征在于包含以下步骤：S1：将超声波换能器探头安装在距离吸水井底部固定高度的位置，探头距离池底高度为h，使用超声波换能器和信号采集系统，获取超声波经吸水井底部泥水界面反射后的回波信号，在不同工况下收集p组回波信号，p满足50≤p≤500；S2：对于超声波回波信号进行9层小波分解，获得各层小波系数，使用Hesusure阈值选取方法进行降噪处理，使用软阈值函数进行系数处理，使用处理后的小波系数进行信号重构，提高信号信噪比；S3：对于降噪后的信号进行特征提取，然后建立随机森林分类器进行分类，将信号分为三种类型，分别对应三种不同的泥水界面情形，分别为单回波信号、双回波信号及三回波信号；S4：对于不同类别的回波信号，采用不同的高斯回波数学模型，分别为单回波模型、双回波模型和三回波模型；S5：对信号做希尔伯特变换，取绝对值得到信号包络；S6：对于信号包络，使用简化后的高斯回波模型进行参数估计，高斯单回波模型为

【技术特征摘要】
1.一种基于参数估计的铅锌矿吸水井淤泥厚度测量方法，其特征在于包含以下步骤：S1：将超声波换能器探头安装在距离吸水井底部固定高度的位置，探头距离池底高度为h，使用超声波换能器和信号采集系统，获取超声波经吸水井底部泥水界面反射后的回波信号，在不同工况下收集p组回波信号，p满足50≤p≤500；S2：对于超声波回波信号进行9层小波分解，获得各层小波系数，使用Hesusure阈值选取方法进行降噪处理，使用软阈值函数进行系数处理，使用处理后的小波系数进行信号重构，提高信号信噪比；S3：对于降噪后的信号进行特征提取，然后建立随机森林分类器进行分类，将信号分为三种类型，分别对应三种不同的泥水界面情形，分别为单回波信号、双回波信号及三回波信号；S4：对于不同类别的回波信号，采用不同的高斯回波数学模型，分别为单回波模型、双回波模型和三回波模型；S5：对信号做希尔伯特变换，取绝对值得到信号包络；S6：对于信号包络，使用简化后的高斯回波模型进行参数估计，高斯单回波模型为其中β为幅度，a为带宽因子，τ为传播时间，p(t)为噪声，t为时间变量，双回波模型为其中β21、β22为两个回波分别的幅度，a21、a22为两个回波分别的带宽因子，τ21、τ22为两个回波分别的传播时间，三回波模型为其中β31、β32、β33为三个回波分别的幅度，a31、a32、a33为三个回波分别的带宽因子，τ31、τ32、τ33为三个回波分别的传播时间；S7：使用人工蜂群算法对包络信号进行参数估计，分别估计出单回波飞行时间τ，双回波飞行时间τ21，τ22，三回波飞行时间τ31，τ32，τ33；S8：对于单回波信号，计算超声波传输距离其中d1为单回波超声波传输距离，v为超声波在水中的传输速度，设hk为淤泥厚度，则hx＝h-d1，对于双回波信号，超声波传输距离为或其中d21、d22分别为第一回波和第二回波的传输距离，淤泥厚度为hx＝h-d21或hx＝h-d22，对于三回波信号，舍弃第三回波对应的参数，超声波传输距离为或淤泥厚度为hx＝h-d31或hx＝h-d32，其中d31、d32分别为第一回波、第二回波的传输距离。2.根据权利要求1所述的基于参数估计的铅锌矿吸水井淤泥厚度测量方法，其特征在于：S1将超声波换能器探头安装在水下固定高度位置，所选探头发射角为6°，探头距离池底高度为h，历史最高泥位为hm，h满足h-hm≥30cm，以探头为中心点，探头面为水平面，向下至池底，半径为h*tan3°+h′的圆柱体空间内为测量空间，h′为冗余量，满足3.根据权利要求...

【专利技术属性】
技术研发人员：唐朝晖，史伟东，王紫勋，曾思迪，唐励雍，
申请(专利权)人：中南大学，
类型：发明
国别省市：湖南,43