一种海底浅部沉积层剪切波速度探测装置及方法制造方法及图纸

技术编号:27200308 阅读:32 留言:0更新日期:2021-01-31 12:07
本发明专利技术公开了一种海底浅部沉积层剪切波速度探测装置及方法,属于海洋地震勘探技术领域。方法实现包括如下步骤。首先进行作业区人工地震探测系统设计,主要是给定合理的海底地震仪器设备投放间距和海面气枪激发间距;按照设计进行海上人工源地震作业。其次开展海底面波的数据预处理,包括格式转换、钟差校正、数据截取等。然后进行海底面波频散能量成像、频散曲线提取。最后开展频散曲线反演。装置基于上述方法实现。通过本发明专利技术可实现海底浅部沉积层剪切波速度成像,为我国海洋地质工程调查、海洋地球物理调查提供新的技术支撑。洋地球物理调查提供新的技术支撑。洋地球物理调查提供新的技术支撑。

【技术实现步骤摘要】
一种海底浅部沉积层剪切波速度探测装置及方法


[0001]本专利技术涉及属于海洋地震勘探
,具体涉及一种海底浅部沉积层剪切波速度探测装置及方法。主要应用于海洋地质工程调查,同时应用于海洋地球物理调查。

技术介绍

[0002]海洋浅部沉积层剪切波速度结构探测在海洋地质工程、海洋地球物理调查中至关重要。在许多海洋工程建设活动中,如海上石油平台、海洋风场、海底管道和海底隧道等,需要评价海底沉积层是否存在“软”的介质、断裂等,除了高成本的海底钻井这种直接探测方式外,浅海沉积物的剪切波速度结构是一种有效的间接评价手段。在海洋地球物理勘探中,为了开展多波多分量海底地震勘探,需要先验的浅部剪切波速度模型。在海洋深部剪切波速度成像中,需要先验的浅部沉积层剪切波速度模型进行剪切波到时差校准。
[0003]利用海底的体波和面波均可开展浅部沉积层剪切波速度研究。相对于海底体波探测技术,海底面波具有横向衰减慢、传播距离更远的优势。因浅部沉积层往往剪切波速度特别低,体波的压缩波-剪切波比值较低,很难在海底采集到可用于研究的高信噪比转换波。但面波具有频散特性,该特性使得海底面波相速度对浅部沉积层的剪切波速度最为敏感,有利于海底面波频散曲线反演剪切波速度结构。
[0004]海底面波是指在海水和沉积物界面之间传播的界面波,因外国学者Scholte于1947年率先完成数学证明,垂向或径向的海底面波命名为舒尔特波(Scholte-wave;Scholte,1947)。切向的海底面波与陆地的切向面波名称一致,均称为勒夫波(Love-wave;Love,1911)。
[0005]在海洋地质工程和地球物理调查中,浅部沉积层剪切波速度探测是基础内容。但因海洋特殊环境,对探测仪器、作业方式和数据处理提出了新的要求,传统海洋钻井方式经济成本高。因此,有必要研发一种海底浅部沉积层剪切波速度探测装置及方法,可行性高、算法稳定、施工成本低、反演结果可靠,将为海洋地质工程、海洋地球物理调查提供可靠技术支撑。

技术实现思路

