【技术实现步骤摘要】
一种基于进退法的横波分裂分析方法
[0001]本专利技术涉及多分量地震资料处理
，尤其是一种针对裂缝型储层预测的横波分裂分析方法。
技术介绍
[0002]裂缝作为检测非常规油气含量的重要参数，它不仅可以连接孤立的孔隙，增加储层的有效孔隙度，还可以作为重要的油气运移空间，提高储层的渗透率。20世纪90年代开始，随着多波多分量勘探技术开始登上舞台，横波分裂技术开始应用转换横波，近年来，随着3D3C多波多分量勘探技术的飞速发展，横波分裂技术的发展又到了一个新的高潮。当横波的偏振方向以斜交地方式穿过含裂缝的地层时，会发生横波分裂现象(SWS)。初始横波分裂为一个快横波和一个慢横波，分裂的快横波偏振方向与裂缝走向平行，快慢波时差与裂缝的强度密切相关，因此，横波分裂裂缝检测技术成为研究裂缝方向和其发育程度最直接最可靠的方法之一。当今常规的横波分裂分析方法如旋转相关法、切向能量法、协方差矩阵法都是一种枚举式的网格扫描算法。旋转相关法通过快慢横波相关性能很好的寻找最佳的裂缝参数(快横波偏振方向和快慢横波时差)。但是，在面对庞大的三分量...
【技术保护点】
【技术特征摘要】
1.一种基于进退法的横波分裂分析方法，其特征在于，步骤如下：S1、根据快慢横波相关性构造目标函数如下：式中，θ和Δt为横波分裂参数，S1(t)和S2(t)分别为一个时窗内的快横波和慢横波的时间序列；和分别是S1(t)和S2(t)在一个时窗内的平均值；S2、通过确定的真实值或实际数据得到R分量数据R(t)和T分量数据T(t)，进一步计算得到S1(t)和S2(t)；S3、将测试区间划分为4个子区间，取每个子区间的中心点作为初始点；S4、在每个子区间分别采用二维进退法确定该区间的结果值，共计得到4个结果值；S5、比较4个结果值，取得到最大COV值对应的时差和角度作为输出值。2.如权利要求1所述的基于进退法的横波分裂分析方法，其特征在于，所述步骤S2中，总区间为180
×
60，180为方位角范围，60为时差范围，分为四个子区间后，每个子区域的大小为90
×
30。3.如权利要求1所述的基于进退法的横波分裂分析方法，其特征在于，所述步骤S3中，由R分量数据R(t)和T分量数据T(t)，根据以下公式计算得到S1(t)和S2(t)：4.如权利要求1所述的基于进退法的横波分裂分析方法，其特征在于，所述步骤S4包括以下子步骤：(1)对Δt方向零点的搜索；第一步、确定范围，范围的上下界由初始点(θ
k
,Δt
k
)和比较点(θ
k
,Δt
k+1
)确定；Δt
k+1
等于Δt
k
加上步长h；如果COV'
Δt
(θ
k
,Δt
k
)
×
COV'
Δt
(θ
k
,Δt
k+1
)小于零，得到范围[Δt
k
,Δt
k+1
]，进行下一步；如果COV'
Δt
(θ
k
,Δt
k
)
×
COV'
Δt
(θ
k
,Δt
k+1
)大于零，得到新的Δt
k+1
，新Δt
k+1
等于新的Δt
k
加上新的h，新的Δt
k
等于上一步的Δt
k+1
，新的h是上一步h的两倍；然后，再次判断COV'
Δt
(θ
k
,Δt
k
)
×
COV'
Δt
(θ
k
,Δt
k+1
)的值是大于零还是小于零，确实是执行下一步或者继续更新这两个点；第二步、确定另外两个比较点(θ
k
,l
t
)和(θ
k
,u
t
)缩小范围；其中，l
t
和u
t
的计算公式如下：l
t
＝Δt
k
+0.382
×
(Δt
k+1
‑
Δt
k
),u
t
＝Δt
k
+0.618
×
(Δt
k+1
‑
Δt
k
).通过比较|COV'
Δt
(θ
k
,l
t
)|与|COV'
Δt
(θ
k
,u
t
)|的大小，更新l
t
或者u
t
：如果|COV'
Δt
(θ
k
,l
t
)|>|COV'
Δt
(θ
k
,u
t
)|，则Δt
k
＝l
t
,l
t+1
＝Δt
k
+0.382
×
(Δt
k+1
‑
Δt
k
)如果|COV'
Δt
(θ
k
,l
t
)|<|COV'
Δt
(θ
k
,u
t
)|，则Δt
k+1
＝u
t
,
u
t+1
＝Δt
k
+0.618
×
(Δt
k+1<...
【专利技术属性】
技术研发人员：周怀来，杨宇勇，何正涛，王元君，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。