井中地震偏移成像方法技术

本发明专利技术提供一种井中地震偏移成像方法，包括：步骤1，获取用于偏移成像的井中地震共炮点道集

【技术实现步骤摘要】
井中地震偏移成像方法


[0001]本专利技术涉及勘探地震
，特别是涉及到一种井中地震偏移成像方法
。

技术介绍

[0002]地震勘探技术作为寻找油气藏的重要手段，已经广泛应用于生产实际中
。
但随着地震勘探的需求不断提高，构造简单
、
易于探测的油气藏逐步减少，取而代之的是一些深层特殊地质体下的复杂
、
精细构造油气藏，勘探目标日趋复杂
。
复杂构造油气藏具有储层薄
、
分布广
、
赋存状态隐蔽等特点，对储层内部小构造的高精度成像有重大需求，这对地面地震勘探技术提出了挑战
。
井中地震勘探技术是在传统的地面地震勘探和
VSP
勘探的基础上发展起来的一门新技术
。
不同于地面地震技术和
VSP
技术，井中地震采用井中激发地面接收的方式进行观测
。
震源位置靠近目的层，且波场减少了一次经过浅层降速带，采集信号信噪比高，对储层识别能力强，可以对井周构造进行高精度成像，弥补了常规地面地震在深层复杂油气藏勘探及井周小构造成像效果不理想的现象，是一种新型物探方法
。
因此，发展井中地震对当前深层复杂油气藏勘探开发至关重要
。
[0003]井中地震技术具有独特的观测方式，其震源和检波器的布置与地面地震不同，这使得地震波的传播路径与地面地震存在较大差异，因此对于井中地震，不能直接套用相对成熟的地面地震成像流程/>。
此外，井中激发地震波需要考虑实际情况，为保护井壁，震源能量较小，使得信号能量弱
。
地层的粘滞效应造成地震信号的高频成分被严重吸收，相比于地面地震，井中地震吸收衰减的影响更为严重，有效信号更加微弱，因此对于井中地震应需研究新的针对性成像方法
。
井中地震技术主要目的是对深层复杂构造进行井周区域精细成像
。
地震资料逆时偏移成像技术基于双程波动理论，不受倾角和横向变速的限制，克服了传统射线类偏移和单程波偏移方法的局限性，在复杂构造油气藏成像方面具有显著优势
。
鉴于上述原因，需要发展一种针对性的井中地震偏移成像方法
。
[0004]在申请号：
CN201610701979.4
的中国专利申请中，涉及到一种地震偏移成像方法及装置，所述方法包括：获取地震数据，根据所述地震数据确定成像速度场数据和品质因子场数据；根据所述品质因子场数据，计算射线路径补偿旅行时间；获取所述地震数据对应的走时表
、
补偿旅行时间表和叠前地震道数据；根据所述走时表和补偿旅行时间表，对所述叠前地震道数据进行补偿；根据所述补偿后的叠前地震道数据，生成成像结果
。
该申请实施例提供的地震偏移成像方法及装置，可以提高偏移成像结果的准确性
。
[0005]在申请号：
CN201110029735.3
的中国专利申请中，涉及到提升波动方程叠前深度偏移成像效果的方法，应用于地震勘探中反射地震资料处理，改善波动方程叠前深度偏移应用效果
。
通过建立深度
‑
角度域剩余动校正关系和层速度反演方法，直接基于角道集的剩余动校量更新波动方程叠前深度偏移的偏移速度模型，以此提升波动方程叠前深度偏移的成像效果
。
利用深度
‑
角度域剩余动校正关系，该方法实现了波动方程叠前深度偏移的偏移道集的剩余动校正和噪音及拉伸切除，改善了偏移叠加剖面的信噪比和分辨率；改进了应用于叠前反演的角道集的质量，能更好地实现直接识别地下构造的含油
、
气或水；能应用于
反射地震资料的二维和三维波动方程叠前深度偏移，对油气
、
矿产资源勘探有重要应用价值
。
[0006]在申请号：
CN201110112909.2
的中国专利申请中，涉及到一种地震多分量联合弹性逆时偏移成像方法，属于勘探地球物理学领域，特征是通过直接利用多分量地震资料作为输入，在对输入数据不进行波场分离的情况下，基于弹性波方程进行正向和逆时外推，联合多分量共同构建地下弹性矢量地震波场，应用照明补偿的弹性波互相关成像条件得到物理意义明确的纵波
、
横波
、
转换纵波和转换横波等四种弹性波成像结果，并进行低通角度域滤波压制弹性逆时偏移中的低波数噪音，得到最终的弹性波偏移成像结果
。
该专利技术通过直接输入多分量地震资料，联合使用多分量共同构建地下弹性矢量地震波场，用于复杂地球介质准确成像
。
[0007]以上现有技术均与本专利技术有较大区别，未能解决我们想要解决的技术问题，为此我们专利技术了一种新的井中地震偏移成像方法
。

