轨道式可控震源极化处理方法技术

本发明专利技术属于地震资料处理方法领域，是一种井间地震资料的前期处理方法。由于轨道式可控震源激发的地震波为环形极化波，与线性极化波存在较大区别，为了与常规的线性极化地震资料进行对比分析，需要对环形极化波向线性极化波进行转换，而提出的轨道式可控震源极化处理方法。采用的技术方案主要包括两个步骤：极化分解和震源坐标系统转换。极化分解把正、反转两套不同方向的资料转化为沿Ｘ、Ｙ两个方向传播的线性极化波场。该技术可以完成轨道式可控震源的环形极化方式转换成线性极化方式，得到的地震波场为线性极化地震波，为轨道式震源地震处理提供良好的资料基础。

【技术实现步骤摘要】

本专利技术属于地震资料处理方法领域，是一种井间地震资料的前期处理方法。
技术介绍
在地震野外资料采集中，常规使用的可控震源产生的地震波都是线性极化波。轨道式可控震源是利用电机使离心质体绕某一轴做旋转，电机以不断增大的频率旋转来撞击井壁进行激发地震波，这样震源产生的地震波是环形极化波。该种震源的优点是既能够产生纵波，也能够产生横波，同时较少产生管波等噪音。
技术实现思路
本专利技术的目的是由于轨道式可控震源激发的地震波为环形极化波，与线性极化波存在较大区别，为了与常规的线性极化地震资料进行对比分析，需要对环形极化波向线性极化波进行转换，而提出的。 采用的技术方案主要包括两个步骤极化分解和震源坐标系统转换。 3.1极化分解 设定质量中心的径向位移为u，则正转转动在X、Y轴的投影分别为μcwx和μcwy，翻转转动在X、Y轴的投影分别为μccwx和μccwy，μcwx、μcwy、μccwx和μccwy都与θ角有关。θ角是从x轴转到质量中心的角度，对于产生的震源函数s(f)，离心质体以f＝ω/2π频率转动，则角度θ＝ωt，假设位移u对所有θ具有相同的振幅A，可得到下式 μcwx＝Acos(θ)(1) μcwy＝Asin(θ) (2) μccwx＝Acos(-θ)＝Acos(θ) (3) μccwy＝Asin(-θ)＝-Asin(θ) (4) 通过正转和翻转记录相加和相减，得到线形极化后的μx(ωt)和μy(ωt)的等同数据。即 μx(ωt)＝μcwx+μcwy+μccwx+μccwy ＝2Acos(ωt)(5) μy(ωt)＝μcwx+μcwy-μccwx-μccwy ＝2Asin(ωt)(6) 这样就完成了把正、反转两套不同方向的资料转化为沿X、Y两个方向传播的线性极化波场。但这样分解出的两个线性分量μx(ωt)和μy(ωt)运动存在相位差，必须采用希尔伯特(Hilbert)变换(90°相位转换)来校正相位差。至此，环形极化的正转和反转记录就被转化为等同的线性极化的X和Y方向的震源记录。 3.2震源坐标系统转换 在震源观测中，假设轨道式可控震源信号在X和Y轴上投影分别为sx(f)和sy(f)，则有 Rij(f)＝Tij(f)Sj(f)(7) Tij表示震源方向为i时在j方向观测到的地球介质的反射系数函数，表示j方向激发，i方向接受，其中i＝x，y，z j＝x，y。如果震源传感器正好位于X方向，则记录下的震源信号是Sx(f)。然而，由于震源在井下水平面内的状态是随机的，震源传感器传出的Sx(f)信号与实际有一定的差别。实际工作中不能检测到Sx(f)，而是检测到S(θ)。 S(θ)＝Sx(f)eiθ(8) 根据3.1中计算可以得到，在X和Y方向上接受的地震波Rx(f)和Ry(f)， Rx(f)＝Txx(f)sx(f)+Txy(f)sy(f)(9) Ry(f)＝Tyx(f)Sx(f)+Tyy(f)sy(f)(10) 考虑到sx(f)和sy(f)的关系，有 Rx(f)＝(Txx(f)+iTxy(f))Sx(f) (11) Ry(f)＝(Tyx(f)+iTyy(f))Sy(f) (12) 作自相关，并假设则有 Cx+(f)＝Txx(f)+iTxy(f)(13) Cy+(f)＝Tyx(f)+iTyy(f)(14) 式中x+和y+中的+表示逆时针观测。 同理，对于顺时针观测有 从而有Cx-(f)＝Txx(f)-iTxy(f)(16) Cy-(f)＝Tyx(f)-iTyy(f)(17) 式中x-和y-中的-表示顺时针观测。 总之，有 Txx(f)＝/2(18) Txy(f)＝-i(Cx+(f)-Cx-(f))/2 (19) Tyx(f)＝/2(20) Tyy(f)＝-i(Cy+(f)-Cy-(f))/2 (21) Tzx(f)＝/2(22) Tzy(f)＝-i(Cz+(f)-Cz-(f))/2 (23) 相当于震源为(x′，y′)，接收源为(x，y)，而且下式成立 S(θ)＝Sxcosθ+Sysin(θ)(24) 同时对其他分量也成立，故可得到 满足下式 因而有 关于θ的求取，可根据以下原理我们知道当在Y方向激发地震波时，进行的观测是SH波观测，此时不存在SV波和P波，观测数据中Z分量不应该有直达P波能量，那么使 和 中直达P波能量完全转移到Tzx上的θ角就是震源坐标系统调整的角度。 专利技术的效果该技术可以完成轨道式可控震源的环形极化方式转换成线性极化方式，得到的地震波场为线性极化地震波，为轨道式震源地震处理提供良好的资料基础。 附图说明 图1震源激发方向与能量传播示意图 图2震源传感器处于某一位置示意图 图3震源转到某一位置时信号传感器接收的记录示意图 图4顺时针原始记录(左)、逆时针原始记录(中)与处理后记录(右)对比图 具体实施例方式 下面结合附图作进一步说明。 轨道式可控震源是一种比较特殊的震源设备，其工作方式为把震源沉放到井中某一深度，通过电缆给震源电机加电，使震源内的偏心体从静止状态开始顺时针旋转，随着转速线形提升，震源频率呈线形提高(偏心体的转动频率就是震源发射的声波信号频率)，达到所要求的频率时断电，得到一套正转记录(ClockWise)。以同样方式从静止状态开始逆时针旋转，则得到一套反转记录(Counter ClockWise)这样就完成一个深度点上地震波的激发。 轨道式可控震源资料与常规震源资料存在很大差异，主要表现在正转和翻转两套记录资料的直达波分别呈现正“叠瓦”和反“叠瓦”状(如图4所示)，无法拾取到准确的初至时间，给后续处理中利用初至进行的直达波射线追踪和层析反演带来很大影响，增大拾取误差，降低了层析成像的精度；另一方面，由于直达波影响，得到的反射波场变得比较零乱，根据成像理论，这样的资料很难实现同相叠加。由于野外记录得到正转和反转两种资料，相互耦合在一起，因此，必须进行环形极化到线形极化的转换。 主要包括两个步骤极化分解和震源坐标系统转换。 5.1极化分解 如图2所示假设质量中心的径向位移为u，则正转转动在X、Y轴的投影分别为μcwx和μcwy，翻转转动在X、Y轴的投影分别为μccwx和μccwy，μcwx、μcwy、μccwx和μccwy都与θ角有关。θ角是从x轴转到质量中心的角度，对于产生的震源函数s(f)，离心质体以f＝ω/2π频率转动，则角度θ＝ωt，假设位移u对所有θ具有相同的振幅A，可得到下式 μcwx＝Acos(θ)(1) μcwy＝Asin(θ)(2) μccwx＝Acos(-θ)＝Acos(θ)(3) μccwy＝Asin(-θ)＝-Asin(θ) (4) 通过正转和翻转记录相加和相减，得到线形极化后的μx(ωt)和μy(ωt)的等同数据。即 μx(ωt)＝μcwx+μcwy+μccwx+μccwy ＝2Acos(ωt)(5) 本文档来自技高网...

