地层弹性参数的确定方法及装置与流程

本发明涉及计算机领域，具体而言，涉及一种地层弹性参数的确定方法及装置。

背景技术：

叠前地震反演是储层预测和油气检测的有效手段，基于叠前反演的流体因子只包含孔隙岩石里面的流体信息，不包含岩石的岩性性质，对储层中的流体性质区分性强，可区分储层中的油、气、水层，而密度参数与储层的物性紧密相联有助于储层特征描述。因此，流体因子和密度对储层流体预测和特征描述具有重要意义。

叠前三参数avo反演是一种常见的获取地层弹性参数的方法，但受三参数反演稳定性和偏移距的限制，无法得到精确的反演结果。即使道集中存在少量的噪声都会导致反演结果精度变差。三参数avo反演很难得到精确的流体因子和密度参数。目前采用的基于含流体因子，剪切模量和密度的三参数叠前反射系数近似式反演流体因子和密度参数的方法。但三参数avo反演稳定性较差，而实际地震数据中又含有噪声，致使反演的流体因子和密度精度和信噪比都较差，无法进行储层流体预测和特征描述。

针对上述的问题，目前尚未提出有效的解决方案。

技术实现要素：

本发明实施例提供了一种地层弹性参数的确定方法及装置，以至少解决相关技术中地层弹性参数的确定效率较低的问题。

根据本发明的一个实施例，提供了一种地层弹性参数的确定方法，包括：根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗，其中，所述第一叠前反射系数近似式以纵波阻抗和横波阻抗为参数；将所述目标纵波阻抗转换为目标纵波速度；利用所述目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式，其中，所述目标叠前反射系数近似式以流体因子和密度为参数；根据所述目标叠前反射系数近似式获取目标地层弹性参数，其中，所述目标地层弹性参数包括目标流体因子和目标密度。

可选地，根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗包括：获取所述纵横波速度和密度的测井曲线；利用所述纵横波速度和密度的测井曲线获取纵横波速度比；利用所述纵横波速度比、纵波阻抗和横波阻抗构造所述第一叠前反射系数近似式；将所述第一叠前反射系数近似式转换为第一线性方程组，其中，所述第一线性方程组以纵波阻抗和横波阻抗为参数；对所述第一线性方程组进行求解，得到所述目标纵波阻抗。

可选地，利用所述目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式包括：构造以纵波速度、流体因子和密度为参数的所述第二叠前反射系数近似式；利用所述目标纵波速度将所述第二叠前反射系数近似式改写为所述目标叠前反射系数近似式。

可选地，根据所述目标叠前反射系数近似式获取目标地层弹性参数包括：将所述目标叠前反射系数近似式转换为第二线性方程组，其中，所述第二线性方程组以流体因子和密度为参数；对所述第二线性方程组进行求解，得到所述目标流体因子和所述目标密度。

可选地，对所述第二线性方程组进行求解，得到所述目标流体因子和所述目标密度包括：对所述第二线性方程组进行求解，得到求解结果，其中，所述求解结果包括：求解流体因子和求解目标密度；利用获取到的测井数据对所述求解流体因子和所述求解目标密度进行道积分和低频补偿处理，得到第一处理结果，其中，所述第一处理结果包括：所述目标流体因子和所述目标密度。

可选地，在根据所述目标叠前反射系数近似式获取目标地层弹性参数之后，所述方法还包括：利用所述目标地层弹性参数对工区进行储层流体预测和特征描述处理，得到第二处理结果；根据所述第二处理结果获取有利钻井区。

根据本发明的另一个实施例，提供了一种地层弹性参数的确定装置，包括：确定模块，用于根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗，其中，所述第一叠前反射系数近似式以纵波阻抗和横波阻抗为参数；第一转换模块，用于将所述目标纵波阻抗转换为目标纵波速度；第二转换模块，用于利用所述目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式，其中，所述目标叠前反射系数近似式以流体因子和密度为参数；第一获取模块，用于根据所述目标叠前反射系数近似式获取目标地层弹性参数，其中，所述目标地层弹性参数包括目标流体因子和目标密度。

可选地，所述装置还包括：处理模块，用于利用所述目标地层弹性参数对工区进行储层流体预测和特征描述处理，得到第二处理结果；第二获取模块，用于根据所述第二处理结果获取有利钻井区。

