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

申请/专利权人：东北石油大学

摘要：本发明属于油藏开发技术领域，具体涉及一种不依赖子波粘声最小二乘逆时偏移成像方法，在方法中采用常分数阶解耦拉普拉斯粘声波动方程模拟地震波在粘声介质中的传播，建立卷积目标函数，构建具有相同震源的模拟数据和观测数据间的匹配，从而降低对震源子波的敏感性，推导正演算子与伴随算子，在迭代过程中补偿衰减效应，实现不依赖子波的粘声最小二乘逆时偏移，在震源子波不准确的情况下，获取可靠的能够补偿深层能量的成像结果。本方法更能满足油气勘探开发对复杂地下构造成像的需求，提高对油气藏的识别精度，提升我国油气勘探开发核心竞争力。

主权项：1.一种不依赖子波的粘声最小二乘逆时偏移成像方法，其特征在于：不依赖子波的粘声最小二乘逆时偏移成像方法为：首先，常分数阶拉普拉斯粘滞声波方程为： 式中：Q为品质因子，为拉普拉斯算子，c0为地震波速度，ω0为参考角频率，ωd＝2πfd为主角频率，fd为子波主频，st；xs为震源子波，xs为震源位置，该方程能够解耦相位频散与振幅衰减效应，并将变分数阶拉普拉斯算子合理的近似为不受混合域影响的常分数阶算子，在Q≥15的情况下，能够精确近似变分数阶拉普拉斯粘声波方程；根据扰动理论，将模型参数分解为背景模型c0和扰动模型δc两部分，分别对应背景波场u0和扰动波场δu，采用泰勒公式展开： 代入公式1，可得： 并与公式1相减，可得Born正演方程：Aδu＝mDu05其中，为反射率；常规的粘声最小二乘逆时偏移基于L2范数建立目标函数： 式中，Usyn为模拟数据，Uobs为观测数据，地震数据可表达为：U＝G*s7式中，*为卷积运算，U为地震数据，G为格林函数，将式7代入式6，可得 式中，ssyn为模拟数据子波，sobs为观测数据子波，当两者不一致时，会影响成像结果的计算；为降低子波对成像结果的影响，建立卷积型目标函数： 式中，和分别为模拟数据参考道和观测数据参考道，将式7代入式9可得 式中，s'＝sobs*ssyn，因此，式10看作是具有相同震源s'的模拟数据与观测数据间的匹配，能够降低真实子波与偏移子波不一致带来的影响；利用拉格朗日乘子法求解约束最优化问题， 式中，H为计算区域，μ为伴随波场，令可得伴随方程 式中，代表互相关算子，由于参考道中直达波后的反射经卷积与互相关计算会引入噪音，因此需要对参考道应用时窗，保留直达波部分，所以式12可以重写为： 式中，W为时窗， 式中，tc为时窗长度，n为时窗阶数；进而得到梯度： 采用共轭梯度法迭代更新反射率成像，公式如下： 式中，gk代表梯度，dk代表共轭梯度方向，βk满足： 反射率成像迭代更新：mk+1＝mk+αkdk19式中，k为迭代次数，αk为迭代更新的步长。

全文数据：

权利要求：

百度查询： 东北石油大学 一种不依赖子波的粘声最小二乘逆时偏移成像方法