一种水陆检波器数据海底反射系数反演的方法技术

一种水陆检波器数据海底反射系数反演的方法，特点是包括以下步骤：用地震人工震源激发和采集水中检波器与陆地检波器地震数据并做预处理；计算水中检波器数据与陆地检波器数据傅立叶变换；计算上行波场以及上行波场一阶和二阶延迟数据；计算上行波场以及上行波场一阶和二阶延迟时间序列数据；计算海底反射系数特征方程系数；构造并求解海底反射系数特征方程；计算上行波数据能量；确定最佳海底反射系数；绘制海底反射系数剖面和存储海底反射系数于地震数据道头中；本发明专利技术实现水中检波器数据与陆地检波器数据的合并处理，消除地震数据中海水鸣震多次波干扰影响，有效提高地震数据信噪比和分辨率，具有一定的抗噪能力。

【技术实现步骤摘要】
【专利摘要】，特点是包括以下步骤：用地震人工震源激发和采集水中检波器与陆地检波器地震数据并做预处理；计算水中检波器数据与陆地检波器数据傅立叶变换；计算上行波场以及上行波场一阶和二阶延迟数据；计算上行波场以及上行波场一阶和二阶延迟时间序列数据；计算海底反射系数特征方程系数；构造并求解海底反射系数特征方程；计算上行波数据能量；确定最佳海底反射系数；绘制海底反射系数剖面和存储海底反射系数于地震数据道头中；本专利技术实现水中检波器数据与陆地检波器数据的合并处理，消除地震数据中海水鸣震多次波干扰影响，有效提高地震数据信噪比和分辨率，具有一定的抗噪能力。【专利说明】
本专利技术涉及油田的勘探、开发、开采技术，具体是为反映地下地层层位、油藏描述提供高分辨率的地震图形和数据的。
技术介绍
随着地震勘探技术的发展，海上油气勘探的难度和深度也越来越大，对地震资料的信噪比和分辨率要求也越来越高。三维地震可以精确地描述储层。OBC(Ocean BottomCable，海底电缆)，是一种联合海上和陆地地震数据采集技术，把检波器固定在海底，可以获得高分辨率三维地震数据。它使用一个固定的排列，把检波器固定在海底接收地震波，而一艘航船仅仅拖拽震源进行地震波激发。在OBC数据采集中，至少三艘船:一艘震源船，仅仅拖拽着气枪震源排列，进行地震波激发；一艘接收船，是固定不动的，连接着海底电缆，接收地震波；一艘船或者几艘船，铺设海底电缆和回收海底电缆。OBC目前仅仅局限于水深不超过150米的海域。但是随着电缆和电缆回收系统的改进，OBC技术正在向更深的海域推进。在一些特殊的地区，为了方便油藏监控试验，设计试验系统，把接收电缆保留在海底多年，以方便多次数据采集。为了得到地下介质高精度的三维图像，必须精确地知道所有炮点和接收点的位置。GPS (Global Posit1ning System,全球定位系统)数据或者基于海岸的无线电定位，提供炮船精确位置。同时也要精确定位电缆船，把每一个检波器放置到预先设计的位置。 在OBC数据采集时，由于海底和海面都是较强的反射界面。随着一个反射地震子波从地下达到海底，海底电缆中的检波器，感应并记录下这个反射地震子波。这个反射子波继续向上前进达到海面，受到海面的反射，然后改变方向向下传播，达到海底。海底电缆中的检波器，再一次感应并记录下这个地震子波。同时这个地震子波受到海底的反射，然后改变方向向上传播，达到海面，受到海面的反射，然后改变方向向下传播，达到海底。这个循环重复进行。而这些原始反射地震子波不希望的二次和后续达到，就是海水鸣震多次波(交混回响)。海水鸣震多次波是海上地震勘探数据中最大的噪声干扰。消除海水鸣震多次波噪声干扰，是海上地震数据处理中最为重要的步骤。 在浅水区，使用反褶积方法，可以有效消除海水鸣震多次波干扰，以恢复一次反射地震子波。但是对于海水深度超过10米的地区，原始反射地震子波与后续海水鸣震多次波干扰之间的时差很大，使得反褶积算法不能有效去除海水鸣震多次波干扰产生的同相轴。如果不能有效地从地震数据中去除海水鸣震多次波干扰，那么几个子波代表着一个反射界面，这样就会模糊地质断层界面。OBC采集的数据，提供了同一位置水中检波器数据和陆地检波器数据两种数据。这两种数据分别使用水中检波器和陆地检波器记录。水中检波器是一种压力检波器，记录的是地震波产生的压力变化；陆地检波器是一种质点速度检波器，记录的是质点速度变化。由于这两种检波器的记录机理不同，对于同一位置处海水鸣震多次波干扰，表现出不同特征。与水中检波器记录的海水鸣震多次波干扰相比，陆地检波器记录的海水鸣震多次波干扰表现出极性和振幅特征差异。两种检波器记录的海水鸣震多次波干扰，其极性是相反的，振幅是不同的，且相差一个与海底反射系数成比例的常数，这个常数值就是标定因子。因此利用这种振幅和极性特征差异，可以有效消除海水鸣震多次波干扰。这样消除海水鸣震多次波干扰的步骤包括:(I)在每一个接收点位置处，记录水中检波器数据和陆地检波器数据两种数据；(2)利用两种检波器的传感器灵敏度(传导常数)，调整陆地检波器数据的振幅，以匹配水中检波器数据的振幅；(3)计算确定标定因子、海水深度、海底反射系数参数，(4)利用标定因子，标定调整振幅后的陆地检波器数据；(5)把标定后的陆地检波器数据，与对应的水中检波器数据相加，得到水陆地检波器标定数据，(6)利用海水深度、海底反射系数参数，消除海水鸣震多次波干扰。这样，计算海底反射系数，构成了室内海上地震数据处理消除海水鸣震多次波干扰方法的关键步骤。 在包含海水鸣震多次波的OBC数据中，由于数据的能量是由地下反射波的能量和海水鸣震多次波的能量构成。消除海水鸣震多次波后，地震数据中，已经消除了多次波能量，剩下地下反射波的能量。因此为了有效估算海底反射系数，通常使用能量最小作为确定海底反射系数的准则。 目前常规处理方法是，采用扫描方法计算确定海底反射系数。采用预先设定一个海底反射系数范围值和扫描步长，采用扫描方法给出一系列的反射系数值，对一系列海底反射系数，计算一系列Backus海水鸣震多次波逆滤波器算子，对频率域水陆检上行波场数据(即，水陆检波器数据之和)乘以Backus海水鸣震多次波逆滤波器算子，得到滤波处理后的上行波场数据，然后使用傅立叶逆变换，将上行波场数据变换到时间域，在时间域计算一系列波场数据能量，其中最小能量所对应的海底反射系数，就是所求最佳海底反射系数。以上方法采用预先设定一个海底反射系数范围值和扫描步长，采用扫描方法，因此计算费时，且计算的海底反射系数数值必定存在误差。而且计算精度低，计算费时，不能满足实际数据处理需要。 专利技术目的 本专利技术提供一种目的在于快速、准确地进行水中检波器数据和陆地检波器数据合并处理，消除地震数据中海水鸣震多次波干扰，计算精确、效率高的水陆检波器数据海底反射系数反演的方法。
技术实现思路
本专利技术包括以下步骤: I)用地震人工震源激发和采集水中检波器与陆地检波器地震数据并做预处理； 步骤I)所述的预处理是指对地震数据置标签、定义观测系统、水中检波器与陆地检波器地震数据分离、速度分析、叠加预处理工作； 2)计算水中检波器数据与陆地检波器数据傅立叶变换； 按照以下公式计算水中检波器数据傅立叶变换Η」(ω): ",(ω) = ￡>,.产(!) H=I 式中，hj为水检数据；i为虚数单位，且i2 = _1 ; ω为角频率,单位为弧度/秒(Ω/s) ;n表示共检波点道集数据时间样点顺序号，η = I, 2, 3，…，IN, IN表示共检波点道集数据时间样点总个数表示共检波点道集数据道顺序号= 1，2，3，…，IL，IL表示共检波点道集数据总道数。 按照以下公式计算陆地检波器数据傅立叶变换(^(ω): ^/(ω)-Χ-/[/?>/,"Η(2) ?-1 式中，gj 为陆检数据。 3)计算上行波场以及上行波场一阶和二阶延迟数据； 按照以下公式计算上行波场数据Cj ( ω ): Cj ( ω ) = Hj ( ω ) +Gj ( ω )(3) 按照以下公式计算上行波场一阶延迟数据Dj(Co): Dj ( ω ) = -本文档来自技高网...