根据本发明的又一个实施例，还提供了一种存储介质，所述存储介质中存储有计算机程序，其中，所述计算机程序被设置为运行时执行上述任一项方法实施例中的步骤。

根据本发明的又一个实施例，还提供了一种电子装置，包括存储器和处理器，所述存储器中存储有计算机程序，所述处理器被设置为运行所述计算机程序以执行上述任一项方法实施例中的步骤。

通过本发明，通过根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗，其中，第一叠前反射系数近似式以纵波阻抗和横波阻抗为参数；将目标纵波阻抗转换为目标纵波速度；利用目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式，其中，目标叠前反射系数近似式以流体因子和密度为参数；根据目标叠前反射系数近似式获取目标地层弹性参数，其中，目标地层弹性参数包括目标流体因子和目标密度的方式，用两参数近似式(第一叠前反射系数近似式)反演纵波阻抗，把得到的纵波阻抗转换成纵波速度，再利用得到的纵波速度重新推导得到仅含流体因子和密度的两项avo近似式(目标叠前反射系数近似式)，利用该两项avo近似式进行两参数avo反演得到高信噪比和精度的目标流体因子和目标密度。因此，可以解决相关技术中地层弹性参数的确定效率较低的问题，达到提高地层弹性参数的确定效率的效果。

附图说明

此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1是本发明实施例的一种地层弹性参数的确定方法的移动终端的硬件结构框图；

图2是根据本发明实施例的地层弹性参数的确定方法的流程图；

图3是根据本发明实施例的地层弹性参数的确定装置的结构框图；

图4是根据本发明可选实施例的地层弹性参数的确定方法的示意图；

图5是根据本发明可选的实施方式的地震角道集的示意图；

图6是根据本发明可选的实施方式的反演的流体因子的示意图；

图7是根据本发明可选的实施方式的反演的密度的示意图。

具体实施方式

下文中将参考附图并结合实施例来详细说明本发明。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互组合。

需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。

本申请实施例一所提供的方法实施例可以在移动终端、计算机终端或者类似的运算装置中执行。以运行在移动终端上为例，图1是本发明实施例的一种地层弹性参数的确定方法的移动终端的硬件结构框图。如图1所示，移动终端10可以包括一个或多个(图1中仅示出一个)处理器102(处理器102可以包括但不限于微处理器mcu或可编程逻辑器件fpga等的处理装置)和用于存储数据的存储器104，可选地，上述移动终端还可以包括用于通信功能的传输设备106以及输入输出设备108。本领域普通技术人员可以理解，图1所示的结构仅为示意，其并不对上述移动终端的结构造成限定。例如，移动终端10还可包括比图1中所示更多或者更少的组件，或者具有与图1所示不同的配置。

存储器104可用于存储计算机程序，例如，应用软件的软件程序以及模块，如本发明实施例中的地层弹性参数的确定方法对应的计算机程序，处理器102通过运行存储在存储器104内的计算机程序，从而执行各种功能应用以及数据处理，即实现上述的方法。存储器104可包括高速随机存储器，还可包括非易失性存储器，如一个或者多个磁性存储装置、闪存、或者其他非易失性固态存储器。在一些实例中，存储器104可进一步包括相对于处理器102远程设置的存储器，这些远程存储器可以通过网络连接至移动终端10。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。

传输装置106用于经由一个网络接收或者发送数据。上述的网络具体实例可包括移动终端10的通信供应商提供的无线网络。在一个实例中，传输装置106包括一个网络适配器(networkinterfacecontroller，简称为nic)，其可通过基站与其他网络设备相连从而可与互联网进行通讯。在一个实例中，传输装置106可以为射频(radiofrequency，简称为rf)模块，其用于通过无线方式与互联网进行通讯。

在本实施例中提供了一种地层弹性参数的确定方法，图2是根据本发明实施例的地层弹性参数的确定方法的流程图，如图2所示，该流程包括如下步骤：

步骤s202，根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗，其中，第一叠前反射系数近似式以纵波阻抗和横波阻抗为参数；

步骤s204，将目标纵波阻抗转换为目标纵波速度；

