[发明专利]海上风浪一致性影响量板制作方法、系统、介质及设备

权利要求书

1.一种海上风浪一致性影响量板制作方法，其特征在于，包括：

步骤(1)、建立勘探地区的深度域地震速度模型；

步骤(2)、采用波动方程有限元数值解法，模拟平静海面条件下的拖缆采集单炮记录；

步骤(3)、重复步骤(2)，根据海况分级简表中的数据，模拟全部浪高级别的起伏海面条件下的拖缆采集单炮记录；

步骤(4)、以平静海面条件下的炮记录为地震数据的基数据，以起伏海面条件下的炮记录地震数据的监测数据，计算目标层位的NRMS值，根据NRMS值评估两次采集地震数据的一致性；

步骤(5)、重复步骤(4)两两计算全部起伏海面下的炮记录上目标层位的NRMS值；

步骤(6)、根据步骤(4)和步骤(5)计算的全部浪高条件下的NRMS值绘制勘探地区目标层位的风浪一致性影响量板；

步骤(2)中，波动方程有限元数值解法采用四边形单元，双线性插值函数；采用四边形网格剖分深度域地震速度模型，将炮点和接收点放置在海面以下预先设定沉放深度的计算节点上，以确保数值计算的精度；

双线性插值函数的构造方法为：

正方向单元四个顶点的编号及其坐标，构造的插值函数为：

其中，N_i表示第i个点的插值函数；ξ和η分别是自然坐标中的横坐标和纵坐标，ξ_i和η_i表示点i的自然坐标中的横坐标和纵坐标，i＝1,2,3,4；

直角坐标系下任意四边形单元的四个顶点坐标(x₁,y₁),……,(x₄,y₄)和函数值u₁,……,u₄表示成：

其中，x,y,u均表示为ξ和η的线性函数，上述插值函数称作双线性插值函数；

采用四边形网格剖分深度域地震速度模型，将炮点和接收点放置在海面以下预先设定沉放深度的计算节点上，以确保数值计算的精度；震源函数采用延时的Ricker子波，计算方法如下式：

其中，R(t)表示t时刻的震源函数值；f_p表示峰值频率；τ表示延迟时间；除海面采用Dirichlet零边界条件以外，其余边界全部采用吸收边界条件。

2.如权利要求1所述的量板制作方法，其特征在于：所述深度域地震速度模型的建立方法为：根据地震解释成果建立勘探地区的地震速度模型，即在层位解释数据的基础上结合岩性解释数据、测井数据及时深关系建立探区深度域地震速度模型；当没有勘探地区的地震解释成果作为参考时，采用深度域或者时间域的地震成像结果作为探区构造模型，并用偏移速度填充为速度模型。

3.如权利要求1所述的量板制作方法，其特征在于：所述波动方程有限元数值解法采用四边形单元，双线性插值函数。

4.如权利要求3所述的四边形单元，其特征在于：采用四边形网格剖分所述深度域地震速度模型，将炮点和接收点放置在海面以下预先设定沉放深度的计算节点上，以确保数值计算的精度。

5.如权利要求1所述的量板制作方法，其特征在于：所述步骤(3)中，起伏海面条件下的炮记录按正弦函数变化，计算公式为：

其中，z和x表示空间变量；h表示海浪波高；λ表示海浪波长；x_s表示炮点横坐标。

6.如权利要求1所述的量板制作方法，其特征在于：所述步骤(4)中，以NRMS值小于0.1作为评判标准，NRMS值小于0.1表明基数据和监测数据的一致性满足要求。

7.如权利要求1所述的量板制作方法，其特征在于：所述步骤(5)中，两两计算全部起伏海面下的炮记录上目标层位的NRMS值的方法为：把浪高a米的炮记录和浪高b米的炮记录分别当作基数据和监测数据，计算NRMS值，并包含a等于b的情况。