主权项 |
一种粮堆中异物的识别方法,其特征在于,步骤如下:a)用探地雷达探测,对粮堆中异物进行定位;b)根据所测粮堆深度与异物深度,将粮堆沿垂直方向从上到下依次分成至少两个虚拟层;基于反演算法和反射折射定律,得出逐层反演公式:第一层介电常数: <mrow> <msqrt> <msub> <mi>ϵ</mi> <mn>1</mn> </msub> </msqrt> <mo>=</mo> <msub> <mi>k</mi> <mn>1</mn> </msub> <mo>·</mo> <mfrac> <mrow> <mn>1</mn> <mo>+</mo> <msub> <mi>f</mi> <mn>1</mn> </msub> <mfrac> <msub> <mi>A</mi> <mn>0</mn> </msub> <msub> <mi>A</mi> <mi>m</mi> </msub> </mfrac> </mrow> <mrow> <mn>1</mn> <mo>-</mo> <msub> <mi>f</mi> <mn>1</mn> </msub> <mfrac> <msub> <mi>A</mi> <mn>0</mn> </msub> <msub> <mi>A</mi> <mi>m</mi> </msub> </mfrac> </mrow> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>第二层及之后层介电常数: <mrow> <msqrt> <msub> <mi>ϵ</mi> <mi>i</mi> </msub> </msqrt> <mo>=</mo> <msub> <mi>k</mi> <mi>i</mi> </msub> <msqrt> <msub> <mi>ϵ</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </msqrt> <mo>·</mo> <mfrac> <mrow> <mfrac> <mrow> <mn>4</mn> <msqrt> <msub> <mi>ϵ</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>2</mn> </mrow> </msub> <msub> <mi>ϵ</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </msqrt> </mrow> <mrow> <msub> <mi>ϵ</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>2</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>ϵ</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </mrow> </mfrac> <mo>-</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mfrac> <msub> <mi>A</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <msub> <mi>A</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>2</mn> </mrow> </msub> </mfrac> </mrow> <mrow> <mfrac> <mrow> <mn>4</mn> <msqrt> <msub> <mi>ϵ</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>2</mn> </mrow> </msub> <msub> <mi>ϵ</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> </msqrt> </mrow> <mrow> <msub> <mi>ϵ</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>2</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>ϵ</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>2</mn> </mrow> </msub> </mrow> </mfrac> <mo>+</mo> <msub> <mi>f</mi> <mi>i</mi> </msub> <mfrac> <msub> <mi>A</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <msub> <mi>A</mi> <mrow> <mi>i</mi> <mo>-</mo> <mn>2</mn> </mrow> </msub> </mfrac> </mrow> </mfrac> <mo>,</mo> <mi>i</mi> <mo>=</mo> <mn>2,3</mn> <mo>,</mo> <mo>.</mo> <mo>.</mo> <mo>.</mo> <mi>N</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>12</mn> <mo>)</mo> </mrow> </mrow>公式(12)中,εi是各虚拟层的介电常数,i表示从上到下层数的序号,ε0是空气介电常数,Ai是各虚拟层的反射回波幅度,ki是各虚拟层的线性误差校准因子,fi是各虚拟层幅度衰减因子;公式(11)中,Am是发射电磁波的幅度,A0是第一虚拟层反射回波幅度,k1是第一虚拟层的线性误差校准因子,f1是第一虚拟层幅度衰减因子;c)将探地雷达天线放置在粮堆表面,测量探地雷达天线对应位置垂直扫描线上各虚拟层的反射回波幅度Ai;d)采用测井法进行两次测量确定公式(11)、(12)中的参数f1,k1,fi,ki;采用一个发射天线与两个接收天线,第一次测量时,将发射天线和接收天线分别插入第一虚拟层,根据测量的电磁波幅度衰减和相移,并结合经几何衰减校正的电磁波传导模型,确定第一虚拟层介电常数ε1,并根据探地雷达测得的发射电磁波的幅度Am,空气介电常数ε0,第一层反射回波幅度A0,结合公式(11),设f1=1,确定k1;第二次测量时,将发射天线和接收天线分别插入第二虚拟层,根据测量的电磁波幅度衰减和相移,并结合经几何衰减校正的电磁波传导模型,确定第二虚拟层介电常数ε2,考虑到是同一介质,各层衰减相同,可把f1作为各层幅度衰减因子fk,即f1=fk,k=2,…,n;并结合公式(12),设ki=k1,确定幅度衰减校准因子fi;e)利用公式(12)迭代计算探地雷达天线对应位置垂直扫描线上各虚拟层的介电常数;f)根据各虚拟层的介电常数,判断异物的性质和种类。 |