基于立体成像条件的一次波和多次波联合成像方法及系统转让专利
申请号 : CN201911059882.8
文献号 : CN110780348B
文献日 : 2021-07-09
发明人 : 李志娜 , 彭思凯 , 李振春 , 王鹏 , 孙苗苗 , 贺紫林 , 丁艺璇
申请人 : 中国石油大学(华东)
摘要 :
权利要求 :
1.一种基于立体成像条件的一次波和多次波联合成像方法,其特征在于,所述成像方法包括:
获取地表观测点的地震数据,所述地震数据包括一次波和多次波;
利用表层多次波去除方法从所述地震数据中提取震源子波;
根据所述地震数据和所述震源子波,采用单程波波动方程沿深度轴进行正向波场延拓,确定深度为z处的波场,记为延拓后震源波场;
根据所述地震数据,采用单程波波动方程沿深度轴进行反向波场延拓,确定深度为z处的波场,记为延拓后检波波场;
利用局部倾斜叠加方法分别对所述延拓后震源波场和所述延拓后检波波场进行分解,得到不同局部斜率的震源波场和不同局部斜率的检波波场;
根据设定的速度场,利用平面波解构滤波器确定震源波场和检波波场局部斜率的取值范围;
在所述局部斜率的取值范围内,将分解的震源波场和检波波场进行互相关确定成像结果。
2.根据权利要求1所述的基于立体成像条件的一次波和多次波联合成像方法,其特征在于,所述利用表层多次波去除方法从所述地震数据中提取震源子波,具体包括:根据所述地震数据,基于反馈迭代的表层多次波去除方法确定地表加权算子;
根据所述地表加权算子,采用傅里叶反变换确定所述震源子波。
3.根据权利要求2所述的基于立体成像条件的一次波和多次波联合成像方法,其特征在于,所述根据所述地震数据,基于反馈迭代的表层多次波去除方法确定地表加权算子,具体公式为:
(n+1) (n+1) (n)ΔP (z0)=P(z0)‑A (w)ΔP (z0)P(z0)(n+1)
其中,n表示迭代次数,n=0,1,2...,z0表示地表观测点的深度,ΔP (z0)表示第n次(0) (n+1)
迭代后一次波数据,ΔP (z0)=P(z0),P(z0)表示地震数据,A (w)表示第n次迭代后地表加权算子,w表示角频率。
4.根据权利要求1所述的基于立体成像条件的一次波和多次波联合成像方法,其特征在于,所述单程波波动方程具体为:其中, u(x,z,t)表示波场,v(x,z,w)表示地震波速度,+符号表示上行波场,-符号表示下行波场,x表示水平方向坐标,z表示垂直方向坐标,t表示地震波传播时间,w表示角频率。
5.根据权利要求1所述的基于立体成像条件的一次波和多次波联合成像方法,其特征在于,所述在所述局部斜率的取值范围内,将分解的震源波场和检波波场进行互相关确定成像结果,具体成像公式为:
I(x,z)=∫∫WS(x,z,p,t)WD(x,z,p,t)dpdt其中,p表示空间和时间的局部斜率三维矢量,WS(x,z,p,t)表示分解的震源波场,WD(x,z,p,t)表示分解的检波波场,x表示水平方向坐标,z表示垂直方向坐标,t表示地震波传播时间。
6.一种基于立体成像条件的一次波和多次波联合成像系统,其特征在于,所述成像系统包括:
地震数据获取模块,用于获取地表观测点的地震数据,所述地震数据包括一次波和多次波;
震源子波提取模块,用于利用表层多次波去除方法从所述地震数据中提取震源子波;
延拓后震源波场确定模块,用于根据所述地震数据和所述震源子波,采用单程波波动方程沿深度轴进行正向波场延拓,确定深度为z处的波场,记为延拓后震源波场;
延拓后检波波场确定模块,用于根据所述地震数据,采用单程波波动方程沿深度轴进行反向波场延拓,确定深度为z处的波场,记为延拓后检波波场;
波场分解模块,用于利用局部倾斜叠加方法分别对所述延拓后震源波场和所述延拓后检波波场进行分解,得到不同局部斜率的震源波场和不同局部斜率的检波波场;
局部斜率取值范围确定模块,根据设定的速度场,利用平面波解构滤波器确定震源波场和检波波场局部斜率的取值范围;
成像结果确定模块,用于在所述局部斜率的取值范围内,将分解的震源波场和检波波场进行互相关确定成像结果。
7.根据权利要求6所述的基于立体成像条件的一次波和多次波联合成像系统,其特征在于,所述震源子波提取模块,具体包括:地表加权算子确定单元,根据所述地震数据,基于反馈迭代的表层多次波去除方法确定地表加权算子;
震源子波确定单元,用于根据所述地表加权算子,采用傅里叶反变换确定所述震源子波。
8.根据权利要求7所述的基于立体成像条件的一次波和多次波联合成像系统,其特征在于,所述地表加权算子确定单元,具体公式为:(n+1) (n+1) (n)ΔP (z0)=P(z0)‑A (w)ΔP (z0)P(z0)(n+1)
其中,n表示迭代次数,n=0,1,2...,z0表示地表观测点的深度,ΔP (z0)表示第n次(0) (n+1)
迭代后一次波数据,ΔP (z0)=P(z0),P(z0)表示地震数据,A (w)表示第n次迭代后地表加权算子,w表示角频率。
9.根据权利要求6所述的基于立体成像条件的一次波和多次波联合成像系统,其特征在于,所述单程波波动方程具体为:其中, u(x,z,t)表示波场,v(x,z,w)表示地震波速度,+符号表示上行波场,-符号表示下行波场,x表示水平方向坐标,z表示垂直方向坐标,t表示地震波传播时间,w表示角频率。
10.根据权利要求6所述的基于立体成像条件的一次波和多次波联合成像系统,其特征在于,所述成像结果确定模块,具体成像公式为:I(x,z)=∫∫WS(x,z,p,t)WD(x,z,p,t)dpdt其中,p表示空间和时间的局部斜率三维矢量,WS(x,z,p,t)表示分解的震源波场,WD(x,z,p,t)表示分解的检波波场,x表示水平方向坐标,z表示垂直方向坐标,t表示地震波传播时间。
说明书 :
基于立体成像条件的一次波和多次波联合成像方法及系统
技术领域
背景技术
的地球内部精细结构信息。传统的多次波成像方法需要首先进行多次波分离,这是个耗时
耗力的过程。为了节省计算成本,部分学者提出了一种一次波和多次波联合成像的思路。
波的选择存在一定的问题,子波能量过小会影响一次反射波的成像效果,而能量过大则会
造成多次波成像信息被一次波成像结果所掩盖,且子波选取不当会带来成像位置不准确等
问题;此外,成像过程中不同阶多次波引起的不相干波场互相关会引入大量的串扰假象,严
重影响成像的质量,二维反褶积成像条件可以压制一定的串扰噪声,但是对复杂构造的适
应性较差,降低了传统方法的应用价值。
发明内容
第n次迭代后一次波数据,ΔP (z0)=P(z0),P(z0)表示地震数据,A (w)表示第n次迭代
后地表加权算子,w表示角频率。
表示地震波传播时间,w表示角频率。
传播时间。
第n次迭代后一次波数据,ΔP (z0)=P(z0),P(z0)表示地震数据,A (w)表示第n次迭代
后地表加权算子,w表示角频率。
表示地震波传播时间,w表示角频率。
传播时间。
计的地层局部斜率基础上,同时利用延拓波场的传播时间和传播方向来约束成像过程,基
于多参数约束的立体成像条件来实现一次波和多次波的联合成像,提高一次波和多次波联
合成像的质量。
附图说明
例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获
得其他的附图。
具体实施方式
本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他
实施例,都属于本发明保护的范围。
括:
第n次迭代后一次波数据,ΔP (z0)=P(z0),P(z0)表示地震数据,A (w)表示第n次迭代
后地表加权算子,w表示角频率。
乘反演理论迭代实现多次波的预测,同时多次迭代以后可以求得逆震源子波A(w),对其进
行求逆并求取傅里叶反变换即可获取震源子波,表示如下:
提取更为可靠,且适用于含多次波数据的子波提取。
波的选择存在一定的问题,震源子波能量过小会影响一次反射波的成像效果,而能量过大
则会造成多次波成像信息被一次波成像结果所掩盖,且震源子波选取不当会造成错误的成
像位置等问题。
的脉冲震源作为虚拟震源即输入震源波场,将记录的地震数据作为接收波场即检波波场,
震源波场和检波波场表示如下:
的震源子波,d(xr,z0,t)表示距离地面z0深度处位于x=xr处的地震数据。
表示地震波传播时间,w表示角频率。
表示与局部斜率相关的一系列震源波场,WD(x,z,p,t)表示分解的检波波场即表示与局部
斜率相关的一系列检波波场,x表示水平方向坐标,z表示垂直方向坐标,t表示地震波传播
时间。
在应用成像条件时利用多参数(时间和局部斜率)同时约束成像过程来进行成像串扰的压
制。
加。
为了提高成像精度,参与叠加的局部斜率的范围应足够大,采样间隔应足够小,而与之伴随
的则是高计算成本。为了兼顾成像精度和计算成本,本发明引入平面解构滤波器从速度场
中估计反射层的斜率,以此作为成像过程中局部斜率p的约束。通过在反射层斜率附近选取
一定范围的数值作为局部斜率的取值范围,该取值范围选取时应视地下构造的复杂情况进
行选取,构造越复杂选取范围越大,简单近似平层的情况下根据反射层角度在附近(‑π/18,
π/18)角度范围内选取即可,这样可以避免了大量局部斜率p循环所导致的高计算成本,同
时也保证了成像的精度。
上,同时利用延拓波场的传播时间和传播方向来约束成像过程,基于多参数约束的立体成
像条件来实现一次波和多次波的联合成像,具体成像公式为:
传播时间。
实现最终的成像。该过程中不同斜率成分的分解波场互相关增加了传播方向的限制,有效
的避免了串扰噪声的产生,同时在估计的较为准确的地层倾角约束下保证了成像的精度,
避免了大量p循环造成的高计算成本。
统,101道接收,道间距为20m。图3为本发明实施例提取的震源子波与真实子波对比图,参见
图3可知,采用本发明方法基于波动理论利用多次波提取的子波与实际子波基本完全吻合,
这一结果验证了本发明中子波提取方法的可行性。图4为本发明实施例简单层状模型单炮
地震成像结果图,参见图4,是采用不同方法对单炮地震记录成像结果的对比,图4中(a)部
分表示任选脉冲震源利用传统互相关成像方法得到的一次波和多次波联合成像结果,选取
子波与实际子波的相位差异造成了成像界面不准确,子波能量过强造成多次波成像信息能
量较弱,成像结果无法体现利用多次波获取更丰富地下构造信息的优势。同时,不相关波场
互相关还引入了大量的串扰噪声。图4中(b)部分应用本发明方法提取震源子波并利用传统
互相关成像条件得到的一次波和多次波联合成像结果,一次波和多次波联合成像的优势显
而易见,多次波在成像过程中提供了更为丰富的地下构造信息,这也是联合成像相比传统
一次波成像的优势,但是成像假象严重降低了成像的质量。图4中(c)部分应用本发明方法
基于立体成像条件得到的一次波和多次波联合成像结果中不仅一次波和多次波联合成像
优势得到了体现,同时通过时间和局部斜率多参数约束,成像噪声得到了有效的压制,有效
的验证了本发明的可行性。
子波设定为雷克子波,主频为25Hz,采用中间放炮两边接收观测系统,101道接收,道间距为
20m。图6为本发明实施例复杂模型多炮地震成像结果图,参见图6,是采用不同方法对所有
炮地震数据成像结果的对比,图6中(a)部分表示任选脉冲震源子波利用传统互相关成像方
法得到的一次波和多次波联合成像结果,图6中(b)部分应用本发明方法提取子波并利用传
统互相关成像条件得到的一次波和多次波联合成像结果,图6中(c)部分应用本发明方法基
于立体成像条件得到的一次波和多次波联合成像结果。图6(b)的成像结果和图6(a)成像结
果相比,采用真实子波能获取相比传统一次波和多次波联合成像更多的地下构造信息(如
图白色竖线右侧成像范围所示),同时传统互相关成像条件仅通过时间约束成像过程,两图
中都存在严重的串扰噪声(如图中白色箭头处所示),这些假象严重降低了成像的可靠性。
图6(c)的成像结果不仅一次波和多次波联合成像优势得到了体现,同时通过在成像过程中
增加了局部斜率约束,成像噪声得到了有效的压制(白色箭头所示),本发明的有效性再次
得到了验证。
统包括:
第n次迭代后一次波数据,ΔP (z0)=P(z0),P(z0)表示地震数据,A (w)表示第n次迭代
后地表加权算子,w表示角频率。
表示地震波传播时间,w表示角频率。
传播时间。
波保证了成像位置的准确性,能更好的平衡一次波和多次波成像的贡献,然后将提取的震
源子波与原始地震数据同时作为一次波和多次波联合成像的虚拟震源,最后在成像过程中
引入立体成像条件通过时间和局部斜率多参数控制实现成像串扰的压制,提高了一次波和
多次波联合成像的质量,达到了同时利用一次波和多次波获取更为可靠、丰富的构造信息
的目的。
而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说
明即可。
本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不
应理解为对本发明的限制。