[发明专利]一种裂缝柔度参数获取方法、装置，电子设备及存储介质

权利要求书

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,f_zz,f_xz分别表示加载在v_x上的震源，加载在v_z上的震源，加载在σ_xx上的震源，加载在σ_zz上的震源，加载σ_xz上的震源；A₃₁＝-(λ+2μ)；A₄₁＝-λ；A₅₂＝-μ；B₃₂＝-λ；B₃₃＝-(λ+2μ)；B₅₁＝-μ；表示正则化项。