[发明专利]基于全波形反演技术确定地质裂缝柔度参数的方法及装置

权利要求书

1.一种基于全波形反演技术确定地质裂缝柔度参数的方法，其特征在于，所述方法包括：

获取真实的地震波信号的传播记录；其中，所述传播记录包括：所述地震波信号的初始传播速度，所述地震波信号产生的震源位置，及表征水平地质裂缝对所述地震波信号传播影响的真实速度值和真实应力值；

获取速度-应力的表达式；其中，所述速度-应力的表达式中包括：所述水平地质裂缝处的柔度参数；

基于预先确定的所述水平地质裂缝处的柔度参数的初始值、预先确定的第一地震参数，所述地震波信号的初始传播速度，所述震源位置，和所述速度-应力的表达式，得到用于模拟所述水平地质裂缝对所述地震波信号传播影响的第一速度值和第一应力值；

基于所述速度-应力的表达式，确定出反传波波动方程；

基于所述真实速度值、所述真实应力值、所述第一速度值、所述第一应力值和所述水平地质裂缝处的柔度参数的初始值对所述波动方程求偏导，得到所述水平地质裂缝处的柔度参数的梯度值；

在确定所述梯度值小于等于预设阈值时，确定所述水平地质裂缝处的柔度参数的初始值为所述水平地质裂缝处的柔度参数的值。

2.根据权利要求1所述的方法，其特征在于，在所述得到所述水平地质裂缝处的柔度参数的梯度值之后，所述方法还包括：

在确定所述梯度值大于所述预设阈值时，对所述水平地质裂缝处的柔度参数的初始值进行更新，得到更新后的地质裂缝处的柔度参数值；

基于所述更新后的地质裂缝处的柔度参数值和所述速度-应力的表达式，得到用于模拟所述水平地质裂缝对所述地震波信号传播影响的第二速度值和第二应力值；

基于所述真实速度值、所述真实应力值、所述第二速度值、所述第二应力值、所述更新后的地质裂缝处的柔度参数值，和所述波动方程，得到所述水平地质裂缝处的柔度参数的新的梯度值；

在确定所述新的梯度值小于所述预设阈值时，确定所述更新后的地质裂缝处的柔度参数值为所述水平地质裂缝处的柔度参数的值。

3.根据权利要求1所述的方法，其特征在于，所述反传波波动方程为其中，w为所述反传波的波场参数向量，u^syn包括所述第一速度值和所述第一应力值；u^obs包括所述真实速度值和所述真实应力值；G表征网格模型，T表征所述地震波信号的接收时长；Lu-f为所述速度-应力的表达式的简写；Ω表征整个时间空间；,代表标量积。

4.根据权利要求3所述的方法，其特征在于，所述基于所述真实速度值、所述真实应力值、所述第一速度值、所述第一应力值和所述水平地质裂缝处的柔度参数的初始值对所述波动方程求偏导，得到所述水平地质裂缝处的柔度参数的梯度值包括：

对所述波动方程求关于u的偏导，并令其结果为0，得到伴随方程；

基于所述真实速度值、所述真实应力值、所述第一速度值、所述第一应力值、所述水平地质裂缝处的柔度参数的初始值，和所述伴随方程，得到所述反传波的波场参数向量的值；

对所述波动方程求关于L的偏导，并令其结果为0，得到所述水平地质裂缝处的柔度参数的梯度方程；

将所述反传波的波场参数向量的值和所述水平地质裂缝处的柔度参数的初始值，输入所述梯度方程，得到所述梯度值。

5.根据权利要求4所述的方法，其特征在于，对所述波动方程求关于L的偏导，并令其结果为0，得到所述水平地质裂缝处的柔度参数的梯度方程，包括：

对所述波动方程求关于L的偏导，并令其结果为0，得到所述水平地质裂缝处的柔度参数的初始梯度方程；

利用所述速度-应力的表达式，替换所述初始梯度方程中的速度，得到所述梯度方程。