CN102621581A - 一种vsp波场分离方法 - Google Patents
一种vsp波场分离方法 Download PDFInfo
- Publication number
- CN102621581A CN102621581A CN2012100998146A CN201210099814A CN102621581A CN 102621581 A CN102621581 A CN 102621581A CN 2012100998146 A CN2012100998146 A CN 2012100998146A CN 201210099814 A CN201210099814 A CN 201210099814A CN 102621581 A CN102621581 A CN 102621581A
- Authority
- CN
- China
- Prior art keywords
- wave
- ripple
- cos
- sin
- wave field
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000926 separation method Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 36
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 claims abstract description 15
- 230000010287 polarization Effects 0.000 claims description 42
- 230000001052 transient effect Effects 0.000 claims description 5
- 238000011144 upstream manufacturing Methods 0.000 claims description 5
- 230000015572 biosynthetic process Effects 0.000 abstract 1
- 230000007547 defect Effects 0.000 abstract 1
- 238000003786 synthesis reaction Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000002950 deficient Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 230000009466 transformation Effects 0.000 description 2
- 229920000535 Tan II Polymers 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种地震资料处理处理技术,尤其是一种基于仿射坐标系的VSP资料波场分离处理领域,具体方法过程是:(A)进行水平方向的x、y分量合成;(A)由水平分量x′方向与垂直分量z方向分离P波与SV波波场;(C)在新坐标系中假定对于不同深度H的上行P-P波、上行P-SV波偏振角α(H)。本发明同现有技术相比,能对地下任意一个反射点,根据该点的上行反射波方向建立仿射坐标系进行波场分离,不同点的坐标系不同。该方法克服了传统直角坐标系波场不能完全分离的缺陷,完全以所需要波场为基础,能够实现波场的最优分离。
Description
[技术领域]
本发明涉及一种地震资料处理处理技术,尤其是一种基于仿射坐标系的VSP资料波场分离处理领域。
[背景技术]
VSP资料接收的波场非常丰富,包括直达波、上行P-P波、下行P-P波、上行P-SV波、下行P-SV波等波场,这些波场通过三分量检波器能够得到完整的记录,但是如何将他们分离成为VSP资料处理的关键。通常最有效的方法就是三分量极化滤波实现波场初步分离。传统的极化滤波方法偏振角的计算是基于直达波的传播方向而确定,因此对于受观测角度影响,除了正交方向外,其他方向上的波不能实现完全分离,特别是上行P-P反射波、上行P-SV反射波等,存在剩余能量,分离效果较差。
[发明内容]
本发明为了克服现有技术的不足,提供了一种基于仿射坐标来实现VSP波场分离的方法。
为了实现上述目的,本发明设计了一种VSP波场分离方法,其特征在于:根据上行反射波的方向构建新的坐标系进行波场分离;
具体方法过程是:
(A)进行水平方向的x、y分量合成,由于不同深度地层反射波到达检波器时在水平方向的偏振角不受影响,其计算公式为:
(B)由水平分量x′方向与垂直分量z方向分离P波与SV波波场,具体操作为,在上行波场几何关系图中设一点O点为反射点,上行P波的偏振角为α,上行SV波的偏振角为β;根据反射波的偏振角构建仿射坐标系,设上行P-P波方向单位矢量为u,上行P-SV波方向单位矢量为v,则建立新坐标系(O,u,v);
(C)计算不同深度H的上行P-P波、上行P-SV波偏振角α(H),β(H)。具体操作为,在一个记录点的水平分量与垂直分量沿层选取上行P-P波或上行P-SV波的时窗,用传统方法计算偏振角α(H)与β(H)。
(D)在新坐标系中对于不同深度H的上行P-P波、上行P-SV波偏振角α(H),β(H),有如下公式:
(E)解上述方程组可得:
即实现了上行P-P波、上行P-SV波的波场分离。
所述的步骤(C)中,其计算偏振角α(H)与β(H)的方法包括三角关系法、矢端曲线和能量准则法或能量加权的瞬时方位直方图法。
本发明同现有技术相比,能对地下任意一个反射点,根据该点的上行反射波方向建立仿射坐标系进行波场分离,不同点的坐标系不同。该方法克服了传统直角坐标系波场不能完全分离的缺陷,完全以所需要波场为基础,能够实现波场的最优分离。
[附图说明]
图1为现有技术基于直角坐标系的三分量VSP观测示意图。
图2为本发明基于空间仿射坐标系VSP上行波反射波几何关系示意图。
图3为采用传统方法波场分离的上行P-P波。
图4为采用传统方法波场分离的上行P-SV波。
图5为采用仿射坐标系波场分离的上行P-P波。
图6为采用仿射坐标系波场分离的上行P-SV波。
[具体实施方式]
下面结合附图对本发明做进一步描述。
如图1所示的三分量VSP观测示意图,传统方法基于直角坐标系,利用直达波的偏振角进行波场分离;
假设从P波震源传到井下第一个深度点上的直达P波,其质点运动方向和波的传播方向一致,都在由震源和井确定的平面内,则这种直达P波的偏振方向是线性的,它在水平面内的投影也是直线。
具体方法过程是:
(A)用直达波偏振方向在水平面内的投影作为参考,测出三分量检波器观测时水平分量的相对方位并可将观测到的水平分量信号转换到以直达P波偏振方向在水平面内投影方向一致,记为Hp方向的参考坐标系;
其三分量检波器的水平分量(x,y)转换到以Hp方向为参考的一致坐标系(x′,y′)的转换公式为:
式中θ为x与Hp的夹角,即直达P波的偏振角,通过上述公式确定线性偏振信号方位。
通常在水平面内一般采用三角关系法、矢端曲线和能量准则法或能量加权的瞬时方位直方图法来计算偏振角。
(1)三角关系法。
利用简单的三角关系进行计算,偏振角θ的计算公式为
θ=arctg(y/x)
(2)矢端曲线和能量准则法。
矢端曲线是一种表示直达波水平分量取向的直观图示方法。将x,y两个方向上一定时窗内的样点值按坐标点在图上,连成曲线,就画出该矢量曲线。幅度最大的方向即为偏振角的方向,即可估计出偏振角θ。
偏振角θ也可由能量准则解析的求出。设能量的表达式为:
能量最大时的方位角即为偏振角,使能量取极大地必要条件是:
从而可以求出:
可解出偏振角θ=θmax
(3)能量加权的瞬时方位直方图法
矢端曲线也可以用瞬时能量Ri和瞬时方位θi写出,即
(B)利用旋转后的水平分量Hp和垂直分量z计算P波偏振方向与铅垂线的夹角θv,再通过坐标变换求出沿P波偏振方向的分量x″方向,以及与P波偏振方向垂直的分量y″方向,即公式为:
即完成了VSP的P-P波、P-SV波波场分离工作。
传统方法是在直角坐标系中实现,坐标变换为直角坐标变换。而实际资料中不同深度地层的反射波到达检波器时θ是不同的,因而在直角坐标变换后的投影也不同,所以无法实现完全的波场分离。
如图2所示的VSP上行反射波场几何关系图,基于仿射坐标系的波场分离方法根据上行反射波的方向构建新的坐标系,从而实现上行波场的最优分离。
具体方法过程是:
(A)首先进行水平方向的x、y分量合成,由于不同深度地层反射波到达检波器时在水平方向的偏振角不受影响,其方法与传统方法一致,其计算公式为:
(B)由水平分量x′方向与垂直分量z方向分离P波与SV波波场,具体操作为,在上行波场几何关系图中设一点O点为反射点,上行P波的偏振角为α,上行SV波的偏振角为β;根据反射波的偏振角构建仿射坐标系,设上行P-P波方向单位矢量为u,上行P-SV波方向单位矢量为v,则建立新坐标系(O,u,v);
(C)计算不同深度H的上行P-P波、上行P-SV波偏振角α(H),β(H)。具体操作为,在一个记录点的水平分量与垂直分量沿层选取上行P-P波或上行P-SV波的时窗,采用与水平分量定位相同方法计算偏振角α(H)与β(H)。
(D)在新坐标系中对于不同深度H的上行P-P波、上行P-SV波偏振角α(H),β(H),有如下公式:
(E)解上述方程组可得:
即实现了上行P-P波、上行P-SV波的波场分离。
传统波场分离方法分离的上行P-P波、上行P-SV波如图3、图4所示。基于仿射坐标系的波场分离方法分离的上行P-P波、上行P-SV波如图5、图6所示。可以看出,传统波场分离方法的上行P-P波与上行P-SV波无法完全分开,分离效果较差,而采用基于仿射坐标系的波场分离方法分离后的上行波场较为干净,分离效果较好。
本发明能够对地下任意一个反射点,根据该点的上行反射波方向建立仿射坐标系,不同点的坐标系不同。该方法克服了传统直角坐标系波场不能完全分离的缺陷,完全以所需要波场为基础,能够实现波场的最优分离。
Claims (2)
1.一种VSP波场分离方法,其特征在于:根据上行反射波的方向构建新的坐标系进行波场分离;
具体方法过程是:
(A)首先进行水平方向的x、y分量合成,由于不同深度地层反射波到达检波器时在水平方向的偏振角不受影响,其计算公式为:
(B)由水平分量x′方向与垂直分量z方向分离P波与SV波波场,具体操作为,在上行波场几何关系图中设一点O点为反射点,上行P波的偏振角为α,上行SV波的偏振角为β;根据反射波的偏振角构建仿射坐标系,设上行P-P波方向单位矢量为u,上行P-SV波方向单位矢量为v,则建立新坐标系(O,u,v);
(C)计算不同深度H的上行P-P波、上行P-SV波偏振角α(H),β(H)。具体操作为,在一个记录点的水平分量与垂直分量沿层选取上行P-P波或上行P-SV波的时窗,并计算偏振角α(H)与β(H);
(D)在新坐标系中对于不同深度H的上行P-P波、上行P-SV波偏振角α(H),β(H),有如下公式:
(E)解上述方程组可得:
即实现了上行P-P波、上行P-SV波的波场分离。
2.根据权利要求1所述的一种VSP波场分离方法,其特征在于:所述的步骤(C)中,其计算偏振角α(H)与β(H)的方法包括三角关系法、矢端曲线和能量准则法或能量加权的瞬时方位直方图法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210099814.6A CN102621581B (zh) | 2012-04-06 | 2012-04-06 | 一种vsp波场分离方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210099814.6A CN102621581B (zh) | 2012-04-06 | 2012-04-06 | 一种vsp波场分离方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102621581A true CN102621581A (zh) | 2012-08-01 |
CN102621581B CN102621581B (zh) | 2014-08-27 |
Family
ID=46561606
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210099814.6A Expired - Fee Related CN102621581B (zh) | 2012-04-06 | 2012-04-06 | 一种vsp波场分离方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102621581B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104316960A (zh) * | 2014-10-29 | 2015-01-28 | 中国石油天然气集团公司 | 一种基于vsp的储层裂缝预测方法及系统 |
CN106353815A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 一种垂直地震剖面剩余极化角能量补偿方法 |
CN107664770A (zh) * | 2016-11-09 | 2018-02-06 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 垂直地震剖面数据逐点矢量合成方法 |
CN109655899A (zh) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | Vsp地震数据瞬时角度谱计算方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4935903A (en) * | 1989-05-30 | 1990-06-19 | Halliburton Geophysical Services, Inc. | Reinforcement of surface seismic wavefields |
CN101246218B (zh) * | 2007-02-15 | 2010-06-02 | 中国石油化工股份有限公司 | 一种三分量vsp波场分离方法 |
-
2012
- 2012-04-06 CN CN201210099814.6A patent/CN102621581B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4935903A (en) * | 1989-05-30 | 1990-06-19 | Halliburton Geophysical Services, Inc. | Reinforcement of surface seismic wavefields |
CN101246218B (zh) * | 2007-02-15 | 2010-06-02 | 中国石油化工股份有限公司 | 一种三分量vsp波场分离方法 |
Non-Patent Citations (1)
Title |
---|
王成礼等: "VSP波场极化特征分析", 《油气地球物理》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104316960A (zh) * | 2014-10-29 | 2015-01-28 | 中国石油天然气集团公司 | 一种基于vsp的储层裂缝预测方法及系统 |
CN106353815A (zh) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | 一种垂直地震剖面剩余极化角能量补偿方法 |
CN106353815B (zh) * | 2015-07-17 | 2020-03-31 | 中国石油化工股份有限公司 | 一种垂直地震剖面剩余极化角能量补偿方法 |
CN107664770A (zh) * | 2016-11-09 | 2018-02-06 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 垂直地震剖面数据逐点矢量合成方法 |
CN107664770B (zh) * | 2016-11-09 | 2019-02-15 | 中国石油天然气集团有限公司 | 垂直地震剖面数据逐点矢量合成方法 |
CN109655899A (zh) * | 2017-10-11 | 2019-04-19 | 中国石油化工股份有限公司 | Vsp地震数据瞬时角度谱计算方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN102621581B (zh) | 2014-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102012521B (zh) | 一种地震储层预测中叠前裂缝的检测方法 | |
CN102692646B (zh) | 一种三维三分量矢量波场分离的方法和系统 | |
Levin et al. | Crust and upper mantle of Kamchatka from teleseismic receiver functions | |
CN103064111B (zh) | 一种基于形态滤波的微地震事件识别方法 | |
CN102053261B (zh) | 一种地震数据处理方法 | |
CN103048655B (zh) | 天波超视距雷达频域超分辨微多径测高方法 | |
CN107132575B (zh) | 基于横波极化分析预测裂缝方位角的方法 | |
CN105093292A (zh) | 一种地震成像的数据处理方法和装置 | |
CN105785440B (zh) | 一种矿井槽波双分量地震信号频散曲线提取方法 | |
CN102621581B (zh) | 一种vsp波场分离方法 | |
CN107290722B (zh) | 微震源的定位方法和装置 | |
CN101893720A (zh) | 一种多波的波场分离与合成的方法和系统 | |
CN104865598A (zh) | 一种精细岩溶古地貌恢复方法 | |
CN105093274A (zh) | 一种水力压裂裂缝震源机制的反演方法及系统 | |
CN102073064B (zh) | 一种利用相位信息提高速度谱分辨率的方法 | |
CN103777245B (zh) | 基于地震资料的油气成藏条件定量评价方法 | |
CN102879817A (zh) | 基于地面地震数据获取地下裂缝信息的控制方法 | |
HAO et al. | A study on the extension of fault zones in Yellow Sea and its adjacent areas based on gravity data | |
CN102305940B (zh) | 流体因子提取方法 | |
CN103941287A (zh) | 一种基于平面导航的快速三维断层解释方法 | |
CN104570086B (zh) | 一种在共偏移距共方位角域进行叠前裂缝预测的方法 | |
CN104360384B (zh) | 基于纵横波能量自动扫描的微地震事件定位方法及装置 | |
CN105425287A (zh) | 地震波的叠前分离方法 | |
CN103105603A (zh) | 一种x波段测波雷达海流反演预处理方法 | |
CN109143362A (zh) | 基于共散射角道集的散射波分离方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20140827 |
|
CF01 | Termination of patent right due to non-payment of annual fee |