【技术实现步骤摘要】
一种裂缝柔度参数获取方法、装置,电子设备及存储介质
[0001]本申请涉及地质勘测领域,具体而言,涉及一种裂缝柔度参数获取方法、装置,电子设备及存储介质。
技术介绍
[0002]在地震波传播过程中,裂缝的存在会引起地震波的响应,其影响可以用裂缝的柔度来表示。
[0003]目前,采用线性滑动理论来描述裂缝,其具有柔度参数变化较大,裂缝厚度趋近于零的特点。该方式使得在全波形反演过程中裂缝的边界不够突出,进而造成裂缝位置以及裂缝柔度参数的反演结果不正确。
技术实现思路
[0004]本申请实施例的目的在于提供一种裂缝柔度参数获取方法、装置,电子设备及存储介质,以在全波形反演中突出裂缝边界、更加准确地确定裂缝位置,以及更加准确地获取裂缝柔度参数的问题。
[0005]本专利技术是这样实现的:
[0006]第一方面,本申请实施例提供一种裂缝柔度参数获取方法,包括:获取裂缝的法向柔度参数初始设定值及切向柔度参数初始设定值;将所述法向柔度参数初始设定值及所述切向柔度参数初始设定值输入至预先构建的全波形反 ...
【技术保护点】
【技术特征摘要】
1.一种裂缝柔度参数获取方法,其特征在于,包括:获取裂缝的法向柔度参数初始设定值及切向柔度参数初始设定值;将所述法向柔度参数初始设定值及所述切向柔度参数初始设定值输入至预先构建的全波形反演模型中进行迭代,当迭代至预设条件时,输出所述裂缝的法向柔度参数目标值及切向柔度参数目标值;其中,所述全波形反演模型包括预设的柔度梯度函数;所述柔度梯度函数基于正则化约束项所构建;所述柔度梯度函数的表达式为:所述柔度梯度函数的表达式为:其中,表示裂缝法向柔度的梯度值;通过所述裂缝法向柔度的梯度值可计算下一次迭代的法向柔度参数;w
xx
表示反传波场中沿x轴方向的应力;w
zz
表示反传波场中沿z轴方向的应力;w
xz
表示反传波场中沿xz轴方向的应力;x表示切向方向;z表示法向方向;λ、μ表示背景介质的拉梅系数,ρ表示所述背景介质的密度;Δz表示有限差分中单位网格的法向方向的长度;S
N
表示全波形反演中的所述裂缝的法向柔度参数值,初始时S
N
为所述法向柔度参数初始设定值;σ
zz
表示正传波场中沿z轴方向的应力;σ
xz
表示正传波长中沿xz轴方向的应力;k表示约束项;β为一设定参数;T表示终止时间;G
ST
表示裂缝切向柔度的梯度值;通过所述裂缝切向柔度的梯度值可计算下一次迭代的切向柔度参数;S
T
表示全波形反演中的所述裂缝
的切向柔度参数值,初始时S
T
为所述切向柔度参数初始设定值;Ω表示整个模型空间。2.根据权利要求1所述的方法,其特征在于,通过如下步骤构建所述柔度梯度函数,包括:构建反演目标函数;其中,所述反演目标函数包括二范数目标函数及所述正则化约束项;基于拉格朗日算法,获取所述二范数目标函数的裂缝柔度梯度表达式;基于所述拉格朗日算法,获取所述正则化约束项对应的裂缝柔度梯度表达式;基于所述二范数目标函数的裂缝柔度梯度表达式及所述正则化约束项对应的裂缝柔度梯度表达式,生成所述柔度梯度函数。3.根据权利要求2所述的方法,其特征在于,所述反演目标函数的表达式为:其中,表示反演目标函数值;m是包含S
N
和S
T
的参数集合;G表示整个目标空间;u
syn
表示合成地震记录;u
obs
表示真实的地震记录;w=(w
x
,w
z
,w
xx
,w
zz
,w
xz
)
T
,表示反传波场的参数;w
x
表示反传波场中速度的切向分量;w
z
表示反传波场中速度的法向分量;u=(v
x
,v
z
,σ
xx
,σ
zz
,σ
xz
)
T
,表示正传波场的参数;v
x
表示正传波场中速度的切向分量;v
z
表示正传波场中速度的法向分量;σ
xx
表示正传波场中沿x轴方向的应力;f=(f
x
,f
z
,f
xx
,f
zz
,f
xz
)
T
;f
x
,f
z
,f
xx
,...
【专利技术属性】
技术研发人员:王康,彭苏萍,卢勇旭,崔晓芹,
申请(专利权)人:中国矿业大学北京,
类型:发明
国别省市:
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。