[发明专利]一种提高测量地震检波器精度的方法

权利要求书

1.一种提高测量地震检波器精度的方法，其特征在于采用如下步骤：

a.确定检波器线圈质量m、选配激励电流I₀；

b.选定初始值：响应信号初始幅值A₀＝50～300mV，

              检波器阻尼B_l＝0.6，

              检波器自然角频率w₀＝2πf，其中f＝24Hz，

              时间t；

c.根据公式：令b₀＝A₀，b₁＝-B_lw₀，得到：

d.对式(1)中t选择8～30个不同的值；

e.用修正的阻尼系数最小二乘法(Modified Damped Least Squares，MDLS)进行数据处理，得到b₀，b₁，b₂的真值。

f.通过(3)求出f；

g.通过(4)求出Bt；

h.由于即通过求出S。

2.根据权利要求1所述的提高测量地震检波器精度的方法，其特征在于步骤步骤d选择的8～30个不同的测量数据值分布均匀并在式(1)前3/4周期内。

3.根据权利要求1所述的提高测量地震检波器精度的方法，其特征在于步骤e的最小二乘法可采用的数学程式为：

(1)设置常系数：ε₁＝0.5e-6、ε₂＝0.5e-14、ε₃＝0.5e-3；阻尼因子λ＝1e-2，修正因子β＝1e-14，参比量g_e＝0.3e-5，待求参数个数n＝(n＝3)，样点数m＝(8≤m≤30)；矩阵Γ＝{0.1，-0.1，0.01，-0.01}，并设待求矩阵初值b＝{500，-30，50}。

(2)将b代入响应方程，结合测量数据求响应方程的残差向量r和最小二乘误差Q。

(3)计算矩阵，e_i为第i个样点电压幅变，b_g为b阵中第I个元素。

(4)计算矩阵g＝A^Tr，H＝A^TA。给H阵对角线上乘一常数(1+λ)，令其等于B。

(5)判断阵B的奇异性；B非奇异，转入(8)；

(6)挑选矩阵g最大元素，确定相应的l；

(7)若|g_l|≤ε_l，转入(13)，否则构造矩阵P_lxn，其中第1个元素P_l＝-g_l，其余为0，转入线性搜索部分；

(8)令矩阵P＝-B^-1g。若g^TP≥0转入(6)，否则令α＝1；

(9)取矩阵，计算相应的残差向量和最小二乘残差

(10)判断，成立取，不成立取，转入(9)，

(11)令矩阵，若中至少有一个能满足，转入(13)，不满足；

(12)取转入(3)；

(13)每次将矩阵b中任一元素b_i(1≤i≤n)修正成(1+γ_κ)b_j，γ_κ为Γ阵中的第κ个元素，保持其它元素不变，可得到4n个矩阵b的修正矩阵，确定出最小二乘残乘最小的

(14)若转入(12)，否则输出b₀，b₁，b₂真值。

4.根据权利要求1、3所述的提高测量地震检波器精度的方法，其特征在于数学程式中的线性搜索程序如下：

(51)取中间量(52)令，计算相应的(53)若令h＝2h，转入(52)，否则；

(54)取，并判断成立与否；

(55)若成立，令，并计算相应的。若，取h＝α，转入(54)，若，取λ＝4λ，转入(11)；

(56)若不成立，取，并计算相应的，若，取λ＝4λ，转入(14)；否则，转入(54)。