步骤s206，利用目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式，其中，目标叠前反射系数近似式以流体因子和密度为参数；

步骤s208，根据目标叠前反射系数近似式获取目标地层弹性参数，其中，目标地层弹性参数包括目标流体因子和目标密度。

可选地，在本实施例中，第一叠前反射系数近似式可以但不限于为以纵波阻抗和横波阻抗为参数的两参数avo近似式。

可选地，在本实施例中，目标叠前反射系数近似式可以但不限于为以流体因子和密度为参数的两参数avo近似式。

可选地，在本实施例中，将目标纵波阻抗转换为目标纵波速度的方式可以但不限于包括利用gardner公式把得到的目标纵波阻抗转换成目标纵波速度。

在一个可选的实施方式中，为了克服叠前三参数avo反演流体因子和密度稳定性差的问题，提供了一种两参数avo反演流体因子和密度的方法，充分利用两参数avo反演的稳定性优势，通过两次两参数avo反演来得到流体因子和密度为参数，具体做法为第一次用fatti两参数近似式反演纵横波阻抗，利用gardner公式把得到的纵波阻抗转换成纵波速度，再利用得到的纵波速度重新推导得到仅含流体因子和密度的两项avo近似式，利用该两项avo近似式进行第二次两参数avo反演得到高信噪比和精度的流体因子和密度，此时可利用反演结果进行工区的储层流体预测和特征描述。

通过上述步骤，通过根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗，其中，第一叠前反射系数近似式以纵波阻抗和横波阻抗为参数；将目标纵波阻抗转换为目标纵波速度；利用目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式，其中，目标叠前反射系数近似式以流体因子和密度为参数；根据目标叠前反射系数近似式获取目标地层弹性参数，其中，目标地层弹性参数包括目标流体因子和目标密度的方式，用两参数近似式(第一叠前反射系数近似式)反演纵波阻抗，把得到的纵波阻抗转换成纵波速度，再利用得到的纵波速度重新推导得到仅含流体因子和密度的两项avo近似式(目标叠前反射系数近似式)，利用该两项avo近似式进行两参数avo反演得到高信噪比和精度的目标流体因子和目标密度。因此，可以解决相关技术中地层弹性参数的确定效率较低的问题，达到提高地层弹性参数的确定效率的效果。

可选地，在上述步骤s202中，获取纵横波速度和密度的测井曲线；利用纵横波速度和密度的测井曲线获取纵横波速度比；利用纵横波速度比、纵波阻抗和横波阻抗构造第一叠前反射系数近似式；将第一叠前反射系数近似式转换为第一线性方程组，其中，第一线性方程组以纵波阻抗和横波阻抗为参数；对第一线性方程组进行求解，得到目标纵波阻抗。

可选地，在本实施例中，可以但不限于通过以下方式获取纵横波速度和密度的测井曲线：获得角度域道集，利用反褶积得到角度域反射系数，同时获得纵横波速度和密度的测井曲线；可选方法如下：

基于野外采集的地震资料利用地震成像方法提取叠前地震角道集数据；

利用叠前地震地质标定方法获得角道集中的子波，利用最小平方脉冲反褶积方法去除角度域道集的子波，获得角度域反射系数。

利用声波测井曲线、偶极声波测井曲线和密度测井曲线获得纵横波速度和密度的测井曲线。

在一个可选的实施方式中，可以利用fatti两项avo近似式反演纵横波阻抗；可选方法如下：

在aki-richards近似式的基础上，当入射角小于35°，只包含纵波阻抗和横波阻抗的avo近似式为：

其中，ip和is分别表示平均纵波阻抗和横波阻抗，θ是入射角度，δip和δis表示两层介质边界处的纵横波阻抗变化值，γsat为饱和岩石的纵横波速度比。

上式(1)可构造出n个线性方程：

其中，

采用svd(奇异值分解)算法对公式(2)进行求解，从而得到纵波阻抗和横波阻抗。

可选地，在本实施例中，在上述步骤s204中，可以通过以下方式将目标纵波阻抗转换成目标纵波速度：

密度和速度之间的数值关系式为：

上式(3)经变换可写为：

公式(4)左右两边同时加纵波速度差，则可给出：

