-
公开(公告)号:CN115685322A
公开(公告)日:2023-02-03
申请号:CN202211409557.1
申请日:2022-11-11
Applicant: 东北石油大学
Abstract: 本发明属于数据处理技术领域,具体涉及一种基于虚拟检波点染色算法的粘滞介质靶区成像方法。1、设置染色位置2、设置震源、检波点、地震记录信息;3、求取震源地震波场值;4、保存计算出的真实波场;5、计算真实波场;6、判断时间循环是否达到最大时间步T,否则重复3至5,是则结束时间循环;7、将真实波场到达染色区域的波场值赋予染色波场;8、计算真实波场;9、判断时间循环是否达到最大时间步T,否则重复7至8,是则结束时间循环;10、加载检波点波场数据;11、求解检波点地震波场值;12、计算检波点地震波场值,判断是否到达时间步长0,否则重复7至9,是则结束循环;13、利用成像条件处理后得到地下构造成像结果。
-
公开(公告)号:CN115079264A
公开(公告)日:2022-09-20
申请号:CN202210746533.9
申请日:2022-06-28
Applicant: 东北石油大学
Abstract: 本发明属于油藏开发技术领域,具体涉及一种不依赖子波粘声最小二乘逆时偏移成像方法,在方法中采用常分数阶解耦拉普拉斯粘声波动方程模拟地震波在粘声介质中的传播,建立卷积目标函数,构建具有相同震源的模拟数据和观测数据间的匹配,从而降低对震源子波的敏感性,推导正演算子与伴随算子,在迭代过程中补偿衰减效应,实现不依赖子波的粘声最小二乘逆时偏移,在震源子波不准确的情况下,获取可靠的能够补偿深层能量的成像结果。本方法更能满足油气勘探开发对复杂地下构造成像的需求,提高对油气藏的识别精度,提升我国油气勘探开发核心竞争力。
-
公开(公告)号:CN107589443B
公开(公告)日:2019-05-28
申请号:CN201710701093.4
申请日:2017-08-16
Applicant: 东北石油大学
Abstract: 本发明提供一种基于弹性波最小二乘逆时偏移成像的方法及系统,其中,方法包括,数据接收步骤,接收多分量地震数据;梯度构建步骤,根据预先构建的弹性波方程构建弹性波最小二乘逆时偏移的梯度;数据模拟步骤,根据梯度,采用逆时反偏移方式获取弹性波波场的模拟数据;残差计算步骤,根据多分量地震数据和模拟数据,以及根据预先构建的目标函数,获取数据残差;成像获取步骤,根据数据残差,采用预先构建的共轭梯度,以迭代计算方式获取成像数据;成像输出步骤,输出最小二乘逆时偏移成像剖面。
-
公开(公告)号:CN115079264B
公开(公告)日:2024-08-30
申请号:CN202210746533.9
申请日:2022-06-28
Applicant: 东北石油大学
Abstract: 本发明属于油藏开发技术领域,具体涉及一种不依赖子波粘声最小二乘逆时偏移成像方法,在方法中采用常分数阶解耦拉普拉斯粘声波动方程模拟地震波在粘声介质中的传播,建立卷积目标函数,构建具有相同震源的模拟数据和观测数据间的匹配,从而降低对震源子波的敏感性,推导正演算子与伴随算子,在迭代过程中补偿衰减效应,实现不依赖子波的粘声最小二乘逆时偏移,在震源子波不准确的情况下,获取可靠的能够补偿深层能量的成像结果。本方法更能满足油气勘探开发对复杂地下构造成像的需求,提高对油气藏的识别精度,提升我国油气勘探开发核心竞争力。
-
公开(公告)号:CN115826040A
公开(公告)日:2023-03-21
申请号:CN202211410788.4
申请日:2022-11-11
Applicant: 东北石油大学
Abstract: 本发明属于数据处理技术领域,具体涉及一种基于广义染色算法的粘滞介质靶区成像方法。1、设置染色位置;2、设置震源、检波点、地震记录信息;3、求取时间步为t的震源地震波场值;4、计算时间步长为t的染色波场值;5、计算时间步长为t+1的真实波场与染色波场;6、判断时间循环是否到达最大时间步长T,否则重复3至5,是则结束时间循环;7、加载时间步长为T‑t的检波点波场数据;8、求解时间步为T‑t的检波点地震波场值;9、计算时间步长为T‑t‑1的检波点地震波场值,判断是否到达时间步长0;否则重复7至9,是则结束循环;10、利用成像条件处理后得到地下构造成像结果。
-
公开(公告)号:CN115826040B
公开(公告)日:2025-05-13
申请号:CN202211410788.4
申请日:2022-11-11
Applicant: 东北石油大学
Abstract: 本发明属于数据处理技术领域,具体涉及一种基于广义染色算法的粘滞介质靶区成像方法。1、设置染色位置;2、设置震源、检波点、地震记录信息;3、求取时间步为t的震源地震波场值;4、计算时间步长为t的染色波场值;5、计算时间步长为t+1的真实波场与染色波场;6、判断时间循环是否到达最大时间步长T,否则重复3至5,是则结束时间循环;7、加载时间步长为T‑t的检波点波场数据;8、求解时间步为T‑t的检波点地震波场值;9、计算时间步长为T‑t‑1的检波点地震波场值,判断是否到达时间步长0;否则重复7至9,是则结束循环;10、利用成像条件处理后得到地下构造成像结果。
-
公开(公告)号:CN116882214B
公开(公告)日:2023-12-26
申请号:CN202311146141.X
申请日:2023-09-07
Applicant: 东北石油大学三亚海洋油气研究院
IPC: G06F30/20 , G06F17/11 , G06F111/10 , G06F119/14
Abstract: 本申请公开了基于DFL粘弹性方程的瑞雷波数值模拟方法及系统,属于地震波模拟技术领域,方法包括以下步骤:根据地震波在衰减介质中的传播特性,建立DFL粘弹性方程;根据面波产生条件和消除人工边界反射方法,设置真空自由表面条件和CPML边界条件;以所述DFL粘弹性方程、真空形式和CPML作为模型,采用SGPS法数值实现瑞雷波正演模拟。本申请以DFL粘弹性方程描述体波的传播,真空形式作为自由表面边界条件,CPML为吸收边界条件,采用SGPS法对浅层各向同性衰减介质进行了包括瑞雷波和体波在内的全波场模拟,该方法既继承了真空形式的简单易实现特点,又避免了使用真空形式时困扰基于GSLS方程的数值频散问题。
-
公开(公告)号:CN116882214A
公开(公告)日:2023-10-13
申请号:CN202311146141.X
申请日:2023-09-07
Applicant: 东北石油大学三亚海洋油气研究院
IPC: G06F30/20 , G06F17/11 , G06F111/10 , G06F119/14
Abstract: 本申请公开了基于DFL粘弹性方程的瑞雷波数值模拟方法及系统,属于地震波模拟技术领域,方法包括以下步骤:根据地震波在衰减介质中的传播特性,建立DFL粘弹性方程;根据面波产生条件和消除人工边界反射方法,设置真空自由表面条件和CPML边界条件;以所述DFL粘弹性方程、真空形式和CPML作为模型,采用SGPS法数值实现瑞雷波正演模拟。本申请以DFL粘弹性方程描述体波的传播,真空形式作为自由表面边界条件,CPML为吸收边界条件,采用SGPS法对浅层各向同性衰减介质进行了包括瑞雷波和体波在内的全波场模拟,该方法既继承了真空形式的简单易实现特点,又避免了使用真空形式时困扰基于GSLS方程的数值频散问题。
-
公开(公告)号:CN108181653A
公开(公告)日:2018-06-19
申请号:CN201810039657.7
申请日:2018-01-16
Applicant: 东北石油大学
IPC: G01V1/28
Abstract: 本发明提供的针对VTI介质逆时偏移方法、设备及介质,方法为:在炮集数据中提取震源坐标,并生成震源子波;根据正演方案,通过qP波波动方程对震源子波沿时间正向延拓,在延拓过程中选择若干部分波场进行存储;通过qP波波动方程利用若干部分波场沿时间逆向延拓,得到震源波场;在炮集数据中提取地震数据,通过qP波波动方程对插值后的地震数据沿时间逆向延拓,得到检波点波场;将震源波场和检波点波场进行互相关成像,得到一炮偏移剖面,重复上述步骤,得到多炮偏移剖面;将多个炮偏移剖面进行叠加,得到最终成像剖面。本发明在大幅减小计算量的情况下,能够获得相同运动学特征的纵波分量,使成像剖面的反射同相轴准确归位,提升复杂介质的成像精度。
-
公开(公告)号:CN107589443A
公开(公告)日:2018-01-16
申请号:CN201710701093.4
申请日:2017-08-16
Applicant: 东北石油大学
Abstract: 本发明提供一种基于弹性波最小二乘逆时偏移成像的方法及系统,其中,方法包括,数据接收步骤,接收多分量地震数据;梯度构建步骤,根据预先构建的弹性波方程构建弹性波最小二乘逆时偏移的梯度;数据模拟步骤,根据梯度,采用逆时反偏移方式获取弹性波波场的模拟数据;残差计算步骤,根据多分量地震数据和模拟数据,以及根据预先构建的目标函数,获取数据残差;成像获取步骤,根据数据残差,采用预先构建的共轭梯度,以迭代计算方式获取成像数据;成像输出步骤,输出最小二乘逆时偏移成像剖面。
-
-
-
-
-
-
-
-
-