【技术保护点】
轨道式可控震源极化处理方法，其特征是包括两个步骤：极化分解和震源坐标系统转换，　１．１极化分解　设定质量中心的径向位移为ｕ，则正转转动在Ｘ、Ｙ轴的投影分别为μｃｗｘ和μｃｗｙ，翻转转动在Ｘ、Ｙ轴的投影分别为μｃｃｗｘ和μｃｃｗｙ ，μｃｗｘ、μｃｗｙ、μｃｃｗｘ和μｃｃｗｙ都与θ角有关，θ角是从ｘ轴转到质量中心的角度，对于产生的震源函数ｓ（ｆ），离心质体以ｆ＝ω／２π频率转动，则角度θ＝ωｔ，假设位移ｕ对所有θ具有相同的振幅Ａ，可得到下式　μｃｗｘ＝Ａｃｏｓ（ θ）　（１）　μｃｗｙ＝Ａｓｉｎ（θ）　（２）　μｃｃｗｘ＝Ａｃｏｓ（－θ）＝Ａｃｏｓ（θ）　（３）　μｃｃｗｙ＝Ａｓｉｎ（－θ）＝－Ａｓｉｎ（θ）　（４）　通过正转和翻转记录相加和相减，得到线形极化后的μｘ （ωｔ）和μｙ（ωｔ）的等同数据，即　μｘ（ωｔ）＝μｃｗｘ＋μｃｗｙ＋μｃｃｗｘ＋μｃｃｗｙ　＝２Ａｃｏｓ（ωｔ）　（５）　μｙ（ωｔ）＝μｃｗｘ＋μｃｗｙ－μｃｃｗｘ－μｃｃｗｙ　＝２Ａｓｉｎ（ωｔ）　（６） 　这样就完成了把正、反转两套不同方向的资料转化为沿Ｘ、Ｙ两个方向传播的线性极化波场，但这样分解出的两个线性分量μｘ（ωｔ）和μｙ（ωｔ）运动存在相位差，必须采用希尔伯特变换来校正相位差，至此，环形极化的正转和反转记录就被转化为等同的线 性极化的Ｘ和Ｙ方向的震源记录，　１．２震源坐标系统转换　在震源观测中，假设轨道式可控震源信号在Ｘ和Ｙ轴上投影分别为：ｓ↓［ｘ］（ｆ）和ｓ↓［ｙ］（ｆ），则有　Ｒ↓［ｉｊ］（ｆ）＝Ｔ↓［ｉｊ］（ｆ）Ｓ↓［ｊ］（ｆ）　（７） 　Ｔ↓［ｉｊ］表示震源方向为ｉ时在ｊ方向观测到的地球介质的反射系数函数，表示ｊ方向激发，ｉ方向接受，其中ｉ＝ｘ，ｙ，ｚ　ｊ＝ｘ，ｙ。如果震源传感器正好位于Ｘ方向，则记录下的震源信号是Ｓ↓［ｘ］（ｆ）。然而，由于震源在井下水平面内的状态 是随机的，震源传感器传出的Ｓ↓［ｘ］（ｆ）信号与实际有一定的差别，实际工作中不能检测到Ｓ↓［ｘ］（ｆ），而是检测到Ｓ（θ），　Ｓ（θ）＝Ｓ↓［ｘ］（ｆ）ｅ↑［ｉθ］　（８）　根据１．１中计算可以得到，在Ｘ和Ｙ方向上接受的地震波 Ｒ↓［ｘ］（ｆ）和Ｒ↓［ｙ］（ｆ），　Ｒ↓［ｘ］（ｆ）＝Ｔ↓［ｘｘ］（ｆ）ｓ↓［ｘ］（ｆ）＋Ｔ↓［ｘｙ］（ｆ）ｓ↓［ｙ］（ｆ）...

【技术特征摘要】

【专利技术属性】
技术研发人员：孔庆丰，乔玉雷，姚忠瑞，
申请(专利权)人：孔庆丰，
类型：发明
国别省市：37[中国|山东]