CN111060964B - 地层弹性参数的确定方法及装置 - Google Patents

地层弹性参数的确定方法及装置 Download PDF

Info

Publication number
CN111060964B
CN111060964B CN201811204319.0A CN201811204319A CN111060964B CN 111060964 B CN111060964 B CN 111060964B CN 201811204319 A CN201811204319 A CN 201811204319A CN 111060964 B CN111060964 B CN 111060964B
Authority
CN
China
Prior art keywords
target
reflection coefficient
longitudinal wave
density
parameters
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201811204319.0A
Other languages
English (en)
Other versions
CN111060964A (zh
Inventor
陈刚
王小军
赵春雪
徐倩茹
张田
罗鑫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201811204319.0A priority Critical patent/CN111060964B/zh
Publication of CN111060964A publication Critical patent/CN111060964A/zh
Application granted granted Critical
Publication of CN111060964B publication Critical patent/CN111060964B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48—Processing data
    • G01V1/50—Analysing data
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/62—Physical property of subsurface
    • G01V2210/624—Reservoir parameters
    • G01V2210/6242—Elastic parameters, e.g. Young, Lamé or Poisson

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种地层弹性参数的确定方法及装置，该方法包括：根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗，其中，第一叠前反射系数近似式以纵波阻抗和横波阻抗为参数；将目标纵波阻抗转换为目标纵波速度；利用目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式，其中，目标叠前反射系数近似式以流体因子和密度为参数；根据目标叠前反射系数近似式获取目标地层弹性参数，其中，目标地层弹性参数包括目标流体因子和目标密度。通过本发明，解决了相关技术中地层弹性参数的确定效率较低的问题，进而达到了提高地层弹性参数的确定效率的效果。

Description

地层弹性参数的确定方法及装置
技术领域
本发明涉及计算机领域，具体而言，涉及一种地层弹性参数的确定方法及装置。
背景技术
叠前地震反演是储层预测和油气检测的有效手段，基于叠前反演的流体因子只包含孔隙岩石里面的流体信息，不包含岩石的岩性性质，对储层中的流体性质区分性强，可区分储层中的油、气、水层，而密度参数与储层的物性紧密相联有助于储层特征描述。因此，流体因子和密度对储层流体预测和特征描述具有重要意义。
叠前三参数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包括一个网络适配器(Network Interface Controller，简称为NIC)，其可通过基站与其他网络设备相连从而可与互联网进行通讯。在一个实例中，传输装置106可以为射频(Radio Frequency，简称为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近似式为：
Figure BDA0001830831020000081
其中，Ip和Is分别表示平均纵波阻抗和横波阻抗，θ是入射角度，ΔIp和ΔIs表示两层介质边界处的纵横波阻抗变化值，γsat为饱和岩石的纵横波速度比。
上式(1)可构造出n个线性方程：
Figure BDA0001830831020000082
其中，
Figure BDA0001830831020000083
采用SVD(奇异值分解)算法对公式(2)进行求解，从而得到纵波阻抗和横波阻抗。
可选地，在本实施例中，在上述步骤S204中，可以通过以下方式将目标纵波阻抗转换成目标纵波速度：
密度和速度之间的数值关系式为：
Figure BDA0001830831020000084
上式(3)经变换可写为：
Figure BDA0001830831020000091
公式(4)左右两边同时加纵波速度差，则可给出：
Figure BDA0001830831020000092
其中，a可将密度和纵波速度测井数据经最小二乘拟合获得，利用公式(5)可将纵波阻抗转换为纵波速度。
可选地，在上述步骤S206中，构造以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式；利用目标纵波速度将第二叠前反射系数近似式改写为目标叠前反射系数近似式。
可选地，在本实施例中，利用得到的纵波速度，重新改写推导的以纵波速度、流体因子和密度为参数的AVO近似式为仅含流体因子和密度的两项AVO近似式，用该近似式反演流体因子和密度；可选方法如下：
为了从叠前地震数据中提取流体因子，以流体因子、剪切模量和密度为参数的AVO近似式：
Figure BDA0001830831020000093
其中，f、ρ、μ分别表示平均流体因子、密度和剪切模量，Δf、Δρ、Δμ表示两层介质边界处流体因子、密度和剪切模量的变化量，γsat、γdry分别为饱和岩石和干岩石的纵横波速度比。
但由于三参数AVO反演具有稳定性差的问题。因此，重新推导了新的以纵波速度、流体因子和密度为参数的AVO近似式：
Figure BDA0001830831020000101
利用上述得到的纵波速度，公式(7)可写为：
Figure BDA0001830831020000102
纵波速度作为已知量，含有两个未知量的线性方程组可表示为：
Figure BDA0001830831020000103
其中，
Figure BDA0001830831020000104
可选地，在上述步骤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近似式：
Figure BDA0001830831020000141
其中，Ip和Is分别表示平均纵波阻抗和横波阻抗，θ是入射角度，ΔIp和ΔIs表示两层介质边界处的纵横波阻抗变化值，γsat为饱和岩石的纵横波速度比。
上式可构造出n个线性方程：
Figure BDA0001830831020000142
其中，
Figure BDA0001830831020000143
采用SVD(奇异值分解)算法对公式(b)进行求解，从而得到纵波阻抗和横波阻抗。
可选地，在上述步骤3中，可以但不限于通过以下方式将纵波阻抗转换为纵波速度：
根据Gardner公式
Figure BDA0001830831020000151
可得到纵波速度与纵波阻抗间的关系式：
Figure BDA0001830831020000152
其中，a可将密度和纵波速度测井数据经最小二乘拟合获得，公式(c)可将纵波阻抗转换为纵波速度。
可选地，在上述步骤4中，可以但不限于通过以下方式反演流体因子和密度：
本实施例中重新推导了新的以纵波速度、流体因子和密度为参数的AVO近似式：
Figure BDA0001830831020000153
其中，f、ρ、Vp分别表示平均流体因子、密度和纵波速度，Δf、Δρ、ΔVp表示两层介质边界处流体因子、密度和纵波速度的变化量，γsat、γdry分别为饱和岩石和干岩石的纵横波速度比。
利用步骤3中得到的纵波速度，公式(4)可写为：
Figure BDA0001830831020000154
公式(e)可写为含有两个未知量的线性方程组：
Figure BDA0001830831020000161
其中，
Figure BDA0001830831020000162
采用SVD(奇异值分解)算法对方程组(f)进行求解，从而反演得到流体因子和密度参数。
可选地，在上述步骤5中，可以但不限于通过以下方式利用流体因子和密度进行工区的储层描述和流体预测：
将反演得到的结果进行道积分，以测井资料计算的流体因子和密度的低频成分作为带限反演结果的低频成份，最终获取补偿后的流体因子和密度参数。依据井上油、气、水层在该流体因子剖面和密度剖面上的显示特征，可对工区进行储层流体和特征描述，从而获得有利钻井区。
在一个可选的实施方式中，通过实际工区的某口测井纵横波速度和密度数据合成地震叠前角道集对上述方法进行测试，图5是根据本发明可选的实施方式的地震角道集的示意图，如图5所示，左侧为合成的地震角道集，右侧为合成的加噪后的角道集。图6是根据本发明可选的实施方式的反演的流体因子的示意图，如图6所示，左侧是加噪前Russell和上述反演方法反演的流体因子，两种方法得到的结果都与真实值一致。右侧是在叠前道集中加噪后的反演结果对比情况，基于Russell近似式反演的流体因子与偏离真实值较大且信噪比较差，而上述方法反演的流体因子的精度和信噪比未受到明显的影响。图7是根据本发明可选的实施方式的反演的密度的示意图，如图7所示，左侧是加噪前Russell和上述反演方法反演的密度，两种方法得到的结果都与真实值一致。右侧是在叠前道集中加噪后的反演结果对比情况，上述方法反演的密度精度和信噪比未受到影响、基于Russell近似式反演的密度与密度实际值偏差较大。可见，上述方法具有较高的反演稳定性，能得到高精度的流体因子和密度反演结果。
综上所述，本可选实施例提供上述反演方法充分利用了两参数AVO反演的稳定性特性，推导了仅含流体因子和密度的两项AVO近似式，用该近似式进行后续反演，即使用含噪的叠前道集来反演，结果依然可靠，保证了从叠前地震资料中提取的流体因子和密度的可靠性，为储层流体检测和特性描述方面提供了可靠的地震资料，提高了储层流体预测、特征描述的准确性和有效性。
本发明的实施例还提供了一种存储介质，该存储介质中存储有计算机程序，其中，该计算机程序被设置为运行时执行上述任一项方法实施例中的步骤。
可选地，在本实施例中，上述存储介质可以被设置为存储用于执行以下步骤的计算机程序：
S1，根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗，其中，第一叠前反射系数近似式以纵波阻抗和横波阻抗为参数；
S2，将目标纵波阻抗转换为目标纵波速度；
S3，利用目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式，其中，目标叠前反射系数近似式以流体因子和密度为参数；
S4，根据目标叠前反射系数近似式获取目标地层弹性参数，其中，目标地层弹性参数包括目标流体因子和目标密度。
可选地，在本实施例中，上述存储介质可以包括但不限于：U盘、只读存储器(Read-Only Memory，简称为ROM)、随机存取存储器(Random Access Memory，简称为RAM)、移动硬盘、磁碟或者光盘等各种可以存储计算机程序的介质。
本发明的实施例还提供了一种电子装置，包括存储器和处理器，该存储器中存储有计算机程序，该处理器被设置为运行计算机程序以执行上述任一项方法实施例中的步骤。
可选地，上述电子装置还可以包括传输设备以及输入输出设备，其中，该传输设备和上述处理器连接，该输入输出设备和上述处理器连接。
可选地，在本实施例中，上述处理器可以被设置为通过计算机程序执行以下步骤：
S1，根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗，其中，第一叠前反射系数近似式以纵波阻抗和横波阻抗为参数；
S2，将目标纵波阻抗转换为目标纵波速度；
S3，利用目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式，其中，目标叠前反射系数近似式以流体因子和密度为参数；
S4，根据目标叠前反射系数近似式获取目标地层弹性参数，其中，目标地层弹性参数包括目标流体因子和目标密度。
可选地，本实施例中的可选示例可以参考上述实施例及可选实施方式中所描述的示例，本实施例在此不再赘述。
显然，本领域的技术人员应该明白，上述的本发明的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个的计算装置上，或者分布在多个计算装置所组成的网络上，可选地，它们可以用计算装置可执行的程序代码来实现，从而，可以将它们存储在存储装置中由计算装置来执行，并且在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明不限制于任何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种地层弹性参数的确定方法，其特征在于，包括：
根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗，其中，所述第一叠前反射系数近似式是以纵波阻抗和横波阻抗为参数的AVO近似式；
将所述目标纵波阻抗转换为目标纵波速度；
利用所述目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式，其中，所述第二叠前反射系数近似式为AVO近似式，所述目标叠前反射系数近似式是以流体因子和密度为参数的AVO近似式；
根据所述目标叠前反射系数近似式获取目标地层弹性参数，其中，所述目标地层弹性参数包括目标流体因子和目标密度。
2.根据权利要求1所述的方法，其特征在于，根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗包括：
获取所述纵横波速度和密度的测井曲线；
利用所述纵横波速度和密度的测井曲线获取纵横波速度比；
利用所述纵横波速度比、纵波阻抗和横波阻抗构造所述第一叠前反射系数近似式；
将所述第一叠前反射系数近似式转换为第一线性方程组，其中，所述第一线性方程组以纵波阻抗和横波阻抗为参数；
对所述第一线性方程组进行求解，得到所述目标纵波阻抗。
3.根据权利要求1所述的方法，其特征在于，利用所述目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式包括：
构造以纵波速度、流体因子和密度为参数的所述第二叠前反射系数近似式；
利用所述目标纵波速度将所述第二叠前反射系数近似式改写为所述目标叠前反射系数近似式。
4.根据权利要求1所述的方法，其特征在于，根据所述目标叠前反射系数近似式获取目标地层弹性参数包括：
将所述目标叠前反射系数近似式转换为第二线性方程组，其中，所述第二线性方程组以流体因子和密度为参数；
对所述第二线性方程组进行求解，得到所述目标流体因子和所述目标密度。
5.根据权利要求4所述的方法，其特征在于，对所述第二线性方程组进行求解，得到所述目标流体因子和所述目标密度包括：
对所述第二线性方程组进行求解，得到求解结果，其中，所述求解结果包括：求解流体因子和求解目标密度；
利用获取到的测井数据对所述求解流体因子和所述求解目标密度进行道积分和低频补偿处理，得到第一处理结果，其中，所述第一处理结果包括：所述目标流体因子和所述目标密度。
6.根据权利要求1至5中任一项所述的方法，其特征在于，在根据所述目标叠前反射系数近似式获取目标地层弹性参数之后，所述方法还包括：
利用所述目标地层弹性参数对工区进行储层流体预测和特征描述处理，得到第二处理结果；
根据所述第二处理结果获取有利钻井区。
7.一种地层弹性参数的确定装置，其特征在于，包括：
确定模块，用于根据获取的纵横波速度和密度的测井曲线以及第一叠前反射系数近似式确定目标纵波阻抗，其中，所述第一叠前反射系数近似式是以纵波阻抗和横波阻抗为参数的AVO近似式；
第一转换模块，用于将所述目标纵波阻抗转换为目标纵波速度；
第二转换模块，用于利用所述目标纵波速度将以纵波速度、流体因子和密度为参数的第二叠前反射系数近似式转换为目标叠前反射系数近似式，其中，所述第二叠前反射系数近似式为AVO近似式，所述目标叠前反射系数近似式是以流体因子和密度为参数的AVO近似式；
第一获取模块，用于根据所述目标叠前反射系数近似式获取目标地层弹性参数，其中，所述目标地层弹性参数包括目标流体因子和目标密度。
8.根据权利要求7所述的装置，其特征在于，所述装置还包括：
处理模块，用于利用所述目标地层弹性参数对工区进行储层流体预测和特征描述处理，得到第二处理结果；
第二获取模块，用于根据所述第二处理结果获取有利钻井区。
9.一种存储介质，其特征在于，所述存储介质中存储有计算机程序，其中，所述计算机程序被设置为运行时执行所述权利要求1至6任一项中所述的方法。
10.一种电子装置，包括存储器和处理器，其特征在于，所述存储器中存储有计算机程序，所述处理器被设置为运行所述计算机程序以执行所述权利要求1至6任一项中所述的方法。
CN201811204319.0A 2018-10-16 2018-10-16 地层弹性参数的确定方法及装置 Active CN111060964B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811204319.0A CN111060964B (zh) 2018-10-16 2018-10-16 地层弹性参数的确定方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811204319.0A CN111060964B (zh) 2018-10-16 2018-10-16 地层弹性参数的确定方法及装置

Publications (2)

Publication Number Publication Date
CN111060964A CN111060964A (zh) 2020-04-24
CN111060964B true CN111060964B (zh) 2022-05-10

Family

ID=70296606

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811204319.0A Active CN111060964B (zh) 2018-10-16 2018-10-16 地层弹性参数的确定方法及装置

Country Status (1)

Country Link
CN (1) CN111060964B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118011485A (zh) * 2022-11-09 2024-05-10 中国石油天然气股份有限公司 页岩油双甜点关键参数一体化反演方法及装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10505669A (ja) * 1994-09-02 1998-06-02 エクソン プロダクション リサーチ カンパニー 地震データのプリスタックインバージョンから貯留層岩質及び流体内容を導出する方法
CN101446645A (zh) * 2007-11-27 2009-06-03 中国石油天然气股份有限公司 一种利用地震流体阻抗进行流体确定的方法
CN102466816A (zh) * 2010-11-04 2012-05-23 中国石油天然气集团公司 一种叠前地震数据地层弹性常数参数反演的方法
CN103713319A (zh) * 2013-12-31 2014-04-09 张远银 一种基于地震约束建模的叠前反演方法
CN106842310A (zh) * 2015-12-04 2017-06-13 中国石油化工股份有限公司 叠前地震四参数同步反演方法
CN106842311A (zh) * 2015-12-04 2017-06-13 中国石油化工股份有限公司 粘弹性流体因子直接反演方法
CN107367762A (zh) * 2017-06-27 2017-11-21 中国石油天然气股份有限公司 一种确定储层参数的方法及装置

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10505669A (ja) * 1994-09-02 1998-06-02 エクソン プロダクション リサーチ カンパニー 地震データのプリスタックインバージョンから貯留層岩質及び流体内容を導出する方法
EP0786098A4 (en) * 1994-09-02 1998-10-07 Exxon Production Research Co METHOD FOR DEDUCING THE LITHOLOGICAL CHARACTERISTICS AND THE FLUID CONTENT OF A TANK OF THE INVERSION BEFORE SUMMATION OF SEISMIC DATA
CN101446645A (zh) * 2007-11-27 2009-06-03 中国石油天然气股份有限公司 一种利用地震流体阻抗进行流体确定的方法
CN102466816A (zh) * 2010-11-04 2012-05-23 中国石油天然气集团公司 一种叠前地震数据地层弹性常数参数反演的方法
CN103713319A (zh) * 2013-12-31 2014-04-09 张远银 一种基于地震约束建模的叠前反演方法
CN106842310A (zh) * 2015-12-04 2017-06-13 中国石油化工股份有限公司 叠前地震四参数同步反演方法
CN106842311A (zh) * 2015-12-04 2017-06-13 中国石油化工股份有限公司 粘弹性流体因子直接反演方法
CN107367762A (zh) * 2017-06-27 2017-11-21 中国石油天然气股份有限公司 一种确定储层参数的方法及装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
AVO流体反演技术应用研究;陈刚等;《石油天然气学报》;20111031;第33卷(第10期);第70-75页 *
PP and PS waves joint AVO inversion of density and shear modulus;Q.Z. Du et.al;《76rd EAGE International Annual Meeting》;20141231;第1-5页 *
Two-step joint PP- and PS-wave three-term amplitude-variation with offset inversion;Qizhen Du;《Interpretation》;20161231;T613-T625 *
纵波与转换横波联合AVO反演轨迹横波速度和密度;陈刚等;《断块油气田》;20151130;第22卷(第6期);第729-734页 *

Also Published As

Publication number Publication date
CN111060964A (zh) 2020-04-24

Similar Documents

Publication Publication Date Title
US8868348B2 (en) Well constrained horizontal variable H-V curve constructing method for seismic wave velocity field construction
CN109581501B (zh) 用于沙漠区深度域速度建模的方法
CN111123354B (zh) 基于频变反射振幅衰减预测致密气层的方法及设备
FR3058180A1 (fr) Determination en temps reel de la lenteur de la boue, du type de formation, et des pointes de lenteur monopolaire dans des applications en fond de puits
CN103245972B (zh) 一种在二维空间内确定复杂地质构造的方法
CN110887772B (zh) 一种碳酸盐岩储层渗透率识别方法、系统及装置
CN103850679A (zh) 一种利用多种测井曲线对声波时差曲线进行重构的方法
Fahimuddin et al. 4D seismic history matching of a real field case with EnKF: use of local analysis for model updating
CN104316966A (zh) 一种流体识别方法及系统
CN107576985B (zh) 一种地震反演的方法和装置
CN111679318B (zh) 储层中断裂系统的识别方法及系统
CN102830170B (zh) 一种基于超声测试获取煤样横波信号的控制方法及装置
CN111060964B (zh) 地层弹性参数的确定方法及装置
WO2017120088A1 (en) Effective permeability upscaling for a discrete fracture network
CN106707345B (zh) 一种角度弹性参数岩性识别方法和装置
CN107991705B (zh) 基于二维统计特征的测井曲线校正方法和装置
CN111502647B (zh) 一种钻井地质环境因素确定方法、装置及存储介质
CN112649848B (zh) 利用波动方程求解地震波阻抗的方法和装置
CN116146202B (zh) 一种礁滩碳酸盐岩的圈闭划分方法和装置
CN111366974B (zh) 基于频变固液解耦流体因子的叠前地震反演方法和装置
CN118226537B (zh) 大洋钻探套管直达波压制方法、装置、设备和存储介质
CN111830562A (zh) 一种油气储层渗透率预测方法及装置
CN103760598A (zh) 一种射线弹性参数反演方法
CN113219533B (zh) 一种叠前时间偏移速度建模方法、装置、介质及电子设备
CN112363244B (zh) 波阻抗反演方法与碳酸盐岩非均质储层预测方法及系统

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant