发明名称 一种用于石油勘探中确定地下速度结构的方法
摘要 本发明涉及石油勘探中确定地下速度结构的方法,是能够预测地下速度结构,确定油气藏的埋藏位置和深度的方法:是以在各检波点拾取的由各炮点所激发产生的地震记录为基准的前提下,依次进行扫描步骤、数据处理步骤、反射波速度参数的拾取步骤和拟合外推输出步骤的方法。本发明可以提高地震波速度精度及速度分析效率,减少人为因素的干扰,克服现有速度分析方法只适用炮检距与深度之比较小的条件限制,达到逐点拾取、快速、人工干预少、精度高的目的。
申请公布号 CN1257412C 申请公布日期 2006.05.24
申请号 CN200310103437.X 申请日期 2003.10.31
申请人 中国石油化工股份有限公司;中国石油化工股份有限公司石油勘探开发研究院 发明人 孙开峰;管路平;郭建;林华根;叶勇;高志凌
分类号 G01V1/28(2006.01);G01V1/00(2006.01) 主分类号 G01V1/28(2006.01)
代理机构 北京思创毕升专利事务所 代理人 韦庆文
主权项 1.一种用于石油勘探中确定地下速度结构的方法,所述方法是以在各检波点记录的由各炮点所激发产生的地震记录为基准,进行确定地下速度结构的方法,其特征在于:A 扫描步骤:输入地震记录后,对小范围偏移距地震资料以不同斜率进行线性扫描:即从0秒到地震记录结束时间,在炮点道集及检波点道集中各取5-7道地震数据,以地震记录采样间隔为步长向下滑动,扫描地震记录,其中滑动步长为采样间隔的1-3倍;B 数据处理步骤:在时间t处,炮点及检波点相邻道的时间增量Δs、Δg分别取为地震记录采样间隔的整数倍,同相轴斜率从-n到+n变化,其中n=0.8毫秒/米,增量取为0.01毫秒/米;并将相对应的5-7道地震数据求和,并比较求和结果,求和值的最大值对应的斜率即为炮点射线参数和检波点射线参数;其中的求和公式:<math> <mrow> <mi>ag</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mo>-</mo> <mi>n</mi> </mrow> <mi>n</mi> </munderover> <mi>S</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mi>k</mi> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>t</mi> <mo>+</mo> <mi>mk&Delta;g</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </math> <math> <mrow> <mi>as</mi> <mrow> <mo>(</mo> <mi>m</mi> <mo>,</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mo>-</mo> <mi>n</mi> </mrow> <mi>n</mi> </munderover> <mi>S</mi> <mrow> <mo>(</mo> <mi>i</mi> <mo>+</mo> <mi>k</mi> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mi>t</mi> <mo>+</mo> <mi>mk&Delta;s</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow> </math> C 反射波速度参数的拾取步骤:通过所述的扫描步骤和数据处理步骤后,利用炮点射线参数、检波点射线参数、旅行时间、炮点位置和检波点位置参数获得每个扫描点的反射波速度参数,其中速度公式为:<math> <mrow> <mi>v</mi> <mo>=</mo> <msqrt> <mfrac> <mrow> <mn>2</mn> <mi>ag</mi> <mo>&CenterDot;</mo> <mi>as</mi> <mo>-</mo> <mi>r</mi> <mrow> <mo>(</mo> <msub> <mi>v</mi> <msub> <mi>t</mi> <mn>1</mn> </msub> </msub> <mo>+</mo> <msub> <mi>v</mi> <msub> <mi>t</mi> <mn>2</mn> </msub> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mi>ag</mi> <mo>+</mo> <mi>as</mi> <mo>-</mo> <mn>2</mn> <mi>I</mi> <mo>/</mo> <mi>t</mi> </mrow> </mfrac> </msqrt> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow> </math> D 拟合外推步骤:对通过步骤C获得的每个扫描点的反射波速度参数进行拟合外推操作,进行速度拟合得到零偏地震波速度参数;所述方法中的拟合外推步骤为对于一反射时间t,以每道速度参数为y坐标值,以每道的位置为x坐标值,进行速度拟合外推操作,得到x=0时的速度值;根据拟合方差的大小,选择使用线性公式或更高阶多项式进行速度拟合得到零偏地震波速度参数;E 输出步骤:即将零偏地震波速度参数绘图输出的过程。
地址 100029北京市朝阳区惠新东街甲6号