CN108363101A - 一种斜井井间地震高斯束叠前深度偏移成像方法 - Google Patents
一种斜井井间地震高斯束叠前深度偏移成像方法 Download PDFInfo
- Publication number
- CN108363101A CN108363101A CN201810106904.0A CN201810106904A CN108363101A CN 108363101 A CN108363101 A CN 108363101A CN 201810106904 A CN201810106904 A CN 201810106904A CN 108363101 A CN108363101 A CN 108363101A
- Authority
- CN
- China
- Prior art keywords
- ray
- imaging
- seismic
- gaussian beam
- point
- 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.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 77
- 238000013508 migration Methods 0.000 title claims abstract description 37
- 230000005012 migration Effects 0.000 title claims abstract description 37
- 230000001186 cumulative effect Effects 0.000 claims abstract description 5
- 230000000644 propagated effect Effects 0.000 claims abstract description 4
- 230000005284 excitation Effects 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000005755 formation reaction Methods 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 5
- 238000000354 decomposition reaction Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims description 2
- 238000010304 firing Methods 0.000 claims description 2
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 238000000034 method Methods 0.000 description 12
- 238000004422 calculation algorithm Methods 0.000 description 4
- 238000010276 construction Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 239000003518 caustics Substances 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004587 chromatography analysis Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
- G01V2210/512—Pre-stack
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
本发明公开了一种斜井井间地震高斯束叠前深度偏移成像方法,包括以下步骤:(1)读入井间地震反射波场记录等文件;(2)将地震波场分解成一系列的局部平面部,得到相应的数据体;(3)从炮点进行射线追踪,利用程函方程计算任意射线传播的时间τ;(4)获取中心射线,进行高斯束动力学追踪,使用龙格库塔方程求解动力学射线追踪方程;(5)计算射线束覆盖范围内网格节点的地震波振幅;(6)从检波点进行运动学和动力学射线追踪,计算并存储每条射线对应射线束范围内网格节点的属性信息;(7)选取炮点和检波点相应的射线束对进行成像计算;(8)累加所有射线束对的成像结果。本发明的有益之处在于:该成像方法具有高效率和高精度。
Description
技术领域
本发明涉及一种成像方法,具体涉及一种斜井井间地震高斯束叠前深度偏移成像方法,属于物质或物体的探测技术领域。
背景技术
井间地震勘探时,分别在井中放置震源和检波器,可有效的避开地表低速带对地震信号高频成分的吸收,获得具有较宽频带、极高分辨率的地震信号,可对井间地层、构造、储层等地质目标进行精细成像。
常规的井间地震反射波偏移成像方法分为:波动方程类偏移和射线类偏移。
波动方程偏移方法求解波动方程的数值解,使用递归波场来成像,一般包括:以单程波动方程为基础的共炮集偏移算法和炮检距域偏移算法以及双程波动方程的逆时偏移算法。这些方法要求精确的成像条件并且单程波动方程偏移不能对陡倾构造准确成像,运算效率较低。
以几何射线理论为基础的偏移算法通过计算地震波场的振幅、相位等实现波场延拓成像,例如Kirchhoff偏移方法。但Kirchhoff偏移方法基于零阶射线理论,不能对焦散区、阴影区和多波至区准确成像,因此在复杂构造区域该成像方法往往不能达到精细勘探的需求。
第85届SEG年会上公开了《The inverse Gaussian beam common-reflection-point stack imaging in crosswell seismic》,介绍了斜井井间地震反射波逆高斯束叠加成像方案。该方案将每一个检波点处的能量逆高斯分解到反射点的实际位置,对反射波资料进行归位成像,有效改善了常规的VSP-CDP成像方案。
CN104391327A公开了《一种海上斜井井间地震叠前逆时深度偏移成像方法》,介绍了斜井条件下井间地震资料深度域偏移成像方案。该方案采用TTI介质层析反演方法获得偏移模型,并对其进行网格划分,然后采用初至波射线追踪方法计算激发时间成像条件,最后对井间地震反射波场进行波场逆时延拓获得不同时刻波场值,并应用成像条件对每一个网格进行成像获得井间地震深度偏移成像剖面。
可见,现有成像方法能够在一定程度上解决斜井井间地震反射波场偏移成像,但是基于射线类成像方法计算精度较低,基于波动方程类成像方法实现过程复杂,计算效率也比较低。
发明内容
为解决现有技术的不足,本发明的目的在于提供一种具有高效率和高精度的斜井井间地震高斯束叠前深度偏移成像方法。
一种斜井井间地震高斯束叠前深度偏移成像方法,其特征在于,包括以下步骤:
Step1:读入井间地震反射波场记录、观测***文件、偏移速度文件及相关参数文件,该相关参数文件包括:参与计算地震工区的网格点数、网格间距、高斯束初始宽度、参考频率、最大频率和地震记录采样点数;
Step2:在震源位置高斯束的波前特征为平面波,将地震波场分解成一系列的局部平面部,得到相应的数据体;
Step3:从炮点沿着激发井向下方向到激发井向上方向进行射线追踪,利用程函方程计算任意射线传播的时间(即地震波走时)τ:
其中,v为离散点处的速度值;
Step4:通过运动学追踪获取中心射线后进行高斯束动力学追踪,使用龙格库塔方程求解动力学射线追踪方程:
其中,n为相邻射线上的点到中心射线的垂直距离,p、q的初值分别为垂直于中心射线方向上慢度矢量的分量和相邻射线离开中心射线的距离;
Step5:根据中心射线计算射线束覆盖范围内网格节点的地震波振幅:
其中,A0为炮点振幅,q(R)为检波点处相邻射线离开中心射线的距离,N为射线到达检波点时穿过的地层数,Ri为第i个界面的反射系数或透射系数,αi、βi分别为第i个界面处的入射角与透射角,ρi(R)、vi(R)分别为射线穿过第i个界面前的地层密度与速度, 分别为射线穿过第i个界面后的地层密度与速度;
Step6:从检波点沿着接收井向下方向到接收井向上方向进行运动学和动力学射线追踪,计算并存储每条射线对应射线束范围内网格节点的属性信息,该属性信息包括:地震波的振幅和走时;
Step7:选取炮点和检波点相应的射线束对进行成像计算,计算方程如下:
其中,Is(x)为点x处的成像值,ps为炮点发出射线的慢度值,pbc为检波点发出射线的慢度值,A为权函数,Ds为局部平面波分解结果,L为单炮记录中划分的不同窗,p′为炮记录用于局部倾斜叠加的慢度参数,τ′为炮点发出射线束到达成像点处的走时与检波点发出射线束到达成像点处的走时之和;
Step8:累加所有射线束对的成像结果,得到最终的偏移成像结果。
本发明的有益之处在于:
1、采用井间地震高斯射线束传播算子,解决了波动方程法成像效率低问题,克服了地面地震成像方法在深部构造成像微弱、分辨率不够等不足;
2、使用动力学射线追踪,克服了射线类偏移成像方法在复杂构造区域(如焦散区、奇异区等)无法精细成像的困难;
3、高斯束具有一定的有效宽度,可以选择对计算点有贡献的高斯束的叠加来计算最终的波场,有效提高了计算效率,是一种兼顾效率和精度的成像方法。
附图说明
图1是本发明的斜井井间地震高斯束叠前深度偏移成像方法的流程图;
图2(a)是井间地震高斯束偏移激发点射线束传播范围示意图;
图2(b)是井间地震高斯束偏移接收点射线束传播范围示意图;
图3是速度模型图;
图4是井间地震高斯束叠前深度偏移成像剖面;
图5是中海油某区块实际资料成像剖面中的过井地面地震剖面;
图6是中海油某区块实际资料斜井井间地震高斯束叠前深度偏移成像剖面。
具体实施方式
以下结合附图和具体实施例对本发明作具体的介绍。
参照图1,本发明的斜井井间地震高斯束叠前深度偏移成像方法包括以下步骤:
Step1:读入文件
读入井间地震反射波场记录、观测***文件、偏移速度文件及相关参数文件,其中,该相关参数文件包括:参与计算地震工区的网格点数、网格间距、高斯束初始宽度、参考频率、最大频率和地震记录采样点数。
Step2:分解地震波场
在震源位置高斯束的波前特征为平面波,将地震波场分解成一系列的局部平面部,得到相应的数据体。
Step3:从炮点进行射线追踪并计算地震波旅行时
从炮点沿着激发井向下方向到激发井向上方向进行射线追踪,利用程函方程计算任意射线传播的时间(即地震波走时)τ:
其中,v为离散点处的速度值。
斜井井间地震高斯束叠前深度偏移炮点(激发点)按照一定角度出射的射线束传播范围如图2(a)所示。
Step4:求解动力学射线追踪方程
通过运动学追踪获取中心射线后进行高斯束动力学追踪,使用龙格库塔方程求解动力学射线追踪方程:
其中,n为相邻射线上的点到中心射线的垂直距离,p、q的初值分别为垂直于中心射线方向上慢度矢量的分量和相邻射线离开中心射线的距离。
Step5:求取射线束覆盖范围内网格节点的振幅信息
根据中心射线计算射线束覆盖范围内网格节点的地震波振幅:
其中,A0为炮点振幅,q(R)为检波点处相邻射线离开中心射线的距离,N为射线到达检波点时穿过的地层数,Ri为第i个界面的反射系数或透射系数,αi、βi分别为第i个界面处的入射角与透射角,ρi(R)、vi(R)分别为射线穿过第i个界面前的地层密度与速度,(R)、vi(R)分别为射线穿过第i个界面后的地层密度与速度。
Step6:从检波点进行射线追踪并计算网格节点的属性信息
从检波点沿着接收井向下方向到接收井向上方向进行运动学与动力学射线追踪,计算并存储每条射线束范围内网格节点的属性信息,该属性信息包括:地震波的振幅和走时。
斜井井间地震高斯束叠前深度偏移检波点(接收点)按照一定角度出射的射线束传播范围如图2(b)所示。
Step7:进行成像计算
选取炮点和检波点相应的射线束对进行互相关成像计算,计算方程如下:
其中,Is(x)为点x处的成像值,ps为炮点发出射线的慢度值,pbc为检波点发出射线的慢度值,A为权函数,Ds为局部平面波分解结果,L为单炮记录中划分的不同窗,p′为炮记录用于局部倾斜叠加的慢度参数,τ′为炮点发出射线束到达成像点处的走时与检波点发出射线束到达成像点处的走时之和。
Step8:累加成像结果
累加所有射线束对的成像结果,得到最终的偏移成像结果。
为了验证本发明提供的成像方法具有很好的成像效果,我们分别在正演模型和中海油某区块进行了井间地震资料成像试验。
1、在正演模型上进行井间地震资料成像试验
速度模型(正演模型)如图3所示。
从图3中可以看出:两井之间存在逆断层、砂体尖灭和倾斜产状的地层等。
我们采用本发明提供的方法进行成像,其中:
(1)读入的参数:井间地震反射波场记录(reflectwavefield.cds)、观测***文件(geom.lge)、平滑后的速度文件(velosmooth.dat)。
(2)设定的相关参数文件:参与计算地震工区的网格点数nx、网格间距dx、参考频率f、最大频率fmax、地震记录采样点数nt以及高斯束初始宽度w。
(3)从炮点沿着激发井向下方向到激发井向上方向进行射线追踪,入射角度范围:-20°~160°。
(4)从检波点沿着接收井向下方向到接收井向上方向进行运动学与动力学射线追踪,入射角度:30°~165°。
斜井井间地震高斯束叠前深度偏移成像剖面如图4所示。
从图4中可以看出:本发明的成像方法对图3中所示的逆断层以及倾斜产状的地层成像效果与模型对应非常好。
由此可见,应用本发明提供的成像方法对典型的井间复杂构造模型进行了成功的高斯束叠前深度偏移成像,检验了方法的正确性、有效性和稳定性。
2、在中海油某区块进行井间地震资料成像试验
中海油某区块地面地震勘探过井剖面如图5所示。
从图5可以看出,目的层在0.4s~0.6s之间,目的层中微小断裂构造发育,地面地震成像方法对目标区微小构造分辨程度较低。
我们采用本发明提供的方法进行成像,其中:
(1)读入的文件:反射波场记录(oriwavefield.cds)、观测***文件(geometry.lge)、平滑后的速度文件(velo_smooth.dat)。
(2)设定的相关参数文件:参与计算地震工区的网格点数nx、网格间距dx、参考频率f、最大频率fmax、地震记录采样点数nt以及高斯束初始宽度w。
(3)从炮点沿着激发井向下方向到激发井向上方向进行射线追踪,入射角度范围:-5°~85°。
(4)从检波点沿着接收井向下方向到接收井向上方向进行运动学与动力学射线追踪,入射角度:-160°~20°。
中海油该区块实际资料斜井井间地震高斯束叠前深度偏移成像剖面如图6所示。
从图6可以看出:斜井井间地震高斯束叠前深度偏移成像剖面目的层间构造特征刻画的更清楚,地层产状与地面地震成像结果基本一致,地面地震剖面中CDP间隔为12.5m,井间地震偏移剖面中CDP间隔为3.125m。
此外,从图6还可以看出:本发明提供的斜井井间地震高斯束叠前深度偏移成像方法能够精细刻画地质构造形态,成像结果分辨率更高。
由此可见,应用本发明提供的成像方法对具有较复杂构造的中海油某井间地震勘探区块实际资料进行的高斯束叠前深度偏移成像处理,获得了较好的地质成像效果。
需要说明的是,上述实施例不以任何形式限制本发明,凡采用等同替换或等效变换的方式所获得的技术方案,均落在本发明的保护范围内。
Claims (3)
1.一种斜井井间地震高斯束叠前深度偏移成像方法,其特征在于,包括以下步骤:
Step1:读入井间地震反射波场记录、观测***文件、偏移速度文件及相关参数文件;
Step2:在震源位置高斯束的波前特征为平面波,将地震波场分解成一系列的局部平面部,得到相应的数据体;
Step3:从炮点沿着激发井向下方向到激发井向上方向进行射线追踪,利用程函方程计算任意射线传播的时间τ:
其中,v为离散点处的速度值;
Step4:通过运动学追踪获取中心射线后进行高斯束动力学追踪,使用龙格库塔方程求解运动学射线追踪方程:
其中,n为相邻射线上的点到中心射线的垂直距离,p、q的初值分别为垂直于中心射线方向上慢度矢量的分量和相邻射线离开中心射线的距离;
Step5:根据中心射线计算射线束覆盖范围内网格节点的地震波振幅:
其中,A0为炮点振幅,q(R)为检波点处相邻射线离开中心射线的距离,N为射线到达检波点时穿过的地层数,Ri为第i个界面的反射系数或透射系数,αi、βi分别为第i个界面处的入射角与透射角,ρi(R)、vi(R)分别为射线穿过第i个界面前的地层密度与速度, 分别为射线穿过第i个界面后的地层密度与速度;
Step6:从检波点沿着接收井向下方向到接收井向上方向进行射线追踪,计算并存储每条射线对应射线束范围内网格节点的属性信息;
Step7:选取炮点和检波点相应的射线束对进行成像计算,计算方程如下:
其中,Is(x)为点x处的成像值,ps为炮点发出射线的慢度值,pbc为检波点发出射线的慢度值,A为权函数,Ds为局部平面波分解结果,L为单炮记录中划分的不同窗,p′为炮记录用于局部倾斜叠加的慢度参数,τ′为炮点发出射线束到达成像点处的走时与检波点发出射线束到达成像点处的走时之和;
Step8:累加所有射线束对的成像结果,得到最终的偏移成像结果。
2.根据权利要求1所述的斜井井间地震高斯束叠前深度偏移成像方法,其特征在于,在Step1中,所述相关参数文件包括:参与计算地震工区的网格点数、网格间距、高斯束初始宽度、参考频率、最大频率和地震记录采样点数。
3.根据权利要求1所述的斜井井间地震高斯束叠前深度偏移成像方法,其特征在于,在Step6中,所述属性信息包括:地震波的振幅和走时。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810106904.0A CN108363101B (zh) | 2018-02-02 | 2018-02-02 | 一种斜井井间地震高斯束叠前深度偏移成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810106904.0A CN108363101B (zh) | 2018-02-02 | 2018-02-02 | 一种斜井井间地震高斯束叠前深度偏移成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108363101A true CN108363101A (zh) | 2018-08-03 |
CN108363101B CN108363101B (zh) | 2019-07-02 |
Family
ID=63004200
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810106904.0A Active CN108363101B (zh) | 2018-02-02 | 2018-02-02 | 一种斜井井间地震高斯束叠前深度偏移成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108363101B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110568496A (zh) * | 2019-09-26 | 2019-12-13 | 核工业北京地质研究院 | 一种复杂介质条件下射线追踪方法 |
CN111624647A (zh) * | 2020-06-05 | 2020-09-04 | 中油奥博(成都)科技有限公司 | 一种变偏移距vsp射线追踪积分叠前时间偏移方法和装置 |
CN111624648A (zh) * | 2020-06-05 | 2020-09-04 | 中油奥博(成都)科技有限公司 | 一种起伏地表变偏移距vsp积分叠前深度偏移方法 |
CN111781635A (zh) * | 2019-04-04 | 2020-10-16 | 中国石油天然气集团有限公司 | 海底四分量弹性波高斯束深度偏移方法和装置 |
CN111999770A (zh) * | 2020-09-03 | 2020-11-27 | 中国地质科学院 | 一种tti介质转换ps波精确束偏移成像方法及*** |
BE1027342B1 (nl) * | 2019-10-15 | 2021-01-07 | Univ Southwest Jiaotong | Werkwijze voor anisotrope seismische beeldvorming |
CN112485826A (zh) * | 2020-11-12 | 2021-03-12 | 中国地质大学(武汉) | 绝对波阻抗反演成像方法、装置、设备及存储介质 |
CN112505765A (zh) * | 2020-11-18 | 2021-03-16 | 东华理工大学 | Lax Friedrichs扫描地震波旅行时间的方法 |
CN112904418A (zh) * | 2021-01-22 | 2021-06-04 | 西南交通大学 | 一种自适应射线加密式克希霍夫型束偏移地震波成像方法 |
CN113805237A (zh) * | 2020-06-11 | 2021-12-17 | 中国石油化工股份有限公司 | 偏移陆地交叉排列地震的使用压缩感测模型的方法和*** |
CN114002742A (zh) * | 2021-08-18 | 2022-02-01 | 南方海洋科学与工程广东省实验室(广州) | 一种欧拉高斯束偏移成像方法及装置 |
CN114428327A (zh) * | 2020-09-21 | 2022-05-03 | 中国石油化工股份有限公司 | 高斯束深度偏移成像方法、装置、设备和存储介质 |
CN114924312A (zh) * | 2022-05-10 | 2022-08-19 | 吉林大学 | 一种基于波前初始化射线追踪技术的高斯波束偏移方法及装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005010797A2 (en) * | 2003-07-23 | 2005-02-03 | Lee Wook B | Improved 3d veloctiy modeling, with calibration and trend fitting using geostatistical techniques, particularly advantageous for curved-ray prestack time migration and for such migration followed by prestack depth migration |
US20100256916A1 (en) * | 2009-04-03 | 2010-10-07 | Chevron U.S.A. Inc. | Method for target-oriented reverse time migration for prestack depth imaging |
CN102116870A (zh) * | 2011-02-12 | 2011-07-06 | 中国石油大学(华东) | 弹性波高斯束叠前深度偏移技术 |
CN102183786A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 双复杂条件下保真振幅高斯束叠前深度偏移方法 |
CN103995288A (zh) * | 2014-05-13 | 2014-08-20 | 中国石油天然气集团公司 | 一种高斯束叠前深度偏移方法及装置 |
CN104391327A (zh) * | 2014-12-04 | 2015-03-04 | 中国海洋石油总公司 | 一种海上斜井井间地震叠前逆时深度偏移成像方法 |
CN106291687A (zh) * | 2016-07-21 | 2017-01-04 | 中国地质科学院地质研究所 | 各向异性多波高斯束叠前深度偏移成像方法 |
CN107561583A (zh) * | 2017-08-01 | 2018-01-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 用于高斯束叠前深度偏移的局部角度计算方法及成像方法 |
-
2018
- 2018-02-02 CN CN201810106904.0A patent/CN108363101B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005010797A2 (en) * | 2003-07-23 | 2005-02-03 | Lee Wook B | Improved 3d veloctiy modeling, with calibration and trend fitting using geostatistical techniques, particularly advantageous for curved-ray prestack time migration and for such migration followed by prestack depth migration |
US20100256916A1 (en) * | 2009-04-03 | 2010-10-07 | Chevron U.S.A. Inc. | Method for target-oriented reverse time migration for prestack depth imaging |
CN102116870A (zh) * | 2011-02-12 | 2011-07-06 | 中国石油大学(华东) | 弹性波高斯束叠前深度偏移技术 |
CN102183786A (zh) * | 2011-02-12 | 2011-09-14 | 中国石油大学(华东) | 双复杂条件下保真振幅高斯束叠前深度偏移方法 |
CN103995288A (zh) * | 2014-05-13 | 2014-08-20 | 中国石油天然气集团公司 | 一种高斯束叠前深度偏移方法及装置 |
CN104391327A (zh) * | 2014-12-04 | 2015-03-04 | 中国海洋石油总公司 | 一种海上斜井井间地震叠前逆时深度偏移成像方法 |
CN106291687A (zh) * | 2016-07-21 | 2017-01-04 | 中国地质科学院地质研究所 | 各向异性多波高斯束叠前深度偏移成像方法 |
CN107561583A (zh) * | 2017-08-01 | 2018-01-09 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 用于高斯束叠前深度偏移的局部角度计算方法及成像方法 |
Non-Patent Citations (5)
Title |
---|
张林等: "《VTI介质井间地震高斯束正演数值模拟》", 《物探化探计算技术》 * |
杨飞龙等: "《VSP高斯射线束法叠加成像》", 《物探与化探》 * |
邵荣峰: "《高斯束层析成像反演方法研究》", 《中国石油大学(华东)硕士学位论文》 * |
韩建光等: "《高斯束叠前深度偏移影响因素分析》", 《物探化探计算技术》 * |
魏峥嵘等: "《斜井井间地震高斯射线束正演方法》", 《西安石油大学学报(自然科学版)》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111781635A (zh) * | 2019-04-04 | 2020-10-16 | 中国石油天然气集团有限公司 | 海底四分量弹性波高斯束深度偏移方法和装置 |
CN110568496A (zh) * | 2019-09-26 | 2019-12-13 | 核工业北京地质研究院 | 一种复杂介质条件下射线追踪方法 |
BE1027342B1 (nl) * | 2019-10-15 | 2021-01-07 | Univ Southwest Jiaotong | Werkwijze voor anisotrope seismische beeldvorming |
NL2024231B1 (en) * | 2019-10-15 | 2021-01-27 | Univ Southwest Jiaotong | Anisotropic seismic imaging method |
CN111624647A (zh) * | 2020-06-05 | 2020-09-04 | 中油奥博(成都)科技有限公司 | 一种变偏移距vsp射线追踪积分叠前时间偏移方法和装置 |
CN111624648A (zh) * | 2020-06-05 | 2020-09-04 | 中油奥博(成都)科技有限公司 | 一种起伏地表变偏移距vsp积分叠前深度偏移方法 |
CN111624647B (zh) * | 2020-06-05 | 2022-06-24 | 中油奥博(成都)科技有限公司 | 一种变偏移距vsp射线追踪积分叠前时间偏移方法和装置 |
CN111624648B (zh) * | 2020-06-05 | 2022-04-01 | 中油奥博(成都)科技有限公司 | 一种起伏地表变偏移距vsp积分叠前深度偏移方法 |
CN113805237A (zh) * | 2020-06-11 | 2021-12-17 | 中国石油化工股份有限公司 | 偏移陆地交叉排列地震的使用压缩感测模型的方法和*** |
CN111999770A (zh) * | 2020-09-03 | 2020-11-27 | 中国地质科学院 | 一种tti介质转换ps波精确束偏移成像方法及*** |
CN111999770B (zh) * | 2020-09-03 | 2024-01-16 | 中国地质科学院 | 一种tti介质转换ps波精确束偏移成像方法及*** |
CN114428327A (zh) * | 2020-09-21 | 2022-05-03 | 中国石油化工股份有限公司 | 高斯束深度偏移成像方法、装置、设备和存储介质 |
CN112485826B (zh) * | 2020-11-12 | 2022-04-26 | 中国地质大学(武汉) | 绝对波阻抗反演成像方法、装置、设备及存储介质 |
CN112485826A (zh) * | 2020-11-12 | 2021-03-12 | 中国地质大学(武汉) | 绝对波阻抗反演成像方法、装置、设备及存储介质 |
CN112505765A (zh) * | 2020-11-18 | 2021-03-16 | 东华理工大学 | Lax Friedrichs扫描地震波旅行时间的方法 |
CN112505765B (zh) * | 2020-11-18 | 2023-05-09 | 东华理工大学 | Lax Friedrichs扫描地震波旅行时间的方法 |
CN112904418B (zh) * | 2021-01-22 | 2021-08-17 | 西南交通大学 | 一种自适应射线加密式克希霍夫型束偏移地震波成像方法 |
CN112904418A (zh) * | 2021-01-22 | 2021-06-04 | 西南交通大学 | 一种自适应射线加密式克希霍夫型束偏移地震波成像方法 |
CN114002742A (zh) * | 2021-08-18 | 2022-02-01 | 南方海洋科学与工程广东省实验室(广州) | 一种欧拉高斯束偏移成像方法及装置 |
CN114002742B (zh) * | 2021-08-18 | 2024-02-23 | 南方海洋科学与工程广东省实验室(广州) | 一种欧拉高斯束偏移成像方法及装置 |
CN114924312A (zh) * | 2022-05-10 | 2022-08-19 | 吉林大学 | 一种基于波前初始化射线追踪技术的高斯波束偏移方法及装置 |
CN114924312B (zh) * | 2022-05-10 | 2024-03-12 | 吉林大学 | 一种基于波前初始化射线追踪技术的高斯波束偏移方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN108363101B (zh) | 2019-07-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108363101B (zh) | 一种斜井井间地震高斯束叠前深度偏移成像方法 | |
CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
Xie et al. | Wave-equation-based seismic illumination analysis | |
CN106094032B (zh) | 一种构建地层速度模型的方法 | |
CN105549081A (zh) | 各向异性介质共炮域高斯束偏移成像方法 | |
CN102879819B (zh) | 保持地震波场运动学特征的地震数据处理方法及装置 | |
CN102841379B (zh) | 一种基于共散射点道集的叠前时间偏移与速度分析方法 | |
CN102841375A (zh) | 一种复杂条件下基于角度域共成像点道集的层析速度反演方法 | |
CN102841376A (zh) | 一种基于起伏地表的层析速度反演方法 | |
CN105093292A (zh) | 一种地震成像的数据处理方法和装置 | |
CN105093319B (zh) | 基于三维地震数据的地面微地震静校正方法 | |
CN107870355B (zh) | 一种复杂地形条件下的克希霍夫型波束偏移方法 | |
CN103995288A (zh) | 一种高斯束叠前深度偏移方法及装置 | |
CN102901984B (zh) | 真地表地震数据倾角道集构建方法 | |
CN107515420A (zh) | 一种用于局部相关同相轴的走时与梯度精确拾取方法 | |
CN106199704A (zh) | 一种三维三分量海底电缆地震资料速度建模方法 | |
CN104749623B (zh) | 一种地震资料成像处理方法 | |
Hu et al. | Ray-illumination compensation for adjoint-state first-arrival traveltime tomography | |
CN102053269A (zh) | 一种对地震资料中速度分析方法 | |
Guo et al. | Topography-dependent eikonal tomography based on the fast-sweeping scheme and the adjoint-state technique | |
CN106054252A (zh) | 一种叠前时间偏移的方法及装置 | |
CN105510958A (zh) | 一种适用复杂介质的三维vsp观测***设计方法 | |
CN103336302B (zh) | 基于高次余弦幅度加权的地震波束形成方法 | |
CN109085644A (zh) | 基于双射线走时的真地表成像方法 | |
CN104101901A (zh) | 一种转换波弯曲射线保幅各向异性叠前时间偏移方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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 |