发明名称 一种曲面地表地震资料速度分析的方法
摘要 本发明为一种曲面地表地震资料速度分析的方法,属于地震勘探资料处理方法技术领域,本发明涉及针对地形起伏剧烈地区地震资料的叠加速度分析。为了突破常规基于水平层状介质和浮动基准面的速度分析技术在地形起伏剧烈地区应用条件、应用效果都受到一定的限制,本发明给出了一种基于曲地表的叠加速度分析方法,避免了常规静校正对常规速度分析结果的影响,解决了起伏地表速度分析时大的地表高差问题及反射点弥散问题,旅行时计算更为准确。本发明为地表起伏剧烈地区地震资料的速度分析提供了一种实用的技术手段,为基于曲地表的叠加与偏移成像提供了更为精确的基础资料。
申请公布号 CN102043167B 申请公布日期 2013.04.24
申请号 CN200910236263.1 申请日期 2009.10.23
申请人 中国石油化工股份有限公司;中国石油化工股份有限公司石油勘探开发研究院 发明人 潘宏勋;方伍宝;杨勤勇;林伯香;邵志东;傅星菊
分类号 G01V1/30(2006.01)I 主分类号 G01V1/30(2006.01)I
代理机构 北京思创毕升专利事务所 11218 代理人 刘明华
主权项 1.一种曲面地表地震资料速度分析的方法,其特征在于,本方法为一种基于曲面地表的叠加速度分析方法,基于曲面地表地震资料,运用曲面地表CMP道集反射波的旅行时方程得到反射波旅行时,并建立叠加速度模型;所述方法包括,①.拾取曲面地表地震资料,并进行预处理;②.建立曲面地表CMP道集反射波的旅行时方程:对于曲面地表,CMP点到反射界面的垂直距离为h<sub>m</sub>,炮点、检波点到CMP点的高差分别记作h<sub>sm</sub>、h<sub>rm</sub>,曲面地表CMP道集反射波的旅行时方程可近似地表示为<maths num="0001"><![CDATA[<math><mrow><msup><mi>t</mi><mn>2</mn></msup><mo>&ap;</mo><msup><mrow><mo>(</mo><msub><mi>t</mi><mn>0</mn></msub><mo>+</mo><mfrac><msub><mi>h</mi><mi>sm</mi></msub><msub><mi>v</mi><mi>s</mi></msub></mfrac><mo>+</mo><mfrac><msub><mi>h</mi><mi>rm</mi></msub><msub><mi>v</mi><mi>s</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><mfrac><msup><mi>x</mi><mn>2</mn></msup><msubsup><mi>v</mi><mi>s</mi><mn>2</mn></msubsup></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></math>]]></maths>式中,<img file="FDA00002244777000012.GIF" wi="202" he="123" />x为炮检距,v<sub>s</sub>为叠加速度;当只有一个水平反射界面时,(1)式是精确的表达式,v<sub>s</sub>则变为层速度;③.建立速度谱模型对于CMP道集,第i道、时窗长度τ=M·Δ,Δ为时间采样间隔,内第j个采样处的振幅值A<sub>i,j</sub>(v<sub>1</sub>),时窗内的样点数目为M,CMP点上非零值的总道数为N,用(1)式得到曲面地表反射波的旅行时间,用(2)式计算速度谱;<maths num="0002"><![CDATA[<math><mrow><msub><mi>S</mi><mi>n</mi></msub><mrow><mo>(</mo><msub><mi>v</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mi>M</mi></munderover><msup><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>A</mi><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><msub><mi>v</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>)</mo></mrow><mn>2</mn></msup></mrow><mrow><mi>N</mi><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mi>M</mi></munderover><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><msup><mi>A</mi><mn>2</mn></msup><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><mo>(</mo><msub><mi>v</mi><mn>1</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></math>]]></maths>④.输出速度数据。
地址 100728 北京市朝阳区朝阳门北大街22号