其中，a可将密度和纵波速度测井数据经最小二乘拟合获得，利用公式(5)可将纵波阻抗转换为纵波速度。

可选地，在上述步骤s206中，构造以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式；利用目标纵波速度将第二叠前反射系数近似式改写为目标叠前反射系数近似式。

可选地，在本实施例中，利用得到的纵波速度，重新改写推导的以纵波速度、流体因子和密度为参数的avo近似式为仅含流体因子和密度的两项avo近似式，用该近似式反演流体因子和密度；可选方法如下：

为了从叠前地震数据中提取流体因子，以流体因子、剪切模量和密度为参数的avo近似式：

其中，f、ρ、μ分别表示平均流体因子、密度和剪切模量，δf、δρ、δμ表示两层介质边界处流体因子、密度和剪切模量的变化量，γsat、γdry分别为饱和岩石和干岩石的纵横波速度比。

但由于三参数avo反演具有稳定性差的问题。因此，重新推导了新的以纵波速度、流体因子和密度为参数的avo近似式：

利用上述得到的纵波速度，公式(7)可写为：

纵波速度作为已知量，含有两个未知量的线性方程组可表示为：

其中，

可选地，在上述步骤s208中，将目标叠前反射系数近似式转换为第二线性方程组，其中，第二线性方程组以流体因子和密度为参数；对第二线性方程组进行求解，得到目标流体因子和目标密度。

可选地，在本实施例中，采用svd(奇异值分解)算法对公式(9)进行求解，从而反演得到流体因子和密度参数。

可选地，对第二线性方程组进行求解，得到目标流体因子和目标密度包括：对第二线性方程组进行求解，得到求解结果，其中，求解结果包括：求解流体因子和求解目标密度；利用获取到的测井数据对求解流体因子和求解目标密度进行道积分和低频补偿处理，得到第一处理结果，其中，第一处理结果包括：目标流体因子和目标密度。

可选地，在本实施例中，可以利用测井数据对反演结果进行道积分和低频补偿，得到最终的流体因子和密度。可选方法如下：

将反演得到的结果进行道积分，以测井资料计算的流体因子和密度的低频成分作为带限反演结果的低频成份，最终获取补偿后的流体因子和密度参数。

可选地，在上述步骤s208之后，利用目标地层弹性参数对工区进行储层流体预测和特征描述处理，得到第二处理结果；根据第二处理结果获取有利钻井区。

可选地，在本实施例中，利用流体因子和密度进行工区的储层描述和流体预测。可选方法如下：

将反演得到的结果进行道积分，以测井资料计算的流体因子和密度的低频成分作为带限反演结果的低频成份，最终获取补偿后的流体因子和密度参数。依据井上油、气、水层在该流体因子剖面和密度剖面上的显示特征，可对工区进行储层流体和特征描述，从而获得有利钻井区。

通过以上的实施方式的描述，本领域的技术人员可以清楚地了解到根据上述实施例的方法可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质(如rom/ram、磁碟、光盘)中，包括若干指令用以使得一台终端设备(可以是手机，计算机，服务器，或者网络设备等)执行本发明各个实施例所述的方法。

在本实施例中还提供了一种地层弹性参数的确定装置，该装置用于实现上述实施例及优选实施方式，已经进行过说明的不再赘述。如以下所使用的，术语“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

图3是根据本发明实施例的地层弹性参数的确定装置的结构框图，如图3所示，该装置包括：

确定模块32，用于根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗，其中，第一叠前反射系数近似式以纵波阻抗和横波阻抗为参数；

第一转换模块34，用于将目标纵波阻抗转换为目标纵波速度；

第二转换模块36，用于利用目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式，其中，目标叠前反射系数近似式以流体因子和密度为参数；

第一获取模块38，用于根据目标叠前反射系数近似式获取目标地层弹性参数，其中，目标地层弹性参数包括目标流体因子和目标密度。

可选地，确定模块包括：第一获取单元，用于获取纵横波速度和密度的测井曲线；第二获取单元，用于利用纵横波速度和密度的测井曲线获取纵横波速度比；第一构造单元，用于利用纵横波速度比、纵波阻抗和横波阻抗构造第一叠前反射系数近似式；第一转换单元，用于将第一叠前反射系数近似式转换为第一线性方程组，其中，第一线性方程组以纵波阻抗和横波阻抗为参数；第一求解单元，用于对第一线性方程组进行求解，得到目标纵波阻抗。

