买专利卖专利找龙图腾,真高效! 查专利查商标用IPTOP,全免费!专利年费监控用IP管家,真方便!
申请/专利权人:中国石油大学(北京)
摘要:本发明实施例提供了一种裂缝介质中横波激发Krauklis波的数值模拟方法及设备。其中,所述方法包括:构建有限元离散方程,对所述有限元离散方程加载横波入射边界条件,得到加载边界条件的有限元离散方程;构建边界辅助矩阵、震源辅助矩阵及震源辅助向量,并结合迭代算法,得到迭代公式,采用所述迭代公式,求解所述加载边界条件的有限元离散方程,实现对裂缝介质中横波激发Krauklis波的数值模拟。本发明实施例提供的裂缝介质中横波激发Krauklis波的数值模拟方法及设备,可以有效实现对裂缝介质中横波激发Krauklis波的数值模拟。
主权项:1.一种裂缝介质中横波激发Krauklis波的数值模拟方法,其特征在于,包括:构建有限元离散方程,对所述有限元离散方程加载横波入射边界条件,得到加载边界条件的有限元离散方程;构建边界辅助矩阵、震源辅助矩阵及震源辅助向量,并结合迭代算法,得到迭代公式,采用所述迭代公式,求解所述加载边界条件的有限元离散方程,实现对裂缝介质中横波激发Krauklis波的数值模拟;所述迭代算法采用β-纽曼算法,所述迭代公式包括: 其中,uk为k时刻质点振动的位移;uk+1为k+1时刻质点振动的位移;BC1为所述边界辅助矩阵;BC2为震源辅助矩阵;b0为震源辅助向量;Δt为时长间隔;为扰动后的阻尼矩阵;为扰动后的刚度矩阵;M为质量矩阵;ak为k时刻质点振动的加速度;ak+1为k+1时刻质点振动的加速度;w为地震子波;vk为k时刻质点振动的速度;vk+1为k+1时刻质点振动的速度;所述边界辅助矩阵的边界节点处的y分量对应的对角线元素为0,其余对角线元素为1;所述震源辅助矩阵的震源节点y分量对应的对角线元素为0,其余对角线元素为1;所述震源辅助向量是由0和1构成的列向量,震源节点x分量对应的元素为1,其余元素为0。
全文数据:
权利要求:
百度查询: 中国石油大学(北京) 裂缝介质中横波激发Krauklis波的数值模拟方法及设备
免责声明
1、本报告根据公开、合法渠道获得相关数据和信息,力求客观、公正,但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解,仅供参考使用,不能作为本公司承担任何法律责任的依据或者凭证。