发明名称 一种三分量检波器埋置方向高精度校正的方法
摘要 本发明是一种利用地面多分量地震数据的偏振特性三分量检波器埋置方向高精度校正的方法。拾取初至时间,计算每一道三分量地震记录给定初至时窗的协方差矩阵,找出绝对值最大的非对角线元素,确定的旋转角度,构造一次旋转矩阵,对协方差矩阵的三个特征值按照从大到小的顺序进行排序,循环计算每一道的主特征向量在XY平面内的投影与X轴的夹角以及主特征向量与Z轴的夹角、方位角,计算检波器偏离角度,根据偏离角度对采集的三分量地震数据进行校正。本发明有很高的精度,有效地减小了随机误差,提高了计算结果的可靠性和校正的计算效率。
申请公布号 CN103675916A 申请公布日期 2014.03.26
申请号 CN201210323186.5 申请日期 2012.09.04
申请人 中国石油天然气集团公司;中国石油集团东方地球物理勘探有限责任公司 发明人 张文波;耿伟峰;李建峰;侯爱源
分类号 G01V1/36(2006.01)I 主分类号 G01V1/36(2006.01)I
代理机构 北京市中实友知识产权代理有限责任公司 11013 代理人 刘天语;武尚修
主权项 1.一种三分量检波器埋置方向高精度校正的方法,特点是采用以下步骤:1)利用正交三分量检波器采集三分量地震记录;2)拾取地震记录初至时间;3)计算每一道三分量地震记录给定初至时窗的协方差矩阵<img file="FDA00002094457200011.GIF" wi="103" he="63" />4)在协方差矩阵<img file="FDA00002094457200012.GIF" wi="77" he="58" />中找出绝对值最大的非对角线元素<img file="FDA00002094457200013.GIF" wi="90" he="63" />p,q是协方差矩阵的行列序号,1≤p,q≤3,然后利用相应的矩阵元素计算出sinθ和cosθ;其中,θ是由协方差矩阵元素<img file="FDA00002094457200014.GIF" wi="89" he="63" /><img file="FDA00002094457200015.GIF" wi="63" he="63" />和<img file="FDA00002094457200016.GIF" wi="63" he="63" />确定的旋转角度,sinθ和cosθ相应的计算过程变量;5)利用步骤4)得到的sinθ和cosθ以及步骤3)得到的协方差矩阵<img file="FDA00002094457200017.GIF" wi="103" he="60" />计算一次旋转后的协方差矩阵<img file="FDA00002094457200018.GIF" wi="72" he="59" />的元素<img file="FDA00002094457200019.GIF" wi="101" he="77" />i,j=1,2,3;6)构造一次旋转矩阵R<sup>(1)</sup>(p,q,θ);7)计算一次旋转后的协方差矩阵<img file="FDA000020944572000110.GIF" wi="73" he="58" />的非对角线元素的平方和E<sup>(1)</sup>(A);8)旋转变换迭代计算;9)构造协方差矩阵<img file="FDA000020944572000111.GIF" wi="77" he="58" />的特征值和特征向量;10)对步骤9)得到的协方差矩阵的三个特征值按照从大到小的顺序进行排序,使λ<sub>1</sub>≥λ<sub>2</sub>≥λ<sub>3</sub>,然后,得到排序后的特征值对应的特征向量V<sub>1</sub>,V<sub>2</sub>,V<sub>3</sub>,其中,V<sub>1</sub>=(v<sub>11</sub>v<sub>21</sub>v<sub>31</sub>)<sup>T</sup>,V<sub>2</sub>=(v<sub>12</sub>v<sub>22</sub>v<sub>32</sub>)<sup>T</sup>,V<sub>3</sub>=(v<sub>13</sub>v<sub>23</sub>v<sub>33</sub>)<sup>T</sup>,v<sub>ij</sub>是特征向量的元素,最大特征值λ<sub>1</sub>对应的特征向量V<sub>1</sub>即为质点振动的主方向;11)特征向量规则化;12)计算特征向量W<sub>1</sub>在XY平面内的投影与X轴的夹角<img file="FDA000020944572000112.GIF" wi="30" he="36" />以及特征向量W<sub>1</sub>与Z轴的夹角α;13)循环计算每一道的主特征向量在XY平面内的投影与X轴的夹角以及主特征向量与Z轴的夹角;14)循环计算每一道地震数据的方位角;15)计算每一个接收点处的三分量检波器的X分量与X坐标轴的偏离角度;16)计算每一个接收点处的三分量检波器的Z分量与垂直方向的偏离角度;17)根据由步骤15)和步骤16)计算得到的三分量检波器的X分量与X坐标轴的偏离角度Ψ和Z分量与垂直方向的偏离角度Φ,对采集的三分量地震数据进行校正。
地址 100007 北京市东城区东直门北大街9号中国石油大厦