发明名称 用于地震资料解释的地质小层数据深时转换方法及装置
摘要 本发明提供一种用于地震资料解释的地质小层数据深时转换方法及装置,所述方法包括:将位于目标层系上部最近的地震标志层位作为计算的起始位置,提取时间和深度初始值;从油田地质分层数据表中获取要计算的地质小层深度值,作为要进行深时转换的目标深度数值;采用声波时差曲线,从起始深度向下逐点递推进行深时转换计算,直至目标深度,求取目标地质小层对应的时间值;以目标层系下方地震解释层位为参考,比较采取测井速度计算的时间值和实际地震波旅行时间值存在的误差,采取误差平均方法开展误差校正,将校正后的结果输出,获得最终的地质小层时间值。
申请公布号 CN102967881B 申请公布日期 2015.07.08
申请号 CN201210450988.2 申请日期 2012.11.12
申请人 中国石油天然气股份有限公司 发明人 刘文岭
分类号 G01V1/28(2006.01)I;G01V1/36(2006.01)I 主分类号 G01V1/28(2006.01)I
代理机构 北京三友知识产权代理有限公司 11127 代理人 李鹤松
主权项 一种用于地震资料解释的地质小层数据深时转换方法,其特征在于,所述方法包括:将位于目标层系上部最近的地震标志层位作为计算的起始位置,提取时间和深度初始值;从油田地质分层数据表中获取要计算的地质小层深度值,作为要进行深时转换的目标深度数值;采用声波时差曲线,从起始深度向下逐点递推进行深时转换计算,直至目标深度,求取目标地质小层对应的时间值;以目标层系下方地震解释层位为参考,比较采取测井速度计算的时间值和实际地震波旅行时间值存在的误差,采取误差平均方法开展误差校正,将校正后的结果输出,获得最终的地质小层时间值;所述提取时间和深度初始值,包括:将位于目标层系上部最近的地震标志层位作为计算的起始位置,将起始层位地震解释结果插值,形成网格数据,求取该网格数据中与待求井位最近点的时间值,与对应地质层位深度值一起,构成起始层位时间值T<sub>0</sub>和深度初始值H<sub>0</sub>;所述获取要计算的地质小层深度值,包括:从油田地质分层数据表中获取要计算的地质小层深度值H<sub>n</sub>,作为要进行深时转换的目标数值;所述采用声波时差曲线从起始深度向下逐点递推进行深时转换计算,包括:以声波时差测井曲线提供的深度采样点间隔和声波时差平均速度,从起始深度H<sub>0</sub>采用如下公式向下逐点递推计算,直至深度值H<sub>n</sub>处,求取目标地质小层对应的时间值T<sub>n</sub>;深时转换的公式为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>T</mi><mi>n</mi></msub><mo>=</mo><msub><mi>T</mi><mn>0</mn></msub><mo>+</mo><mn>2</mn><mi>&Sigma;</mi><mrow><mo>(</mo><mi>&Delta;h</mi><mo>&times;</mo><mfrac><mrow><mi>&Delta;</mi><msub><mi>t</mi><mi>i</mi></msub><mo>+</mo><mi>&Delta;</mi><msub><mi>t</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow><mn>2</mn></mfrac><mo>)</mo></mrow><mo>/</mo><mn>1000</mn></mrow>]]></math><img file="FDA0000667805930000021.GIF" wi="742" he="133" /></maths>其中T<sub>n</sub>:目标地质小层对应的时间值;T<sub>0</sub>:起始层位时间值;Δh:声波时差曲线深度采样间隔,为一常数0.05m;Δt<sub>i</sub>:第i个采样点处的声波时差值;Δt<sub>i+1</sub>:第i+1个采样点处的声波时差值;所述采取误差平均方法开展误差校正,包括:采取深时转换公式计算至目标层系下方的地震解释层位,比较采取测井速度计算的时间值和实际地震波旅行时间值存在的误差,采取采样点误差平均分担的方法,开展误差校正,将校正后的结果输出,获得最终的地质小层时间值;采样点误差平均分担校正公式为:<img file="FDA0000667805930000022.GIF" wi="714" he="156" />其中<img file="FDA0000667805930000023.GIF" wi="144" he="86" />校正后的时间值;T<sub>实际</sub>:目标层系下方最近的地震标志层位实际解释时间值;T<sub>计算</sub>:采取深时转换公式计算的目标层系下方地震标志层位时间值。
地址 100007 北京市东城区东直门北大街9号