CN103984024A - 一种三分量检波器水平分量数据自动校正方法 - Google Patents
一种三分量检波器水平分量数据自动校正方法 Download PDFInfo
- Publication number
- CN103984024A CN103984024A CN201410226691.7A CN201410226691A CN103984024A CN 103984024 A CN103984024 A CN 103984024A CN 201410226691 A CN201410226691 A CN 201410226691A CN 103984024 A CN103984024 A CN 103984024A
- Authority
- CN
- China
- Prior art keywords
- component
- horizontal
- data
- seismometer
- orientation
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种三分量检波器水平分量数据自动校正方法,属于地震勘探数据处理领域,其校正步骤为:设置锤击震源的方位角,作为实际水平X分量;实时测量检波器的姿态角,并自动记录三分量检波器水平X分量和水平Y分量的方位角;自动计算所测量的三分量检波器水平X分量方位到锤击方位的角度差;利用该角度差将测量的水平X分量和水平Y分量地震数据重新投影到锤击方位上;将投影后的水平X分量和水平Y分量地震数据重新叠加合成,获得实际水平X分量的地震数据。同现有技术相比,本发明可以实现自动实时校正,且校正质量高,操作效率高,能够从原始数据本身提高地震数据后期处理精度。
Description
技术领域
本发明涉及三分量地震勘探的数据处理领域,特别涉及一种自动实时、校正质量和效率高的三分量检波器水平分量数据自动校正方法。
背景技术
作为复杂地形及海上油气勘探的重要物探方法之一,三分量地震勘探由于单次的施工即可获得丰富的全波信息而被广泛应用,相对成本低、破坏性小且效率较高。三分量地震勘探主要采用三分量检波器来获取包括横波、纵波、转换波等在内多种类型地震数据,从而提高地震勘探的效率,同时保证后期数据处理有更多选择组合,以提高原始地震数据的质量。由于三分量地震勘探是利用三分量检波器来获取全波信息,在实际施工过程中,并不能保证将检波器完全垂直放置,且各个检波器之间的水平分量方位也存在较大差异,因此相对于同一震源激发位置,各个检波器接收到的数据同样会存在差异,同时也势必导致各个分量之间的各种波形存在相互干扰的现象,即各检波器分量不可避免地会存在不同类型的干扰,这种现象也会进一步影响后期地震数据的成果分析。
总体来说,由于实际现场施工限制,通常各个三分量检波器相互之间的参考方向不一致,因此相对于同一震源激发位置,各个三分量检波器接收的振动信号也存在较大差异,由此导致地震剖面波形图初至混乱,各分量兴趣子波能量强弱分布不均,且对震源系统的适应性不高,特别是用于横波勘探时,与震源敲击方向同向或垂直方向上的横波分量分布规律不明显,甚至导致其不具有横波波形特征,正反向激发震源极性不相反,横纵波波组速度差异不明显,极大影响地震数据处理精度。
针对这一问题,当前所采用的主流方法是多次激发测试炮的方式。即通过合理选择震源激发点随机放炮产生振动信号,根据各个检波器接收到的数据差异逐渐递推来确定检波器之间的方位角信息,这种方式在一定程度上提高了原始数据的精度。但是由于测试炮的数据本身就存在干扰误差,因此利用有误差的数据所递推的方位角信息同样存在误差,且这一误差是原始数据本身决定的,因此是随机性的。
发明内容
针对上述问题,本发明的目的是提供一种精度高、处理效率高的三分量检波器水平分量数据自动校正方法。
为了实现上述目的,本发明采用如下技术方案:一种三分量检波器水平分量数据自动校正方法,包括以下步骤:
1)设置震源激发点的锤击方位,将其作为最终校正完成后的实际水平X分量方位;
2)实时获取三分量检波器的姿态角,并自动记录三分量检波器水平X分量和水平Y分量的方位角;
3)计算所测量的三分量检波器水平X分量方位到锤击方位的角度差;
4)利用所述角度差将测量的水平X分量和水平Y分量地震数据重新投影到锤击方位上;
5)将投影后的水平X分量和水平Y分量地震数据重新叠加合成,获得实际水平X分量的地震数据;
6)由于水平X分量和水平Y分量是垂直的关系,利用上述同样步骤1)~5)获得实际水平Y分量的地震数据。
本方法校正对象为三分量检波器水平X分量和水平Y分量数据。
其中校正对象即地震数据均为是严格正交的三分量模拟检波器或数字检波器获取的三分量地震数据;校正前的地震数据分别由水平X分量和水平Y分量检波器测量来获取,校正后的X分量数据为锤击方位的数据,校正后的Y分量数据为与锤击方位垂直方位的数据,水平X分量和水平Y分量的数据均是由测量的数据投影后合成而来的;用于测量检波器方位角的姿态传感器内置安装于检波器内部,实时准确获取水平X分量和水平Y分量的方位角。
作为优选的技术方案:利用内置安装于三分量检波器内部的姿态传感器实时获取三分量检波器的姿态角。
作为优选的技术方案:所述实际水平X分量和实际水平Y分量的数据均由两组数据合成,分别是测量的水平X分量和水平Y分量地震数据进行投影后的数据。
作为优选的技术方案:所述校正前的数据分别为水平X分量和水平Y分量所测量的数据,校正后的水平X分量数据为锤击方位的数据,校正后的水平Y分量数据为与锤击方位垂直方位的数据。
作为优选的技术方案:所述地震数据为三个严格正交的三分量地震数据。
作为优选的技术方案:所述三分量检波器为三分量模拟检波器或三分量数字检波器。
本发明的有益效果是:本发明通过自动获取倾角方位信息,主动结合震源激发方位对水平分量数据进行重新旋转投影合成,实现对实际水平分量数据重新校正。本发明通过以上过程实现三分量地震数据水平分量数据的自动实时校正,改进水平分量地震波形同相轴混乱及初至拾取困难的问题,且具有校正质量高、操作效率高的特点,能够从原始数据本身提高地震数据后期处理精度。
附图说明
图1为本发明校正流程图;
图2为三分量检波器各分量分布图;
图3为数据校正原理图;
图4为数据校正效果对比图。
具体实施方式
为了使本发明更容易被理解,下面结合附图和具体实施方式对本发明做进一步说明。
参阅图1,一种三分量检波器水平分量数据自动校正方法,主要包括以下步骤:设置锤击方位,即设置震源激发点的锤击方位,将其作为最终校正完成后的实际水平X分量方位;自动获取检波器姿态角,即利用内置安装于三分量检波器内部的姿态传感器实时准确获取检波器的姿态角,并自动记录三分量检波器水平X分量和水平Y分量的方位角;旋转角差计算,即自动计算所测量的三分量检波器水平X分量方位到锤击方位的角度差;水平分量数据重新投影,即利用该角度差将测量的水平X分量和水平Y分量地震数据重新投影到锤击方位上,即实际水平X分量;水平分量数据重新叠加,即将投影后的水平X分量和水平Y分量地震数据重新叠加合成,获得实际水平X分量的地震数据;通过以上过程来完成水平X分量地震数据的校正,由于水平X分量和水平Y分量是垂直的关系,可以利用上述同样步骤获得实际水平Y分量的地震数据。
参阅图1,其中所设置的锤击方位,是在选取震源激发位置时已经确定好的,即震源激发位置与正北方向的夹角,因此在进行自动校正之前可以实现设置;自动获取检波器姿态角,是由集成在检波器内部的姿态传感器自动检测,并由软件自动记录检波器的方位角信息,为数据校正提供参考角;其中的旋转角差计算,主要是计算自动获取的方位角与锤击方位角的角度差;其中的水平分量数据重新投影,是将测量的水平X分量和水平Y分量数据结合旋转角度差进行重新投影,最终投影到实际的水平X分量和水平Y分量上;其中的水平分量数据重新叠加,由于重新投影后的水平X分量和水平Y分量数据均是分别来自于测量的水平X分量和水平Y分量数据,因此需要将重新投影后的实际水平X分量和水平Y分量数据进行重新叠加。上述过程均通过编程软件自动执行完成。
参阅图2,本发明涉及的检波器为三分量检波器,既可以是模拟检波器也可以是数字检波器。其三个分量分别是水平X分量、水平Y分量、Z分量,三个分量为正交关系,即三个分量之间均相互垂直。这种情况下,在随机任意的震源激发位置,对于水平分量来说,无论在任何方位激发震源,水平X分量和水平Y分量均能接收到信号,但是两者的能量会存在差异,甚至会出现其中一个分量没有信号的情况(锤击方位与其中一个分量的方位相同),这样就会出现原本属于水平X分量的数据被分配到水平Y分量的情况,而且对于水平Y分量来说,也引入了干扰。同样这种情况对于水平Y分量也是存在的。因此在校正时实际水平X分量的数据是来自于所测量的水平X分量和水平Y分量,实际水平Y分量同样如此。
参阅图3,水平分量数据投影叠加的具体过程如下:设姿态传感器测量的当前检波器姿态值为S测(X测,Y测,Z测),即所测量检波器方位,设X测与正北方向夹角为β。且由于检波器的三个分量是严格正交的,因此水平X分量与水平Y分量是相互垂直的。
这里假设实际敲击方向如图3中的X实,其与正北方向的夹角为α,α可以在0~360度范围内任意选取,而且是可以提前测量的,即是已知的。
在进行校正时,针对水平X分量,需要将所测量的水平X测分量与Y测分量的数据进行一定角度的旋转,由于校正的数据来自于测量的水平X分量和测量的水平Y分量,因此这一过程分为两部分:即从测量的水平X分量X测方位旋转投影到X实方位(锤击方位),以及测量的水平Y分量旋转到X实方位,进而对数据进行重新投影。Y分量也同样按此方法处理。
如图3中从测量的水平X分量X测方位旋转投影到X实方位,需要旋转τ=α-β的角度。
在X实方位进行投影时,针对X测需取X测*cos(α-β),而对Y测由于需要多旋转90度(水平X分量和水平Y分量是相互垂直的),因此需取Y测*cos(α-90-β),即-Y测*sin(α-β)。
则X实=X测*cos(α-β)+Y测*cos(α-β-90);
即X实=X测*cos(α-β)-Y测*sin(α-β);
以此类推,则Y实=X测*cos(α-β+90)+Y测*cos(α-β);
即Y实=X测*sin(α-β)+Y测*cos(α-β);
如锤击方位为α=245,则按如下公式进行计算:
X实=X测*cos(245-β)-Y测*sin(245-β);
Y实=X测*sin(245-β)+Y测*cos(245-β)。
其中β为自动测量的检波器水平X分量方位(与正北方向夹角)。
水平X分量校正前后的对比效果如图4所示,从图4中可以看出,未校正前的各道地震波形分布杂乱,并没有按照由浅至深(从左向右)的波形走时规律分布,各道波形的起跳点也没有按照规律逐级分布,造成初至无法拾取;而校正后的波形分布非常均匀,且严格按照由浅至深(从左向右)的波形走时规律分布,各道波形的起跳点也是由浅至深(从左向右)的波形走时规律分布,因此能清晰的辨认和拾取初至。
Claims (6)
1.一种三分量检波器水平分量数据自动校正方法,其特征在于:包括以下步骤:
1)设置震源激发点的锤击方位,将其作为最终校正完成后的实际水平X分量方位;
2)实时获取三分量检波器的姿态角,并自动记录三分量检波器水平X分量和水平Y分量的方位角;
3)计算所测量的三分量检波器水平X分量方位到锤击方位的角度差;
4)利用所述角度差将测量的水平X分量和水平Y分量地震数据重新投影到锤击方位上;
5)将投影后的水平X分量和水平Y分量地震数据重新叠加合成,获得实际水平X分量的地震数据;
6)利用上述同样步骤1)~5)获得实际水平Y分量的地震数据。
2.根据权利要求1所述一种三分量检波器水平分量数据自动校正方法,其特征在于:利用内置安装于三分量检波器内部的姿态传感器实时获取三分量检波器的姿态角。
3.根据权利要求1所述一种三分量检波器水平分量数据自动校正方法,其特征在于:所述实际水平X分量和实际水平Y分量的数据均由两组数据合成,分别是测量的水平X分量和水平Y分量地震数据进行投影后的数据。
4.根据权利要求1所述一种三分量检波器水平分量数据自动校正方法,其特征在于:所述校正前的数据分别为水平X分量和水平Y分量所测量的数据,校正后的水平X分量数据为锤击方位的数据,校正后的水平Y分量数据为与锤击方位垂直方位的数据。
5.根据权利要求1所述一种三分量检波器水平分量数据自动校正方法,其特征在于:所述地震数据为三个严格正交的三分量地震数据。
6.根据权利要求1至5中任意一项所述的一种三分量检波器水平分量数据自动校正方法,其特征在于:所述三分量检波器为三分量模拟检波器或三分量数字检波器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410226691.7A CN103984024B (zh) | 2014-05-26 | 2014-05-26 | 一种三分量检波器水平分量数据自动校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410226691.7A CN103984024B (zh) | 2014-05-26 | 2014-05-26 | 一种三分量检波器水平分量数据自动校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103984024A true CN103984024A (zh) | 2014-08-13 |
CN103984024B CN103984024B (zh) | 2017-03-01 |
Family
ID=51276066
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410226691.7A Expired - Fee Related CN103984024B (zh) | 2014-05-26 | 2014-05-26 | 一种三分量检波器水平分量数据自动校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103984024B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155685A (zh) * | 2014-08-15 | 2014-11-19 | 中国科学院地质与地球物理研究所 | 航空母机型地震采集站布设和数据采集系统 |
CN104182651A (zh) * | 2014-09-12 | 2014-12-03 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 用于三分量检波器接收的微地震事件方位角自动质控方法 |
CN104570111A (zh) * | 2015-01-21 | 2015-04-29 | 中国矿业大学(北京) | 共姿态道集方位角分析和校正方法及装置 |
CN106226558A (zh) * | 2016-09-23 | 2016-12-14 | 中国地震局工程力学研究所 | 一种检测力平衡加速度传感器极性的方法及装置 |
CN107300716A (zh) * | 2017-06-30 | 2017-10-27 | 黄河勘测规划设计有限公司 | 用于三维sh波或多波地震勘探的智能检波器 |
CN107561590A (zh) * | 2017-08-11 | 2018-01-09 | 武汉中岩科技有限公司 | 一种土层剪切波测试装置及其方法 |
CN107607990A (zh) * | 2017-08-07 | 2018-01-19 | 中国石油天然气集团公司 | 三分量检波器水平分量的方向检测方法和装置 |
CN110687606A (zh) * | 2019-10-25 | 2020-01-14 | 长安大学 | 一种海底节点地震仪三分量定向校正方法 |
CN111257940A (zh) * | 2020-02-14 | 2020-06-09 | 中国科学院地质与地球物理研究所 | 一种组合式海底地震仪姿态自动校正装置及方法 |
CN111257941A (zh) * | 2020-02-18 | 2020-06-09 | 中国科学院地质与地球物理研究所 | 一种组合式海底地震仪方位角自动识别装置及方法 |
CN111856617A (zh) * | 2020-05-12 | 2020-10-30 | 吉林大学 | 地震检波器的自调校方法、装置、介质及电子设备 |
CN112444842A (zh) * | 2019-08-27 | 2021-03-05 | 中国石油天然气集团有限公司 | 横波可控震源激发横波的角度确定方法、装置及系统 |
CN112464438A (zh) * | 2020-10-30 | 2021-03-09 | 中国石油天然气集团有限公司 | 横波极性判断方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4803669A (en) * | 1986-12-19 | 1989-02-07 | Atlantic Richfield Company | Three-component three-dimensional seismic data acquisition |
US6205403B1 (en) * | 1999-04-30 | 2001-03-20 | Baker Hughes Incorporated | Method for determining horizontal geophone orientation in ocean bottom cables |
-
2014
- 2014-05-26 CN CN201410226691.7A patent/CN103984024B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4803669A (en) * | 1986-12-19 | 1989-02-07 | Atlantic Richfield Company | Three-component three-dimensional seismic data acquisition |
US6205403B1 (en) * | 1999-04-30 | 2001-03-20 | Baker Hughes Incorporated | Method for determining horizontal geophone orientation in ocean bottom cables |
Non-Patent Citations (6)
Title |
---|
HENRY C.BLAND ET AL.: "《Geophone orientation, location, and polarity checking for 3-C seismic surveys》", 《CREWES RESEARCH REPORT》 * |
傅旦丹等: "《海上多分量地震勘探水平分量方位校正》", 《石油物探》 * |
周建新等: "《海上多波地震勘探检波器方位误差校正》", 《中国海上油气(地质)》 * |
李怀良: "《复杂山地多波宽频带地震数据采集关键技术研究》", 《中国博士学位论文全文数据库 基础科学辑》 * |
王立娟等: "《斜井三维VSP方位校正、波场分离处理方法及应用》", 《内蒙古石油化工》 * |
郑儒等: "《VSP 三分量检波器定向》", 《油气地球物理》 * |
Cited By (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155685A (zh) * | 2014-08-15 | 2014-11-19 | 中国科学院地质与地球物理研究所 | 航空母机型地震采集站布设和数据采集系统 |
CN104182651A (zh) * | 2014-09-12 | 2014-12-03 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 用于三分量检波器接收的微地震事件方位角自动质控方法 |
CN104182651B (zh) * | 2014-09-12 | 2017-03-15 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 用于三分量检波器接收的微地震事件方位角自动质控方法 |
CN104570111A (zh) * | 2015-01-21 | 2015-04-29 | 中国矿业大学(北京) | 共姿态道集方位角分析和校正方法及装置 |
CN104570111B (zh) * | 2015-01-21 | 2016-03-02 | 中国矿业大学(北京) | 共姿态道集方位角分析和校正方法及装置 |
WO2016115763A1 (zh) * | 2015-01-21 | 2016-07-28 | 中国矿业大学(北京) | 共姿态道集方位角分析和校正方法及装置 |
CN106226558A (zh) * | 2016-09-23 | 2016-12-14 | 中国地震局工程力学研究所 | 一种检测力平衡加速度传感器极性的方法及装置 |
CN107300716B (zh) * | 2017-06-30 | 2018-11-23 | 黄河勘测规划设计有限公司 | 用于三维sh波或多波地震勘探的智能检波器 |
CN107300716A (zh) * | 2017-06-30 | 2017-10-27 | 黄河勘测规划设计有限公司 | 用于三维sh波或多波地震勘探的智能检波器 |
CN107607990A (zh) * | 2017-08-07 | 2018-01-19 | 中国石油天然气集团公司 | 三分量检波器水平分量的方向检测方法和装置 |
CN107607990B (zh) * | 2017-08-07 | 2019-09-10 | 中国石油天然气集团公司 | 三分量检波器水平分量的方向检测方法和装置 |
CN107561590A (zh) * | 2017-08-11 | 2018-01-09 | 武汉中岩科技有限公司 | 一种土层剪切波测试装置及其方法 |
CN112444842A (zh) * | 2019-08-27 | 2021-03-05 | 中国石油天然气集团有限公司 | 横波可控震源激发横波的角度确定方法、装置及系统 |
CN110687606A (zh) * | 2019-10-25 | 2020-01-14 | 长安大学 | 一种海底节点地震仪三分量定向校正方法 |
CN110687606B (zh) * | 2019-10-25 | 2021-04-20 | 长安大学 | 一种海底节点地震仪三分量定向校正方法 |
CN111257940A (zh) * | 2020-02-14 | 2020-06-09 | 中国科学院地质与地球物理研究所 | 一种组合式海底地震仪姿态自动校正装置及方法 |
CN111257940B (zh) * | 2020-02-14 | 2020-09-04 | 中国科学院地质与地球物理研究所 | 一种组合式海底地震仪姿态自动校正装置及方法 |
CN111257941A (zh) * | 2020-02-18 | 2020-06-09 | 中国科学院地质与地球物理研究所 | 一种组合式海底地震仪方位角自动识别装置及方法 |
CN111257941B (zh) * | 2020-02-18 | 2020-09-04 | 中国科学院地质与地球物理研究所 | 一种组合式海底地震仪方位角自动识别装置及方法 |
CN111856617A (zh) * | 2020-05-12 | 2020-10-30 | 吉林大学 | 地震检波器的自调校方法、装置、介质及电子设备 |
CN112464438A (zh) * | 2020-10-30 | 2021-03-09 | 中国石油天然气集团有限公司 | 横波极性判断方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103984024B (zh) | 2017-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103984024A (zh) | 一种三分量检波器水平分量数据自动校正方法 | |
CN102967882B (zh) | 地层的层速度模型的建模方法 | |
CN104459797B (zh) | 一种井中微地震事件识别和拾取的方法 | |
CN101630016B (zh) | 一种提高垂直地震剖面成像质量的方法 | |
CN102565853A (zh) | 对地质构造的速度模型进行建模的方法 | |
Ao et al. | The correction of shot and OBS position in the 3D seismic experiment of the SW Indian Ocean Ridge | |
CN104570111B (zh) | 共姿态道集方位角分析和校正方法及装置 | |
CN101609160B (zh) | 一种抗背景噪音强的地震数据初至波自动拾取方法 | |
EP2920615A1 (en) | Process for creating image gathers | |
CN104570102B (zh) | 一种近地表速度模型与中深层速度模型的融合方法 | |
CN101915938B (zh) | 一种转换波的偏移成像方法及装置 | |
CN104483705A (zh) | 一种三维剩余静校正方法 | |
CN111257941B (zh) | 一种组合式海底地震仪方位角自动识别装置及方法 | |
CN110146924B (zh) | 基于水波初至偏振方位的海底地震仪位置及方位反演方法 | |
Trad et al. | Fast and robust deblending using apex shifted Radon transform | |
CN105093319A (zh) | 基于三维地震数据的地面微地震静校正方法 | |
CN104199103A (zh) | 转换波静校正方法及装置 | |
CN104182651B (zh) | 用于三分量检波器接收的微地震事件方位角自动质控方法 | |
CN104199098B (zh) | 利用反射转换横波的偏振特性预测地下裂缝的方法 | |
CN112415571B (zh) | 一种微地震定位方法、存储介质及设备 | |
CN108375789A (zh) | 联合采集地震数据的同步匹配方法 | |
CN109085642A (zh) | 一种各向异性介质微地震事件定位方法 | |
CN102778691B (zh) | 一种计算检波器组内静校正时差的方法 | |
CN110579799A (zh) | 一种等旅行时间间隔的地震采集观测方法及系统 | |
CN114545495A (zh) | 一种ps波地震道数据的ovt道集处理方法、装置及电子设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170301 Termination date: 20210526 |