买专利卖专利找龙图腾，真高效！ 查专利查商标用IPTOP,全免费！专利年费监控用IP管家,真方便！

摘要：本说明书实施例提供了一种频率域粘弹性波正演模拟方法、装置、设备及存储介质，该方法包括：获取频率波数域一维粘弹性波动方程，在将其转换为等效积分弱形式后，离散化为有限元控制方程；加载目标工区的三维点力源和吸收边界条件至有限元控制方程；对加载有三维点力源和吸收边界条件的有限元控制方程进行求解，获得目标工区的频率波数域波场；根据自适应波数采样策略确定频率波数域波场对应的波数集合；对波数集合中每个波数作为输入时对应的频率波数域波场，进行非等间隔二维傅里叶反变换，获得对应的频率域三维弹性波波场。本说明书实施例可以提高频率域粘弹性波正演模拟的精度。

主权项：1.一种频率域粘弹性波正演模拟方法，其特征在于，包括：获取频率波数域一维粘弹性波动方程；将所述频率波数域一维粘弹性波动方程转换为等效积分弱形式；将频率波数域一维粘弹性波动方程的等效积分弱形式离散化为有限元控制方程；加载目标工区的三维点力源和吸收边界条件至所述有限元控制方程；对加载有所述三维点力源和所述吸收边界条件的有限元控制方程进行求解，获得所述目标工区的频率波数域波场；根据自适应波数采样策略确定所述频率波数域波场对应的波数集合；对所述波数集合中每个波数作为输入时对应的频率波数域波场，进行非等间隔二维傅里叶反变换，获得对应的频率域三维弹性波波场；所述频率波数域一维粘弹性波动方程，包括： 其中，分别为σ,e,f,L,u关于x和y方向的二维傅里叶变换，L为偏微分算子矩阵，且σ为应力向量；x,y,z分别为空间三维方向；f为震源矢量；ρ为介质密度；ω为角频率，u为频率域位移矢量；e为应变向量；C为粘弹性介质的复刚度矩阵；T为矩阵转置；分别为关于x,y,z方向的偏导；所述频率域三维弹性波波场包括： 其中，uξ,ζ为ξ,ζ下的频率域三维弹性波波场；n为二维波数域中第n个离散节点的编号，n＝1,2,…,Nd；Nd为离散节点总数；Wn为第n个离散节点处的积分权系数；Ne为二维波数域中划分出的离散单元数量；C＝12π；m＝1,2,…,6分别为第e个离散单元的6个节点编号；为第e个三角单元中第m个离散节点的权系数；为第e个三角单元中第m个离散节点的位移分量，为第n个离散节点的位移分量。

全文数据：

权利要求：

百度查询： 中国石油大学(北京) 频率域粘弹性波正演模拟方法、装置、设备及存储介质