[发明专利]基于波形驱动下的叠前纵横波速度比模拟方法

权利要求书

1.一种基于波形驱动下的叠前纵横波速度比模拟方法，其特征在于，所述模拟方法包括以下步骤：

计算得到叠前纵横波速度比模拟的输入数据体；

对所述输入数据体进行相控选样；

根据所述相控选样结果，模拟出伪井数据；

根据所述伪井数据和已知井数据建立初始模型，并统计已知测井样本的概率分布作为先验概率分布；

根据所述初始模型和所述输入数据体，计算得到似然函数；

根据所述初始模型、所述先验概率分布和所述似然函数，在贝叶斯框架下进行马尔科夫链-蒙特卡洛随机模拟，得到符合地震中频阻抗和井曲线结构特征模拟结果。

2.根据权利要求1所述的基于波形驱动下的叠前纵横波速度比模拟方法，其特征在于，通过叠前弹性参数纵横波速度比反演体与综合子波褶积得到所述输入数据体。

3.根据权利要求2所述的基于波形驱动下的叠前纵横波速度比模拟方法，其特征在于，通过叠前同时反演得到所述叠前弹性参数纵横波速度比反演体。

4.根据权利要求3所述的基于波形驱动下的叠前纵横波速度比模拟方法，其特征在于，所述叠前同时反演根据佐伊普里兹方程或其近似式，运用井约束稀疏脉冲算法求取。

5.根据权利要求2所述的基于波形驱动下的叠前纵横波速度比模拟方法，其特征在于，根据叠后地震数据和井数据提取所述综合子波。

6.根据权利要求1所述的基于波形驱动下的叠前纵横波速度比模拟方法，其特征在于，所述相控选样利用了地震波形相似性特征和最大空间距离。

7.根据权利要求6所述的基于波形驱动下的叠前纵横波速度比模拟方法，其特征在于，所述最大空间距离为通过空间函数求取的矢量距离。

8.根据权利要求1所述的基于波形驱动下的叠前纵横波速度比模拟方法，其特征在于，所述伪井数据的模拟包括：

根据式1计算N口已知井的所述输入数据体井旁道与一口未知井的所述输入数据体井旁道的道间关系函数，

式1为：S_new＝∑w_iS_i，i＝1,2,3...,N，

其中，w_i为所述输入数据体的道间关系函数，S_new为未知井的所述输入数据体井旁道，S_i为已知井的所述输入数据体井旁道；

根据式2求取不同波形参数的相关匹配关系函数，

式2为：P_b＝∑λ_abP_a，

其中,λ_ab为不同波形参数的相关匹配关系函数，P_b为所述输入数据体的波形参数，P_a为所述已知井数据波形参数；

根据式3求取不同波形参数之间的空间关系函数，

式3为：

其中，为不同波形参数之间的空间关系函数，w_i为所述输入数据体的道间关系函数，λ_ab为所述不同波形参数的相关匹配关系函数；

根据式4求取所述伪井数据，

式4为：

其中，L_new为所述伪井数据，为所述不同波形参数之间的空间关系函数，L_i为所述已知井数据。

9.根据权利要求1所述的基于波形驱动下的叠前纵横波速度比模拟方法，其特征在于，所述伪井数据为波阻抗数据。

10.根据权利要求1所述的基于波形驱动下的叠前纵横波速度比模拟方法，其特征在于，将所述初始模型和所述输入数据体进行匹配滤波计算得到似然函数。