基于椭球展开的三维地震数据叠加速度求取方法
【专利摘要】本发明涉及一种基于椭球展开的三维地震数据叠加速度求取方法，包括以下步骤：在工区内划分共法向出射点面元，确定各面元位置；以面元中心所在位置为分析点，抽取该分析点处的共法向出射点道集；基于地震数据空间和叠加数据空间三维关系方程进行椭球展开并计算速度谱；在速度谱上拾取均方根速度；对所有面元分析点重复上述过程，得到各分析点处均方根速度；通过将各分析点处均方根速度进行内插外推，建立研究区叠加速度场。本发明得到的叠加速度场准确可靠，可为偏移速度分析方法提供高精度初始速度场，获得的零偏移距叠加剖面主要反射同相轴层次齐全、波组特征清楚、信噪比高。
【专利说明】基于椭球展开的三维地震数据叠加速度求取方法
【技术领域】
[0001] 本发明属于勘探地球物理学领域，具体地，涉及一种基于椭球展开的三维地震数 据叠加速度求取方法。
【背景技术】
[0002] 地球物理勘探主要涉及到两个空间：地震数据空间D和地质体空间S。其中，地震 数据空间涉及到地震观测系统、运动学特征（旅行时等）、动力学特征（振幅、频率、相位、波 形等）等信息；地质体空间描述的则是目标体几何形状、介质弹性性质、速度分布等信息。 地震勘探的一般问题即为通过研究地震数据空间信息来获得地质体空间信息。然而，该问 题的解决过程因地震信号的分辨率而复杂化。1976年Kondrashkov等基于三角映射原理 提出的椭圆展开法（EDR法）为地震数据处理提供了一种无需使用先验信息的处理方法，此 方法可以识别和获得与某个波场分量有关的信号，从而实现从空间S和空间D获取有效信 息。为实现该方法，引入第三个空间--叠加数据空间\，作为以上两个空间的补充。在地 震数据空间内地震资料展开到叠加数据空间的过程中，当且仅当所用速度正确时叠加效果 最佳，因此可以在不同扫描速度下对地震数据进行椭圆展开叠加并提取速度谱，并在能量 最大处拾取速度值构建速度场，此即椭圆法速度分析方法的基本原理。该速度分析方法未 对地下介质作水平层假设，其速度分析过程以椭圆展开和参数展开方法为基础，能够处理 不规则观测系统。与传统方法相比，该方法的主要优点在于不受地层倾角影响。需要指出的 是，在三维地震勘探情况下，现有椭圆展开速度分析方法仍基于二维反射波椭圆展开理论， 难以实现对三维地震数据的有效处理，无法满足研究与实际生产工作的要求。
【发明内容】
[0003] 为了克服现有椭圆法速度分析方法无法实现对三维地震数据进行有效处理的不 足，本发明提出一种基于椭球展开的三维地震数据叠加速度求取方法。该方法通过首先抽 取分析点处共法向出射点道集，然后进行椭球展开并求取速度谱，最后进行速度拾取从而 获得目标区域的均方根速度场。
[0004] 为实现上述目的，本发明的技术方案如下：
[0005] 基于椭球展开的三维地震数据叠加速度求取方法，其特征在于，包括以下步骤： [0006] 步骤1 :划分共法向出射点面兀
[0007] 步骤2 :抽取共法向出射点面元处的共法向出射点道集
[0008] 步骤3 :利用地震数据空间和叠加数据空间三维关系方程进行椭球展开并求取速 度谱
[0009] 步骤4 :在速度谱上进行速度拾取以确定分析点处均方根速度 [0010] 步骤5 :对所有分析点重复步骤2-4得到各分析点处均方根速度
[0011] 步骤6 :通过数据内插、外推得到最终均方根速度场
[0012] 相对于现有技术，本发明的有益效果如下：基于共法向出射点道集开展速度分析 过程，有效降低了反射点弥散的问题，提高了速度分析精度；在常规处理已抽取的CMP道集 基础上开展共法向出射点道集抽取工作，在确保得到有效信息的前提下，提高了计算效率， 适于交互速度分析和拾取；以三维椭球展开法理论为基础进行速度分析，抛开地下水平层 状假设，真正适用于三维条件下的地震资料处理，可有效获得高精度的均方根速度场，切实 提高资料使用效率和处理精度。
【专利附图】
【附图说明】
[0013] 图1是基于椭球展开的三维地震数据叠加速度求取方法的流程图；
[0014] 图2是某三维工区内某两相邻分析点的共法向出射点道集示意图；
[0015] 图3 (a)是利用基于椭球展开的三维地震数据叠加速度求取方法所得的某分析点 处共法向出射点速度谱；
[0016] 图3 (b)是使用某商业软件进行速度分析所得同一分析点处的速度谱；
[0017] 图4(a)是利用基于椭球展开的三维地震数据叠加速度求取方法所得的某分析点 处椭圆拉平校正剖面；
[0018] 图4(b)是使用某商业软件进行速度分析所得同一分析点处的拉平校正剖面；
[0019] 图5 (a)是基于椭球展开的三维地震数据叠加速度求取方法所得叠加剖面；
[0020] 图5 (b)是使用某商业软件速度分析所得叠加剖面；
[0021] 图6 (a)是基于椭球展开的三维地震数据叠加速度求取方法所得三维速度场示意 图；
[0022] 图6(b)是使用某商业软件速度分析所得三维速度场示意图。
【具体实施方式】
[0023] 如图1所示，基于椭球展开的三维地震数据叠加速度求取方法，包括以下步骤：
[0024] 步骤1 :划分共法向出射点面元
[0025] 考虑工区研究对于精度和分辨率的要求并参考CMP点相应面元的划分，在工区内 划分共法向出射点面元，确定各面元位置。具体方法如下：
[0026] 根据工区特点，确定所需目标区域速度场的分辨率与精度，计算出各目标分析点 所在三维法向出射点面元，完成共法向出射点面元划分。简单起见，可令共法向出射点面元 划分情况与CMP面元划分情况保持一致。
[0027] 步骤2 :抽取共法向出射点面元处共法向出射点道集
[0028] 以步骤1划分的各共法向出射点面元的中心作为分析点，计算出各法向出射点位 于该面元内的地震信号对应CMP点的分布范围。根据计算出的CMP点分布范围，将所有可 能具有同一法向出射点的各地震道记录按照炮检距大小排列起来形成的道集称为共法向 出射点道集，。具体方法如下：
[0029] 根据法向出射点与CMP点的几何关系，CMP道集中某一检波点处法向出射点相对 于CMP道集中法向出射点的发散距离rg为：
[0030]
【权利要求】
1. 一种基于椭球展开的三维地震数据叠加速度求取方法，其特征在于，包括以下步 骤： 步骤1 :划分共法向出射点面元 步骤2 :抽取共法向出射点面元处的共法向出射点道集 步骤3 :利用地震数据空间和叠加数据空间三维关系方程进行椭球展开并计算得到均 方根速度谱 步骤4 :在速度谱上进行速度拾取主要能量团的均方根速度 步骤5 :对所有分析点重复步骤2-4,得到各分析点处均方根速度 步骤6 :对各分析点的均方根速度进行内插外推，得到最终均方根速度场。
2. 根据权利要求1所述的基于椭球展开的三维地震数据叠加速度求取方法，其特征在 于，步骤1为：根据工区特点，确定目标区域速度场的分辨率与精度，计算出各目标分析点 所在三维法向出射点面元，完成共法向出射点面元划分。为简单起见，可令共法向出射点面 元划分情况与CMP面元划分情况保持一致。
3. 根据权利要求1-2所述的基于椭球展开的三维地震数据叠加速度求取方法，其特征 在于，步骤2为：以步骤1划分的各共法向出射点面元的中心作为分析点，计算出各法向出 射点位于该面元内的地震信号对应CMP点的分布范围。根据计算出的CMP点分布范围，对 所有CMP点的道集记录按按照炮检距大小重新排列起来形成新的道集称为共法向出射点 道集。具体方法如下： 根据法向出射点与CMP点的几何关系，CMP道集中某一检波点处法向出射点相对于CMP 道集中法向出射点的发散距离rg为： rg =- -sin 沒， CMP道集法向出射点最大散布范围半径rg(l为： I 2 Winax 法向出射点处对应CMP点的分布范围为 X _ range = rgi} -f 0.5gridx 11 5 y _range = rg{) t- 0.5gridy 其中法向出射点是指地下反射点处界面法线与地表的交点。以上各式中，hB为地层法 向深度，1为地震信号对应的偏移距，Θ为地下界面倾角，hmax为目地层的最大埋深，为 最大偏移距，Θ ^为深层地下界面倾角估计值，X_range、y_range分别表示所有法向出射点 分散范围与分析点处法向出射点面元有交点的CMP道集对应CMP点的分布范围在工区内平 面直角坐标系L Qx-LQy中LQx、LQy方向的网格边长，gridx、gridy为L Qx、LQy方向的分析点所在 法向出射点面兀网格边长。 以目标区域内的面元中心作为分析点，计算出隶属于该面元的所有CMP道集所对应的 CMP点分布范围，提取该范围内所有CMP道集，将所有CMP道集记录按照炮检距大小排列起 来组成共法向出射点道集。
4. 根据权利要求1-3所述的基于椭球展开的三维地震数据叠加速度求取方法，其特征 在于，步骤3为对步骤2抽取得到的分析点处共法向出射点道集进行三维椭球展开，计算得 到零偏移距域均方根速度谱。具体方法如下： 三维地震资料反射波椭球展开方程：
其中，I^-Lw为一般工区所用的固定平面坐标系，〇^，1^)为法向出射点平面坐标，在 椭球展开速度分析过程中代表当前分析点的平面坐标；（Lto，L_)分别表示平 面坐标系Ι^χ-Ι^中炮点、检波点的平面坐标；L为零偏移距域内沿零偏移距射线的双程旅 行时（亦即自激自收时间），v为地震波传播速度，t为进行椭球展开的反射波的旅行时；变 量L、Θ的定义则如上式所示。 选定速度扫描范围，利用椭球展开方程式对当前分析点处共法向出射点道集中各地震 道依次进行椭球展开并进行相干叠加，计算得到零偏移距域速度谱。
5. 根据权利要求1-4所述的基于椭球展开的三维地震数据叠加速度求取方法，其特征 在于，步骤4将步骤3所得速度谱进行交互拾取，获得当前分析点处主要能量团的均方根速 度。具体方法如下： 依据共法向出射点道集在零偏移距域内的拉平情况和叠加效果，在椭球展开过程提 取的零偏移距域速度谱上，由浅至深拾取能量聚焦点处的速度值作为当前分析点处地下介 质的均方根速度，从而确定当前速度分析点的速度分布情况。
6. 根据权利要求1-5所述的基于椭球展开的三维地震数据叠加速度求取方法，其特征 在于，步骤5为对其它分析点重复步骤2-4的速度分析过程，直至处理完所有分析点为止， 拾取所有分析点各主要能量团的均方根速度。
7. 根据权利要求1-6所述的基于椭球展开的三维地震数据叠加速度求取方法，其特征 在于，步骤6为：以步骤5所得的均方根速度为输入，在工区目标区域范围内，利用三维插值 方法进行内插、外推，求取工区范围内各离散网格点处的均方根速度，最终获得工区内目标 区域的叠加速度场。
【文档编号】G01V1/30GK104155690SQ201410391204
【公开日】2014年11月19日 申请日期:2014年8月8日 优先权日:2014年8月8日
【发明者】杜启振, 李文滨, 韦亮 申请人:中国石油大学(华东)