发明名称 一种无表层调查资料约束的大炮初至综合建模静校正方法
摘要 本发明是地球物理勘探的地震勘探数据的表层建模及静校正处理方法,利用野外地震采集资料生成文件,对炮点进行井深深度校正,得到炮点激发面的高程,先对所有炮点记录拾取一定排列范围内的初至时间,然后在共中心点道集上进行折射分层,计算折射层速度和检波点接收高程位置处、炮点激发高程位置处的延迟时间;同时利用该数据库进行层析反演,得到该测线的空间速度模型,从该模型中提取与分析折射层对应的高速层速度界面,可得到高速层界面以上的表层厚度,然后利用折射分析的速度、延迟时以及表层厚度,计算出折射层以上的平均表层速度,定义静校正的基准面和速度参数,计算出基准面静校正,本发明静校正量用于地震资料的处理消除表层对地震资料成像的影响。
申请公布号 CN103869368B 申请公布日期 2016.05.11
申请号 CN201210529724.6 申请日期 2012.12.10
申请人 中国石油天然气集团公司;中国石油集团东方地球物理勘探有限责任公司 发明人 冯发全;张宇生;王海立
分类号 G01V1/36(2006.01)I 主分类号 G01V1/36(2006.01)I
代理机构 北京三友知识产权代理有限公司 11127 代理人 汤在彦
主权项 一种无表层调查资料约束的大炮初至综合建模静校正方法,特点是通过以下具体步骤实现: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>&times;</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>&times;</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>&times;</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>&times;</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)将计算得到的检波点基准面静校正和炮点基准面静校正用于地震资料的处理。
地址 100007 北京市东城区东直门北大街9号中国石油大厦