CN107831533A - 基于能量统计的采集方案检波器排列长度定量分析方法 - Google Patents
基于能量统计的采集方案检波器排列长度定量分析方法 Download PDFInfo
- Publication number
- CN107831533A CN107831533A CN201711143018.7A CN201711143018A CN107831533A CN 107831533 A CN107831533 A CN 107831533A CN 201711143018 A CN201711143018 A CN 201711143018A CN 107831533 A CN107831533 A CN 107831533A
- Authority
- CN
- China
- Prior art keywords
- energy
- geophone
- energy value
- focal point
- acquisition scheme
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/16—Receiving elements for seismic signals; Arrangements or adaptations of receiving elements
- G01V1/20—Arrangements of receiving elements, e.g. geophone pattern
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种基于能量统计的采集方案检波器排列长度定量分析方法。所述方法包括:选取勘探目的层,布设震源点和检波点;从震源点位置进行波场能量向下传播延拓,经目的层反射,形成二次震源,将此时目的层的能量作为二次震源的能量值向地表进行反向传播延拓,记录检波点接收到的能量值,记为炮检对能量值;得到每个炮检对的能量值,绘制炮检对能量值图,以图的对角线为零偏移距线,根据炮检对能量值的分布范围,得平行于零偏移距线的参考线,该参考线与检波点线的交点得到的值则是针对当前目的层的合理的检波器排列长度。能够实现了对检波器排列长度的定量评估,进而优化采集观测系统,节约施工成本,提高采集资料质量。
Description
技术领域
本发明石油地震勘探采集方案优化设计领域,更具体地讲,涉及一种一种采集方案检波器排列长度定量分析方法。
背景技术
西部复杂山地地区是是勘探程度较低的油气富集地区。由于山地地区地表地震条件和地下地质条件均很复杂,对地震勘探设备和地震勘探技术的要求高,因此许多常规的地震勘探技术无法适用于山地地区。随着科学技术的发展和国民生活水平的不断提高,人们对能源的需求量越来越大,同时在平原地区资源日益减少的形势下,开发利用山地油气资源是解决能源接替的途径之一。
地震观测系统设计和采集参数优化是整个石油勘探中非常重要的环节,是地震资料品质好坏的先决条件,也是地震资料处理和解释的基础。常规地震观测系统设计的覆盖次数只依赖于炮点和检波点排列范围,而与地下结构、目的层深度无关。然而在构造复杂的地区,覆盖次数除了与炮点和检波点排列有关外,还与目的层的深度以及上覆地质结构有关,此时,常规的基于水平层状假设的针对CMP点进行设计的观测系统精度降低。由于目的层各CRP点上的覆盖次数和地震波照射能量共同决定了该点的成像质量,但是复杂的上覆地质结构会造成下伏地层地震照明阴影区,从而引起下伏勘探目的层照明强度显著下降,最终使得这些目的层界面成像的效果不理想。因此研究基于波动方程的观测系统优化设计方法势在必行。
基于能量统计的采集方案分析方法目前主要是照明度分析,而地震波照明又可以分为基于射线追踪的地震波照明方法和基于波动方程的地震波照明方法两类。射线追踪地震波照明方法目前应用最为普遍,该方法简单直观、计算效率高,但由于射线理论本身的缺陷(如高频近似、射线盲区等)和复杂介质中的奇异性问题,基于射线的照明分析方法在复杂构造区会出现很大的误差,照明精度较低。而基于波动方程的地震波照明方法能有效地克服射线照明分析方法在复杂构造区的不足,尤其是利用波动方程中的单程波方程的地下照明分析是近几年的一个研究热点,它克服了射线照明方法的缺陷,能适应强横向变速介质,使得目标照明分析更加合理、准确。
基于波动方程的照明度计算均是分析能量在地下介质分布的情况,而没有从目标反射层和检波点接收能量的角度去分析观测系统的优劣。使得照明度分析在指导观测系统设计更多的是定性的角度,无法从定量的角度指导观测系统的优化设计,最大炮检距分析是二维观测系统设计中的重要参数,实际生产正急需这种定量分析。
发明内容
针对现有技术中存在的问题,本发明的目的在于解决现有技术存在的上述不足中的至少一项。例如,本发明的目的在于提出一种针对山地复杂构造的采集方案检波器排列长度的定量评估,进而优化采集观测系统,节约施工成本,提高采集资料质量。
为了实现上述目的,本发明提供了一种基于能量统计的采集方案检波器排列长度定量分析方法。所述方法包括:选取勘探目的层,在模型的地表布设震源点和检波点;从每个震源点位置进行波场能量向下传播延拓,经过目的层的反射,在目的层形成二次震源,将此时目的层的能量作为二次震源的能量值向地表进行反向传播延拓,记录每个检波点接收到的对应于每个震源点的能量值;根据每个检波点接收到的对应于每个震源点的能量值绘制出当前模型的能量值分布图,所述能量值分布图的对角线为零偏移距线,根据图中能量值的分布范围,得到一条平行于零偏移距线的参考线,该参考线与检波点线的交点得到的值则是针对当前目的层的合理的检波器排列长度。
在本发明的一个示例性实施例中,所述布设震源点和检波点可以是在模型的地表上的每个点都有震源点和检波点,即为满布设。
在本发明的一个示例性实施例中,所述延拓可以利用高精度的有限差分算子。
在本发明的一个示例性实施例中,所述延拓的计算步骤可以包括:
式1为在频率域中,密度固定的声波方程。
式1为:(k2+Δ)U(rs,r,ω)=-f(ω)δ(r-rs)。
其中:k=ω/c0(r)为波数,Δ为拉普拉斯算子,c0(r)为点r处的背景场速度值,f(ω)为震源函数,rs为震源点位置,w为频率,r为地下介质中某一点的位置,U为波场函数,δ为狄拉克函数。
式2为通过向下传播延拓和反向传播延拓后,在接收点处的波场表达式。
式2为:P(rs,rg,ω)=2k2∫V(r)G(r,rs;ω)m(r)G(r,rg;ω)d3r。
其中:m(r)=δc(r)/c0(r)是在地下r=(x,y,z)处速度模型扰动,G(r,rs;ω)和G(r,rg;ω)分别是频率域中震源和检波点位置的格林函数,V(r)为目标反射层位置,rg为检波点位置。
式3为任一个检波点接收到的对应于任一个震源点的能量值。
式3为:Psr(rs,rg,ω)=|P(rs,rg,ω)|。
在本发明的一个示例性实施例中,所述分布范围可以包括炮检对能量值的80%至95%。
在本发明的一个示例性实施例中,所述合理的检波器排列长度可以是山地复杂构造的采集方案检波器排列长度。
与现有技术相比,本发明的有益效果包括:能够克服常规的基于波动方程的地震采集方案评价方法无法定量分析采集方案检波器的合理摆放个数的问题。从基于能量统计的角度实现了对检波器排列长度的定量评估,进而优化采集观测系统,节约施工成本,提高采集资料质量。对提高山地地震采集设计效率及优化提供了保障,具有重要应用价值。
附图说明
图1示出了根据本发明的一个示例性实施例的基于能量统计的采集方案检波器排列长度定量分析方法的起伏地表的示意图。
图2示出了根据本发明的一个示例性实施例的基于能量统计的采集方案检波器排列长度定量分析方法的地震波传播和反射过程的示意图。
图3示出了根据本发明的一个示例性实施例的基于能量统计的采集方案检波器排列长度定量分析方法的所有震源和检波点的能量分布图。
具体实施方式
在下文中,将结合示例性实施例和附图来详细说明本发明的基于能量统计的采集方案检波器排列长度定量分析方法,即基于能量统计的方式来定量分析采集方案中检波器排列长度的方法。
根据本发明示例性实施例的基于能量统计的采集方案检波器排列长度定量分析方法可以包括如下步骤:
步骤(1),选取目标反射层。对于地震勘探,有主要的勘探目的层,施工设计时采集参数应该主要侧重于适应主要勘探目的层,通过选取目标反射层如图1所示,选择地震波的反射点如图2所示。
步骤(2),设计采集方案。由于本技术论证的是采集方案的检波器排列长度,因此将震源点和检波器在模型的地表全部布设,即布设满模型的整个表面。
步骤(3),对选定的一个勘探目的层位置,从震源点位置利用高精度有限差分算子进行波场能量向下传播延拓,经过勘探目的层反射,在目的层形成二次震源,将此时勘探目的层的能量作为二次震源能量值,利用高精度有限差分算子向地表进行反向传播延拓,然后记录检波器接收到的能量值,需要说明的是对有限差分算子的选择不限于此。对所有的震源点进行同样操作,因此就可以得到整个观测系统的炮检对在检波点处的能量分布图(也可称为能量值分布图)。
计算方法如下:
在频率域中,密度固定的声波方程可以表示为:
(k2+Δ)U(rs,r,ω)=-f(ω)δ(r-rs) (1)
其中:k=ω/c0(r)为波数,Δ为拉普拉斯算子,c0(r)为点r处的背景场速度值,f(ω)为震源函数,rs为震源点位置,w为频率,r为地下介质中某一点的位置,U为波场函数,δ为狄拉克函数。通过正向延拓和反向延拓,在接收点处的波场可以表示为:
P(rs,rg,ω)=2k2∫V(r)G(r,rs;ω)m(r)G(r,rg;ω)d3r (2)
其中:m(r)=δc(r)/c0(r)是在地下r=(x,y,z)处速度模型扰动,G(r,rs;ω)和G(r,rg;ω)分别是频率域中震源和检波点位置的格林函数,V(r)为目标反射层位置,rg为检波点位置。对于一个单独的炮检对,公式(2)的振幅值可以认为是这个炮点所对应的那个检波点的响应,即为炮检对能量值,公式如下:
Psr(rs,rg,ω)=|P(rs,rg,ω)| (3)
炮检对能见度由炮点所对应的所有检波点的能量值组合而成,计算公式如下:
步骤(4),得到所有震源在检波点处的能量分布图。如图3所示,图中横纵坐标的单位是km或m,分别对应炮点和检波点的坐标,图中箭头所指的A区域(梭形黑色区域)为能量最高的区域,B区域(A区域外的环形白色区域)的能量相对于A区域减弱,C区域(B区域外的环形黑色区域)的能量相对于B区域减弱,D区域(C区域外的环形白色区域及向外颜色逐渐变黑的环形区域)的能量相对于C区域减弱,从A区域到B区域到C区域到D区域能量逐渐减弱,D区域以外的黑色背景区域为能量零值区域。图形的对角线则为零偏移距参考线,选择一条与零偏移距参考线平行,且包括了炮检对主要能量的一个范围(一般选择90%或者85%)的线(如图3中的虚线),这条线与检波点线的交点(如图3中的两个黑色原点E和F)得到的值则是针对当前目标层的合理的检波器排列长度。
如图3所示为利用该技术可实现对山地复杂构造的采集方案检波器排列长度的定量评估,进而优化采集观测系统,节约施工成本,提高采集资料质量。
尽管上面已经结合附图和示例性实施例描述了本发明,但是本领域普通技术人员应该清楚,在不脱离权利要求的精神和范围的情况下,可以对上述实施例进行各种修改。
Claims (6)
1.一种基于能量统计的采集方案检波器排列长度定量分析方法,其特征在于,所述方法包括以下步骤:
选取勘探目的层,在模型的地表布设震源点和检波点;
从每个震源点位置进行波场能量向下传播延拓,经过目的层的反射,在目的层形成二次震源,将此时目的层的能量作为二次震源的能量值向地表进行反向传播延拓,记录每个检波点接收到的对应于每个震源点的能量值;
根据每个检波点接收到的对应于每个震源点的能量值绘制出当前模型的能量值分布图,所述能量值分布图的对角线为零偏移距线,根据图中能量值的分布范围,得到一条平行于零偏移距线的参考线,该参考线与检波点线的交点得到的值则是针对当前目的层的合理的检波器排列长度。
2.根据权利要求1所述的基于能量统计的采集方案检波器排列长度定量分析方法,其特征在于,所述布设震源点和检波点是在模型的地表上的每个点都设置震源点和检波点。
3.根据权利要求1所述的基于能量统计的采集方案检波器排列长度定量分析方法,其特征在于,所述延拓利用了高精度的有限差分算子。
4.根据权利要求1所述的基于能量统计的采集方案检波器排列长度定量分析方法,其特征在于,所述延拓的计算步骤包括:
式1为在频率域中,密度固定的声波方程,
式1为:(k2+Δ)U(rs,r,ω)=-f(ω)δ(r-rs),
其中:k=ω/c0(r)为波数,Δ为拉普拉斯算子,c0(r)为点r处的背景场速度值,f(ω)为震源函数,rs为震源点位置,w为频率,r为地下介质中某一点的位置,U为波场函数,δ为狄拉克δ函数;
式2为通过向下传播延拓和反向传播延拓后,在接收点处的波场表达式,
式2为:P(rs,rg,ω)=2k2∫V(r)G(r,rs;ω)m(r)G(r,rg;ω)d3r,
其中:m(r)=δc(r)/c0(r)是在地下r=(x,y,z)处速度模型扰动,G(r,rs;ω)和G(r,rg;ω)分别是频率域中震源和检波点位置的格林函数,V(r)为目标反射层位置,rg为检波点位置;
式3为任一个检波点接收到的对应于任一个震源点的能量值,
式3为:Psr(rs,rg,ω)=|P(rs,rg,ω)|。
5.根据权利要求1所述的基于能量统计的采集方案检波器排列长度定量分析方法,其特征在于,所述分布范围包括了炮检对能量值的80%至95%。
6.根据权利要求1所述的基于能量统计的采集方案检波器排列长度定量分析方法,其特征在于,所述合理的检波器排列长度为山地复杂构造的采集方案检波器排列长度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711143018.7A CN107831533A (zh) | 2017-11-17 | 2017-11-17 | 基于能量统计的采集方案检波器排列长度定量分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711143018.7A CN107831533A (zh) | 2017-11-17 | 2017-11-17 | 基于能量统计的采集方案检波器排列长度定量分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107831533A true CN107831533A (zh) | 2018-03-23 |
Family
ID=61651890
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711143018.7A Pending CN107831533A (zh) | 2017-11-17 | 2017-11-17 | 基于能量统计的采集方案检波器排列长度定量分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107831533A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110988995A (zh) * | 2019-12-13 | 2020-04-10 | 中国海洋石油集团有限公司 | 一种基于长基距的采集参数确定方法和装置 |
CN113031061A (zh) * | 2021-04-02 | 2021-06-25 | 中油奥博(成都)科技有限公司 | 一种气化腔边界识别方法 |
WO2023201865A1 (zh) * | 2022-04-23 | 2023-10-26 | 中国海洋大学 | 一种观测系统排列长度选择方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103645460A (zh) * | 2013-11-21 | 2014-03-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于目标层可视性分析的激发位置确定方法 |
CN103645501A (zh) * | 2013-11-21 | 2014-03-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 最大炮检距确定方法 |
CN104155703A (zh) * | 2014-09-02 | 2014-11-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 评价三维观测系统的方法和设备 |
-
2017
- 2017-11-17 CN CN201711143018.7A patent/CN107831533A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103645460A (zh) * | 2013-11-21 | 2014-03-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于目标层可视性分析的激发位置确定方法 |
CN103645501A (zh) * | 2013-11-21 | 2014-03-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 最大炮检距确定方法 |
CN104155703A (zh) * | 2014-09-02 | 2014-11-19 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 评价三维观测系统的方法和设备 |
Non-Patent Citations (1)
Title |
---|
徐峰 等: "基于目的层反射能量二次导数的炮检距设计", 《西南石油大学学报(自然科学版)》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110988995A (zh) * | 2019-12-13 | 2020-04-10 | 中国海洋石油集团有限公司 | 一种基于长基距的采集参数确定方法和装置 |
CN110988995B (zh) * | 2019-12-13 | 2022-07-08 | 中国海洋石油集团有限公司 | 一种基于长基距的采集参数确定方法和装置 |
CN113031061A (zh) * | 2021-04-02 | 2021-06-25 | 中油奥博(成都)科技有限公司 | 一种气化腔边界识别方法 |
CN113031061B (zh) * | 2021-04-02 | 2022-09-27 | 中油奥博(成都)科技有限公司 | 一种气化腔边界识别方法 |
WO2023201865A1 (zh) * | 2022-04-23 | 2023-10-26 | 中国海洋大学 | 一种观测系统排列长度选择方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101625417B (zh) | 一种优化垂直地震剖面观测系统设计的方法 | |
CN106094032B (zh) | 一种构建地层速度模型的方法 | |
CN103454680B (zh) | Walk-away VSP观测系统垂向覆盖次数的计算方法 | |
CN107831533A (zh) | 基于能量统计的采集方案检波器排列长度定量分析方法 | |
CN103605157B (zh) | 衰减近地表散射波的方法 | |
CN108983284A (zh) | 一种适用于海上斜缆数据的f-p域鬼波压制方法 | |
CN106501841B (zh) | 一种二维宽线地震观测系统的优化设计方法及装置 | |
CN103645460B (zh) | 基于目标层可视性分析的激发位置确定方法 | |
CN106154315A (zh) | 一种对称均匀的高密度地震采集方法 | |
CN105242311B (zh) | 预估初至拾取的时窗的方法和装置 | |
CN107678057B (zh) | 三维地震中确定或优化检波器或设备投入量及排布的方法 | |
CN102540254B (zh) | 一种基于照明分析确定有效激发区域的方法 | |
CN103645501A (zh) | 最大炮检距确定方法 | |
CN104977615B (zh) | 一种基于模型统计拾取的深水obc资料多次波压制方法 | |
CN110501751B (zh) | 一种基于多分量梯度数据联合和深度加权的相关成像方法 | |
CN110907995A (zh) | 井中vsp地震数据的逆时偏移方法及装置 | |
CN103605158A (zh) | 一种最大炮检距确定方法及装置 | |
CN107957594A (zh) | 地震数据的椭圆校正方法、动校正方法及速度谱计算方法 | |
CN103592684A (zh) | 一种保持空间属性信息的海量地震数据压缩方法及装置 | |
CN110579799B (zh) | 一种等旅行时间间隔的地震采集观测方法及系统 | |
CN104267434A (zh) | 一种三维多分量地震观测系统炮检距分布获取方法及装置 | |
CN104155703A (zh) | 评价三维观测系统的方法和设备 | |
CN106842334A (zh) | 一种电磁勘探方法和装置 | |
CN111158050A (zh) | 数据采集系统、方法及隧道地震波超前预报方法 | |
CN109655879B (zh) | 目的层信噪比一致性能量照明的观测系统优化方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180323 |