CN112540407B - 一种叠前深度偏移各向异性场建立方法 - Google Patents

一种叠前深度偏移各向异性场建立方法 Download PDF

Info

Publication number
CN112540407B
CN112540407B CN202011389428.1A CN202011389428A CN112540407B CN 112540407 B CN112540407 B CN 112540407B CN 202011389428 A CN202011389428 A CN 202011389428A CN 112540407 B CN112540407 B CN 112540407B
Authority
CN
China
Prior art keywords
depth
target layer
seismic
depth migration
borehole
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.)
Active
Application number
CN202011389428.1A
Other languages
English (en)
Other versions
CN112540407A (zh
Inventor
王建军
杜百灵
郎玉泉
牛鹏程
孟建盛
窦国兴
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hydrogeology Bureau of China National Administration of Coal Geology
Original Assignee
Hydrogeology Bureau of China National Administration of Coal Geology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hydrogeology Bureau of China National Administration of Coal Geology filed Critical Hydrogeology Bureau of China National Administration of Coal Geology
Priority to CN202011389428.1A priority Critical patent/CN112540407B/zh
Publication of CN112540407A publication Critical patent/CN112540407A/zh
Application granted granted Critical
Publication of CN112540407B publication Critical patent/CN112540407B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/301—Analysis for determining seismic cross-sections or geostructures
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/51—Migration
    • G01V2210/512—Pre-stack
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
    • G01V2210/632—Amplitude variation versus offset or angle of incidence [AVA, AVO, AVI]

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种叠前深度偏移各向异性场建立方法，包括以下步骤：对三维地震各向同性叠前深度偏移成果进行目的层层位解释，得到工区内每个点目的层深度D；在三维地震各向同性叠前深度偏移成果剖面上读取钻孔处目的层的深度D1；根据D1和钻孔目的层的深度D2计算钻孔处各向异性值δ1；对D1和δ1进行线性回归拟合，求取D1和δ1的关系式Y；利用关系式Y和目的层深度D计算工区内每个点的各向异性值δ；利用各向异性值δ建立各向异性场进行三维地震各向异性叠前深度偏移；在三维地震各向异性叠前深度偏移成果剖面上读取钻孔处目的层的深度D2，线性回归拟合得到的关系式Y可用于目的层深度与各向异性值的转化计算。

Description

一种叠前深度偏移各向异性场建立方法
技术领域
本发明专利涉及一种叠前深度偏移各向异性场建立方法，适用于对深度精度要求高但是钻孔数据少的工区建立各向异性场进行三维地震各向异性叠前深度偏移使用。
背景技术
目前叠前深度偏移各向异性场建立只是依靠钻孔处的各向异性值建立各向异性场，方法就是求取每个钻孔处的各向异性值，然后进行插值生成各向异性场，利用各向异性场进行各向异性叠前深度偏移。
问题和缺点：以往方法在钻孔少的工区，通过插值得到的钻孔之间和边界的个向异性值，往往与真实的各向异性值存在较大的误差，造成叠前深度偏移深度精度低。
发明内容
为了有效解决背景技术的问题，有必要提供一种可提高叠前深度偏移深度精度的叠前深度偏移各向异性场建立方法。
本发明的目的是通过以下技术方案实现的：
一种叠前深度偏移各向异性场建立方法，其包括以下步骤：
S1、对三维地震各向同性叠前深度偏移成果进行目的层层位解释，得到工区内每个点目的层深度D；
S2、在三维地震各向同性叠前深度偏移成果剖面上读取钻孔处目的层的深度D1；
S3、根据D1和钻孔目的层的深度D2计算钻孔处各向异性值δ1；
S4、对D1和δ1进行线性回归拟合，求取D1和δ1的关系式Y；
S5、利用关系式Y和目的层深度D计算工区内每个点的各向异性值δ；
S6、利用各向异性值δ建立各向异性场进行三维地震各向异性叠前深度偏移；
S7、在三维地震各向异性叠前深度偏移成果剖面上读取钻孔处目的层的深度D2，并和D1对比，如果两者误差小，则线性回归拟合得到的关系式Y可用于目的层深度与各向异性值的转化计算；如果误差大，对关系式Y进行调整，直到误差满足要求。
所述S3所述的钻孔处各向异性值δ1的计算方法为：
Figure BDA0002811769390000021
所述S4具体为：Y的计算公式为：Y＝aX2+bX+c,Y为各向异性值δ1，X为在三维地震各向同性叠前深度偏移成果剖面上的深度D1，a、b、c为通过线性回归拟合得到的系数项。
所述S5中工区内每个点的各向异性值δ的计算方法为：
δ＝a*D2+b*D+c。
相对于现有技术，本发明的优点在于：
本发明专利主要以一种基于目的层深度建立各向异性场的方法，可以提高各向异性场的准确度，提高叠前深度偏移深度精度，提高勘探目的层的深度预测精度。
附图说明
图1为本发明实施例中对利用工区内的27口钻孔处的D1和δ1值，进行线性回归拟合，求取D1和δ1的关系式Y；
图2为不同的叠前深度偏移的示意图。
具体实施方式
为了使本发明的方法流程、特征、成效易于了解，下面结合具体实施方式，进一步阐述本发明。
本实施例提供一种叠前深度偏移各向异性场建立方法，其包括以下步骤：
1、收集工区三维地震各向同性叠前深度偏移数据，在解释工作站上对三维地震各向同性叠前深度偏移数据剖面上对目的层进行平面间隔20m*20m解释，得到工区内平面间隔20m*20m的目的层深度D；
2、收集工区内的钻孔数据，在过钻孔的三维地震各向同性叠前深度偏移成果剖面上读取每个钻孔处目的层的深度D1；
3、利用各向异性值计算公式，根据D1和钻孔目的层的深度D2计算钻孔处的各向异性值δ1；
Figure BDA0002811769390000031
4、对钻孔处的D1和δ1值，进行线性回归拟合，求取D1和δ1的关系式Y：
Y＝aX2+bX+c,Y为各向异性值δ1，X为在三维地震各向同性叠前深度偏移成果剖面上的深度D1，a、b、c为通过线性回归拟合得到的系数项；通过所述a、b、c，可以算出关系式Y的相关系数R，当相关系数R大于0.8的时候，认为系数项a、b、c满足要求；
相关系数R的计算公式为：
Figure BDA0002811769390000032
其中Cov(X,Y)为X与Y的协方差，Var[X]为X的方差，Var[Y]为Y的方差。
5、利用关系式Y和目的层深度D计算工区内每个点的各向异性值δ；
δ＝a*D2+b*D+c；
6、利用各向异性值δ建立各向异性场进行三维地震各向异性叠前深度偏移：此步骤属于现有技术，利用已有的插值软件将上一步之前计算的目的层的各向异性值δ插值成一个三维的各向异性值δ体，完成利用各向异性值δ建立各向异性场；
7、在三维地震各向异性叠前深度偏移成果剖面上读取钻孔处目的层的深度D2，将D2和D1对比，如果两者误差小，则线性回归拟合得到的关系式Y可用于目的层深度与各向异性值的转化计算；如果D2和D1两者误差大，则需对关系式Y进行调整，具体可以求取三阶Y＝aX3+bX2+cX+d或者四阶的关系式Y＝aX4+bX3+c X2+dX+e，用重新求取的各向异性值δ建立各向异性场进行三维地震各向异性叠前深度偏移，直到误差满足要求。
图1是对利用工区内的27口钻孔处的D1和δ1值，进行线性回归拟合，求取D1和δ1的关系式Y，从图中可以看出，D1和δ1有很好的相关性，可以利用深度D1来计算各向异性值δ1。
表1为井震误差统计表，各向同性叠前深度偏移后，井震误差最大为47.67m；以往方法建立各向异性场进行各向异性叠前深度偏移后，井震误差最大为25.1m，误差较各向同性叠前深度偏移有所减小，但依然较大；完全满足了煤炭等领域对解释成果精度高的要求；新方法建立各向异性场进行各向异性叠前深度偏移后，误差较以往方法建立各向异性场进行各向异性叠前深度偏移有所减小，井震误差最大为2.49m，完全满足了煤炭等领域对解释成果精度高的要求。
附表1井震误差统计表
Figure BDA0002811769390000041
Figure BDA0002811769390000051
Figure BDA0002811769390000061
Figure BDA0002811769390000071
Figure BDA0002811769390000081
图2显示不同的叠前深度偏移，图2中上部的图是三维地震各向同性的叠前深度偏移连井剖面，剖面深度与钻孔的深度差值大；中间的图是以往根据钻孔插值生成各向异性场的三维地震各向异性的叠前深度偏移连井剖面，剖面深度与钻孔的深度差值较上变的图有所减小，但依然较大；下部的图是利用关系式和目的层深度求取各向异性值，建立各向异性场进行性进行三维地震各向异性的叠前深度偏移连井剖面，得到速度2，剖面深度与钻孔的深度差值很小，深度精度较中间的图有明显提高。

Claims (3)

1.一种叠前深度偏移各向异性场建立方法，其特征在于，包括以下步骤：
S1、对三维地震各向同性叠前深度偏移成果进行目的层层位解释，得到工区内每个点目的层深度D；
S2、在三维地震各向同性叠前深度偏移成果剖面上读取钻孔处目的层的深度D1；
S3、根据D1和钻孔目的层的深度D2计算钻孔处各向异性值δ1；
S4、对D1和δ1进行线性回归拟合，求取D1和δ1的关系式Y：Y的计算公式为：Y＝aX2+bX+c,Y为各向异性值δ1，X为在三维地震各向同性叠前深度偏移成果剖面上的深度D1，a、b、c为通过线性回归拟合得到的系数项；如果D2和D1两者之间的误差较大，则需对关系式Y进行调整：求取三阶Y＝aX3+bX2+cX+d或者四阶的关系式Y＝aX4+bX3+cX2+dX+e，用重新求取的各向异性值δ建立各向异性场进行三维地震各向异性叠前深度偏移，直到误差满足要求；
S5、利用关系式Y和目的层深度D计算工区内每个点的各向异性值δ，δ的计算方法为：
δ＝a*D2+b*D+c；
S6、利用各向异性值δ建立各向异性场进行三维地震各向异性叠前深度偏移；
S7、在三维地震各向异性叠前深度偏移成果剖面上读取钻孔处目的层的深度D2，将D2和D1对比，如果两者误差小，则线性回归拟合得到的关系式Y可用于目的层深度与各向异性值的转化计算；如果误差大，对关系式Y进行调整，直到误差满足要求。
2.根据权利要求1所述的叠前深度偏移各向异性场建立方法，其特征在于，所述S3所述的钻孔处各向异性值δ1的计算方法为：
Figure FDA0004137225340000011
3.根据权利要求1所述的叠前深度偏移各向异性场建立方法，其特征在于，
所述S4还包括：通过所述a、b、c，算出关系式Y的相关系数R，当相关系数R大于0.8的时候，认为系数项a、b、c满足要求；相关系数R的计算公式为：
Figure FDA0004137225340000021
其中Cov(X,Y)为X与Y的协方差，Var[X]为X的方差，Var[Y]为Y的方差。
CN202011389428.1A 2020-12-01 2020-12-01 一种叠前深度偏移各向异性场建立方法 Active CN112540407B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011389428.1A CN112540407B (zh) 2020-12-01 2020-12-01 一种叠前深度偏移各向异性场建立方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011389428.1A CN112540407B (zh) 2020-12-01 2020-12-01 一种叠前深度偏移各向异性场建立方法

Publications (2)

Publication Number Publication Date
CN112540407A CN112540407A (zh) 2021-03-23
CN112540407B true CN112540407B (zh) 2023-04-25

Family

ID=75017034

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011389428.1A Active CN112540407B (zh) 2020-12-01 2020-12-01 一种叠前深度偏移各向异性场建立方法

Country Status (1)

Country Link
CN (1) CN112540407B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114296137B (zh) * 2021-12-21 2024-04-26 中国煤炭地质总局地球物理勘探研究院 一种地震q场建立方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6785612B1 (en) * 2003-05-29 2004-08-31 Pgs Americas, Inc. Seismic velocity update for anisotropic depth migration
CN102025531A (zh) * 2010-08-16 2011-04-20 北京亿阳信通软件研究院有限公司 一种性能数据的填补方法及其装置
CN103149588A (zh) * 2013-02-20 2013-06-12 中国石油天然气股份有限公司 一种利用井震标定计算vti各向异性参数的方法及系统
CN109581499A (zh) * 2018-11-08 2019-04-05 成都捷科思石油天然气技术发展有限公司 一种利用各向异性叠前深度偏移生成构造图的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100135115A1 (en) * 2008-12-03 2010-06-03 Chevron U.S.A. Inc. Multiple anisotropic parameter inversion for a tti earth model

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6785612B1 (en) * 2003-05-29 2004-08-31 Pgs Americas, Inc. Seismic velocity update for anisotropic depth migration
CN102025531A (zh) * 2010-08-16 2011-04-20 北京亿阳信通软件研究院有限公司 一种性能数据的填补方法及其装置
CN103149588A (zh) * 2013-02-20 2013-06-12 中国石油天然气股份有限公司 一种利用井震标定计算vti各向异性参数的方法及系统
CN109581499A (zh) * 2018-11-08 2019-04-05 成都捷科思石油天然气技术发展有限公司 一种利用各向异性叠前深度偏移生成构造图的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
任婷 等.基于趋势速度约束的各向异性深度偏移建模与应用.《 CPS/SEG北京2018国际地球物理会议暨展览电子论文集》.2018,479-482. *

Also Published As

Publication number Publication date
CN112540407A (zh) 2021-03-23

Similar Documents

Publication Publication Date Title
CN104730579B (zh) 一种基于表层横波速度反演的纵横波联合静校正方法
CN109541685B (zh) 一种河道砂体识别方法
CN106646613B (zh) 深度域多尺度井控建模与成像联合处理方法
WO2010118624A1 (zh) 地震波速度场构建中井约束横向可变h-v曲线构建方法
CN113031068B (zh) 一种基于反射系数精确式的基追踪叠前地震反演方法
CN108254780A (zh) 一种微地震定位及各向异性速度结构层析成像方法
CN109188522B (zh) 速度场构建方法及装置
CN111175825B (zh) 深度域速度建模方法
CN116184493B (zh) 一种面向水平井开发的井震联合高精度速度建场方法
CN112540407B (zh) 一种叠前深度偏移各向异性场建立方法
CN112925022B (zh) 页岩vti介质各向异性参数的预测方法
CN111239818A (zh) 一种基于三维倾角属性体校正的古地貌分析方法
CN106324678B (zh) 一种基于测井数据约束的全波形反演方法及系统
CN114428275A (zh) 复杂地表的二维拟三维层析建模方法
CN112529981B (zh) 一种深度域煤层底板等高线成图的方法
CN111190224B (zh) 一种基于三维地震波反向照明的动态采样全波形反演系统及方法
CN107576987A (zh) 一种碳酸盐串珠储层水平井轨迹的设计方法
CN111596348B (zh) 基于声波时差测井资料的地层层速度校正方法及装置
CN108508479B (zh) 一种空地井立体重磁数据协同目标位置反演方法
CN105095634B (zh) 一种基于地质块体的偏移速度模型建立方法
CN115561831A (zh) 一种地震合成记录制作方法
CN116106970A (zh) 山前带砾石层近地表速度的建模方法
CN103984021A (zh) 一种基于系列速度误差剖面的偏移速度分析方法
CN113534257A (zh) 深度域初始速度模型建立方法及处理终端
CN119148224A (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
GR01 Patent grant
GR01 Patent grant