正演数值模拟三维钻井噪音的方法及装置、设备、介质制造方法及图纸

本申请的实施例揭示了一种正演数值模拟三维钻井噪音的方法及装置、设备、介质。方法包括：根据钻机震源的特点，确定模拟钻机震源，估算模拟炸药震源延迟时间，并确定模拟炸药震源，在波动方程中加入模拟钻机震源和模拟炸药震源获得一个新方程，利用有限差分法对新方程进行数值模拟得到三维钻井噪音。本申请可以快速并准确地模拟出三维钻井噪音，为后续机器学习去除钻井噪音提供可靠训练集。习去除钻井噪音提供可靠训练集。习去除钻井噪音提供可靠训练集。

【技术实现步骤摘要】
正演数值模拟三维钻井噪音的方法及装置、设备、介质


[0001]本申请涉及声学处理
，具体涉及一种正演数值模拟三维钻井噪音的方法及装置、设备、介质。

技术介绍

[0002]随着对地震勘探精度要求的提高，我们对地震资料的采集要求也逐渐提升。但是，在地震勘探采集过程中，勘探区往往会存在或多或少的钻井工作，由此钻机振动产生的噪音与地震勘探震源产生的信号相互干扰，严重影响地震勘探资料。
[0003]由于钻井作业的连续性特点，若是时常中途停止钻井，易造成卡钻、塌井等情况，协调难度大，因此去除地震勘探资料的钻井噪音非常重要。如图1所示，图1中a表示钻机噪音，b表示钻机噪音振幅谱，已有技术提出钻井噪音在地震记录上呈现能量强，双曲线特征，频谱上频带很窄，近似一单频信号；且采用了叠加去噪音和模型去噪法对钻井噪音进行去除(详参[AN Ling
‑
Fang，安玲芳，ZHANG Jin，等.钻井干扰的特征分析及去噪方法浅析[C]//国家安全地球物理丛书.中国地球物理学会国家安全地球物理专业委员会；陕西省地球物理学会，2014.])。
[0004]随着计算机处理数据能力的提高和机器学习在数值信号处理方面能力的应用，采取机器学习方法，去除钻井噪音已是一种可行的方法。但是机器学习，需要大量的样本作为训练集，依靠实测来获取大量的样本的方式成本高，时间周期长，因此，高精度有效的三维钻井噪音数据的模拟就十分重要。

技术实现思路

[0005]为解决上述技术问题，本申请的实施例分别提供了正演数值模拟三维钻井噪音的方法及装置、设备、计算机可读存储介质，通过计算机模拟可以快速精准地获取到大量三维钻井噪音数据以满足机器学习的需要。
[0006]本申请的其他特性和优点将通过下面的详细描述变得显然，或部分地通过本申请的实践而习得。
[0007]根据本申请实施例的一个方面，提供了一种正演数值模拟三维钻井噪音的方法，所述方法包括：根据钻机震源的特点，确定如下公式(1)所示的模拟钻机震源：
[0008]A
i
(AX
i
，AY
i
，AZ
i
，t)＝a*sin(2πf
c
t)
ꢀꢀꢀꢀꢀꢀꢀ
(1)
[0009]式中，A
i
(AX
i
，AY
i
，AZ
i
，t)为第i个模拟钻机处的模拟钻机震源；(AX
i
，AY
i
，AZ
i
) 表示第i个模拟钻机位置，f
c
为模拟钻机机械振动频率；a为模拟钻机震源振幅；
[0010]通过如下公式(2)估算模拟炸药震源延迟时间t
D
：
[0011][0012]式中，(RX
i
，RY
i
，RZ
i
)为距离(AX
i
，AY
i
，AZ
i
)最远的接受点的坐标；v
i
为 (AX
i
，AY
i
，
AZ
i
)到(RX
i
，RY
i
，RZ
i
)之间的平均速度；
[0013]基于所述模拟炸药震源延迟时间t
D
，确定模拟炸药震源如下公式(3)所示：
[0014][0015]式中，S(SX，SY，SZ，t)为模拟炸药震源，(SX，SY，SZ)为模拟炸药震源位置， f
m
为雷克子波主频；
[0016]将波动方程中加入模拟炸药震源S(SX，SY，SZ，t)和模拟钻机震源 A
i
(AX
i
，AY
i
，AZ
i
，t)后得到如下公式(4)所示的方程：
[0017][0018]式中，P为波场值；V为波在介质中传播的速度；N为钻机个数；
[0019]通过有限差分法对公式(4)中的方程进行数值模拟。
[0020]根据本申请实施例的一个方面，提供了一种正演数值模拟三维钻井噪音的装置，包括：模拟钻机震源确定模块，被配置为根据钻机震源的特点，确定如下公式(1)所示的模拟钻机震源：
[0021]A
i
(AX
i
，AY
i
，AZ
i
，t)＝a*sin(2πf
c
t)
ꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0022]式中，A
i
(AX
i
，AY
i
，AZ
i
，t)为第i个模拟钻机处的模拟钻机震源；(AX
i
，AY
i
，AZ
i
) 表示第i个模拟钻机位置，f
c
为模拟钻机机械振动频率；a为模拟钻机震源振幅；
[0023]模拟炸药震源延迟时间估算模块，被配置为通过如下公式(2)估算模拟炸药震源延迟时间t
D
：
[0024][0025]式中，(RX
i
，RY
i
，RZ
i
)为距离(AX
i
，AY
i
，AZ
i
)最远的接受点的坐标；v
i
为 (AX
i
，AY
i
，AZ
i
)到(RX
i
，RY
i
，RZ
i
)之间的平均速度；
[0026]模拟炸药震源确定模块，被配置为基于所述模拟炸药震源延迟时间t
D
，确定模拟炸药震源如下公式(3)所示：
[0027][0028]式中，S(SX，SY，SZ，t)为模拟炸药震源，(SX，SY，SZ)为模拟炸药震源位置， f
m
为雷克子波主频；
[0029]第一计算模块，被配置为将波动方程中加入模拟炸药震源S(SX，SY，SZ，t)和模拟钻机震源A
i
(AX
i
，AY
i
，AZ
i
，t)后得到如下公式(4)所示的方程：
[0030][0031]式中，P为波场值；V为波在介质中传播的速度；N为钻机个数；
[0032]数值模拟模块，被配置为通过有限差分法对公式(4)中的方程进行数值模拟。
[0033]根据本申请实施例的一个方面，提供了一种电子设备，包括：控制器；存储器，用于存储一个或多个程序，当所述一个或多个程序被所述控制器执行时，使得所述控制器实现上所述的正演数值模拟三维钻井噪音的方法。
[0034]根据本申请实施例的一个方面，还提供了一种计算机可读存储介质，其上存储有计算机可读本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种正演数值模拟三维钻井噪音的方法，其特征在于，所述方法包括：根据钻机震源的特点，确定如下公式(1)所示的模拟钻机震源：A
i
(AX
i
，AY
i
，AZ
i
，t)＝a*sin(2πf
c
t)
ꢀꢀꢀꢀ
(1)式中，A
i
(AX
i
，AY
i
，AZ
i
，t)为第i个模拟钻机处的模拟钻机震源；(AX
i
，AY
i
，AZ
i
)表示第i个模拟钻机位置，f
c
为模拟钻机机械振动频率；a为模拟钻机震源振幅；通过如下公式(2)估算模拟炸药震源延迟时间t
D
：式中，(RX
i
，RY
i
，RZ
i
)为距离(AX
i
，AY
i
，AZ
i
)最远的接受点的坐标；v
i
为(AX
i
，AY
i
，AZ
i
)到(RX
i
，RY
i
，RZ
i
)之间的平均速度；基于所述模拟炸药震源延迟时间t
D
，确定模拟炸药震源如下公式(3)所示：式中，S(SX，SY，SZ，t)为模拟炸药震源，(SX，SY，SZ)为模拟炸药震源位置，f
m
为雷克子波主频；将波动方程中加入模拟炸药震源S(SX，SY，SZ，t)和模拟钻机震源A
i
(AX
i
，AY
i
，AZ
i
，t)后得到如下公式(4)所示的方程：式中，P为波场值；V为波在介质中传播的速度；N为钻机个数；通过有限差分法对公式(4)中的方程进行数值模拟。2.一种正演数值模拟三维钻井噪音的装置，其特征在于，包括：模拟钻机震源确定模块，被配置为根据钻机震源的特点，确定如下公式(1)所示的模拟钻机震源：A
i
(AX
i
，AY
i
，...

【专利技术属性】
技术研发人员：文晓涛，唐超，刘炀，李波，任红萍，
申请(专利权)人：成都理工大学，
类型：发明
国别省市：