技术特征：
1.一种基于广义染色算法的粘滞介质靶区成像方法，其特征在于：基于广义染色算法的粘滞介质靶区成像方法包括以下步骤：步骤一、设置染色位置，选择与复数域染色算法标记染色位置相同的方法，复数域染色算法方程：复数域染色算法方程：方程中：所代表的原波场，所代表的染色波场，为真实速度场，为染色区域速度场，当且仅当时该式成立，当标记的散射点与目标位置重合时，染色波场的波前将与原波场同步传播；标记区域的值只取1或0，1代表标记位置，0代表未标记位置；步骤二、设置震源、检波点、地震记录信息；步骤三、利用分数阶拉普拉斯算子常q解耦粘声波方程和边界条件求取时间步为t的震源地震波场值，并保存波场值；分数阶拉普拉斯算子常q解耦粘声波方程：算子常q解耦粘声波方程：方程中:c为速度，t为时间，p为粘声波波场；q为品质因子；c0为参考频率ω0处定义的相速度；s为震源；步骤四、将时间步长为t由分数阶拉普拉斯算子常q解耦粘声波方程所计算出的波场值应用于粘滞介质的广义染色算法方程，来计算时间步长为t的染色波场值，保存波场值；粘滞介质的广义染色算法方程：粘滞介质的广义染色算法方程：
方程中步骤五、同时利用方程(2)和方程(3)计算时间步长为t+1的真实波场与染色波场；步骤六、判断时间循环是否到达最大时间步长t，否则重复步骤三至步骤五，是则结束时间循环；步骤七、加载时间步长为t-t的检波点波场数据；步骤八、利用方程(2)的波动方程和边界条件求解时间步为t-t的检波点地震波场值；步骤九、计算时间步长为t-t-1的检波点地震波场值，判断是否到达时间步长0；否则重复步骤七至步骤九，是则结束循环；步骤十、利用成像条件处理所有炮的染色波场、检波点反传波场信息、去噪后得到地下构造成像结果；成像条件：其中r为检波点反传波场，为染色波场。2.根据权利要1所述的一种微基于广义染色算法的粘滞介质靶区成像方法，其特征在于：采用有限差分法求解等号左边的时间导数，利用伪谱法求取等号右端的空间算子；因此，时间精度为二阶精度，空间精度为谱精度，其计算方程为：此，时间精度为二阶精度，空间精度为谱精度，其计算方程为：其中f和f-1
分别为一维正反傅里叶变化，k为离散波数。
技术总结
本发明属于数据处理技术领域，具体涉及一种基于广义染色算法的粘滞介质靶区成像方法。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、利用成像条件处理后得到地下构造成像结果。构造成像结果。构造成像结果。
技术研发人员：王宁 石颖 王子墨 王维红 包培楠 李松龄
受保护的技术使用者：东北石油大学
技术研发日：2022.11.11
技术公布日：2023/3/21