可选地，第二转换模块包括：第二构造单元，用于构造以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式；改写单元，用于利用目标纵波速度将第二叠前反射系数近似式改写为目标叠前反射系数近似式。

可选地，第一获取模块包括：第二转换单元，用于将目标叠前反射系数近似式转换为第二线性方程组，其中，第二线性方程组以流体因子和密度为参数；第二求解单元，用于对第二线性方程组进行求解，得到目标流体因子和目标密度。

可选地，第二求解单元包括：求解子单元，用于对第二线性方程组进行求解，得到求解结果，其中，求解结果包括：求解流体因子和求解目标密度；处理子单元，用于利用获取到的测井数据对求解流体因子和求解目标密度进行道积分和低频补偿处理，得到第一处理结果，其中，第一处理结果包括：目标流体因子和目标密度。

可选地，上述装置还包括：处理模块，用于利用目标地层弹性参数对工区进行储层流体预测和特征描述处理，得到第二处理结果；第二获取模块，用于根据第二处理结果获取有利钻井区。

需要说明的是，上述各个模块是可以通过软件或硬件来实现的，对于后者，可以通过以下方式实现，但不限于此：上述模块均位于同一处理器中；或者，上述各个模块以任意组合的形式分别位于不同的处理器中。

下面结合本发明可选实施例进行详细说明。

本发明可选实施例提供了一种地层弹性参数的确定方法，图4是根据本发明可选实施例的地层弹性参数的确定方法的示意图，如图4所示，该方法包括如下流程：

步骤1，获得角度域道集，利用反褶积得到角度域反射系数，同时获得纵横波速度和密度的测井曲线；

步骤2，利用fatti两项avo近似式反演纵横波阻抗；

步骤3，利用gardner公式将纵波阻抗转换为纵波速度；

步骤4，利用得到的纵波速度，重新改写推导的以纵波速度、流体因子和密度为参数的avo近似式为仅含流体因子和密度的两项avo近似式，用该近似式反演流体因子和密度；

步骤5，利用测井数据对反演结果进行道积分和低频补偿，得到最终的流体因子和密度，利用流体因子和密度进行工区的储层描述和流体预测。

可选地，在上述步骤1中，可以但不限于通过以下方式获得纵横波速度和密度的测井曲线：

基于野外采集的地震资料利用地震成像方法提取叠前地震角道集数据；

利用叠前地震地质标定方法获得角道集中的子波，利用最小平方脉冲反褶积方法去除角度域道集的子波，获得角度域反射系数。

利用声波测井曲线、偶极声波测井曲线和密度测井曲线获得纵横波速度和密度的测井曲线

可选地，在上述步骤2中，可以但不限于通过以下方式反演纵横波阻抗：

在aki-richards近似式的基础上，fatti推导了当入射角小于35°，只包含纵波阻抗和横波阻抗的avo近似式：

其中，ip和is分别表示平均纵波阻抗和横波阻抗，θ是入射角度，δip和δis表示两层介质边界处的纵横波阻抗变化值，γsat为饱和岩石的纵横波速度比。

上式可构造出n个线性方程：

其中，

采用svd(奇异值分解)算法对公式(b)进行求解，从而得到纵波阻抗和横波阻抗。

可选地，在上述步骤3中，可以但不限于通过以下方式将纵波阻抗转换为纵波速度：

根据gardner公式可得到纵波速度与纵波阻抗间的关系式：

其中，a可将密度和纵波速度测井数据经最小二乘拟合获得，公式(c)可将纵波阻抗转换为纵波速度。

可选地，在上述步骤4中，可以但不限于通过以下方式反演流体因子和密度：

本实施例中重新推导了新的以纵波速度、流体因子和密度为参数的avo近似式：

其中，f、ρ、vp分别表示平均流体因子、密度和纵波速度，δf、δρ、δvp表示两层介质边界处流体因子、密度和纵波速度的变化量，γsat、γdry分别为饱和岩石和干岩石的纵横波速度比。

利用步骤3中得到的纵波速度，公式(4)可写为：

公式(e)可写为含有两个未知量的线性方程组：

