主权项 |
一种无表层调查资料约束的大炮初至综合建模静校正方法,特点是通过以下具体步骤实现:1)根据野外实际地震采集资料和测量成果整理生成测线的炮点、检波点和关系文件;2)对炮点地表高程进行井深激发深度校正,采用公式Z'<sub>s</sub>=Z<sub>s</sub>‑H<sub>s</sub>计算,其中:Z'<sub>s</sub>为炮点激发面高程,Z<sub>s</sub>为炮点原始地表高程,H<sub>s</sub>为炮点激发深度;3)对全部生产炮记录拾取偏移距范围内的初至时间;4)在共中心点道集内进行折射分层,然后在检波点地表高程Z<sub>R</sub>和炮点激发面高程Z′<sub>s</sub>构成的数据体上,用互换速度分析方法计算折射速度,然后采用高斯—赛德尔方法计算延迟时,得到每一个检波点和炮点的折射速度和延迟时;5)对测线的大炮初至数据采用层析反演方法进行层析反演,得到测线近地表的速度空间模型,然后根据折射分析得到的折射速度,从速度空间模型数据中逐点提取各检波点和炮点位置处与折射速度相等时对应的高程,作为折射层界面高程;6)把提取的折射层界面高程作为折射反演的约束参数,利用延迟时和折射速度计算出折射层以上的表层平均速度,得到每一个检波点地表和炮点激发面处的表层平均速度、高速顶界面高程、高速层速度,建立整个测线的等效表层模型,采用以下公式计算表层平均速度:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>V</mi><mrow><mi>s</mi><mi>w</mi><mo>,</mo><mi>r</mi></mrow></msub><mo>=</mo><msqrt><mrow><mo>(</mo><msup><msub><mi>h</mi><mi>r</mi></msub><mn>2</mn></msup><mo>×</mo><msubsup><mi>v</mi><mrow><mn>1</mn><mo>,</mo><mi>r</mi></mrow><mn>2</mn></msubsup><mo>)</mo><mo>/</mo><mo>(</mo><msubsup><mi>d</mi><mrow><mi>t</mi><mo>,</mo><mi>r</mi></mrow><mn>2</mn></msubsup><mo>×</mo><msubsup><mi>v</mi><mrow><mn>1</mn><mo>,</mo><mi>r</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msup><msub><mi>h</mi><mi>r</mi></msub><mn>2</mn></msup><mo>)</mo></mrow></msqrt></mrow>]]></math><img file="FDA0000912485730000011.GIF" wi="838" he="93" /></maths>式中:V<sub>sw,r</sub>为炮点r处的表层平均速度;h<sub>r</sub>为炮点r处的表层厚度;d<sub>t,r</sub>为炮点r处的延迟时;v<sub>1,r</sub>为炮点r处的折射速度;<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>V</mi><mrow><mi>r</mi><mi>w</mi><mo>,</mo><mi>j</mi></mrow></msub><mo>=</mo><msqrt><mrow><mo>(</mo><msup><msub><mi>h</mi><mi>j</mi></msub><mn>2</mn></msup><mo>×</mo><msubsup><mi>v</mi><mrow><mn>1</mn><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup><mo>)</mo><mo>/</mo><mo>(</mo><msubsup><mi>d</mi><mrow><mi>t</mi><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup><mo>×</mo><msubsup><mi>v</mi><mrow><mn>1</mn><mo>,</mo><mi>j</mi></mrow><mn>2</mn></msubsup><mo>+</mo><msup><msub><mi>h</mi><mi>j</mi></msub><mn>2</mn></msup><mo>)</mo></mrow></msqrt></mrow>]]></math><img file="FDA0000912485730000012.GIF" wi="845" he="95" /></maths>式中:V<sub>rw,j</sub>为检波点j处的表层平均速度;h<sub>j</sub>为检波点j处的表层厚度;d<sub>t,j</sub>为检波点j处的延迟时;v<sub>1,j</sub>为检波点j处的折射速度;7)采用以下公式对检波点基准面静校正计算:<img file="FDA0000912485730000021.GIF" wi="702" he="126" />式中:TR<sub>j</sub>为检波点基准面静校正,j为检波点编号;E<sub>D</sub>为全区定义的水平基准面高程;Z<sub>RG,j</sub>为检波点j处的高速顶界面高程;Z<sub>R,j</sub>为检波点j处的地表高程;V<sub>R</sub>为全区定义的填充速度;V<sub>rw,j</sub>为检波点j处的表层平均速度;8)采用以下公式计算炮点基准面静校正量:<img file="FDA0000912485730000022.GIF" wi="774" he="167" />式中:TS<sub>r</sub>为炮点r点的基准面静校正;r表示测线内炮点的位置编号,范围为第一炮1到最后一炮n;E<sub>D</sub>为全区定义的水平基准面高程;Z<sub>SG,r</sub>为炮点r处的高速顶界面高程;Z'<sub>S,r</sub>为炮点r处的激发面高程;V<sub>R</sub>为全区定义的填充速度;V<sub>sw,r</sub>为炮点r处的表层平均速度,即折射层界面到炮点激发面之间的平均速度;9)将计算得到的检波点基准面静校正和炮点基准面静校正用于地震资料的处理。 |