[0006]本专利技术的目的在于提供一种海底浅部沉积层剪切波速度探测方法,以解决海洋地质工程调查中浅部沉积层剪切波速度的成像问题,通过海底面波数据采集、主动源道集数据预处理、频散能量成像与多模频散曲线反演,实现海底浅部沉积层剪切波速度成像。这将大大降低海洋地质工程中大面积浅部沉积层地质普查的成本,为快速高效实现区域海洋地质工程调查提供可能,更为海洋多波多分量地震勘探、海底深部剪切波成像提供了重要先验参考模型。
[0007]具体的,海底浅部沉积层剪切波速度探测方法,所述方法包括以下步骤:
[0008]步骤101:针对探测区域设计观测系统,包括海底地震仪投放点位与气枪激发间距的确定;按照设计的观测系统进行海上地震数据采集,包括海底地震仪投放、气枪放炮、海
底地震仪回收、海底地震仪数据回收;
[0009]步骤102:利用人工震源激发的时间与位置导航信息,从海底地震仪连续记录数据序列中截取海底地震仪所在台站对应的共接收点道集,对共接收点道集先后完成时间漂移校正、道均衡处理与带通滤波处理,将结果输出保存为SU或SEGY格式;
[0010]步骤103:读取海底地震仪垂向分量的SU或SEGY格式数据,利用相移法进行频散能量成像;从频散能量结果图中,按照能量振幅极大值,提取多阶舒尔特波频散曲线;
[0011]步骤104:根据半波长理论选取最大反演深度,建立一维初始地球物理模型,建立多模频散曲线反演目标函数,基于阻尼最小二乘优化算法求取该目标函数,迭代至预设目标函数收敛,获得最优反演模型;
[0012]步骤105:保存输出最优反演模型,作为最终的海底浅部沉积层剪切波速度结果。
[0013]另一方面,本专利技术还提供了一种用于实现上述方法的海底浅部沉积层剪切波速度探测装置,所述装置包括:海底人工源地震数据预处理模块,用于转换原始数据格式、提取共接收点道集、进行姿态校准、时钟漂移校正、道均衡及带通滤波;舒尔特波频散能量成像模块,用于计算共接收点道集的频散能量谱;舒尔特波频散曲线提取模块,用于从频散能量谱中通过振幅极值提取频散曲线;舒尔特波频散曲线反演模块,用于反演拟合最优浅部沉积层剪切波速度模型;剪切波模型保存输出模块,用于保存和输出最优剪切波反演迭代模型。
[0014]本专利技术的有益效果在于:本专利技术设计一种海底浅部沉积层剪切波速度探测装置及方法,对于海底浅部沉积层剪切波速度探测难题,利用海底地震观测仪器与人工震源激发的海上作业方式,经后期面波频散分析成像获得浅部沉积层剪切波速度结构。通过本专利技术可实现海底浅部沉积层剪切波速度成像,为我国海洋地质工程调查、海洋地球物理调查提供新的技术支撑。
附图说明
[0015]图1是本专利技术方法的流程图。
[0016]图2是本专利技术装置的结构框图。
[0017]图3是本专利技术观测系统设计示意图。
具体实施方式
[0018]下面结合附图对本专利技术的具体实施方式进行说明:
[0019]如图1所示,本专利技术提出的一种海底浅部沉积层剪切波速度探测方法,包括如下步骤:
[0020]步骤101:针对探测区域设计观测系统,包括海底地震仪投放点位与气枪激发间距。按照设计观测系统进行海上地震数据采集,包括海底地震仪投放、气枪放炮、海底地震仪回收、海底地震仪数据回收。
[0021]步骤102:利用人工震源激发的时间与位置导航信息,从海底地震仪连续记录数据序列中截取该台站对应的共接收点道集,并先后完成时间漂移校正、道均衡处理与带通滤波处理,并输出保存为SU或SEGY格式。
[0022]步骤103:读取海底地震仪垂向分量的SU或SEGY格式数据,利用相移法进行频散能
量成像。从频散能量结果图中,按照能量振幅极大值,提取多阶舒尔特波频散曲线。
[0023]步骤104:根据半波长理论选取最大反演深度,建立一维初始地球物理模型,建立多模频散曲线反演目标函数,基于阻尼最小二乘优化算法求取该目标函数,迭代至预设目标函数收敛。
[0024]步骤105:保存输出最优反演模型,作为最终的海底浅部沉积层剪切波速度结果。
[0025]上述步骤101中,以图3为例进行说明,图3中,S点位海底地震观测设备仪器投放点位,O
m
为气枪激发震源最大偏移距位置,Δ为最大偏移距,A点位观测点位(即剪切波速度反演结果对应位置)。图3的附图标记包括海地地震仪1、海面2、海底3、水深4。
[0026]为获得图中点A点处海底浅部沉积层剪切波速度结构,需要先设定气枪激发震源的最大偏移距Δ,一般取值最小为2500m,更大偏移距会提高频散能量成像分辨率;在距离最大偏移距Δ一半的位置设计海底地震仪器观测点位(图3中S点);最后设计气枪的投放深度与激发点距,气枪投放深度范围为5-10m,视作业船只、海况和气枪设备而定,气枪激发点距Δ
x
需满足下式:
[0027][0028]其中,V
sch本文档来自技高网
...

【技术保护点】

【技术特征摘要】
1.一种海底浅部沉积层剪切波速度探测方法,其特征在于,所述方法包括以下步骤:步骤101:针对探测区域设计观测系统,包括海底地震仪投放点位与气枪激发间距的确定;按照设计的观测系统进行海上地震数据采集,包括海底地震仪投放、气枪放炮、海底地震仪回收、海底地震仪数据回收;步骤102:利用人工震源激发的时间与位置导航信息,从海底地震仪连续记录数据序列中截取海底地震仪所在台站对应的共接收点道集,对共接收点道集先后完成时间漂移校正、道均衡处理与带通滤波处理,将结果输出保存为SU或SEGY格式;步骤103:读取海底地震仪垂向分量的SU或SEGY格式数据,利用相移法进行频散能量成像;从频散能量结果图中,按照能量振幅极大值,提取多阶舒尔特波频散曲线;步骤104:根据半波长理论选取最大反演深度,建立一维初始地球物理模型,建立多模频散曲线反演目标函数,基于阻尼最小二乘优化算法求取该目标函数,迭代至预设目标函数收敛,获得最优反演模型;步骤105:保存输出最优反演模型,作为最终的海底浅部沉积层剪切波速度结果。2.如权利要求1所述的海底浅部沉积层剪切波速度探测方法,其特征在于,所述方法步骤101中,先设定气枪激发震源的最大偏移距Δ,取值最小为2500m,在距离最大偏移距Δ一半的位置设计海底地震仪器观测点位;最后设计气枪的投放深度与激发点距,气枪投放深度范围为5-10m,气枪激发点距Δ
x
满足以下公式:其中,V
sch
为探测区域的最小舒尔特波速度,f
max
为本次地震观测记录到的舒尔特波的最大频率。3.如权利要求1所述的海底浅部沉积层剪切波速度探测方法,其特征在于,所述方法步骤103中,用相移法计算舒尔特波的频散能量谱,方法如下;对时空域垂向分量海底舒尔特波共接收点道集r
i
进行傅里叶变换:R
i
(x
k
,ω)=FFT(r
i
)其中,i为道数序号,x
k
为偏移距,R
i
(x
k
,ω)为频谱。然后对频谱进行振幅归一化处理获得归一化频谱T
i
(x
k
,ω):其中,|R
i
(x
k
,ω)|表示取模运算;给定任一点的舒尔特波相速度v
s
和频率ω,可获得该点的频散能量值为:其中,j为复数单位,选取一定的舒尔特波相速度...

【专利技术属性】
技术研发人员:王元游庆瑜郝天珧胡耀星赵春蕾张妍徐锡强
申请(专利权)人:中国科学院地质与地球物理研究所
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1