其中，

采用svd(奇异值分解)算法对方程组(f)进行求解，从而反演得到流体因子和密度参数。

可选地，在上述步骤5中，可以但不限于通过以下方式利用流体因子和密度进行工区的储层描述和流体预测：

将反演得到的结果进行道积分，以测井资料计算的流体因子和密度的低频成分作为带限反演结果的低频成份，最终获取补偿后的流体因子和密度参数。依据井上油、气、水层在该流体因子剖面和密度剖面上的显示特征，可对工区进行储层流体和特征描述，从而获得有利钻井区。

在一个可选的实施方式中，通过实际工区的某口测井纵横波速度和密度数据合成地震叠前角道集对上述方法进行测试，图5是根据本发明可选的实施方式的地震角道集的示意图，如图5所示，左侧为合成的地震角道集，右侧为合成的加噪后的角道集。图6是根据本发明可选的实施方式的反演的流体因子的示意图，如图6所示，左侧是加噪前russell和上述反演方法反演的流体因子，两种方法得到的结果都与真实值一致。右侧是在叠前道集中加噪后的反演结果对比情况，基于russell近似式反演的流体因子与偏离真实值较大且信噪比较差，而上述方法反演的流体因子的精度和信噪比未受到明显的影响。图7是根据本发明可选的实施方式的反演的密度的示意图，如图7所示，左侧是加噪前russell和上述反演方法反演的密度，两种方法得到的结果都与真实值一致。右侧是在叠前道集中加噪后的反演结果对比情况，上述方法反演的密度精度和信噪比未受到影响、基于russell近似式反演的密度与密度实际值偏差较大。可见，上述方法具有较高的反演稳定性，能得到高精度的流体因子和密度反演结果。

综上所述，本可选实施例提供上述反演方法充分利用了两参数avo反演的稳定性特性，推导了仅含流体因子和密度的两项avo近似式，用该近似式进行后续反演，即使用含噪的叠前道集来反演，结果依然可靠，保证了从叠前地震资料中提取的流体因子和密度的可靠性，为储层流体检测和特性描述方面提供了可靠的地震资料，提高了储层流体预测、特征描述的准确性和有效性。

本发明的实施例还提供了一种存储介质，该存储介质中存储有计算机程序，其中，该计算机程序被设置为运行时执行上述任一项方法实施例中的步骤。

可选地，在本实施例中，上述存储介质可以被设置为存储用于执行以下步骤的计算机程序：

s1，根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗，其中，第一叠前反射系数近似式以纵波阻抗和横波阻抗为参数；

s2，将目标纵波阻抗转换为目标纵波速度；

s3，利用目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式，其中，目标叠前反射系数近似式以流体因子和密度为参数；

s4，根据目标叠前反射系数近似式获取目标地层弹性参数，其中，目标地层弹性参数包括目标流体因子和目标密度。

可选地，在本实施例中，上述存储介质可以包括但不限于：u盘、只读存储器(read-onlymemory，简称为rom)、随机存取存储器(randomaccessmemory，简称为ram)、移动硬盘、磁碟或者光盘等各种可以存储计算机程序的介质。

本发明的实施例还提供了一种电子装置，包括存储器和处理器，该存储器中存储有计算机程序，该处理器被设置为运行计算机程序以执行上述任一项方法实施例中的步骤。

可选地，上述电子装置还可以包括传输设备以及输入输出设备，其中，该传输设备和上述处理器连接，该输入输出设备和上述处理器连接。

可选地，在本实施例中，上述处理器可以被设置为通过计算机程序执行以下步骤：

s1，根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗，其中，第一叠前反射系数近似式以纵波阻抗和横波阻抗为参数；

s2，将目标纵波阻抗转换为目标纵波速度；

s3，利用目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式，其中，目标叠前反射系数近似式以流体因子和密度为参数；

s4，根据目标叠前反射系数近似式获取目标地层弹性参数，其中，目标地层弹性参数包括目标流体因子和目标密度。

可选地，本实施例中的可选示例可以参考上述实施例及可选实施方式中所描述的示例，本实施例在此不再赘述。

显然，本领域的技术人员应该明白，上述的本发明的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明不限制于任何特定的硬件和软件结合。

以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。