[发明专利]一种水陆检波器地震数据合并处理的方法及装置

权利要求书

1.一种水陆检波器地震数据合并处理的方法，其特征在于，提供有目标区域的水中检波器地震数据和陆地检波器地震数据，所述方法包括：

对所述陆地检波器地震数据进行相位处理，得到相位处理后的陆地检波器地震数据；

对所述水中检波器地震数据进行第一振幅处理，得到振幅处理后的水中检波器地震数据，以及对所述相位处理后的陆地检波器地震数据进行第二振幅处理，得到振幅处理后的陆地检波器地震数据；

对所述振幅处理后的水中检波器地震数据和陆地检波器地震数据分别进行傅里叶变换，得到傅里叶变换后的水中检波器地震数据和傅里叶变换后的陆地检波器地震数据；

根据所述振幅处理后的水中检波器地震数据和所述振幅处理后的陆地检波器地震数据计算标定因子，并根据所述傅里叶变换后的水中检波器地震数据和所述傅里叶变换后的陆地检波器地震数据，计算所述目标区域的海水深度和海底反射系数；

基于所述海水深度，确定镜像组合因子；

基于所述标定因子、所述海水深度和所述海底反射系数，以及所述镜像组合因子，对所述目标区域的水中检波器地震数据和所述目标区域的陆地检波器地震数据进行镜像组合处理，得到所述目标区域的水陆检波器合并数据；其中，采用下述公式进行所述镜像组合处理：

S＝O_H[n]*H+O_G[n]*G

其中，S表示所述水陆检波器合并数据，H表示所述目标区域的水中检波器地震数据，G表示所述目标区域的陆地检波器地震数据；*表示褶积运算符号；O_H[n]表示水检数据组合算子，O_G[n]表示陆检数据组合算子，n表示镜像组合算子顺序号，n＝-NOp,-NOp+1,-NOp+2,...,-1,0,1,...,NOp-2,NOp-1,NOp，(2NOp+1)表示镜像组合算子长度；Δt表示地震数据时间采样间隔；α表示最佳标定因子；A表示合并后数据期望绝对值振幅；A_H表示水中检波器道集数据平均绝对值振幅，A_G表示陆地检波器道集数据平均绝对值振幅；R_s表示水面反射系数；R表示最佳海底反射系数；Z表示延迟算子，(1-Z)和(1+Z)表示所述镜像组合因子；l为虚数单位，且l²＝-1，ω表示角频率，D表示海水深度，V表示海水速度，Φ表示对陆地检波器地震数据进行相位处理时得到的相位校正因子。

2.根据权利要求1所述的一种水陆检波器地震数据合并处理的方法，其特征在于，采用下述公式对所述水中检波器地震数据进行第一振幅处理：

其中，表示所述振幅处理后的水中检波器地震数据，H_i,j表示水中检波器地震数据，A_H表示水中检波器道集数据平均绝对值振幅，i表示时窗参考道的顺序号，i＝1,2,...,L，L表示时窗参数总道数，j表示时窗参数时间样点的顺序号，j＝1,2,...,N，N表示时窗参数时间样点数，符号||·||为取绝对值运算符。

3.根据权利要求1所述的一种水陆检波器地震数据合并处理的方法，其特征在于，采用下述公式对所述相位处理后的陆地检波器地震数据进行第二振幅处理：

其中，表示所述振幅处理后的陆地检波器地震数据，G_i,j表示所述相位处理后的陆地检波器地震数据，A_G表示陆地检波器道集数据平均绝对值振幅，i表示时窗参考道的顺序号，i＝1,2,...,L，L表示时窗参数总道数，j表示时窗参数时间样点的顺序号，j＝1,2,...,N，N表示时窗参数时间样点数，符号||·||为取绝对值运算符。