CN111077574A - 一种确定地层弹性参数的方法、装置及*** - Google Patents
一种确定地层弹性参数的方法、装置及*** Download PDFInfo
- Publication number
- CN111077574A CN111077574A CN201911398538.1A CN201911398538A CN111077574A CN 111077574 A CN111077574 A CN 111077574A CN 201911398538 A CN201911398538 A CN 201911398538A CN 111077574 A CN111077574 A CN 111077574A
- Authority
- CN
- China
- Prior art keywords
- parameter
- data
- parameters
- fluid
- wave
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 78
- 239000012530 fluid Substances 0.000 claims abstract description 143
- 238000012545 processing Methods 0.000 claims abstract description 30
- 230000006870 function Effects 0.000 claims description 42
- 239000011159 matrix material Substances 0.000 claims description 30
- 230000008859 change Effects 0.000 claims description 27
- 230000015572 biosynthetic process Effects 0.000 claims description 22
- 239000011435 rock Substances 0.000 claims description 18
- 230000005540 biological transmission Effects 0.000 claims description 12
- 230000015654 memory Effects 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 8
- 229920006395 saturated elastomer Polymers 0.000 claims description 6
- 238000005755 formation reaction Methods 0.000 description 19
- 238000010586 diagram Methods 0.000 description 15
- 238000004364 calculation method Methods 0.000 description 14
- 230000008569 process Effects 0.000 description 12
- 230000035945 sensitivity Effects 0.000 description 12
- 239000011148 porous material Substances 0.000 description 8
- 238000011156 evaluation Methods 0.000 description 7
- 238000004590 computer program Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000012512 characterization method Methods 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000010835 comparative analysis Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000036039 immunity Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 239000000047 product Substances 0.000 description 2
- 230000000750 progressive effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 229910021389 graphene Inorganic materials 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 230000001568 sexual effect Effects 0.000 description 1
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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- 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
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
技术领域
本发明涉及地震勘探处理技术领域，特别地，涉及一种确定地层弹性参数的方法、装置及***。
背景技术
地震储层和流体识别具体指通过地震数据所携带的信息，对地下储层和储层内部的流体进行检测和描述的过程。由地震波数据提取的弹性参数中蕴含了丰富的有助于揭示地下地质体以及储层内部孔隙流体特征的信息。随着勘探难度的加大以及勘探技术的进步，基于弹性参数的地下岩层岩性以及孔隙流体识别已成为现阶段油气勘探的关键技术。利用地震信息进行储层和其孔隙内填充流体性质的识别研究，需要结合岩石物理理论，找出可以突出储层和流体的表征手段，然后进行进一步的综合预测和描述。在高品质的地震数据基础上，地下储层和孔隙流体识别的精度主要依赖于选取的表征参数对于储层或流体的敏感程度，以及提取相关表征参数方法的可靠性。通过对弹性参数的综合分析是探索地下储层和孔隙流体特征，对地下岩层的情况进行定性乃至定量的描述的一个有效途径。
国内外学者基于弹性参数针对储层和流体检测的研究也做了大量的研究。其中，Russell等结合了孔隙弹性理论提出了新型流体因子Pfluid，可以用来直接进行流体类型识别。剪切模量μ代表了储层骨架的特性，可以用来描述岩石的岩性特征，密度ρ也能显示储层的岩性特征，Goodway等阐述了两者的乘积μρ＝Preservoir在识别储层岩性方面有着较好的表现。而目前提取参数Pfluid、Preservoir的方法基本均依靠于纵波反演间接计算的方法，使得参数提取的精度受到限制。
发明内容
本说明书实施例的目的在于提供一种确定地层弹性参数的方法、装置及***，可以大幅度提高反演获得弹性参数的准确性。
本说明书提供一种确定地层弹性参数的方法、装置及***是包括如下方式实现的：
一种确定地层弹性参数的方法，包括：
获取目标工区的实测的纵波分角度地震数据以及转换波分角度地震数据；
获取用于表征实测地震数据与基于弹性参数计算的理论地震数据之间接近关系的目标函数，所述目标函数根据用于直接表征储层参数和流体参数分布特征的纵波近似模型以及转换波近似模型构建，所述弹性参数包括储层参数和流体参数；
利用所述目标函数对所述纵波分角度地震数据、转换波分角度地震数据进行反演处理，获得所述目标工区的弹性参数数据。
本说明书所述的方法的另一个实施例中，所述纵波近似模型包括：
Rresevoir＝2(μ2ρ2-μ1ρ1)/(μ2ρ2+μ1ρ1)
Rfluid＝2(L2-L1)/(L2+L2)
其中，r(α)表示纵波数据，α是纵波数据的反射角度和透射角度的平均，Rresevoir为储层参数变化率，μ2为上层介质的剪切模量，μ1为下层介质的剪切模量，ρ2为上层介质的密度参数，ρ1为下层介质的密度参数；Rfluid为流体参数变化率，L1为下层介质的流体参数，L2为下层介质的流体参数，k1为干岩石速度比的二次方，k2为饱和岩石的速度比二次方。
本说明书所述的方法的另一个实施例中，所述转换波近似模型包括：
其中，s(α,β)表示转换波数据，β是转换波数据的反射角度和透射角度的平均。
本说明书所述的方法的另一个实施例中，所述目标函数包括：
其中，r表示根据纵波近似模型确定的纵波数据；s表示根据转换波近似模型确定的转换波数据，W1为纵波分角度叠加数据体的子波组成的矩阵，W2为转换波分角度叠加数据体的子波组成的矩阵；C1为纵波分角度叠加数据体的子波组成的矩阵，C2为转换波分角度叠加数据体的子波组成的矩阵；Ω为相关矩阵；mt为初始的储层参数、流体参数和密度参数的低频分量；G为低频因子，Gm为待求的储层参数、流体参数和密度参数的低频分量；λ1、λ2、λ3为调节权重参数；N为每一个地震道的采样点数。
本说明书所述的方法的另一个实施例中，所述利用所述目标函数对所述纵波分角度地震数据、转换波分角度地震数据进行反演处理，获得所述目标工区的弹性参数数据，包括：
根据所述目标函数对所述纵波分角度地震数据、转换波分角度地震数据进行反演处理，获得储层参数变化率和流体参数变化率；
对所述储层参数变化率和流体参数变化率进行积分，获得储层参数以及流体参数。
另一方面，本说明书实施例还提供一种确定地层弹性参数的装置，包括：
数据获取模块，用于获取目标工区的实测的纵波分角度地震数据以及转换波分角度地震数据；
反演函数获取模块，用于获取用于表征实测地震数据与基于弹性参数计算的理论地震数据之间接近关系的目标函数，所述目标函数根据用于直接表征储层参数和流体参数分布特征的纵波近似模型以及转换波近似模型构建，所述弹性参数包括储层参数和流体参数；
反演模块，用于利用所述目标函数对所述纵波分角度地震数据、转换波分角度地震数据进行反演处理，获得所述目标工区的弹性参数数据。
本说明书所述的装置的另一个实施例中，所述纵波近似模型包括：
Rresevoir＝2(μ2ρ2-μ1ρ1)/(μ2ρ2+μ1ρ1)
Rfluid＝2(L2-L1)/(L2+L2)
其中，r(α)表示纵波数据，α是纵波数据的反射角度和透射角度的平均，Rresevoir为储层参数变化率，μ2为上层介质的剪切模量，μ1为下层介质的剪切模量，ρ2为上层介质的密度参数，ρ1为下层介质的密度参数；Rfluid为流体参数变化率，L1为下层介质的流体参数，L2为下层介质的流体参数，k1为干岩石速度比的二次方，k2为饱和岩石的速度比二次方。
本说明书所述的装置的另一个实施例中，所述转换波近似模型包括：
其中，s(α,β)表示转换波数据，β是转换波数据的反射角度和透射角度的平均。
本说明书所述的装置的另一个实施例中，所述目标函数包括：
其中，r表示根据纵波近似模型确定的纵波数据；s表示根据转换波近似模型确定的转换波数据，W1为纵波分角度叠加数据体的子波组成的矩阵，W2为转换波分角度叠加数据体的子波组成的矩阵；C1为纵波分角度叠加数据体的子波组成的矩阵，C2为转换波分角度叠加数据体的子波组成的矩阵；Ω为相关矩阵；mt为初始的储层参数、流体参数和密度参数的低频分量；G为低频因子，Gm为待求的储层参数、流体参数和密度参数的低频分量；λ1、λ2、λ3为调节权重参数；N为每一个地震道的采样点数。
另一方面，本说明书实施例还提供一种确定地层弹性参数的***，所述***包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现上述任意一个实施例所述方法的步骤。
本说明书一个或多个实施例提供的确定地层弹性参数的方法、装置及***，可以利用直接表征储层参数和流体参数分布特征的反演模型，在分角度叠加数据体的基础上直接计算得到表征储层和流体的弹性参数，不需要间接计算，降低了间接计算引入的累积误差。并且在反演过程中加入了转换波信息，使得反演结果的精度得到较大的提高。
附图说明
为了更清楚地说明本说明书实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本说明书中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本说明书提供的一种确定地层弹性参数的方法实施例的流程示意图；
图2为本说明书提供的一个实施例中的弹性参数确定方法流程示意图；
图3为本说明书提供的另一个实施例中的不同流体指示因子的敏感度示意图；
图4为本说明书提供的另一个实施例中的反演测试理论数据及初始模型示意图；
图5为本说明书提供的另一个实施例中的未含噪声的角道集示意图；
图6为本说明书提供的另一个实施例中的含有噪声的角道集示意图；
图7为本说明书提供的另一个实施例中的储层参数和流体参数反演结果示意图；
图8为本说明书提供的另一个实施例中的储层参数和流体参数反演结果示意图；
图9为本说明书提供的另一个实施例中的储层参数和流体参数反演结果示意图；
图10为本说明书提供的一种确定地层弹性参数的装置实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本说明书中的技术方案，下面将结合本说明书一个或多个实施例中的附图，对本说明书一个或多个实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是说明书一部分实施例，而不是全部的实施例。基于说明书一个或多个实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本说明书实施例方案保护的范围。
图1是本说明书提供的所述确定地层弹性参数的方法实施例流程示意图。虽然本说明书提供了如下述实施例或附图所示的方法操作步骤或装置结构，但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者部分合并后更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中，这些步骤的执行顺序或装置的模块结构不限于本说明书实施例或附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置、服务器或终端产品应用时，可以按照实施例或者附图所示的方法或模块结构进行顺序执行或者并行执行(例如并行处理器或者多线程处理的环境、甚至包括分布式处理、服务器集群的实施环境)。
具体的一个实施例如图1所示，本说明书提供的确定地层弹性参数的方法的一个实施例中，所述方法可以应用于执行地层的弹性参数反演处理的处理设备，所述方法可以包括：
S20：获取目标工区的实测的纵波分角度地震数据以及转换波分角度地震数据。
处理设备可以获取目标工区的纵波分角度地震数据以及转换波分角度地震数据。如可以获取目标工区的实测的纵波叠前地震数据以及转换波叠前地震数据，对地震数据进行包括振幅补偿处理、能量补偿及均衡处理、规则化等处理、噪声的去除、多次波去除、高分辨率、全三维去噪、道集拉平等优化处理。对优化处理后的叠前地震数据进行分角度叠加，获得目标工区的实测纵波分角度地震数据以及转换波分角度地震数据。
可以利用测井数据计算目的层段的理论角度道集数据，计算目的层段内理论道集的振幅随入射角变化(AVO)曲线；计算井点位置处实际地震道集的AVO曲线值。对合成道集的AVO曲线和实际地震道集的AVO曲线进行相关性分析，相关性高的角度范围即为有效角度，将有效入射角度αi(i＝1,2...M)依次等分为M个的分角度数据，其中M大于等于3。然后，分别进行叠加得到M个纵波分角度叠加数据体(r1,r2...rM)以及M个转换波分角度叠加数据体(s1,s2...sM)。
另一些实施例中，还可以根据所述纵波分角度地震数据以及转换波分角度地震数据分别提取纵波地震子波以及转换波地震子波。处理设备可以获取目标工区的实测纵波分角度地震数据以及转换波分角度地震数据，根据所述纵波分角度地震数据以及转换波分角度地震数据分别提取纵波地震子波以及转换波地震子波。
例如，可以利用测井数据制作井旁纵波合成地震道集，基于纵波合成地震道集以及观测得到的纵波分角度叠加数据体估算每个纵波分角度叠加数据体对应的子波(ωp1,ωp2...ωpM)。利用测井数据制作井旁转换横波合成地震道集，基于转换横波合成地震道集以及观测得到的转换横波分角度叠加数据体估算每个转换波分角度叠加数据体对应的子波(ωs1,ωs2...ωsM)。
S22：获取用于表征实测地震数据与基于弹性参数计算的理论地震数据之间接近关系的目标函数，所述目标函数根据用于直接表征储层参数和流体参数分布特征的纵波近似模型以及转换波近似模型构建。
处理设备可以获取目标函数，所述目标函数可以表示实测地震数据与根据弹性参数计算的理论地震数据之间的接近关系。所述接近关系可以是指实测地震数据与根据弹性参数计算的理论地震数据之间的差异性大小，差异性越小，则两个数据越接近。
一些实施例中，所述目标函数可以根据纵波近似模型以及转换波近似模型构建。其中，所述纵波近似模型、转换波近似模型可以为用于指示储层参数和流体参数分布特征的纵波数据反演模型、转换波数据反演模型，所述纵波近似模型以及转换波近似模型如可以根据岩石物理关系构建。
一些实施例中，构建的所述纵波近似模型可以包括：
Rresevoir＝2(μ2ρ2-μ1ρ1)/(μ2ρ2+μ1ρ1)
Rfluid＝2(L2-L1)/(L2+L2)
其中，r(α)表示纵波数据，α是纵波数据的反射角度和透射角度的平均，Rresevoir为储层参数变化率，μ2为上层介质的剪切模量，μ1为下层介质的剪切模量，ρ2为上层介质的密度参数，ρ1为下层介质的密度参数；Rfluid为流体参数变化率，L1为下层介质的流体参数，L2为下层介质的流体参数，L1、L2可以根据Gassmann流体参数与密度的乘积确定，k1为干岩石速度比的二次方，k2为饱和岩石的速度比二次方。
一些实施例中，构建的所述转换波近似模型可以包括：
其中，s(α,β)表示转换波数据，β是转换波数据的反射角度和透射角度的平均。
然后，可以利用测井数据、共反射点地震道集以及上述纵波近似模型、转换波近似模型构建目标函数。一些实施例中，所述目标函数可以表示为：
上式中，r表示纵波数据，可以根据上述纵波近似模型确定；s表示转换波数据，可以根据上述转换波近似模型确定；W1为纵波分角度叠加数据体的子波组成的矩阵，W2为转换波分角度叠加数据体的子波组成的矩阵；C1为纵波数据对应的积分矩阵，C2为转换波数据对应的积分矩阵；Ωj为相关矩阵；mt为初始的储层参数、流体参数和密度参数的低频分量；G为低频因子，Gm为待求的储层参数、流体参数和密度参数的低频分量；λ1、λ2、λ3为调节权重参数；N为每一个地震道的采样点数。
S24：利用所述目标函数对所述纵波分角度地震数据、转换波分角度地震数据进行反演处理，获得所述目标工区的弹性参数数据。
可以将所述纵波分角度地震数据、转换波分角度地震数据代入所述目标函数中，并求解目标函数的最小值，获得所述目标工区的弹性参数数据。一些实施例中，可以根据所述目标函数对所述纵波分角度地震数据、转换波分角度地震数据进行反演处理，获得储层参数变化率和流体参数变化率；对所述储层参数变化率和流体参数变化率进行积分，获得储层参数以及流体参数。
可以在初始模型的基础上对变化率参数积分，求得表征储层和流体的弹性参数，具体公式为：
其中，Presevoir为储层参数，Pfluid为流体参数，Pdensity为密度参数；Presevoir_1为初始模型中第一个采样点的储层参数值，Pfluid为初始模型中第一个采样点的流体参数值，为初始模型中第一个采样点的密度参数值。所述初始模型可以为根据井数据以及地质框架数据构建的低通滤波模型，用于表征反演求解中的弹性参数初始值以及弹性参数真实的低频趋势。
本说明书上述实施例提供的方案，通过利用直接表征储层参数和流体参数与纵波反射系数和转换波反射系数关系的纵波近似模型以及转换波近似模型，在测井曲线和分角度叠加数据体的基础上直接计算得到表征储层和流体的弹性参数，降低了间接计算引入的累积误差。同时，在构建的目标函数中进一步加入了转换波信息以及多项约束信息，使得反演结果的精度得到较大的提高，使得得到的储层参数和流体参数可以用来进行油藏的综合分析，指导后续的解释工作。
基于本说明书上述实施例提供的方案，本说明书还提供一种具体场景示例，以说明本说明书各实施例提供的方案的实用性。图2表示本说明书提供的一个场景示例中的性参数确定方法流程示意图。如图2所示，所述方法可以包括如下步骤：
S1：获取并处理研究区域的测井资料和叠前地震共反射点道集；
S2：根据测井数据和共反射点道集进行AVO属性分析，确定地震数据的分角度方案；
S3：基于处理后的测井数据和井旁地震道数据分别进行纵波以及转换波分角度叠加数据的子波提取；
S4：根据岩石物理关系构建能够指示储层和流体特性的纵波、转换波近似公式；
S5：利用测井数据、共反射点地震道集以及近似公式构建反演目标函数；
S6：根据新构建的目标函数求取储层参数变化率和流体参数变化率；
S7：在初始模型的基础上对反射率参数积分，求得表征储层和流体的弹性参数。
基于上述方案，本说明书还提供了具体的实例对比分析，如图3至图9所示，具体对比分析结果如下。首先，针对不同的地质类型区域，所选取的敏感流体因子都是不同的，为了对流体因子的敏感性进行定量评价，借鉴FIC(Fluid Indicator Coefficient)的定义，建立流体评价因子敏感系数Sfluid可以表征流体评价因子对流体的区分能力，具体表达式为：
上式中，S1和S2分别为岩层孔隙填充不同类型流体时的流体因子值。Sfluid的值在0-200之间，若S1和S2为同一种流体填充时，则Sfluid为0；S1和S2的差值越大，Sfluid的值就越趋近于200。Sfluid值越大表示识别因子对流体的敏感程度越高，即区分流体性质的能力越强，利用该评价参数可以对流体因子的敏感性进行定量评价。
在Hilterman给出的三类典型砂岩模型基础上进行了改编，得到了表1中的基础参数，结合流体评价因子敏感系数公式，可以计算得到常用流体评价因子的敏感系数(表2)。图3为根据表2绘制的不同流体指示因子在不同储层条件中的敏感度示意图。图3中线1代表第一类砂岩模型求得的流体判别因子敏感系数，线2代表第二类砂岩模型求得的流体判别因子敏感系数，线3代表第三类砂岩模型求得的流体判别因子敏感系数。图3中的横坐标号为流体评价因子序号，与表2相对应。从表2和图3可见，本说明书实施例计算的目标流体参数Pfluid相对其他流体因子在不同的储层内均有较敏感的指示，说明本说明书实施例依据的流体参数可以有效地指示孔隙流体类型。
表1三类砂岩模型参数
表2流体评价因子敏感系数
图4给出了反演测试的理论数据以及初始模型数据示意图。其中，1、3、5、7、9实线表示理论值，2、4、6、8、10实线表示反演过程中的初始模型值。1、2线代表纵波速度，3、4线代表横波速度，9、10线代表密度，5、6线以及7、8线分别表示利用纵、横波速度以及密度值计算得到的储层参数和流体参数。
图5为在入射角度分别为5°、15°、25°、35°时基于Zeoppritz公式，利用图4中的理论值与雷克子波褶积得到的未含噪声的角度道集示意图。其中，图5中的(a)图为纵波地震角度道集，图5中的(b)图为转换波地震角度道集。
图6是在图5给出的角度道集基础上增加了噪声干扰后的角度道集的示意图，其中图6中的(a)图为纵波地震角度道集，图6中的(b)图为转换波角度道集，用来验证本说明书实施例提出方法的抗噪性。
图7是基于图6中的(b)图的数据，基于单独纵波数据，分别利用本说明书实施例提出的纵波近似模型以及转换波近似模型(图6中的(a)图)以及间接计算(图6中的(b)图)得到的储层和流体参数结果示意图。其中，实线11、13为反演得到的储层参数曲线，实线12、14为反演得到的流体参数曲线。从图7中可见虽然间接计算可以获得与理论值较为相似的解，但是局部存在较大误差，而利用说明书实施例的方程直接求解得到的结果与理论值更加吻合。
图8是基于图6中的(a)图的数据，分别利用本说明书实施例提出的在纵波和转换波共同约束下的储层参数和流体参数提取方法(图8中的(a)图)以及基于本说明书实施例提出的方程的单独纵波方法(图8中的(b)图)反演示意图。其中，实线15、17为反演的储层参数结果，实线16、18为反演的流体参数结果。通过对比图8中的(a)图和图8中的(b)图可见，在纵波和转换波共同约束下提取的储层和流体参数相对仅在纵波约束下提取的参数精度要高。
为了进一步验证本说明书实施例的抗噪性，对含有噪声的数据(图6中的(a)图)也进行了多波(图9中的(a)图)以及单独纵波(图9中的(b)图)方法直接提取储层参数和流体参数的反演方法对比。其中，实线19、21为反演的储层参数结果，实线20、22为反演的流体参数结果。通过对比可见，在观测数据含有噪声干扰时，本说明书实施例提出的方法相对单独用纵波约束提取的方法更加稳定，抗噪性更强。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。具体的可以参照前述相关处理相关实施例的描述，在此不做一一赘述。
上述对本说明书特定实施例进行了描述。其它实施例在所附权利要求书的范围内。在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描绘的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。
本说明书一个或多个实施例提供的确定地层弹性参数的方法，可以利用直接表征储层参数和流体参数分布特征的反演模型，在分角度叠加数据体的基础上直接计算得到表征储层和流体的弹性参数，不需要间接计算，降低了间接计算引入的累积误差。并且在反演过程中加入了转换波信息，使得反演结果的精度得到较大的提高。
基于上述所述的确定地层弹性参数的方法，本说明书一个或多个实施例还提供一种确定地层弹性参数的装置。所述的装置可以包括使用了本说明书实施例所述方法的***、软件(应用)、模块、组件、服务器等并结合必要的实施硬件的装置。基于同一创新构思，本说明书实施例提供的一个或多个实施例中的装置如下面的实施例所述。由于装置解决问题的实现方案与方法相似，因此本说明书实施例具体的装置的实施可以参见前述方法的实施，重复之处不再赘述。以下所使用的，术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。具体的，图10表示说明书提供的一种确定地层弹性参数的装置实施例的模块结构示意图，如图10所示，所述装置可以包括：
数据获取模块102，可以用于获取目标工区的实测的纵波分角度地震数据以及转换波分角度地震数据；
反演函数获取模块104，可以用于获取用于表征实测地震数据与基于弹性参数计算的理论地震数据之间接近关系的目标函数，所述目标函数根据用于直接表征储层参数和流体参数分布特征的纵波近似模型以及转换波近似模型构建，所述弹性参数包括储层参数和流体参数；
反演模块106，可以用于利用所述目标函数对所述纵波分角度地震数据、转换波分角度地震数据进行反演处理，获得所述目标工区的弹性参数数据。
本说明书的另一个实施例中，所述纵波近似模型可以包括：
Rresevoir＝2(μ2ρ2-μ1ρ1)/(μ2ρ2+μ1ρ1)
Rfluid＝2(L2-L1)/(L2+L2)
其中，r(α)表示纵波数据，α是纵波数据的反射角度和透射角度的平均，Rresevoir为储层参数变化率，μ2为上层介质的剪切模量，μ1为下层介质的剪切模量，ρ2为上层介质的密度参数，ρ1为下层介质的密度参数；Rfluid为流体参数变化率，L1为下层介质的流体参数，L2为下层介质的流体参数，k1为干岩石速度比的二次方，k2为饱和岩石的速度比二次方。
本说明书的另一个实施例中，所述转换波近似模型可以包括：
其中，s(α,β)表示转换波数据，β是转换波数据的反射角度和透射角度的平均。
本说明书的另一个实施例中，所述目标函数可以包括：
其中，r表示根据纵波近似模型确定的纵波数据；s表示根据转换波近似模型确定的转换波数据，W1为纵波分角度叠加数据体的子波组成的矩阵，W2为转换波分角度叠加数据体的子波组成的矩阵；C1为纵波分角度叠加数据体的子波组成的矩阵，C2为转换波分角度叠加数据体的子波组成的矩阵；Ω为相关矩阵；mt为初始的储层参数、流体参数和密度参数的低频分量；G为低频因子，Gm为待求的储层参数、流体参数和密度参数的低频分量；λ1、λ2、λ3为调节权重参数；N为每一个地震道的采样点数。
需要说明的，上述所述的装置根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
本说明书一个或多个实施例提供的确定地层弹性参数的装置，可以利用直接表征储层参数和流体参数分布特征的反演模型，在分角度叠加数据体的基础上直接计算得到表征储层和流体的弹性参数，不需要间接计算，降低了间接计算引入的累积误差。并且在反演过程中加入了转换波信息，使得反演结果的精度得到较大的提高。
本说明书提供的上述实施例所述的方法或装置可以通过计算机程序实现业务逻辑并记录在存储介质上，所述的存储介质可以计算机读取并执行，实现本说明书实施例所描述方案的效果。因此，本说明书还提供一种确定地层弹性参数的设备，包括处理器及存储处理器可执行指令的存储器，所述指令被所述处理器执行时实现包括上述任意一个实施例所述方法的步骤。
所述存储介质可以包括用于存储信息的物理装置，通常是将信息数字化后再以利用电、磁或者光学等方式的媒体加以存储。所述存储介质有可以包括：利用电能方式存储信息的装置如，各式存储器，如RAM、ROM等；利用磁能方式存储信息的装置如，硬盘、软盘、磁带、磁芯存储器、磁泡存储器、U盘；利用光学方式存储信息的装置如，CD或DVD。当然，还有其他方式的可读存储介质，例如量子存储器、石墨烯存储器等等。
需要说明的，上述所述的设备根据方法实施例的描述还可以包括其他的实施方式。具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
上述实施例所述的确定地层弹性参数的设备，可以利用直接表征储层参数和流体参数分布特征的反演模型，在分角度叠加数据体的基础上直接计算得到表征储层和流体的弹性参数，不需要间接计算，降低了间接计算引入的累积误差。并且在反演过程中加入了转换波信息，使得反演结果的精度得到较大的提高。
本说明书还提供一种确定地层弹性参数的***，所述***可以为单独的确定地层弹性参数的***，也可以应用在多种计算机数据处理***中。所述的***可以为单独的服务器，也可以包括使用了本说明书的一个或多个所述方法或一个或多个实施例装置的服务器集群、***(包括分布式***)、软件(应用)、实际操作装置、逻辑门电路装置、量子计算机等并结合必要的实施硬件的终端装置。所述确定地层弹性参数的***可以包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现上述任意一个或者多个实施例中所述方法的步骤。
需要说明的，上述所述的***根据方法或者装置实施例的描述还可以包括其他的实施方式，具体的实现方式可以参照相关方法实施例的描述，在此不作一一赘述。
上述实施例所述的确定地层弹性参数的***，可以利用直接表征储层参数和流体参数分布特征的反演模型，在分角度叠加数据体的基础上直接计算得到表征储层和流体的弹性参数，不需要间接计算，降低了间接计算引入的累积误差。并且在反演过程中加入了转换波信息，使得反演结果的精度得到较大的提高。
本说明书实施例并不局限于必须是符合标准数据模型/模板或本说明书实施例所描述的情况。某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、存储、判断、处理方式等获取的实施例，仍然可以属于本说明书的可选实施方案范围之内。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法或者设备中还存在另外的相同要素。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于***实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本说明书的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述并不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
以上所述仅为本说明书的实施例而已，并不用于限制本说明书。对于本领域技术人员来说，本说明书可以有各种更改和变化。凡在本说明书的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本说明书的权利要求范围之内。
Claims (10)
1.一种确定地层弹性参数的方法，其特征在于，包括：
获取目标工区的实测的纵波分角度地震数据以及转换波分角度地震数据；
获取用于表征实测地震数据与基于弹性参数计算的理论地震数据之间接近关系的目标函数，所述目标函数根据用于直接表征储层参数和流体参数分布特征的纵波近似模型以及转换波近似模型构建，所述弹性参数包括储层参数和流体参数；
利用所述目标函数对所述纵波分角度地震数据、转换波分角度地震数据进行反演处理，获得所述目标工区的弹性参数数据。
2.根据权利要求1所述的方法，其特征在于，所述纵波近似模型包括：
Rresevoir＝2(μ2ρ2-μ1ρ1)/(μ2ρ2+μ1ρ1)
Rfluid＝2(L2-L1)/(L2+L2)
其中，r(α)表示纵波数据，α是纵波数据的反射角度和透射角度的平均，Rresevoir为储层参数变化率，μ2为上层介质的剪切模量，μ1为下层介质的剪切模量，ρ2为上层介质的密度参数，ρ1为下层介质的密度参数；Rfluid为流体参数变化率，L1为下层介质的流体参数，L2为下层介质的流体参数，k1为干岩石速度比的二次方，k2为饱和岩石的速度比二次方。
3.根据权利要求2所述的方法，其特征在于，所述转换波近似模型包括：
其中，s(α,β)表示转换波数据，β是转换波数据的反射角度和透射角度的平均。
4.根据权利要求3所述的方法，其特征在于，所述目标函数包括：
其中，r表示根据纵波近似模型确定的纵波数据；s表示根据转换波近似模型确定的转换波数据，W1为纵波分角度叠加数据体的子波组成的矩阵，W2为转换波分角度叠加数据体的子波组成的矩阵；C1为纵波分角度叠加数据体的子波组成的矩阵，C2为转换波分角度叠加数据体的子波组成的矩阵；Ω为相关矩阵；mt为初始的储层参数、流体参数和密度参数的低频分量；G为低频因子，Gm为待求的储层参数、流体参数和密度参数的低频分量；λ1、λ2、λ3为调节权重参数；N为每一个地震道的采样点数。
5.根据权利要求1所述的方法，其特征在于，所述利用所述目标函数对所述纵波分角度地震数据、转换波分角度地震数据进行反演处理，获得所述目标工区的弹性参数数据，包括：
根据所述目标函数对所述纵波分角度地震数据、转换波分角度地震数据进行反演处理，获得储层参数变化率和流体参数变化率；
对所述储层参数变化率和流体参数变化率进行积分，获得储层参数以及流体参数。
6.一种确定地层弹性参数的装置，其特征在于，包括：
数据获取模块，用于获取目标工区的实测的纵波分角度地震数据以及转换波分角度地震数据；
反演函数获取模块，用于获取用于表征实测地震数据与基于弹性参数计算的理论地震数据之间接近关系的目标函数，所述目标函数根据用于直接表征储层参数和流体参数分布特征的纵波近似模型以及转换波近似模型构建，所述弹性参数包括储层参数和流体参数；
反演模块，用于利用所述目标函数对所述纵波分角度地震数据、转换波分角度地震数据进行反演处理，获得所述目标工区的弹性参数数据。
7.根据权利要求6所述的装置，其特征在于，所述纵波近似模型包括：
Rresevoir＝2(μ2ρ2-μ1ρ1)/(μ2ρ2+μ1ρ1)
Rfluid＝2(L2-L1)/(L2+L2)
其中，r(α)表示纵波数据，α是纵波数据的反射角度和透射角度的平均，Rresevoir为储层参数变化率，μ2为上层介质的剪切模量，μ1为下层介质的剪切模量，ρ2为上层介质的密度参数，ρ1为下层介质的密度参数；Rfluid为流体参数变化率，L1为下层介质的流体参数，L2为下层介质的流体参数，k1为干岩石速度比的二次方，k2为饱和岩石的速度比二次方。
8.根据权利要求7所述的装置，其特征在于，所述转换波近似模型包括：
其中，s(α,β)表示转换波数据，β是转换波数据的反射角度和透射角度的平均。
9.根据权利要求8所述的装置，其特征在于，所述目标函数包括：
其中，r表示根据纵波近似模型确定的纵波数据；s表示根据转换波近似模型确定的转换波数据，W1为纵波分角度叠加数据体的子波组成的矩阵，W2为转换波分角度叠加数据体的子波组成的矩阵；C1为纵波分角度叠加数据体的子波组成的矩阵，C2为转换波分角度叠加数据体的子波组成的矩阵；Ω为相关矩阵；mt为初始的储层参数、流体参数和密度参数的低频分量；G为低频因子，Gm为待求的储层参数、流体参数和密度参数的低频分量；λ1、λ2、λ3为调节权重参数；N为每一个地震道的采样点数。
10.一种确定地层弹性参数的***，其特征在于，所述***包括至少一个处理器以及存储计算机可执行指令的存储器，所述处理器执行所述指令时实现所述权利要求1-5任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911398538.1A CN111077574A (zh) | 2019-12-30 | 2019-12-30 | 一种确定地层弹性参数的方法、装置及*** |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911398538.1A CN111077574A (zh) | 2019-12-30 | 2019-12-30 | 一种确定地层弹性参数的方法、装置及*** |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111077574A true CN111077574A (zh) | 2020-04-28 |
Family
ID=70319846
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911398538.1A Pending CN111077574A (zh) | 2019-12-30 | 2019-12-30 | 一种确定地层弹性参数的方法、装置及*** |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111077574A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111651874A (zh) * | 2020-05-25 | 2020-09-11 | 国网江西省电力有限公司电力科学研究院 | 一种纤维材料吸声性能计算方法 |
| CN112748467A (zh) * | 2020-12-24 | 2021-05-04 | 中国海洋石油集团有限公司 | 利用地震资料确定Russell流体因子的方法及*** |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102692645A (zh) * | 2012-06-01 | 2012-09-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 利用纵波、转换波数据联合反演储层纵横波速度比的方法 |
-
2019
- 2019-12-30 CN CN201911398538.1A patent/CN111077574A/zh active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102692645A (zh) * | 2012-06-01 | 2012-09-26 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 利用纵波、转换波数据联合反演储层纵横波速度比的方法 |
Non-Patent Citations (4)
| Title |
|---|
| 宋永忠 等: ""流体识别因子在扶余油层油气检测中的应用"", 《大庆石油地质与开发》 * |
| 杨文博 等: ""基于叠前地震反演的储层流体识别方法"", 《物探化探计算技术》 * |
| 边会媛 等: ""利用弹性参数识别致密砂岩储层流体性质"", 《地球科学与环境学报》 * |
| 马琦琦 等: ""基于双约束的流体和岩性因子的叠前直接提取方法"", 《石油科学通报》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111651874A (zh) * | 2020-05-25 | 2020-09-11 | 国网江西省电力有限公司电力科学研究院 | 一种纤维材料吸声性能计算方法 |
| CN112748467A (zh) * | 2020-12-24 | 2021-05-04 | 中国海洋石油集团有限公司 | 利用地震资料确定Russell流体因子的方法及*** |
| CN112748467B (zh) * | 2020-12-24 | 2022-10-14 | 中国海洋石油集团有限公司 | 利用地震资料确定Russell流体因子的方法及*** |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7463552B1 (en) | Method for deriving 3D output volumes using filters derived from flat spot direction vectors | |
| Botter et al. | From mechanical modeling to seismic imaging of faults: A synthetic workflow to study the impact of faults on seismic | |
| US7266041B1 (en) | Multi-attribute background relative scanning of 3D geophysical datasets for locally anomaluous data points | |
| EP2915064B1 (en) | Method of indicating the presence of gas hydrate and shallow gas in deepwater environment | |
| US9939506B2 (en) | Methods of investigating formation samples using NMR data | |
| NO346583B1 (no) | Fremgangsmåte for Identifisering av en Kildebergartkandidat | |
| Chehrazi et al. | Pore-facies as a tool for incorporation of small-scale dynamic information in integrated reservoir studies | |
| CN111897011B (zh) | 储层孔隙特征确定方法、装置及设备 | |
| AU2019237361B2 (en) | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic inversions | |
| US10534100B2 (en) | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on time-lapse seismic data | |
| Gholami et al. | Shear wave velocity prediction using seismic attributes and well log data | |
| Faraji et al. | Integration of core data, well logs and seismic attributes for identification of the low reservoir quality units with unswept gas in the carbonate rocks of the world’s largest gas field | |
| CN111077574A (zh) | 一种确定地层弹性参数的方法、装置及*** | |
| CN112114380A (zh) | 一种用于砂砾岩油藏储层的预测方法 | |
| Rotimi et al. | Reservoir characterization and modeling of lateral heterogeneity using multivariate analysis | |
| Kadkhodaie et al. | Acoustic, density, and seismic attribute analysis to aid gas detection and delineation of reservoir properties | |
| US10718876B2 (en) | System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic inversions | |
| CN111077573A (zh) | 一种确定地层弹性参数的方法、装置及*** | |
| Malvić et al. | Reducing variogram uncertainties using the ‘jack-knifing’method, a case study of the Stari Gradac–Barcs-Nyugat field | |
| Refunjol et al. | Fluid discrimination using detrended seismic impedance | |
| Lee et al. | Predicting S-wave sonic logs using machine learning with conventional logs for the Delaware Basin, Texas | |
| Van Riel et al. | Full integration of seismic data into geostatistical reservoir modeling | |
| Bagheripour et al. | Vertical resolution enhancement of petrophysical Nuclear Magnetic Resonance (NMR) log using ordinary kriging | |
| Harris et al. | How Organic Carbon Content and Thermal Maturity Affect Acoustic Properties (and Ultimately Seismic Response) in a Shale Gas/Oil Formation: Woodford Shale, Permian Basin | |
| CN114428293A (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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200428 |