主权项 |
一种针对特定地区的指向性聚束检测方法,其特征在于:包括以下步骤:S1、选取一个合适的指向性聚束监测台阵;S2、计算特定地区相对步骤S1中确定的监测台阵的慢度向量,具体步骤如下:S2.1、读取特定地区历史核爆炸地震事件在该监测台阵各子台上的信号记录,去掉数据的线性趋势项;S2.2、根据公式(5)对步骤S2.1得到的数据进行频率‑波数分析,得到在不同慢度条件下的信号谱P:<maths num="0001"><math><![CDATA[<mrow><mi>P</mi><mrow><mo>(</mo><msub><mi>S</mi><mi>n</mi></msub><mo>,</mo><msub><mi>S</mi><mi>e</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>Σ</mi><mrow><mi>f</mi><mo>=</mo><msub><mi>f</mi><mn>1</mn></msub></mrow><msub><mi>f</mi><mn>2</mn></msub></munderover><mo>|</mo><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><msub><mi>F</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>·</mo><msup><mi>e</mi><mrow><mn>2</mn><mi>π</mi><msqrt><mrow><mo>-</mo><mn>1</mn></mrow></msqrt><mi>f</mi><mrow><mo>(</mo><msub><mi>S</mi><mi>n</mi></msub><mo>·</mo><mi>d</mi><mi>n</mi><mi>o</mi><mi>r</mi><mi>t</mi><mi>h</mi><mo>+</mo><msub><mi>S</mi><mi>e</mi></msub><mo>·</mo><mi>d</mi><mi>e</mi><mi>a</mi><mi>s</mi><mi>t</mi><mo>)</mo></mrow></mrow></msup><mo>|</mo></mrow><mrow><mi>J</mi><mo>·</mo><munderover><mi>Σ</mi><mrow><mi>f</mi><mo>=</mo><msub><mi>f</mi><mn>1</mn></msub></mrow><msub><mi>f</mi><mn>2</mn></msub></munderover><mo>{</mo><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><msub><mi>F</mi><mi>k</mi></msub><msup><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>}</mo></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001109383900000011.GIF" wi="1446" he="391" /></maths>式中S<sub>n</sub>和S<sub>e</sub>是慢度向量在南北向和东西向上的分量,F(f)是信号与自身作卷积后的序列,f<sub>1</sub>和f<sub>2</sub>分别是滤波频带的起止频率,J为该监测台阵的子台个数,dnorth和deast是各子台相对于中心子台在南北向和东西向的偏移量;S2.3、记录信号谱P在最大值时对应的慢度向量的分量S<sub>nm</sub>、S<sub>em</sub>;S3、对特定地区核爆炸地震事件进行检测,具体步骤如下:S3.1、获取监测台阵各子台的实时监测数据;S3.2、结合步骤S3.1的实时监测数据,将步骤S2.3中得到的S<sub>nm</sub>、S<sub>em</sub>代入到公式(6)中计算台阵信号检测的F值:<maths num="0002"><math><![CDATA[<mrow><mi>F</mi><mrow><mo>(</mo><msub><mi>s</mi><mrow><mi>n</mi><mi>m</mi></mrow></msub><mo>,</mo><msub><mi>s</mi><mrow><mi>e</mi><mi>m</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mo>[</mo><mfrac><mrow><mi>J</mi><mo>-</mo><mn>1</mn></mrow><mi>J</mi></mfrac><mo>]</mo><mo>·</mo><mfrac><mrow><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><msub><mi>n</mi><mn>0</mn></msub></mrow><mrow><msub><mi>n</mi><mn>0</mn></msub><mo>+</mo><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></munderover><msup><mrow><mo>|</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>+</mo><msub><mi>lag</mi><mi>i</mi></msub><mo>(</mo><mrow><msub><mi>s</mi><mrow><mi>n</mi><mi>m</mi></mrow></msub><mo>,</mo><msub><mi>s</mi><mrow><mi>e</mi><mi>m</mi></mrow></msub></mrow><mo>)</mo><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow><mrow><munderover><mi>Σ</mi><mrow><mi>n</mi><mo>=</mo><msub><mi>n</mi><mn>0</mn></msub></mrow><mrow><msub><mi>n</mi><mn>0</mn></msub><mo>+</mo><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></munderover><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><msup><mrow><mo>|</mo><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>+</mo><msub><mi>lag</mi><mi>i</mi></msub><mo>(</mo><mrow><msub><mi>s</mi><mrow><mi>n</mi><mi>m</mi></mrow></msub><mo>,</mo><msub><mi>s</mi><mrow><mi>e</mi><mi>m</mi></mrow></msub></mrow><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>[</mo><mfrac><mn>1</mn><mi>J</mi></mfrac><mo>]</mo><mo>·</mo><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>J</mi></munderover><msub><mi>x</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>n</mi><mo>+</mo><msub><mi>lag</mi><mi>j</mi></msub><mo>(</mo><mrow><msub><mi>s</mi><mrow><mi>n</mi><mi>m</mi></mrow></msub><mo>,</mo><msub><mi>s</mi><mrow><mi>e</mi><mi>m</mi></mrow></msub></mrow><mo>)</mo><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001109383900000012.GIF" wi="1902" he="286" /></maths>式中n<sub>0</sub>为开始计算的采样点,N为时域的采样窗长,J为监测台阵的子台个数,lag<sub>i</sub>(s<sub>nm</sub>,s<sub>em</sub>)为给定慢度矢量(s<sub>nm</sub>,s<sub>em</sub>)下聚束的时间延迟,x<sub>i</sub>(n+lag<sub>i</sub>(s<sub>nm</sub>,s<sub>em</sub>))是第i个子台在采样点n处延迟后的零均值滤波数据的幅值;S3.3、将步骤S3.2中计算得到的F值代入服从公式(7)的F分布模型中得到信号的检测概率P:F~F(N<sub>1</sub>,N<sub>2</sub>,λ) (7)式中N<sub>1</sub>=2BT,N<sub>2</sub>=(J‑1)N<sub>1</sub>,λ=N<sub>1</sub>(SNR)<sup>2</sup>,其中B为有效带宽,T为有效时长,J为该监测台阵的子台个数,λ为F分布非中心因子,(SNR)<sup>2</sup>为信号与噪声的平均能量之比;S3.4、当信号检测概率P大于等于设定的检测阈值时,即可确定该特定地区发生核爆炸地震事件;当信号检测概率P小于设定的检测阈值时,重复上述步骤S3.1‑S3.3。 |