[发明专利]一种确定炮检点组合的方法

摘要

本发明涉及地球物理勘探地震数据采集中的确定炮检点组合的方法，采集或收集以往的地震资料，获取确定炮检组合的地球物理参数，确定炮点组内距、检波点组内距的绕射波场最大时间变化率，确定炮检组合对绕射波的最大压制特性曲线，根据炮检组合对绕射波的最大振幅压制特性曲线确定炮检组合。本发明可以通过地下绕射点的深度和上伏地层的均方根速度，以及偏移距矢量和炮检距矢量，结合需要保护绕射波的最高频率（高截频），确定炮点组合基距和检波点组合基距。

主权项

一种确定炮检点组合的方法，特点是通过以下具体步骤实现：1)采集或收集地震资料，获取确定炮检组合的地球物理参数；2)确定炮点组内距、检波点组内距的绕射波场最大时间变化率；(1)通过绕射点与炮点、检波点的关系建立以下炮点组合中心到绕射点再到检波点组合中心的时距方程：$t=\frac{1}{v}(\sqrt{x_s^2+y_s^2+z^2}+\sqrt{x_r^2+y_r^2+z^2})---\left(2\right)$式中x_s、y_s表示炮点组合中心的横坐标和纵坐标，x_r、y_r表示检波点组合中心的横坐标和纵坐标，z为绕射点的深度，v表示绕射点上伏地层的均方根速度；(2)对以上时距方程求导取极值得到检波点组内距的最大时间变化率和炮点组内距的最大时间变化率：${\mathrm{dt}}_{r,max}=\frac{(r_0+h)}{\sqrt{{(r_0+h)}^2+z^2}}\frac{\mathrm{\δ}r}{v}=br\frac{\mathrm{\δ}r}{v}---\left(3\right)$${\mathrm{dt}}_{s\·\max }=\frac{(r_0-h)}{\sqrt{{(r_0-h)}^2+z^2}}\frac{\mathrm{\δ}s}{v}=bs\frac{\mathrm{\δ}s}{v}---\left(4\right)$式中，dt_r,max为检波点组内距的最大时间变化率，dt_s,max为炮点组内距的最大时间变化率，h为半炮检距，r₀为绕射点在地面投影到炮检距中点的连线距离，δs表示炮点组合內距，δr表示检波点组合內距，z为绕射点的深度，br为绕射点至检波点的射线倾角的正弦，bs为炮点至绕射点的射线倾角的正弦；(3)将步骤1)得到的地球物理参数代入公式(3)和公式(4)计算，绘制出炮点组内距和检波点组内距的最大时间变化率随半炮检距变化的关系曲线；3)确定炮检组合对绕射波的最大压制特性曲线；(1)将检波点组内距的最大时间变化率公式和炮点组内距的最大时间变化率公式代入炮检组合的振幅压制特性公式，得到炮检组合对绕射波的最大压制特性曲线公式：炮检组合的振幅压制特性公式是：$\left|H\left(f\right)\right|=\frac{\sin \left({\mathrm{\πfNdt}}_s\right)}{N\sin \left({\mathrm{\πfdt}}_s\right)}\frac{\sin \left({\mathrm{\πfMdt}}_r\right)}{M\sin \left({\mathrm{\πfdt}}_r\right)}---\left(5\right)$式中：H(f)表示炮检组合振幅压制特性值；N为炮点组合个数；M为检波点组合个数；dt_r为检波点组内距的时间变化率，dt_s为炮点组内距的时间变化率；炮检组合对绕射波的最大压制特性曲线公式是：$\begin{array}{c}{\left|H\left(f\right)\right|}_{\max }=\frac{\sin \left({\mathrm{\πfNdt}}_{s\·\max }\right)}{N\sin \left({\mathrm{\πfdt}}_{s\·\max }\right)}\frac{\sin \left({\mathrm{\πfMdt}}_{r\·\max }\right)}{M\sin \left({\mathrm{\πfdt}}_{r\·\max }\right)}\\ =\frac{\sin (\mathrm{\π}f\·bs\·Ls/v)}{N\sin (\mathrm{\π}f\·bs\·Ls/(N\×v\left)\right)}\frac{\sin (\mathrm{\π}f\·br\·Lr/v)}{M\sin (\mathrm{\π}f\·br\·Lr/(M\×v\left)\right)}\end{array}---\left(6\right)$式中：|H(f)|_max表示炮检组合的最大振幅压制特性值；N为炮点组合个数；M为检波点组合个数；Ls表示炮点组合基距，Ls＝N×δs；Lr表示检波点组合基距，Lr＝M×δr；v表示绕射点上伏地层的均方根速度，δs表示炮点组合內距，δr表示检波点组合內距，br为绕射点至检波点的射线倾角的正弦，bs为炮点至绕射点的射线倾角的正弦；(2)根据以上公式计算绘制炮检组合对绕射波的最大压制特性曲线；4)根据炮检组合对绕射波的最大振幅压制特性曲线确定炮检组合。