一种基于高斯射线束目的层能量照明的补炮方法
一种基于高斯射线束目的层能量照明的补炮方法
【技术领域】
[0001] 本发明涉及地球物理勘探技术领域，尤其涉及一种基于高斯射线束目的层能量照 明的补炮方法。
【背景技术】
[0002] 地震采集是地震勘探的第一环节，是决定整个地震资料品质好坏的最关键一步， 而观测系统设计又是地震采集环节中的重点。观测系统设计主要是合理设置炮点和检波点 的分布位置，提高野外数据采集工作的效率和准确程度。
[0003] 常规观测系统设计主要是基于水平层状介质理论，观测系统属性参数都是基于水 平层状介质得到的，没有考虑实际地质构造的影响。由于实际地质构造中的目的层并非呈 严格的水平层状分布，现有观测系统得到的地震资料会出现信号能量较弱的区域。
[0004] 送部分能量较弱的区域中信噪比较低，常规观测系统中不会对能量弱区做任何处 理，因此观测到的地震资料难W达到地震处理和解释的要求。导致在后续的地震处理过程 中不能针对能量弱区还原真实的地下地质构造。
[0005] 因此，亟需一种能够对目的层能量较弱的区域进行补炮的方法。
【发明内容】
[0006] 本发明针对现有技术的不足，提出了一种基于高斯射线束目的层能量照明的补炮 方法，包括W下步骤：
[0007] 根据现有地质资料建立勘探区域的Η维地质模型，基于炮点高斯射线束对地质模 型做照明分析来确定目的层面元偏移能量；
[0008] 统计目的层面元偏移能量分布情况W圈定目的层中能量弱区；
[0009] 根据各炮点对能量弱区的照明参考值确定目的层中能量弱区对应的地表补炮区 域；
[0010] 在地表补炮区域的原有炮点网格中等间隔设置补炮炮点。
[0011] 根据本发明的一个实施例，所述基于炮点高斯射线束对地质模型做照明分析来确 定目的层面元偏移能量的步骤中，进一步包括：
[0012] 基于炮点高斯射线束对地质模型做照明分析，确定入射到每一目的层面元的全部 高斯射线束的入射能量；
[0013] 计算经由目的层反射之后，由每一地表检波器接收的覆盖检波器所在网格范围的 高斯射线束能量；
[0014] 根据入射到一个目的层面元的并且由全部地表检波器接收的高斯线束能量确定 目的层面元偏移能量。
[0015] 根据本发明的一个实施例，所述基于炮点高斯射线束对地质模型做照明分析，确 定入射到每一目的层面元的全部高斯射线束的入射能量包括：
[0016] 统计入射到目的层的全部高斯射线束；
[0017] 分别统计覆盖每一目的层面元的高斯射线束；
[0018] 叠加计算入射到每一目的层面元的高斯射线束的入射能量，得到每一目的层面元 入射能量。
[0019] 根据本发明的一个实施例，所述计算经由目的层反射之后，由每一地表检波器接 收的覆盖检波器所在网格范围的高斯射线束能量包括：
[0020] 统计通过目的层反射至地表的所有高斯射线束；
[0021] 分别统计覆盖每一地表检波器所在网格范围的高斯射线束；
[0022] 叠加计算由每一地表检波器接收的全部高斯射线束的入射能量，得到每一地表检 波器的接收能量。
[0023] 根据本发明的一个实施例，所述根据入射到一个目的层面元的并且由全部地表检 波器接收的高斯线束能量确定目的层面元偏移能量包括：
[0024] 计算每一目的层面元入射能量占炮点高斯射线束总入射能量的比值；
[0025] 统计每一地表检波器的接收能量，计算所有地表检波器接收的由一个目的层面元 反射的来自每一炮点的高斯射线束能量；
[0026] 叠加计算由目的层面元反射的来自所有炮点的高斯射线束能量，获得目的层面元 的有效反射能量；
[0027] 由目的层面元的有效反射能量与目的层面元的入射能量占炮点高斯射线束总入 射能量的比值的乘积得到目的层面元的面元偏移能量。
[0028] 根据本发明的一个实施例，所述根据各炮点对能量弱区的照明参考值确定目的层 中能量弱区对应的补炮区域包括：
[0029] 统计入射到目的层能量弱区的高斯射线束对应的炮点；
[0030] 分别计算入射到目的层能量弱区的来自各炮点的高斯射线束能量与能量弱区总 入射能量的比值，获得各炮点对能量弱区的照明参考值；
[0031] 选择对能量弱区的照明参考值预大于或等于设阔值的炮点组成地表补炮区域。
[0032] 根据本发明的一个实施例，所述在补炮区域的原有炮点网格中等间隔设置补炮炮 点为：
[0033] 在地表补炮区域的原有的两个相邻炮点的中间位置设置一个补炮炮点。
[0034] 根据本发明的一个实施例，所述根据现有地质资料建立勘探区域的Η维地质模型 包括：
[0035] 利用现有地质资料分别构建中深层构造模型和地表模型；
[0036] 设定每一地层的纵波速度、横波速度和密度，得到Η维块状地质模型。
[0037] 根据本发明的一个实施例，所述目的层面元的面元偏移能量表示为：
[0038]
[003引其中，（i，如表示目的层中的水平坐标，Ρ为目的层面元反射到地表的射线束编 号，ER(p)为目的层面元内第P条射线束反射到达地面被所有检波器接收能量叠加总和，N 为目的层面元反射到地面的射线束的总数，PCell (i，j)表示目的层面元入射能量占炮点射 线束总入射能量的比值。
[0040] 根据本发明的一个实施例，所述预设能量阔值为60%。
[0041] 本发明提出目的层面元偏移能量的概念用于表征真正可用于地面成像的地下面 元有效接收能量，从而可W定量评价观测系统设计的优劣，更合理地设计观测系统参数。同 时面元偏移能量的提出使得地震采集与资料处理成像更加一体化，为W成像为目的的地震 采集设计提供依据。
[0042] 本发明还可W准确圈定目的层的能量弱区，定量确定补炮数目，从而使得目的层 中照明能量弱的位置得到补偿，提高地震资料处理后的信噪比，改善地震采集资料品质。
[0043] 本发明的其它特征和优点将在随后的说明书中阐述，并且部分地从说明书中变得 显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要 求书W及附图中所特别指出的结构来实现和获得。
【附图说明】
[0044] 图la是根据本发明实施例的高斯射线束振幅分布示意图；
[0045] 图化是根据本发明实施例的基于高斯射线束能量叠加的原理示意图；
[0046] 图2是根据本发明实施例的射线坐标系示意图；
[0047] 图3是根据本发明实施例的基于高斯射线束目的层能量照明的补炮方法的步骤 流程图；
[0048] 图4是根据本发明实施例的确定目的层面元偏移能量的步骤流程图；
[0049] 图5a是入射到目的层面元的高斯射线束分布示意图；
[0050] 图化是地表检波器接收到的目的层反射线束的分布示意图；
[0051] 图5c是目的层面元偏移能量射线束分布示意图；
[0052] 图6a是多线束入射照明的Η维地质结构示意图；
[0053] 图化是多线束接 收照明的Η维地质结构示意图；
[0054] 图6c是多线束偏移照明的Η维地质结构示意图；
[00巧]图7是某工程示例的Η维地质构造模型示意图；
[0056] 图8是补炮前观测系统对目的层照明的分布图；
[0057] 图9是补炮后观测系统对目的层照明的分布图；
[0058] 图10a是地表补炮区域的分布图；
[0059] 图10b是地表补炮区域的放大显示结果。
【具体实施方式】
[0060] 本发明的实施例中构建基于勘探区域的Η维地质模型，利用基于高斯束射线束的 照明技术获得地下目的层的偏移能量分布。由分布结果找到目的层能量弱的位置适当在地 表相应位置进行补炮，定量确定补炮数目，使得目的层中照明能量弱的位置得到补偿。使得 地震资料处理后目的层能量和信噪比得到提高，最终达到改善地震采集资料品质的目的。
[0061] W下对基于高斯射线束的照明技术进行详细说明。
[0062] 如图la所示，高斯射线束可W看作是一条从震源出发W射线为中必的能量管，射 线束的振幅分布W偏离中必射线的距离平方呈指数衰减。如图化所示，接收点R或地下成 像点处的波场可看作是由多条从震源点sO出发，在R点一定范围内的高斯射线束能量的叠 加。
[0063] 在图2中，P为空间一点，其在中必射线S上的垂直投影点为P'，s为P'点到震源 点sO的弧长，η表示P'点到P点的距离，即P点到中必射线S的距离。图2中射线坐标系 由向量η(沿射线的法矢量）和t (沿射线的切矢量）来定义。
[0064] 高斯射线束公式是建立在射线中必坐标系下的波动方程集中于射线附近的高频 渐近解。中必射线附近高频能量的分布用下式表示：
[0065]
(1)
[0066] 上式中U表7K纵波位移，指数部分一号表W正向延拓（正演），+号表W反向延拓。 Ψ为某一设定常数值，ω表示高斯射线束波场频率。（S，η)为计算点在中必射线坐标系下 的射线坐标。V表示中必射线的速度，τ为中必射线旅行时。
[0067] (1)式中P(s)和q(s)为沿中必射线变化的复值动力学参数，满足如下常微分方程 组：
[0068]
(2)
[0069] 对前面的公式进行变换，可W得到如下更具物理意义的高斯射线束表达形式：
[0075] 式中A为高斯射线束振幅，K表示射线束的波前曲率，L表示射线的有效半宽度，L 决定了高斯射线束振幅在中必射线附近的分布，由于振幅分布类似于高斯分布故而得名为 高斯射线束。送样可W利用上面的公式进行高斯射线束的正演。
[0076] 本发明的实施例提供一种基于高斯射线束目的层能量照明的补炮方法。图3所示 为本实施例的补炮方法的步骤流程图。
[0077] 首先，在步骤S101中根据现有地质资料建立勘探区域的Η维地质模型，基于炮点 高斯射线束对地质模型做照明分析，确定目的层面元偏移能量。
[0078] 其中，Η维地质建模是本实施例的后续步骤中照明和偏移的基础。本实施例主要 针对射线束进行建模，因此需要确定Η维块状地质模型。
[0079] Η维建模的基本方法是利用插值构造层位、断层曲面，通过裁剪、缝合建立拓扑一 致的曲面模型，利用自动Η维块体追踪方法建立Η维块体模型。具体来说，利用过工区的深 度域解释剖面或深度域层位数据/时间域层位数据得到中深层构造模型，其中，时间域数 据需要时深尺进行时深转换；然后，利用GIS高程或SI^检波点坐标生成真实地表模型。进 而设定每一地层的纵波速度、横波速度和密度，得到Η维块状地质模型。
[0080] 利用高斯射线束进行正向照明，可W将虚拟检波器放置在反射界面或地表来计算 虚拟检波器所接收到的能量，确定目的层面元偏移能量。
[0081] W下参考图4详细说明步骤S101中确定目的层面元偏移能量的具体过程。
[0082] 首先，在子步骤S1011中基于炮点高斯射线束对地质模型做照明分析，确定入射 到每一目的层面元的全部高斯射线束的入射能量。
[0083] 具体而言，如图5a所示，统计所有入射到目的层的高斯射线束，从而分别统计覆 盖每一目的层面元的全部高斯射线束。叠加计算入射到每一目的层面元上的全部高斯射线 束的入射能量，得到每一目的层面元入射能量ECell (X，y)，形成图6a所示的多束射线的入 射照明结果。
[0084] 接下来，执行子步骤S1012,计算经由目的层反射之后，由每一地表检波器接收的 覆盖检波器所在网格范围的高斯射线束能量。即计算检波器的接收能量，也就是检波器接 收到来自炮点的能量。
[0085] 如图化所示，统计所有通过目的层反射至地表的高斯射线束，分别统计覆盖每一 地表检波器所在网格范围的高斯射线束。叠加计算由每一地表检波器接收的全部高斯射线 束的入射能量，得到每一地表检波器的接收能量邸ecv(x，y)，形成图化所示的多束射线的 接收照明结果。
[0086] 随后，在步骤S1013中根据入射到一个目的层面元的并且由全部地表检波器接收 的高斯线束能量确定目的层面元偏移能量，即通过面元反射且被地表检波器接收到的可用 于实现偏移成像的有效能量。
[0087] 如图5c中所示，面元偏移能量是指震源激发的地震波通过目的层面元反射后，被 观测系统中多个检波器接收到的能量总和，代表了真正可W用于成像地下面元的有效接收 能量。
[0088] 具体来说，首先计算每一目的层面元入射能量占炮点射线束总入射能量的比值 PCelia，j)。
[0089] 统计每一地表检波器的接收能量，计算所有地表检波器接收的由一个目的层面元 反射的来自每一炮点的高斯射线束能量，然后叠加计算由一个目的层面元反射的来自所有 炮点的高斯射线束能量，获得一个目的层面元的有效反射能量。
[0090] 最后，计算该目的层面元的有效反射能量与该目的层面元的入射能量占炮点高斯 射线束总入射能量的比值的乘积，得到该目的层面元的面元偏移能量：
[0091]
。）
[0092] 其中，（i，j)表示目的层中的水平坐标，P为目的层面元反射到地面的射线束编 号，ER(p)为目的层面元内第P条射线束反射到达地面被所有检波器接收能量叠加总和，N 为目的层面元反射到地面的射线束的总数，
表示目的层面元的有效反射能量。
[0093] ER(p)的数值可W根据表达式（1)至（6)来计算，再根据表达式（7)可获得图6c 所示的多束射线的偏移照明结果。送样W来，可W基于高斯束射线束的照明技术对Η维地 质模型进行分析，获得目的层的偏移能量分布。
[0094] 与现有技术相比，本实施例中提出目的层面元偏移能量的概念用于表征真正可用 于地面成像的地下面元有效接收能量，从而可W定量评价观测系统设计的优劣，更合理地 设计观测系统参数。同时面元偏移能量的提出使得地震采集与资料处理成像更加一体化， 为W成像为目的的地震采集设计提供依据。
[009引随后，执化步骤S102,统计目的层面兀偏移战垂分布情况W圈志目的层中战垂弱 区。
[0096] 具体的，首先选定地震勘探的目的层，即此次地震勘探主要想解决地下哪一目的 层的地质问题，然后针对选定的地下目 的层进行面元偏移能量统计，根据地质区域的实际 构造或者勘探目的设定能量弱区的阔值，圈定目的层中能量弱区。
[0097] 接下来，在步骤S103中根据各炮点对能量弱区的照明参考值确定目的层中能量 弱区对应的补炮区域。
[009引首先设定炮点照射到能量弱区的能量阔值。例如可W根据炮点照射能量占区域总 能量的60 %作为能量阔值。
[0099] 然后统计到达目的层能量弱区的高斯射线束对应的炮点，分别计算到达目的层能 量弱区的来自各炮点的高斯射线束能量与能量弱区总入射能量的比值，获得各炮点对能量 弱区的照明参考值。选择对能量弱区的照明参考值大于或等于预设阔值的炮点组成地表补 炮区域。
[0100] 最后，执行步骤S104,在补炮区域的原有炮点网格中等间隔设置补炮炮点。
[0101] 优选的，在地表补炮区域的原有的两个相邻炮点的中间位置设置一个补炮炮点， 送样设置加密炮点可使得原有炮点网格加密一倍。不限于此，也可在原有的两个相邻炮点 的中间位置设置两个补炮炮点，或者更多，可根据实施成本具体设定。如此W来，可W通过 在补炮区域加密炮点，使得目的层中照明能量弱的位置得到补偿，提高地震资料的信噪比。 引 工招连例
[010引下面W D邸地区为例，详细说明整个补炮过程。观测系统参数如下表1所示。
[0104] 表 1
[0105]
[0106] 通过过工区的Η条深度剖面建立一个Η维构造模型（如图7所示）。图8是补炮 前观测系统对目的层照明的分布，从图8可知，面元偏移能量分布差别很大，因此选中图中 黑圆圈所示位置进行补炮演示。图9为补炮后的目的层能量的分布，可知黑圆圈位置能量 得W补偿（补偿炮数为100炮，W目的层最大能量为标准）。图10a和图1化为地面补炮位 置和补炮数的放大显示，从图中可W看到，补炮的位置都是在原先炮点网格加密后的位置， 使得原有炮点网格成倍加密。对比图8和图9可W看出，本发明可W准确圈定目的层的能 量弱区，定量确定补炮数目，从而使得目的层中照明能量弱的位置得到补偿，提高地震资料 处理后的信噪比，改善地震采集资料品质。
[0107] 虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采 用的实施方式，并非用W限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本 发明所公开的精神和范围的前提下，可W在实施的形式上及细节上作任何的修改与变化， 但本发明的专利保护范围，仍须W所附的权利要求书所界定的范围为准。
【主权项】
1. 一种基于高斯射线束目的层能量照明的补炮方法，其特征在于，包括以下步骤： 根据现有地质资料建立勘探区域的三维地质模型，基于炮点高斯射线束对地质模型做 照明分析来确定目的层面元偏移能量； 统计目的层面元偏移能量分布情况以圈定目的层中能量弱区； 根据各炮点对能量弱区的照明参考值确定目的层中能量弱区对应的地表补炮区域； 在地表补炮区域的原有炮点网格中等间隔设置补炮炮点。2. 根据权利要求1所述的方法，其特征在于，所述基于炮点高斯射线束对地质模型做 照明分析来确定目的层面元偏移能量的步骤包括： 基于炮点高斯射线束对地质模型做照明分析，确定入射到每一目的层面元的全部高斯 射线束的入射能量； 计算经由目的层反射之后，由每一地表检波器接收的覆盖检波器所在网格范围的高斯 射线束能量； 根据入射到一个目的层面元的并且由全部地表检波器接收的高斯线束能量确定目的 层面元偏移能量。3. 根据权利要求2所述的方法，其特征在于，所述基于炮点高斯射线束对地质模型做 照明分析，确定入射到每一目的层面元的全部高斯射线束的入射能量包括： 统计入射到目的层的全部高斯射线束； 分别统计覆盖每一目的层面元的高斯射线束； 叠加计算入射到每一目的层面元的高斯射线束的入射能量，得到每一目的层面元入射 能量。4. 根据权利要求3所述的方法，其特征在于，所述计算经由目的层反射之后，由每一地 表检波器接收的覆盖检波器所在网格范围的高斯射线束能量包括： 统计通过目的层反射至地表的所有高斯射线束； 分别统计覆盖每一地表检波器所在网格范围的高斯射线束； 叠加计算由每一地表检波器接收的全部高斯射线束的入射能量，得到每一地表检波器 的接收能量。5. 根据权利要求4所述的方法，其特征在于，所述根据入射到一个目的层面元的并且 由全部地表检波器接收的高斯线束能量确定目的层面元偏移能量包括： 计算每一目的层面元入射能量占炮点高斯射线束总入射能量的比值； 统计每一地表检波器的接收能量，计算所有地表检波器接收的由一个目的层面元反射 的来自每一炮点的高斯射线束能量； 叠加计算由目的层面元反射的来自所有炮点的高斯射线束能量，获得目的层面元的有 效反射能量； 由目的层面元的有效反射能量与目的层面元的入射能量占炮点高斯射线束总入射能 量的比值的乘积得到目的层面元的面元偏移能量。6. 根据权利要求1-5中任一项所述的方法，其特征在于，所述根据各炮点对能量弱区 的照明参考值确定目的层中能量弱区对应的补炮区域包括： 统计入射到目的层能量弱区的高斯射线束对应的炮点； 分别计算入射到目的层能量弱区的来自各炮点的高斯射线束能量与能量弱区总入射 能量的比值，获得各炮点对能量弱区的照明参考值； 选择对能量弱区的照明参考值大于或等于预设能量阈值的炮点组成地表补炮区域。7. 根据权利要求1-5中任一项所述的方法，其特征在于，所述在补炮区域的原有炮点 网格中等间隔设置补炮炮点为： 在地表补炮区域的原有的两个相邻炮点的中间位置设置一个补炮炮点。8. 根据权利要求1-5中任一项所述的方法，其特征在于，所述根据现有地质资料建立 勘探区域的三维地质模型包括： 利用现有地质资料分别构建中深层构造模型和地表模型； 设定每一地层的纵波速度、横波速度和密度，获得三维块状地质模型。9. 根据权利要求5所述的方法，其特征在于，所述目的层面元的面元偏移能量表示为：其中，（i，j)表示目的层中的水平坐标，P为目的层面元反射到地表的射线束编号， ER(P)为目的层面元内第p条射线束反射到达地面被所有检波器接收能量叠加总和，N为目 的层面元反射到地面的射线束的总数，PCelUi，j)表示目的层面元入射能量占炮点射线束 总入射能量的比值。10. 根据权利要求6所述的方法，其特征在于，所述预设能量阈值为60%。
【专利摘要】本发明公开了一种基于高斯射线束目的层能量照明的补炮方法，包括以下步骤：根据现有地质资料建立勘探区域的三维地质模型，基于炮点高斯射线束对地质模型做照明分析来确定目的层面元偏移能量；统计目的层面元偏移能量分布情况以圈定目的层中能量弱区；根据各炮点对能量弱区的照明参考值确定目的层中能量弱区对应的地表补炮区域；在地表补炮区域的原有炮点网格中等间隔设置补炮炮点。本发明可以准确圈定目的层的能量弱区，定量确定补炮数目，从而使得目的层中照明能量弱的位置得到补偿，提高地震资料处理后的信噪比，改善地震采集资料品质。
【IPC分类】G01V1/00
【公开号】CN105487106
【申请号】CN201410479944
【发明人】肖云飞, 殷厚成
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院
【公开日】2016年4月13日
【申请日】2014年9月18日
【技术领域】
[0001] 本发明涉及地球物理勘探技术领域，尤其涉及一种基于高斯射线束目的层能量照 明的补炮方法。
【背景技术】
[0002] 地震采集是地震勘探的第一环节，是决定整个地震资料品质好坏的最关键一步， 而观测系统设计又是地震采集环节中的重点。观测系统设计主要是合理设置炮点和检波点 的分布位置，提高野外数据采集工作的效率和准确程度。
[0003] 常规观测系统设计主要是基于水平层状介质理论，观测系统属性参数都是基于水 平层状介质得到的，没有考虑实际地质构造的影响。由于实际地质构造中的目的层并非呈 严格的水平层状分布，现有观测系统得到的地震资料会出现信号能量较弱的区域。
[0004] 送部分能量较弱的区域中信噪比较低，常规观测系统中不会对能量弱区做任何处 理，因此观测到的地震资料难W达到地震处理和解释的要求。导致在后续的地震处理过程 中不能针对能量弱区还原真实的地下地质构造。
[0005] 因此，亟需一种能够对目的层能量较弱的区域进行补炮的方法。
【发明内容】
[0006] 本发明针对现有技术的不足，提出了一种基于高斯射线束目的层能量照明的补炮 方法，包括W下步骤：
[0007] 根据现有地质资料建立勘探区域的Η维地质模型，基于炮点高斯射线束对地质模 型做照明分析来确定目的层面元偏移能量；
[0008] 统计目的层面元偏移能量分布情况W圈定目的层中能量弱区；
[0009] 根据各炮点对能量弱区的照明参考值确定目的层中能量弱区对应的地表补炮区 域；
[0010] 在地表补炮区域的原有炮点网格中等间隔设置补炮炮点。
[0011] 根据本发明的一个实施例，所述基于炮点高斯射线束对地质模型做照明分析来确 定目的层面元偏移能量的步骤中，进一步包括：
[0012] 基于炮点高斯射线束对地质模型做照明分析，确定入射到每一目的层面元的全部 高斯射线束的入射能量；
[0013] 计算经由目的层反射之后，由每一地表检波器接收的覆盖检波器所在网格范围的 高斯射线束能量；
[0014] 根据入射到一个目的层面元的并且由全部地表检波器接收的高斯线束能量确定 目的层面元偏移能量。
[0015] 根据本发明的一个实施例，所述基于炮点高斯射线束对地质模型做照明分析，确 定入射到每一目的层面元的全部高斯射线束的入射能量包括：
[0016] 统计入射到目的层的全部高斯射线束；
[0017] 分别统计覆盖每一目的层面元的高斯射线束；
[0018] 叠加计算入射到每一目的层面元的高斯射线束的入射能量，得到每一目的层面元 入射能量。
[0019] 根据本发明的一个实施例，所述计算经由目的层反射之后，由每一地表检波器接 收的覆盖检波器所在网格范围的高斯射线束能量包括：
[0020] 统计通过目的层反射至地表的所有高斯射线束；
[0021] 分别统计覆盖每一地表检波器所在网格范围的高斯射线束；
[0022] 叠加计算由每一地表检波器接收的全部高斯射线束的入射能量，得到每一地表检 波器的接收能量。
[0023] 根据本发明的一个实施例，所述根据入射到一个目的层面元的并且由全部地表检 波器接收的高斯线束能量确定目的层面元偏移能量包括：
[0024] 计算每一目的层面元入射能量占炮点高斯射线束总入射能量的比值；
[0025] 统计每一地表检波器的接收能量，计算所有地表检波器接收的由一个目的层面元 反射的来自每一炮点的高斯射线束能量；
[0026] 叠加计算由目的层面元反射的来自所有炮点的高斯射线束能量，获得目的层面元 的有效反射能量；
[0027] 由目的层面元的有效反射能量与目的层面元的入射能量占炮点高斯射线束总入 射能量的比值的乘积得到目的层面元的面元偏移能量。
[0028] 根据本发明的一个实施例，所述根据各炮点对能量弱区的照明参考值确定目的层 中能量弱区对应的补炮区域包括：
[0029] 统计入射到目的层能量弱区的高斯射线束对应的炮点；
[0030] 分别计算入射到目的层能量弱区的来自各炮点的高斯射线束能量与能量弱区总 入射能量的比值，获得各炮点对能量弱区的照明参考值；
[0031] 选择对能量弱区的照明参考值预大于或等于设阔值的炮点组成地表补炮区域。
[0032] 根据本发明的一个实施例，所述在补炮区域的原有炮点网格中等间隔设置补炮炮 点为：
[0033] 在地表补炮区域的原有的两个相邻炮点的中间位置设置一个补炮炮点。
[0034] 根据本发明的一个实施例，所述根据现有地质资料建立勘探区域的Η维地质模型 包括：
[0035] 利用现有地质资料分别构建中深层构造模型和地表模型；
[0036] 设定每一地层的纵波速度、横波速度和密度，得到Η维块状地质模型。
[0037] 根据本发明的一个实施例，所述目的层面元的面元偏移能量表示为：
[0038]
[003引其中，（i，如表示目的层中的水平坐标，Ρ为目的层面元反射到地表的射线束编 号，ER(p)为目的层面元内第P条射线束反射到达地面被所有检波器接收能量叠加总和，N 为目的层面元反射到地面的射线束的总数，PCell (i，j)表示目的层面元入射能量占炮点射 线束总入射能量的比值。
[0040] 根据本发明的一个实施例，所述预设能量阔值为60%。
[0041] 本发明提出目的层面元偏移能量的概念用于表征真正可用于地面成像的地下面 元有效接收能量，从而可W定量评价观测系统设计的优劣，更合理地设计观测系统参数。同 时面元偏移能量的提出使得地震采集与资料处理成像更加一体化，为W成像为目的的地震 采集设计提供依据。
[0042] 本发明还可W准确圈定目的层的能量弱区，定量确定补炮数目，从而使得目的层 中照明能量弱的位置得到补偿，提高地震资料处理后的信噪比，改善地震采集资料品质。
[0043] 本发明的其它特征和优点将在随后的说明书中阐述，并且部分地从说明书中变得 显而易见，或者通过实施本发明而了解。本发明的目的和其他优点可通过在说明书、权利要 求书W及附图中所特别指出的结构来实现和获得。
【附图说明】
[0044] 图la是根据本发明实施例的高斯射线束振幅分布示意图；
[0045] 图化是根据本发明实施例的基于高斯射线束能量叠加的原理示意图；
[0046] 图2是根据本发明实施例的射线坐标系示意图；
[0047] 图3是根据本发明实施例的基于高斯射线束目的层能量照明的补炮方法的步骤 流程图；
[0048] 图4是根据本发明实施例的确定目的层面元偏移能量的步骤流程图；
[0049] 图5a是入射到目的层面元的高斯射线束分布示意图；
[0050] 图化是地表检波器接收到的目的层反射线束的分布示意图；
[0051] 图5c是目的层面元偏移能量射线束分布示意图；
[0052] 图6a是多线束入射照明的Η维地质结构示意图；
[0053] 图化是多线束接 收照明的Η维地质结构示意图；
[0054] 图6c是多线束偏移照明的Η维地质结构示意图；
[00巧]图7是某工程示例的Η维地质构造模型示意图；
[0056] 图8是补炮前观测系统对目的层照明的分布图；
[0057] 图9是补炮后观测系统对目的层照明的分布图；
[0058] 图10a是地表补炮区域的分布图；
[0059] 图10b是地表补炮区域的放大显示结果。
【具体实施方式】
[0060] 本发明的实施例中构建基于勘探区域的Η维地质模型，利用基于高斯束射线束的 照明技术获得地下目的层的偏移能量分布。由分布结果找到目的层能量弱的位置适当在地 表相应位置进行补炮，定量确定补炮数目，使得目的层中照明能量弱的位置得到补偿。使得 地震资料处理后目的层能量和信噪比得到提高，最终达到改善地震采集资料品质的目的。
[0061] W下对基于高斯射线束的照明技术进行详细说明。
[0062] 如图la所示，高斯射线束可W看作是一条从震源出发W射线为中必的能量管，射 线束的振幅分布W偏离中必射线的距离平方呈指数衰减。如图化所示，接收点R或地下成 像点处的波场可看作是由多条从震源点sO出发，在R点一定范围内的高斯射线束能量的叠 加。
[0063] 在图2中，P为空间一点，其在中必射线S上的垂直投影点为P'，s为P'点到震源 点sO的弧长，η表示P'点到P点的距离，即P点到中必射线S的距离。图2中射线坐标系 由向量η(沿射线的法矢量）和t (沿射线的切矢量）来定义。
[0064] 高斯射线束公式是建立在射线中必坐标系下的波动方程集中于射线附近的高频 渐近解。中必射线附近高频能量的分布用下式表示：
[0065]
(1)
[0066] 上式中U表7K纵波位移，指数部分一号表W正向延拓（正演），+号表W反向延拓。 Ψ为某一设定常数值，ω表示高斯射线束波场频率。（S，η)为计算点在中必射线坐标系下 的射线坐标。V表示中必射线的速度，τ为中必射线旅行时。
[0067] (1)式中P(s)和q(s)为沿中必射线变化的复值动力学参数，满足如下常微分方程 组：
[0068]
(2)
[0069] 对前面的公式进行变换，可W得到如下更具物理意义的高斯射线束表达形式：
[0075] 式中A为高斯射线束振幅，K表示射线束的波前曲率，L表示射线的有效半宽度，L 决定了高斯射线束振幅在中必射线附近的分布，由于振幅分布类似于高斯分布故而得名为 高斯射线束。送样可W利用上面的公式进行高斯射线束的正演。
[0076] 本发明的实施例提供一种基于高斯射线束目的层能量照明的补炮方法。图3所示 为本实施例的补炮方法的步骤流程图。
[0077] 首先，在步骤S101中根据现有地质资料建立勘探区域的Η维地质模型，基于炮点 高斯射线束对地质模型做照明分析，确定目的层面元偏移能量。
[0078] 其中，Η维地质建模是本实施例的后续步骤中照明和偏移的基础。本实施例主要 针对射线束进行建模，因此需要确定Η维块状地质模型。
[0079] Η维建模的基本方法是利用插值构造层位、断层曲面，通过裁剪、缝合建立拓扑一 致的曲面模型，利用自动Η维块体追踪方法建立Η维块体模型。具体来说，利用过工区的深 度域解释剖面或深度域层位数据/时间域层位数据得到中深层构造模型，其中，时间域数 据需要时深尺进行时深转换；然后，利用GIS高程或SI^检波点坐标生成真实地表模型。进 而设定每一地层的纵波速度、横波速度和密度，得到Η维块状地质模型。
[0080] 利用高斯射线束进行正向照明，可W将虚拟检波器放置在反射界面或地表来计算 虚拟检波器所接收到的能量，确定目的层面元偏移能量。
[0081] W下参考图4详细说明步骤S101中确定目的层面元偏移能量的具体过程。
[0082] 首先，在子步骤S1011中基于炮点高斯射线束对地质模型做照明分析，确定入射 到每一目的层面元的全部高斯射线束的入射能量。
[0083] 具体而言，如图5a所示，统计所有入射到目的层的高斯射线束，从而分别统计覆 盖每一目的层面元的全部高斯射线束。叠加计算入射到每一目的层面元上的全部高斯射线 束的入射能量，得到每一目的层面元入射能量ECell (X，y)，形成图6a所示的多束射线的入 射照明结果。
[0084] 接下来，执行子步骤S1012,计算经由目的层反射之后，由每一地表检波器接收的 覆盖检波器所在网格范围的高斯射线束能量。即计算检波器的接收能量，也就是检波器接 收到来自炮点的能量。
[0085] 如图化所示，统计所有通过目的层反射至地表的高斯射线束，分别统计覆盖每一 地表检波器所在网格范围的高斯射线束。叠加计算由每一地表检波器接收的全部高斯射线 束的入射能量，得到每一地表检波器的接收能量邸ecv(x，y)，形成图化所示的多束射线的 接收照明结果。
[0086] 随后，在步骤S1013中根据入射到一个目的层面元的并且由全部地表检波器接收 的高斯线束能量确定目的层面元偏移能量，即通过面元反射且被地表检波器接收到的可用 于实现偏移成像的有效能量。
[0087] 如图5c中所示，面元偏移能量是指震源激发的地震波通过目的层面元反射后，被 观测系统中多个检波器接收到的能量总和，代表了真正可W用于成像地下面元的有效接收 能量。
[0088] 具体来说，首先计算每一目的层面元入射能量占炮点射线束总入射能量的比值 PCelia，j)。
[0089] 统计每一地表检波器的接收能量，计算所有地表检波器接收的由一个目的层面元 反射的来自每一炮点的高斯射线束能量，然后叠加计算由一个目的层面元反射的来自所有 炮点的高斯射线束能量，获得一个目的层面元的有效反射能量。
[0090] 最后，计算该目的层面元的有效反射能量与该目的层面元的入射能量占炮点高斯 射线束总入射能量的比值的乘积，得到该目的层面元的面元偏移能量：
[0091]
。）
[0092] 其中，（i，j)表示目的层中的水平坐标，P为目的层面元反射到地面的射线束编 号，ER(p)为目的层面元内第P条射线束反射到达地面被所有检波器接收能量叠加总和，N 为目的层面元反射到地面的射线束的总数，
表示目的层面元的有效反射能量。
[0093] ER(p)的数值可W根据表达式（1)至（6)来计算，再根据表达式（7)可获得图6c 所示的多束射线的偏移照明结果。送样W来，可W基于高斯束射线束的照明技术对Η维地 质模型进行分析，获得目的层的偏移能量分布。
[0094] 与现有技术相比，本实施例中提出目的层面元偏移能量的概念用于表征真正可用 于地面成像的地下面元有效接收能量，从而可W定量评价观测系统设计的优劣，更合理地 设计观测系统参数。同时面元偏移能量的提出使得地震采集与资料处理成像更加一体化， 为W成像为目的的地震采集设计提供依据。
[009引随后，执化步骤S102,统计目的层面兀偏移战垂分布情况W圈志目的层中战垂弱 区。
[0096] 具体的，首先选定地震勘探的目的层，即此次地震勘探主要想解决地下哪一目的 层的地质问题，然后针对选定的地下目 的层进行面元偏移能量统计，根据地质区域的实际 构造或者勘探目的设定能量弱区的阔值，圈定目的层中能量弱区。
[0097] 接下来，在步骤S103中根据各炮点对能量弱区的照明参考值确定目的层中能量 弱区对应的补炮区域。
[009引首先设定炮点照射到能量弱区的能量阔值。例如可W根据炮点照射能量占区域总 能量的60 %作为能量阔值。
[0099] 然后统计到达目的层能量弱区的高斯射线束对应的炮点，分别计算到达目的层能 量弱区的来自各炮点的高斯射线束能量与能量弱区总入射能量的比值，获得各炮点对能量 弱区的照明参考值。选择对能量弱区的照明参考值大于或等于预设阔值的炮点组成地表补 炮区域。
[0100] 最后，执行步骤S104,在补炮区域的原有炮点网格中等间隔设置补炮炮点。
[0101] 优选的，在地表补炮区域的原有的两个相邻炮点的中间位置设置一个补炮炮点， 送样设置加密炮点可使得原有炮点网格加密一倍。不限于此，也可在原有的两个相邻炮点 的中间位置设置两个补炮炮点，或者更多，可根据实施成本具体设定。如此W来，可W通过 在补炮区域加密炮点，使得目的层中照明能量弱的位置得到补偿，提高地震资料的信噪比。 引 工招连例
[010引下面W D邸地区为例，详细说明整个补炮过程。观测系统参数如下表1所示。
[0104] 表 1
[0105]
[0106] 通过过工区的Η条深度剖面建立一个Η维构造模型（如图7所示）。图8是补炮 前观测系统对目的层照明的分布，从图8可知，面元偏移能量分布差别很大，因此选中图中 黑圆圈所示位置进行补炮演示。图9为补炮后的目的层能量的分布，可知黑圆圈位置能量 得W补偿（补偿炮数为100炮，W目的层最大能量为标准）。图10a和图1化为地面补炮位 置和补炮数的放大显示，从图中可W看到，补炮的位置都是在原先炮点网格加密后的位置， 使得原有炮点网格成倍加密。对比图8和图9可W看出，本发明可W准确圈定目的层的能 量弱区，定量确定补炮数目，从而使得目的层中照明能量弱的位置得到补偿，提高地震资料 处理后的信噪比，改善地震采集资料品质。
[0107] 虽然本发明所公开的实施方式如上，但所述的内容只是为了便于理解本发明而采 用的实施方式，并非用W限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本 发明所公开的精神和范围的前提下，可W在实施的形式上及细节上作任何的修改与变化， 但本发明的专利保护范围，仍须W所附的权利要求书所界定的范围为准。
【主权项】
1. 一种基于高斯射线束目的层能量照明的补炮方法，其特征在于，包括以下步骤： 根据现有地质资料建立勘探区域的三维地质模型，基于炮点高斯射线束对地质模型做 照明分析来确定目的层面元偏移能量； 统计目的层面元偏移能量分布情况以圈定目的层中能量弱区； 根据各炮点对能量弱区的照明参考值确定目的层中能量弱区对应的地表补炮区域； 在地表补炮区域的原有炮点网格中等间隔设置补炮炮点。2. 根据权利要求1所述的方法，其特征在于，所述基于炮点高斯射线束对地质模型做 照明分析来确定目的层面元偏移能量的步骤包括： 基于炮点高斯射线束对地质模型做照明分析，确定入射到每一目的层面元的全部高斯 射线束的入射能量； 计算经由目的层反射之后，由每一地表检波器接收的覆盖检波器所在网格范围的高斯 射线束能量； 根据入射到一个目的层面元的并且由全部地表检波器接收的高斯线束能量确定目的 层面元偏移能量。3. 根据权利要求2所述的方法，其特征在于，所述基于炮点高斯射线束对地质模型做 照明分析，确定入射到每一目的层面元的全部高斯射线束的入射能量包括： 统计入射到目的层的全部高斯射线束； 分别统计覆盖每一目的层面元的高斯射线束； 叠加计算入射到每一目的层面元的高斯射线束的入射能量，得到每一目的层面元入射 能量。4. 根据权利要求3所述的方法，其特征在于，所述计算经由目的层反射之后，由每一地 表检波器接收的覆盖检波器所在网格范围的高斯射线束能量包括： 统计通过目的层反射至地表的所有高斯射线束； 分别统计覆盖每一地表检波器所在网格范围的高斯射线束； 叠加计算由每一地表检波器接收的全部高斯射线束的入射能量，得到每一地表检波器 的接收能量。5. 根据权利要求4所述的方法，其特征在于，所述根据入射到一个目的层面元的并且 由全部地表检波器接收的高斯线束能量确定目的层面元偏移能量包括： 计算每一目的层面元入射能量占炮点高斯射线束总入射能量的比值； 统计每一地表检波器的接收能量，计算所有地表检波器接收的由一个目的层面元反射 的来自每一炮点的高斯射线束能量； 叠加计算由目的层面元反射的来自所有炮点的高斯射线束能量，获得目的层面元的有 效反射能量； 由目的层面元的有效反射能量与目的层面元的入射能量占炮点高斯射线束总入射能 量的比值的乘积得到目的层面元的面元偏移能量。6. 根据权利要求1-5中任一项所述的方法，其特征在于，所述根据各炮点对能量弱区 的照明参考值确定目的层中能量弱区对应的补炮区域包括： 统计入射到目的层能量弱区的高斯射线束对应的炮点； 分别计算入射到目的层能量弱区的来自各炮点的高斯射线束能量与能量弱区总入射 能量的比值，获得各炮点对能量弱区的照明参考值； 选择对能量弱区的照明参考值大于或等于预设能量阈值的炮点组成地表补炮区域。7. 根据权利要求1-5中任一项所述的方法，其特征在于，所述在补炮区域的原有炮点 网格中等间隔设置补炮炮点为： 在地表补炮区域的原有的两个相邻炮点的中间位置设置一个补炮炮点。8. 根据权利要求1-5中任一项所述的方法，其特征在于，所述根据现有地质资料建立 勘探区域的三维地质模型包括： 利用现有地质资料分别构建中深层构造模型和地表模型； 设定每一地层的纵波速度、横波速度和密度，获得三维块状地质模型。9. 根据权利要求5所述的方法，其特征在于，所述目的层面元的面元偏移能量表示为：其中，（i，j)表示目的层中的水平坐标，P为目的层面元反射到地表的射线束编号， ER(P)为目的层面元内第p条射线束反射到达地面被所有检波器接收能量叠加总和，N为目 的层面元反射到地面的射线束的总数，PCelUi，j)表示目的层面元入射能量占炮点射线束 总入射能量的比值。10. 根据权利要求6所述的方法，其特征在于，所述预设能量阈值为60%。
【专利摘要】本发明公开了一种基于高斯射线束目的层能量照明的补炮方法，包括以下步骤：根据现有地质资料建立勘探区域的三维地质模型，基于炮点高斯射线束对地质模型做照明分析来确定目的层面元偏移能量；统计目的层面元偏移能量分布情况以圈定目的层中能量弱区；根据各炮点对能量弱区的照明参考值确定目的层中能量弱区对应的地表补炮区域；在地表补炮区域的原有炮点网格中等间隔设置补炮炮点。本发明可以准确圈定目的层的能量弱区，定量确定补炮数目，从而使得目的层中照明能量弱的位置得到补偿，提高地震资料处理后的信噪比，改善地震采集资料品质。
【IPC分类】G01V1/00
【公开号】CN105487106
【申请号】CN201410479944
【发明人】肖云飞, 殷厚成
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院
【公开日】2016年4月13日
【申请日】2014年9月18日
最新回复(0)