主权项 |
1.一种地震资料品质分析技术,其特征在于采用如下步骤:1)根据分析内容和已知野外班报记录选定控制点共炮集或共检波点集数据;2)根据地震数据的采样率确定最高分析频率,根据目标区块有效反射信号频率的范围选择分析频带,对分析数据的每一炮的各地震道做频率滤波处理,得到每一道数据的分频结果;3)采用以下公式(1)对同一炮同一频带的所有地震道做时间域中值滤波,得到该炮各频带的统计能量曲线结果;<math><mrow><msub><mi>A</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>T</mi><mo>,</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>MedianFilter</mi><mrow><mi>j</mi><mo>=</mo><msub><mi>j</mi><mn>1</mn></msub><mo>·</mo><mo>·</mo><mo>·</mo><msub><mi>j</mi><mn>2</mn></msub></mrow></munder><mo>[</mo><mo>|</mo><munderover><mi>Σ</mi><mrow><mi>t</mi><mo>=</mo><msub><mi>t</mi><mn>1</mn></msub></mrow><msub><mi>t</mi><mn>2</mn></msub></munderover><mrow><mo>(</mo><msub><mi>x</mi><mi>ij</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>,</mo><mi>f</mi><mo>)</mo></mrow></mrow><mo>|</mo><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math> 其中:xij(t,f)-采集的地震道i=1,2,…N炮集号;j=j1,…j2炮集内道号;t=t1…t2计算时窗内样点Ai(T,f)——时频分析的炮集统计能量,T-计算时窗中点,f-为分析频率4)将计算出的同一炮的各频带的统计能量曲线结果绘制到同一张图上,得到炮集内的时频分析图;将不同炮同一频带的统计能量曲线结果绘制到同一张图上,得到炮集间的时频分析图;5)根据地震数据的道长确定最大分析时间,根据目标区块地层情况确定分析时窗大小,对分析数据的每一炮的各地震道在选定的分析时窗内做快速富立叶变换,得到每一道数据各个分析时窗在频率域的变换结果;6)采用以下公式对同一炮同一分析时窗的所有地震道做频率域中值滤波,得到该炮的各分析时窗的统计频谱曲线;<math><mrow><msub><mi>A</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>MedianFilter</mi><mrow><mi>j</mi><mo>=</mo><msub><mi>j</mi><mn>1</mn></msub><mo>·</mo><mo>·</mo><mo>·</mo><msub><mi>j</mi><mn>2</mn></msub></mrow></munder><mo>[</mo><mo>|</mo><mi>FFT</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>ij</mi></msub><mrow><mo>(</mo><msub><mi>t</mi><mi>l</mi></msub><mo>)</mo></mrow><mo>|</mo><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math> 其中:xij(tl)——采集的地震道i=1,2,…N炮集号;j=j1,…j2炮集内道号;tl=t1…tm反射信号时窗数据Ai(f)——第i炮的频率统计振幅,f——为频率;7)将同一炮的各分析时窗的统计频谱曲线结果绘制到同一张图上,得到炮集内的频时分析图;将不同炮同一分析时窗的统计频谱曲线结果绘制到同一张图上,即得到炮集间的频时分析图;8)采用以下公式3获得地震数据的瞬时频率信噪比图,并确定信噪比临界点;<math><mrow><msub><mi>R</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>f</mi><mo>)</mo></mrow><mo>=</mo><munder><mi>MedianFilter</mi><mrow><mi>j</mi><mo>=</mo><msub><mi>j</mi><mn>1</mn></msub><mo>·</mo><mo>·</mo><mo>·</mo><msub><mi>j</mi><mn>2</mn></msub></mrow></munder><mo>[</mo><mo>|</mo><mi>FFT</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>ij</mi></msub><mrow><mo>(</mo><msubsup><mi>t</mi><mi>l</mi><mi>S</mi></msubsup><mo>)</mo></mrow><mo>|</mo><mo>)</mo></mrow><mo>]</mo><mo>/</mo><munder><mi>MedianFilter</mi><mrow><mi>j</mi><mo>=</mo><msub><mi>j</mi><mn>1</mn></msub><mo>·</mo><mo>·</mo><mo>·</mo><msub><mi>j</mi><mn>2</mn></msub></mrow></munder><mo>[</mo><mo>|</mo><mi>FFT</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>ij</mi></msub><mrow><mo>(</mo><msubsup><mi>t</mi><mi>k</mi><mi>N</mi></msubsup><mo>)</mo></mrow><mo>|</mo><mo>)</mo></mrow><mo>]</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></math> 其中:xij(tl S)——采集的地震道i=1,2,…N炮集号;j=j1,…j2炮集内道号;<math><mrow><msubsup><mi>t</mi><mi>l</mi><mi>S</mi></msubsup><mo>=</mo><msubsup><mi>t</mi><mn>1</mn><mi>S</mi></msubsup><mo>·</mo><mo>·</mo><mo>·</mo><msubsup><mi>t</mi><mi>M</mi><mi>S</mi></msubsup></mrow></math>反射信号时窗数据xij(tk N)——采集地震道(<math><mrow><msubsup><mi>t</mi><mi>k</mi><mi>N</mi></msubsup><mo>=</mo><msubsup><mi>t</mi><mn>1</mn><mi>N</mi></msubsup><mo>·</mo><mo>·</mo><mo>·</mo><msubsup><mi>t</mi><mi>P</mi><mi>N</mi></msubsup></mrow></math>初至前干扰时窗数据)Ri(f)-瞬时频率信噪比,f——为频率;计算窗根据道长取200-300毫秒,计算时窗每次向下滑动半个时窗进行计算;9)通过以上步骤4)、7)、8)得到的曲线的交汇图,直接读出能量相对差异值;大地吸收衰减值,瞬时频率信噪比值及信噪比临界点的位置。 |