【技术保护点】
一种水陆检波器数据海底反射系数反演的方法，特点是包括以下步骤：1)用地震人工震源激发和采集水中检波器与陆地检波器地震数据并做预处理；2)计算水中检波器数据与陆地检波器数据傅立叶变换；3)计算上行波场以及上行波场一阶和二阶延迟数据；按照以下公式计算上行波场数据Cj(ω)：Cj(ω)＝Hj(ω)+Gj(ω)               (3)按照以下公式计算上行波场一阶延迟数据Dj(ω)：Dj(ω)＝‑2Rs(Hj(ω)+Gj(ω))eiωτ         (4)式中，Rs为海面反射系数，τ为水层双程旅行时间，单位为秒(s)。按照以下公式计算上行波场二阶延迟数据Ej(ω)：Ej(ω)=Rs2(Hj(ω)+Gj(ω))ei2ωτ---(5);]]>4)计算上行波场以及上行波场一阶和二阶延迟时间序列数据；按照以下公式计算上行波场时间序列数据cj[n]：cj[n]=12π∫-∞∞Cj(ω)e-iωndω---(6)]]>按照以下公式计算上行波场一阶延迟时间序列数据dj[n]：dj[n]=12π∫-∞∞Dj(ω)e-iωndω---(7)]]>按照以下公式计算上行波场二阶延迟时间序列数据ej[n]：ej[n]=12π∫-∞∞Ej(ω)e-iωndω---(8)]]>5)计算海底反射系数特征方程系数；按照以下公式计算海底反射系数特征方程系数a0、a1、a2、a3和a4：a0=Σj=1ILΣn=1INcj2[n]]]>a1=2Σj=1ILΣn=1INcj[n]dj[n]]]>a2=Σj=1ILΣn=1IN(dj2[n]+2cj[n]ej[n])]]>a3=2Σj=1ILΣn=1INdj[n]ej[n]]]>a4=Σj=1ILΣn=1INej2[n]---(9)]]>6)构造并求解海底反射系数特征方程；按照以下公式构造海底反射系数特征方程：a1+2a2R+3a3R2+4a4R3＝0            (10)式中，R为海底反射系数特征方程的变元，求解特征方程(10)，可以得到三个海底反射系数数值R1、R2和R3；7)计算上行波数据能量；按照以下公式计算上行波数据能量Q1、Q2和Q3：Q1=a0+a1R1+a2R12+a3R13+a4R14]]>Q2=a0+a1R2+a2R22+a3R23+a4R24]]>Q3=a0+a1R3+a2R32+a3R33+a4R34---(11)]]>8)确定最佳海底反射系数；按照以下公式确定最佳海底反射系数Rbest：Rbest=MinRbest∈[R1,R2,R3]{Q1,Q2,Q3}---(12)]]>9)绘制海底反射系数剖面和存储海底反射系数于地震数据道头中。...

【技术特征摘要】

【专利技术属性】
技术研发人员：高少武，赵波，罗国安，钱忠平，颉冬莲，张文栋，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，中油油气勘探软件国家工程研究中心有限公司，
类型：发明
国别省市：河北;13