技术实现思路

[0008]本专利技术的目的是提供一种可高效获取高质量井中地震资料叠前深度偏移剖面的井中地震偏移成像方法
。
[0009]本专利技术的目的可通过如下技术措施来实现：井中地震偏移成像方法，该井中地震偏移成像方法包括：
[0010]步骤1，获取用于偏移成像的井中地震共炮点道集
、
初始速度模型
、
初始品质因子模型；
[0011]步骤2，确定偏移所用的横向网格间隔和共炮点道集的有效频带范围及主频；
[0012]步骤3，确定偏移速度模型
v
p
和偏移品质因子模型
Q
p
；
[0013]步骤4，获取每个共炮点道集对应的观测参数，确定偏移孔径及有效成像范围；
[0014]步骤5，基于数值算法求解波动方程，获取该炮的成像时间及相应的震源波场；
[0015]步骤6，利用成像条件实现对该共炮点道集的偏移成像；
[0016]步骤7，将所有炮的单炮偏移剖面进行叠加，获得最终的偏移剖面
。
[0017]本专利技术的目的还可通过如下技术措施来实现：
[0018]在步骤1中，依据探测目标及探区地质条件，获取以
SEGY
或
SGY
格式存储的用于偏移成像的井中地震共炮点道集
D(x
r
,t
；
x
s
)、
初始速度模型
v、
初始品质因子模型
Q
；其中，
x
s
＝
(x
s
,y
s
,z
s
)
表示震源空间位置矢量，
x
r
＝
(x
r
,y
r
,z
r
)
表示检波点空间位置矢量，
t
表示波传播时间
。
[0019]在步骤2中，基于本探区目标面元尺寸，确定偏本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.
井中地震偏移成像方法，其特征在于，该井中地震偏移成像方法包括：步骤1，获取用于偏移成像的井中地震共炮点道集
、
初始速度模型
、
初始品质因子模型；步骤2，确定偏移所用的横向网格间隔和共炮点道集的有效频带范围及主频；步骤3，确定偏移速度模型
v
p
和偏移品质因子模型
Q
p
；步骤4，获取每个共炮点道集对应的观测参数，确定偏移孔径及有效成像范围；步骤5，基于数值算法求解波动方程，获取该炮的成像时间及相应的震源波场；步骤6，利用成像条件实现对该共炮点道集的偏移成像；步骤7，将所有炮的单炮偏移剖面进行叠加，获得最终的偏移剖面
。2.
根据权利要求1所述的井中地震偏移成像方法，其特征在于，在步骤1中，依据探测目标及探区地质条件，获取以
SEGY
或
SGY
格式存储的用于偏移成像的井中地震共炮点道集
D(x
r
,t
；
x
s
)、
初始速度模型
v、
初始品质因子模型
Q
；其中，
x
s
＝
(x
s
,y
s
,z
s
)
表示震源空间位置矢量，
x
r
＝
(x
r
,y
r
,z
r
)
表示检波点空间位置矢量，
t
表示波传播时间
。3.
根据权利要求1所述的井中地震偏移成像方法，其特征在于，在步骤2中，基于本探区目标面元尺寸，确定偏移所用的横向网格间隔
dx
和
dy。4.
根据权利要求1所述的井中地震偏移成像方法，其特征在于，在步骤2中，对共炮点道集
D(x
r
,t
；
x
s
)
进行频谱分析，确定其有效频带范围及主频
f
dom
；其中，
x
s
＝
(x
s
,y
s
,z
s
)
表示震源空间位置矢量，
x
r
＝
(x
r
,y
r
,z
r
)
表示检波点空间位置矢量，
t
表示波传播时间
。5.
根据权利要求1所述的井中地震偏移成像方法，其特征在于，在步骤3中，扫描共炮点道集所涉及的大地坐标，建立工区观测系统，确定偏移速度模型和偏移品质因子模型
。6.
根据权利要求5所述的井中地震偏移成像方法，其特征在于，在步骤3中，对于井中地震共炮点道集
D(x
r
,t
；
x
s
)
，逐道进行道头参数读取，其中，
x
s
＝
(x
s
,y
s
,z
s
)
表示震源空间位置矢量，
x
r
＝
(x
r
,y
r
,z
r
)
表示检波点空间位置矢量，
t
表示波传播时间；获取每一道的四个道头关键字
Sx
，
Sy
，
Gx
，
Gy
；通过对所有道的四个道头关键字分别进行大小比较，获取一下四个参数：
Min
x
，
Min
y
、Max
x
，
Max
y
，其中
Min
x
表示所有道的两个道头关键字
Sx
和
Gx
的最小值，
Max
x
表示所有道的两个道头关键字
Sx
和
Gx
的最大值，
Min
y
表示所有道的两个道头关键字
Sy
和
Gy
的最小值，
Max
y
表示所有道的两个道头关键字
Sy
和
Gy
的最大值；以
Min
x
和
Min
y
为坐标原点，基横向网格间隔
dx
和
dy
建立工区观测系统，确定工区最终成像剖面对应的横向空间网格点数
Nx
和
Ny。7.
根据权利要求6所述的井中地震偏移成像方法，其特征在于，在步骤3中，横向空间网格点数
Nx
和
Ny
具体为：和
8.
根据权利要求7所述的井中地震偏移成像方法，其特征在于，在步骤3中，依据初始速度模型
v
和初始品质因子模型
Q
的大地坐标位置及所建立的工区观测系统，将初始速度模型
v、
初始品质因子模型
Q
投影至所建立的工区观测系统，确定用于偏移的偏移速度模型
v
p
和偏移品质因子模型
Q
p
；将初始速度模型
v、
初始品质因子模型
Q
投影至所建立的工区观测系
统，确定用于偏移的偏移速度模型
v
p
和偏移品质因子模型
Q
p
。9.
根据权利要求8所述的井中地震偏移成像方法，其特征在于，在步骤3中，将工区观测系统覆盖到的初始速度模型
v
和初始品质因子模型
Q
提取出来，以所建立的工区观测系统为坐标系统，形成规则的速度模型和品质因子模型，即为所述的偏移速度模型
v
p
和偏移品质因子模型
Q
p
。10.
根据权利要求1所述的井中地震偏移成像方法，其特征在于，在步骤4中，获取每个共炮点道集对应的观测参数，包括炮点空间位置和检波点空间位置，提取该共炮点道集对应的局部偏移速度模型及局部偏移品质因子模型，确定偏移孔径及有效成像范围
。11.
根据权利要求
10
所述的井中地震偏移成像方法，其特征在于，在步骤4中，针对每一炮，读取该炮的共炮点道集的道头参数，获取该炮的炮点深度
selev
，炮点空间位置
sx
和
sy
，每一道对应的检波点深度
gelev
以及检波点空间位置
gx
和
gy
；对该炮所有道的四个道头关键字
sx
，
sy
，
gx
，
gy
分别进行大小比较，获取一下四个参数：
min
x
，
min
y
、max
x
，
max
y
，其中
min
x
表示该炮的所有道的两个道头关键字
sx
和
gx
的最小值，
max
x
表示该炮的所有道的两个道头关键字
sx
和
gx
的最大值，
min
y
表示该炮的所有道的两个道头关键字
sy
和
...

【专利技术属性】
技术研发人员：孔庆丰，孟阳，孙卫国，杨宏伟，赵茂强，杨海中，徐永辉，杜正，
申请(专利权)人：中国石油化工股份有限公司胜利油田分公司物探研究院，
类型：发明
国别省市：