(19)中华人民共和国国家知识产权局
(12)发明专利说明书 | ||
(10)申请公布号 CN 112987088 A (43)申请公布日 2021.06.18 | ||
(21)申请号 CN202110197468.4
(22)申请日 2021.02.22
(71)申请人 成都理工大学
地址 610059 四川省成都市成华区二仙桥东三路1号
(72)发明人 张祖豪 贺锡雷 刘厚军 贺振华
(74)专利代理机构 50230 重庆市信立达专利代理事务所(普通合伙)
代理人 陈炳萍
(51)Int.CI
G01V1/28(20060101)
G01V1/30(20060101)
G01V1/36(20060101)
权利要求说明书 说明书 幅图 |
(54)发明名称
一种渗流介质地震横波数值模拟和成像方法 | |
(57)摘要
本发明属于多维成像技术领域，公开一种渗流介质地震横波数值模拟和成像方法，包括：按照爆炸反射截面原理把初始模型通过傅氏变换扩展到频率‑空间域；在深度间隔内使用多组参考参数进行波场延拓得到多组波场延拓结果；按地质模型中多组参考参数与实际参数之间的关系，用插值方程求得最终的波场延拓结果；执行频率和深度的运算，完成整个二维剖面的波场延拓运算；对频率域波场做反傅氏变换得到地震记录。本发明通过改进的波场延拓方法，实现非均匀介质的地震横波波场延拓，发展出适用于渗流介质的有效的多物理量、多参数模型的地震横波数值模拟和偏移成像新方法；通过地震横波波场正演和反演成像来预测含流体储层，有助于正确部署勘探和开发方案。 | |
法律状态
法律状态公告日 | 法律状态信息 | 法律状态 |
2021-06-18 | 公开 | 公开 |
2021-07-06 最强地震 | 实质审查的生效 | 实质审查的生效 |
2023-04-18 | 授权 | 发明专利权授予 |
权 利 要 求 说 明 书
1.一种渗流介质地震横波数值模拟和成像方法,其特征在于,所述渗流介质地震横波数值模拟和成像方法包括以下步骤:
按照爆炸反射截面原理,把初始模型u(x,z)通过傅氏变换扩展到频率-空间域U(x,z,ω);
在每一个深度间隔内使用多组参考参数进行波场延拓得到多组波场延拓结果;
按地质模型中多组参考参数与实际参数之间的关系,用适当插值方程求得最终的波场延拓结果;
执行频率ω=1...ω
max
和深度
max
的运算,完成整个二维剖面的波场延拓运算;
对频率域波场U(x,ω)做反傅氏变换得到地震记录u(x,t)。
2.如权利要求1所述的渗流介质地震横波数值模拟和成像方法,其特征在于,将一维方程拓展到二维,将初始模型通过傅氏变换到频率-波数域,纵波和两个横波的方程如下:
其中,U
0
为岩石骨架位移,Q
0
为达西速度场,P
0
为流体压力,φ为孔隙度,μ为剪切模量,k
x
为水平波数,κ为渗透率,η为粘滞系数,ρ
f
为流体密度,ρ
b
为等效体密度,τ为松弛时间,系数
是流体粘滞度和孔隙空间形态有关的函数:
ω为频率,
为无量纲形态因子,其他为中间参数。
3.如权利要求1所述的渗流介质地震横波数值模拟和成像方法,其特征在于,预先取若干组参数,在每一个深度间隔内使用多组参考参数进行波场延拓得到多组波场延拓结果;在ΔZ=Z
i+1
-Z
i
的间隔内的第一波场延拓公式为:
其中,PS-Factor为相移延拓因子,对于横波有:
其中,ρ
b
为地层体密度,μ为地层剪切模量,ω为角频率;
ε=jωρ
f
κ/η,其中
ρ
f
为流体密度,κ为渗透率,η为粘滞系数。
4.如权利要求1所述的渗流介质地震横波数值模拟和成像方法,其特征在于,预先取若干组参数
若N=3,则有3组参考参数P
1
,P
2
,P
3
;对波场U(k
x
,z
m
,ω)延拓到波场
版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。
发表评论