CN112147675B - 一种储层流体检测方法及装置 - Google Patents
一种储层流体检测方法及装置 Download PDFInfo
- Publication number
- CN112147675B CN112147675B CN201910575689.3A CN201910575689A CN112147675B CN 112147675 B CN112147675 B CN 112147675B CN 201910575689 A CN201910575689 A CN 201910575689A CN 112147675 B CN112147675 B CN 112147675B
- Authority
- CN
- China
- Prior art keywords
- frequency
- detection point
- spectrum
- wave data
- target
- 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
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 265
- 239000012530 fluid Substances 0.000 title claims abstract description 88
- 238000001228 spectrum Methods 0.000 claims abstract description 213
- 238000012545 processing Methods 0.000 claims abstract description 33
- 238000004458 analytical method Methods 0.000 claims abstract description 26
- 238000000034 method Methods 0.000 claims abstract description 26
- 238000000926 separation method Methods 0.000 claims abstract description 21
- 238000004590 computer program Methods 0.000 claims description 17
- 230000006870 function Effects 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 11
- 238000012360 testing method Methods 0.000 claims description 4
- 239000002156 adsorbate Substances 0.000 abstract 1
- 238000001179 sorption measurement Methods 0.000 abstract 1
- 239000006185 dispersion Substances 0.000 description 26
- 238000010586 diagram Methods 0.000 description 21
- 230000015572 biosynthetic process Effects 0.000 description 14
- 230000008859 change Effects 0.000 description 14
- 239000011435 rock Substances 0.000 description 12
- 230000007423 decrease Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 238000004891 communication Methods 0.000 description 6
- 238000010183 spectrum analysis Methods 0.000 description 4
- 238000009826 distribution Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000009434 installation Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 240000004371 Panax ginseng Species 0.000 description 1
- 235000002789 Panax ginseng Nutrition 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 238000002310 reflectometry Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/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/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (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
本发明提供一种储层流体检测方法及装置，所述方法包括：对目标储层对应的SVP地震数据进行上下行波分离，获得下行波数据及上行波数据；对下行波数据及上行波数据进行时频分析，获得下行波数据对应的振幅谱和上行波数据对应的振幅谱；根据下行波数据对应的振幅谱参考频率对应的最大振幅值及预设频率范围内各个频率对应的最大振幅值，获得预设频率范围内各个频率对应的谱均衡系数；根据各个频率对应的谱均衡系数对上行波数据在预设频率范围内的振幅谱进行谱均衡处理，获得目标振幅谱；根据目标振幅谱，则确定目标储层内包括流体。所述装置用于执行上述方法。本发明实施例提供的吸附质气体的吸附方式确定方法及装置，提高了储层流体检测的准确性。
Description
技术领域
本发明涉及地质勘探技术领域，尤其涉及一种储层流体检测方法及装置。
背景技术
储层流体识别是油气藏勘探与储层评价的重要环节。随着石油工业的发展，油气勘探开发的目标不再是以构造为主的油气藏，而是更为隐蔽的以岩性圈闭等类型为主的复杂油气藏，研究内容从构造、储层发展到储层流体检测。
近年来实验室模拟、理论研究及实际地震资料均表明含流体填充介质常常导致地震波发生不同程度的频散和衰减，所以利用含烃储层导致的地震波频散，引起反射系数随频率变化这一现象可以进行油气检测。现有技术中，从地震数据中提取与频率有关的地震振幅属性都是基于常规地面地震的频谱分析技术，例如Partyka等人提出了利用短时离散傅里叶变换对地震数据进行频谱分析来描述储层的特征，开创了利用频谱分析技术进行储层描述的先河，但频谱分析的结果通常包含“子波叠印”的影响，使得原有的反射能量受子波的影响在主频位置能量最大，而在主频位置两侧能量减小，各频率成分能量分布不均衡，所以通常要进行谱均衡来消除这种子波效应。对于常规地面地震，谱均衡过程中需要选取弹性反射层作为参考界面对其它层进行处理，而在实际地震资料中却很难找到一个合适的弹性反射层，通常选取的是一个上覆反射性较强的连续性界面作为参考界面，在处理过程中会引起较大的误差，分析结果受到地层构造和弹性特征的影响很大，从而影响地震资料频散特征的提取，导致对储层流体的检测不准确。
因此，如何提出一种储层流体检测方法，避免假设弹性层对储层流体检测的影响，以提高储层流体检测的准确性成为本领域需要解决的重要课题。
发明内容
针对现有技术中的问题，本发明实施例提供一种储层流体检测方法及装置，能够提高储层流体检测的准确性。
一方面，本发明提出一种储层流体检测方法，包括：
获取目标储层对应的垂直地震剖面地震数据；
对所述垂直地震剖面地震数据进行上下行波分离，获得所述目标储层对应的参考检测点的下行波数据以及所述目标储层对应的反射检测点的上行波数据；其中，所述参考检测点和所述反射检测点是预先设置的；
对所述参考检测点的下行波数据以及所述反射检测点的上行波数据进行时频分析，获得所述参考检测点的下行波数据对应的振幅谱和所述反射检测点的上行波数据对应的振幅谱；其中，述振幅谱为采样时间和频率的函数；
根据所述参考检测点的下行波数据对应的振幅谱中参考频率对应的最大振幅值以及预设频率范围内各个频率对应的最大振幅值，分别获得所述预设频率范围内各个所述频率对应的谱均衡系数；其中，所述参考频率是所述垂直地震剖面地震数据对应的主频；
根据所述预设频率范围内各个所述频率对应的谱均衡系数对所述反射检测点的上行波数据在所述预设频率范围内的振幅谱进行谱均衡处理，获得目标振幅谱；
若判断获知所述目标振幅谱中在所述预设频率范围内的各个所述频率对应的振幅不稳定，则确定所述目标储层内包括流体。
另一方面，本发明提供一种储层流体检测装置，包括：
获取单元，用于获取目标储层对应的垂直地震剖面地震数据；
分离单元，用于对所述垂直地震剖面地震数据进行上下行波分离，获得所述目标储层对应的参考检测点的下行波数据以及所述目标储层对应的反射检测点的上行波数据；其中，所述参考检测点和所述反射检测点是预先设置的；
时频分析单元，用于对所述参考检测点的下行波数据以及所述反射检测点的上行波数据进行时频分析，获得所述参考检测点的下行波数据对应的振幅谱和所述反射检测点的上行波数据对应的振幅谱；所述振幅谱为采样时间和频率的函数；
获得单元，用于根据所述参考检测点的下行波数据对应的振幅谱中参考频率对应的最大振幅值以及预设频率范围内各个频率对应的最大振幅值，分别获得所述预设频率范围内各个所述频率对应的谱均衡系数；其中，所述参考频率是所述垂直地震剖面地震数据对应的主频；
处理单元，用于根据所述预设频率范围内各个所述频率对应的谱均衡系数对所述反射检测点的上行波数据在所述预设频率范围内的振幅谱进行谱均衡处理，获得目标振幅谱；
确定单元，用于在判断获知所述目标振幅谱中在所述预设频率范围内的各个所述频率对应的振幅不稳定之后，则确定目标储层内包括流体。
再一方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现上述任一实施例所述储层流体检测方法的步骤。
又一方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现上述任一实施例所述储层流体检测方法的步骤。
本发明实施例提供的储层流体检测方法及装置，能够对目标储层对应的VSP地震数据进行上下行波分离，获得目标储层对应的参考检测点的下行波数据以及所述目标储层对应的反射检测点的上行波数据，并对参考检测点的下行波数据以及反射检测点的上行波数据进行时频分析，获得参考检测点的下行波数据对应的振幅谱和反射检测点的上行波数据对应的振幅谱，基于参考检测点的下行波数据对应的振幅谱对反射检测点的上行波数据在预设频率范围内的振幅谱进行谱均衡处理，获得目标振幅谱，在判断获知目标振幅谱中在预设频率范围内的各个所述频率对应的振幅不稳定之后，确定所述目标储层内包括流体，由于能够使上行波数据只体现与频率有关的反射系数的变化情况，准确表征地震资料的频散特征，提高了储层流体检测的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1是本发明一实施例提供的一种储层流体检测方法的流程示意图。
图2是本发明一实施例提供的储层地质模型的结构示意图。
图3是本发明一实施例提供的VSP地震记录的示意图。
图4a是本发明一实施例提供的目标振幅谱的示意图。
图4b是本发明一实施例提供的目标振幅谱的振幅变化率的示意图。
图5a是本发明另一实施例提供的目标振幅谱的示意图。
图5b是本发明另一实施例提供的目标振幅谱的振幅变化率的示意图。
图6是本发明一实施例提供的储层流体检测装置的结构示意图。
图7是本发明一实施例提供的储层流体检测装置的结构示意图。
图8为本发明一实施例提供的电子设备的实体结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
图1是本发明一实施例提供的一种储层流体检测方法的流程示意图，如图1所示，本发明实施例提供的储层流体检测方法，包括：
S101、获取目标储层对应的垂直地震剖面地震数据；
具体地，采集包括目标储层在内的区域的垂直地震剖面(Vertical SeismicProfiling，简称VSP)地震数据，即在所述区域的地表设置震源激发地震波，预先设置在井内的各个检波器接收地震波，在垂直方向观测一维人工场，然后对所观测得到的资料经过校正、叠加、滤波等处理，可以得到VSP地震数据，得到的VSP地震数据与所述目标储层对应。储层流体检测装置可以获取所述目标储层对应的VSP地震数据。其中，所述目标储层是预先确定的；在本发明实施例中，各个检波器的安装位置称为检测点。
S102、对所述垂直地震剖面地震数据进行上下行波分离，获得所述目标储层对应的参考检测点的下行波数据以及所述目标储层对应的反射检测点的上行波数据；其中，所述参考检测点和所述反射检测点是预先设置的垂直地震剖面地震数据；
具体地，所述VSP地震数据包括上行波数据和下行波数据，所述储层流体检测装置可以对所述VSP地震数据进行上下行波分离，获得每个所述检测点的上行波数据和/或下行波数据，所述检测点的上行波数据为所述检测点对应的检波器接收到的上行波数据，所述检测点的下行波数据为所述检测点对应的检波器接收到的下行波数据。在各个所述检测点中，预先设置所述目标储层的参考检测点和反射检测点，所述储层流体检测装置可以获得所述参考检测点的下行波数据和所述反射检测点的上行波数据。其中，所述VSP地震数据的上下行波的分离，可以通过专业软件实现，专业软件根据实际情况进行选择，本发明实施例不做限定。
例如，所述参考检测点为各个检测点中最靠近所述目标储层的上分界面的检测点，或者为各个所述检测点中最靠近所述目标储层的下分界面的检测点。设置在所述反射检测点的检波器能够接收到所述目标储层的分界面的反射波，即如果所述参考检测点靠近所述目标储层的上分界面，设置在所述反射检测点的检波器能够接收到所述上分界面的发射波，如果所述参考检测点靠近所述目标储层的下分界面，设置在所述反射检测点的检波器能够接收到所述下分界面的发射波。其中，在本发明实施例中所有的上行波都是反射波。
S103、对所述参考检测点的下行波数据以及所述反射检测点的上行波数据进行时频分析，获得所述参考检测点的下行波数据对应的振幅谱和所述反射检测点的上行波数据对应的振幅谱；所述振幅谱为采样时间和频率的函数；
具体地，所述储层流体检测装置在获得所述参考检测点的下行波数据和所述反射检测点的上行波数据之后，可以对所述参考检测点的下行波数据进行时频分析，得到所述参考检测点的下行波数据对应的振幅谱，并对所述反射参考点的上行波数据进行时频分析，得到所述反射参考点的上行波数据对应的振幅谱。其中，所述振幅谱是采样时间和频率的函数。其中，所述储层流体检测装置可以采用维格纳(Wigner-Ville)分布等时频分析方法进行时频分析。
S104、根据所述参考检测点的下行波数据对应的振幅谱中参考频率对应的最大振幅值以及预设频率范围内各个频率对应的最大振幅值，分别获得所述预设频率范围内各个所述频率对应的谱均衡系数；其中，所述参考频率是所述垂直地震剖面地震数据对应的主频；
具体地，所述VSP地震数据对应的地震子波存在主频，所述储层流体检测装置对所述VSP地震数据进行处理，可以获得所述VSP地震数据对应的主频，并将所述VSP地震数据对应的主频作为参考频率。所述储层流体检测装置从所述参考检测点的下行波数据对应的振幅谱中可以获得所述参考频率对应的最大振幅值，所述储层流体检测装置还可以从所述反射检测点的下行波数据对应的振幅谱中获得预设频率范围内每个频率对应的最大振幅值。所述储层流体检测装置根据所述参考频率对应的最大振幅值和所述预设频率范围内每个所述频率对应的最大振幅值，可以获得所述预设频率范围内每个所述频率对应的谱均衡系数。所述预设频率范围根据实际需要进行设置，本发明实施例不做限定。由于利用参考检测点的下行波数据获得谱均衡系数，避免了弹性参考界面的选择，能够减小后续谱均衡处理中由于参考界面的选取造成的误差。
例如，所述储层流体检测装置获得所述参考频率f参对应的最大振幅值为a，所述预设频率范围内的一个频率fi对应的最大振幅值为b，可以得到频率fi对应的谱均衡系数ω＝a/b。
S105、根据所述预设频率范围内各个所述频率对应的谱均衡系数对所述反射检测点的上行波数据在所述预设频率范围内的振幅谱进行谱均衡处理，获得目标振幅谱；
具体地，所述储层流体检测装置在获得所述预设频率范围内的各个所述频率对应的谱均衡系数之后，利用各个所述频率对应的谱均衡系数对所述反射点的上行波数据在所述预设频率范围内的振幅谱进行谱均衡处理，可以获得目标振幅谱。谱均衡处理的目的是消除“子波叠印”对所述反射检测点的上行波数据对应的振幅谱的影响，使目标振幅谱能够真实反映目标储层的频散特征。
S106、若判断获知所述目标振幅谱中在所述预设频率范围内的各个所述频率对应的振幅不稳定，则确定所述目标储层内包括流体。
具体地，所述储层流体检测装置在获得所述目标振幅谱之后，可以判断所述目标振幅谱在所述预设频率范围内的振幅变化情况，如果所述目标储层不包括流体，那么所述目标储层的下界面反射的上行波，发生的频散和衰减较小，最终得到的所述目标振幅谱在所述预设频率范围内的振幅基本保持不变，如果所述目标储层包括流体，会导致所述目标储层的下界面反射的上行波，发生不同程度的频散和衰减，最终得到的所述目标振幅谱在所述预设频率范围内会表现出明显的频散特性，不同频率对应的最大振幅值能够清晰的反映出上述频散特性，例如，所述目标振幅谱在所述预设频率范围内的振幅值会随着频率的增大而逐渐减小。在本发明实施例中，所述储层流体检测装置可以判断所述目标振幅谱在所述预设频率范围内的振幅变化是否在预设区间内，如果所述振幅变化在所述预设区间内，那么所述目标振幅谱中在所述预设频率范围内的各个所述频率对应的振幅稳定，所述目标储层内不包括流体；如果所述振幅变化不在所述预设区间内，那么所述目标振幅谱中在所述预设频率范围内的各个所述频率对应的振幅不稳定，所述目标储层内包括流体。
本发明实施例提供的储层流体检测方法，能够对目标储层对应的VSP地震数据进行上下行波分离，获得目标储层对应的参考检测点的下行波数据以及所述目标储层对应的反射检测点的上行波数据，并对参考检测点的下行波数据以及反射检测点的上行波数据进行时频分析，获得参考检测点的下行波数据对应的振幅谱和反射检测点的上行波数据对应的振幅谱，基于参考检测点的下行波数据对应的振幅谱对反射检测点的上行波数据在预设频率范围内的振幅谱进行谱均衡处理，获得目标振幅谱，在判断获知目标振幅谱中在预设频率范围内的各个所述频率对应的振幅不稳定之后，确定所述目标储层内包括流体，由于能够使上行波数据只体现与频率有关的反射系数的变化情况，准确表征地震资料的频散特征，提高了储层流体检测的准确性。
在上述各实施例的基础上，进一步地，所述参考检测点为各个检测点中最靠近所述目标储层的上分界面的检测点，或者为各个所述检测点中最靠近所述目标储层的下分界面的检测点。
具体地，为了采集VSP地震数据，会在井内设置多个检波器，为了使获得的所述参考检测点的下行波数据受到的噪声干扰较少，以及在谱均衡时能够消除“子波叠印”的影响，将所述多个检波器中最靠近所述目标储层的上分界面检波器对应的检测点，作为所述参考检测点，或者将所述多个检波器中最靠近所述目标储层的下分界面检波器对应的检测点，作为所述参考检测点。可理解的是，所述目标储层的上分界面是位于所述目标储层相邻的上方地层与所述目标储层的分界面，所述目标储层的下分界面是位于所述目标储层相邻的下方地层与所述目标储层的分界面。
在上述各实施例的基础上，进一步地，所述根据所述参考检测点的下行波数据对应的振幅谱中参考频率对应的最大振幅值以及预设频率范围内各个频率对应的最大振幅值，分别获得所述预设频率范围内各个所述频率对应的谱均衡系数振幅谱包括：
将所述参考检测点的下行波数据对应的振幅谱中所述参考频率对应的最大振幅值除以所述预设频率范围内每个所述频率对应的最大振幅值的结果，作为每个所述频率对应的谱均衡系数。
具体地，所述储层流体检测装置在计算所述预设频率范围内每个所述频率对应的谱均衡系数时，可以用所述参考检测点的下行波数据对应的振幅谱中所述参考频率对应的最大振幅值除以所述预设频率范围内每个所述频率对应的最大振幅值。
例如，所述预设频率范围内的第i个频率为fi，fi对应的最大振幅值为Ai，所述参考频率对应的最大振幅值为B，那么fi对应的谱均衡系数ωi＝B/Ai。
在上述各实施例的基础上，进一步地，所述若判断获知所述目标振幅谱中在所述预设频率范围内的各个所述频率对应的振幅不稳定，则确定所述目标储层内包括流体包括：
若判断获知所述预设频率范围内的各个所述频率对应的振幅值随着所述频率的增大而减小，则确定所述目标储层内包括流体。
具体地，所述储层流体检测装置在获得所述目标振幅谱之后，在所述预设频率范围内对比所有相邻两个所述频率的振幅值，如果对比的结果显示，随着所述频率的增大，对应的振幅值减小，说明所述振幅随频率发生了变化，这种变化是由于目标储层内包括的流体导致的，可以确定所述目标储层内包括流体。
在上述各实施例的基础上，进一步地，本发明实施例提供的储层流体检测方法还包括：
将所述目标振幅谱以图形显示。
具体地，所述储层流体检测装置在获得所述目标振幅图之后，可以将所述目标振幅谱以图形的方式显示出来，将所述目标振幅谱以图形的方式显示，便于相关人员查看。
为了验证本发明实施例提供的储层流体检测方法，可以通过正演模拟的方法建立了一个常规的储层地质模型对本发明实施例提供的储层流体检测方法进行测试。
图2是本发明一实施例提供的储层地质模型的结构示意图，如图2所示，由上到下依次设置上围岩层、含水层、含油层和下围岩层，各层的厚度相等，为100m，由上到下均匀设置16个检波器，相邻两个检波器之间的间距为25m。设置上围岩层和下围岩层的地震波传播速度为3700m/s，上围岩层和下围岩层的密度为2.3g/cm3；设置含水层的地震波传播速度为2644m/s，密度为2.32g/cm3；上围岩层、下围岩层和含水层不设置频散和衰减，含油层设置频散和衰减，设置含油层的速度频散范围为2900m/s～3135m/s。所述储层地质模型的参数设置以及含油层的频散衰减特征都是基于实际测井资料及岩石物理实验得到的，能够反映实际储层的情况。
利用所述储层地质模型模拟VSP地震数据采集时，使用的地震子波为雷克子波，所述雷克子波的主频为40Hz，采样时间间隔为1ms。图3是本发明一实施例提供的VSP地震记录的示意图，如图3所示，上行波1为上围岩层与含水层的分界面的反射波，上行波2为含水层与含油层的分界面的反射波，上行波3为含油层与下围岩层的分界面的反射波。
以含油层为目标储层，设置第8个检波器对应的检测点为参考检测点，并设置第5个检波器对应的检测点为反射检测点。采用本发明实施例提供的储层流体检测方法，基于第8个检波器对应的参考检测点的下行波数据对第5个检波器对应的反射检测点的上行波数据进行频谱均衡处理之后，得到的目标振幅谱如图4a所示。在图4a中，上分界面是指含水层与含油层的分界面，下分界面是指含油层与下围岩的分界面，第5个检波器可以接收到上行波2和上行波3，上行波2经过频谱均衡处理之后可以得到上行波2对应的目标振幅谱，在图4a中对应于上分界面的振幅，第5个检波器接收到的上行波3经过频谱均衡处理之后可以得到上行波3对应的目标振幅谱，在图4a中对应于下分界面的振幅。上行波2没有经过含油层，从图4a可以看出，在24～60Hz频率范围内，对应的振幅值变化不大，基本保持稳定；上行波3经过了含油层，从图4a可以看出，在24～60Hz频率范围内，对应的振幅值随着频率的增大而减小，这种反射波振幅随频率的变化特征是由于含油层导致的。对图4a中振幅计算其随频率的变化率，得到图4b显示的结果，上行波的振幅变化率可以反映频散程度的大小，振幅变化率的正负表示反射系数随频率增大或减小的趋势，绝对值的大小表示频散程度的大小，可以看出经过含油层的下分界面反射的反射波的频散程度明显要大于含油层的上分界面反射的反射波的频散程度。
以含油层为目标储层，设置第12个检波器对应的检测点为参考检测点，并设置第9个检波器对应的检测点为反射检测点。采用本发明实施例提供的储层流体检测方法，基于第12个检波器对应的参考检测点的下行波数据对第9个检波器对应的反射检测点的上行波数据进行频谱均衡处理之后，得到的目标振幅谱如图5a所示。在图5a中，下分界面是指含油层与下围岩层的分界面，第9个检波器可以接收到上行波3，第9个检波器接收到的上行波3经过频谱均衡处理之后可以得到对应的目标振幅谱，在图5a中对应于下分界面的振幅。上行波3经过了含油层，从图5a可以看出，在20～60Hz频率范围内，对应的振幅值随着频率的增大而减小，这种反射波振幅随频率的变化特征是由于含油层导致的。对图5a中振幅计算其随频率的变化率，得到图5b显示的结果，上行波的振幅变化率可以反映频散程度的大小，振幅变化率的正负表示反射系数随频率增大或减小的趋势，绝对值的大小表示频散程度的大小。
图6是本发明一实施例提供的储层流体检测装置的结构示意图，如图6所示，本发明实施例提供的储层流体检测装置包括获取单元601、分离单元602、时频分析单元603、获得单元604、处理单元605和确定单元606，其中：
获取单元601用于获取目标储层对应的垂直地震剖面地震数据；分离单元602用于对所述垂直地震剖面地震数据进行上下行波分离，获得所述目标储层对应的参考检测点的下行波数据以及所述目标储层对应的反射检测点的上行波数据；其中，所述参考检测点和所述反射检测点是预先设置的；时频分析单元603用于对所述参考检测点的下行波数据以及所述反射检测点的上行波数据进行时频分析，获得所述参考检测点的下行波数据对应的振幅谱和所述反射检测点的上行波数据对应的振幅谱；所述振幅谱为采样时间和频率的函数；获得单元604用于根据所述参考检测点的下行波数据对应的振幅谱中参考频率对应的最大振幅值以及预设频率范围内各个频率对应的最大振幅值，分别获得所述预设频率范围内各个所述频率对应的谱均衡系数；其中，所述参考频率是所述垂直地震剖面地震数据对应的主频；处理单元605用于根据所述预设频率范围内各个所述频率对应的谱均衡系数对所述反射检测点的上行波数据在所述预设频率范围内的振幅谱进行谱均衡处理，获得目标振幅谱；确定单元606用于在判断获知所述目标振幅谱中在所述预设频率范围内的各个所述频率对应的振幅不稳定之后，则确定目标储层内包括流体。
具体地，采集包括目标储层在内的区域的VSP地震数据，即在所述区域的地表设置震源激发地震波，预先设置在井内的各个检波器接收地震波，在垂直方向观测一维人工场，然后对所观测得到的资料经过校正、叠加、滤波等处理，可以得到VSP地震数据，得到的VSP地震数据与所述目标储层对应。获取单元601可以获取所述目标储层对应的VSP地震数据。其中，所述目标储层是预先确定的；在本发明实施例中，各个检波器的安装位置称为检测点。
所述VSP地震数据包括上行波数据和下行波数据，分离单元602可以对所述VSP地震数据进行上下行波分离，获得每个所述检测点的上行波数据和/或下行波数据，所述检测点的上行波数据为所述检测点对应的检波器接收到的上行波数据，所述检测点的下行波数据为所述检测点对应的检波器接收到的下行波数据。在各个所述检测点中，预先设置所述目标储层的参考检测点和反射检测点，所述储层流体检测装置可以获得所述参考检测点的下行波数据和所述反射检测点的上行波数据。其中，所述VSP地震数据的上下行波的分离，可以通过专业软件实现，专业软件根据实际情况进行选择，本发明实施例不做限定。
在获得所述参考检测点的下行波数据和所述反射检测点的上行波数据之后，时频分析单元603可以对所述参考检测点的下行波数据进行时频分析，得到所述参考检测点的下行波数据对应的振幅谱，并对所述反射参考点的上行波数据进行时频分析，得到所述反射参考点的上行波数据对应的振幅谱。其中，所述振幅谱是采样时间和频率的函数。其中，所述储层流体检测装置可以采用维格纳(Wigner-Ville)分布等时频分析方法进行时频分析。
所述VSP地震数据对应的地震子波存在主频，获得单元604对所述VSP地震数据进行处理，可以获得所述VSP地震数据对应的主频，并将所述VSP地震数据对应的主频作为参考频率。获得单元604从所述参考检测点的下行波数据对应的振幅谱中可以获得所述参考频率对应的最大振幅值，获得单元604还可以从所述反射检测点的下行波数据对应的振幅谱中获得预设频率范围内每个频率对应的最大振幅值。获得单元604根据所述参考频率对应的最大振幅值和所述预设频率范围内每个所述频率对应的最大振幅值，可以获得所述预设频率范围内每个所述频率对应的谱均衡系数。所述预设频率范围根据实际需要进行设置，本发明实施例不做限定。由于利用参考检测点的下行波数据获得谱均衡系数，避免了弹性参考界面的选择，能够减小后续谱均衡处理中由于参考界面的选取造成的误差。
在获得所述预设频率范围内的各个所述频率对应的谱均衡系数之后，处理单元605利用各个所述频率对应的谱均衡系数对所述反射点的上行波数据在所述预设频率范围内的振幅谱进行谱均衡处理，可以获得目标振幅谱。谱均衡处理的目的是消除“子波叠印”对所述反射检测点的上行波数据对应的振幅谱的影响，使目标振幅谱能够真实反映目标储层的频散特征。
在获得所述目标振幅谱之后，确定单元606可以判断所述目标振幅谱在所述预设频率范围内的振幅变化情况，如果所述目标储层不包括流体，那么所述目标储层的下界面反射的上行波，发生的频散和衰减较小，最终得到的所述目标振幅谱在所述预设频率范围内的振幅基本保持不变，如果所述目标储层包括流体，会导致所述目标储层的下界面反射的上行波，发生不同程度的频散和衰减，最终得到的所述目标振幅谱在所述预设频率范围内会表现出明显的频散特性，不同频率对应的最大振幅值能够清晰的反映出上述频散特性，例如，所述目标振幅谱在所述预设频率范围内的振幅值会随着频率的增大而逐渐减小。在本发明实施例中，确定单元606可以判断所述目标振幅谱在所述预设频率范围内的振幅变化是否在预设区间内，如果所述振幅变化在所述预设区间内，那么所述目标振幅谱中在所述预设频率范围内的各个所述频率对应的振幅稳定，所述目标储层内不包括流体；如果所述振幅变化不在所述预设区间内，那么所述目标振幅谱中在所述预设频率范围内的各个所述频率对应的振幅不稳定，所述目标储层内包括流体。
本发明实施例提供的储层流体检测装置，能够对目标储层对应的VSP地震数据进行上下行波分离，获得目标储层对应的参考检测点的下行波数据以及所述目标储层对应的反射检测点的上行波数据，并对参考检测点的下行波数据以及反射检测点的上行波数据进行时频分析，获得参考检测点的下行波数据对应的振幅谱和反射检测点的上行波数据对应的振幅谱，基于参考检测点的下行波数据对应的振幅谱对反射检测点的上行波数据在预设频率范围内的振幅谱进行谱均衡处理，获得目标振幅谱，在判断获知目标振幅谱中在预设频率范围内的各个所述频率对应的振幅不稳定之后，确定所述目标储层内包括流体，由于能够使上行波数据只体现与频率有关的反射系数的变化情况，准确表征地震资料的频散特征，提高了储层流体检测的准确性。
在上述各实施例的基础上，进一步地，所述参考检测点为各个检测点中最靠近所述目标储层的上分界面的检测点，或者为各个所述检测点中最靠近所述目标储层的下分界面的检测点。
具体地，为了采集VSP地震数据，会在井内设置多个检波器，为了使获得的所述参考检测点的下行波数据受到的噪声干扰较少，以及在谱均衡时能够消除“子波叠印”的影响，将所述多个检波器中最靠近所述目标储层的上分界面检波器对应的检测点，作为所述参考检测点，或者将所述多个检波器中最靠近所述目标储层的下分界面检波器对应的检测点，作为所述参考检测点。可理解的是，所述目标储层的上分界面是位于所述目标储层相邻的上方地层与所述目标储层的分界面，所述目标储层的下分界面是位于所述目标储层相邻的下方地层与所述目标储层的分界面。
在上述各实施例的基础上，进一步地，获得单元604具体用于：
将所述参考检测点的下行波数据对应的振幅谱中所述参考频率对应的最大振幅值除以所述预设频率范围内每个所述频率对应的最大振幅值的结果，作为每个所述频率对应的谱均衡系数。
具体地，获得单元604在计算所述预设频率范围内每个所述频率对应的谱均衡系数时，可以用所述参考检测点的下行波数据对应的振幅谱中所述参考频率对应的最大振幅值除以所述预设频率范围内每个所述频率对应的最大振幅值。
在上述各实施例的基础上，进一步地，确定单元606具体用于：
若判断获知所述预设频率范围内的各个所述频率对应的振幅值随着所述频率的增大而减小，则确定所述目标储层内包括流体。
具体地，确定单元606在获得所述目标振幅谱之后，在所述预设频率范围内对比所有相邻两个所述频率的振幅值，如果对比的结果显示，随着所述频率的增大，对应的振幅值减小，说明所述振幅随频率发生了变化，这种变化是由于目标储层内包括的流体导致的，可以确定所述目标储层内包括流体。
图7是本发明一实施例提供的储层流体检测装置的结构示意图，如图7所示，在上述各实施例的基础上，进一步地，本发明实施例提供的储层流体检测装置还包括显示单元607，其中：
所述显示单元用于将所述目标振幅谱以图形显示。
具体地，显示单元607在获得所述目标振幅图之后，可以将所述目标振幅谱以图形的方式显示出来，将所述目标振幅谱以图形的方式显示，便于相关人员查看。
图8为本发明一实施例提供的电子设备的实体结构示意图，如图8所示，该电子设备可以包括：处理器(processor)801、通信接口(Communications Interface)802、存储器(memory)803和通信总线804，其中，处理器801，通信接口802，存储器803通过通信总线804完成相互间的通信。处理器801可以调用存储器803中的逻辑指令，以执行如下方法：获取目标储层对应的垂直地震剖面地震数据；对所述垂直地震剖面地震数据进行上下行波分离，获得所述目标储层对应的参考检测点的下行波数据以及所述目标储层对应的反射检测点的上行波数据；其中，所述参考检测点和所述反射检测点是预先设置的；对所述参考检测点的下行波数据以及所述反射检测点的上行波数据进行时频分析，获得所述参考检测点的下行波数据对应的振幅谱和所述反射检测点的上行波数据对应的振幅谱；所述振幅谱为采样时间和频率的函数；根据所述参考检测点的下行波数据对应的振幅谱中参考频率对应的最大振幅值以及预设频率范围内各个频率对应的最大振幅值，分别获得所述预设频率范围内各个所述频率对应的谱均衡系数；其中，所述参考频率是所述垂直地震剖面地震数据对应的主频；根据所述预设频率范围内各个所述频率对应的谱均衡系数对所述反射检测点的上行波数据在所述预设频率范围内的振幅谱进行谱均衡处理，获得目标振幅谱；若判断获知所述目标振幅谱中在所述预设频率范围内的各个所述频率对应的振幅不稳定，则确定所述目标储层内包括流体。
此外，上述的存储器803中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
本实施例公开一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，计算机能够执行上述各方法实施例所提供的方法，例如包括：获取目标储层对应的垂直地震剖面地震数据；对所述垂直地震剖面地震数据进行上下行波分离，获得所述目标储层对应的参考检测点的下行波数据以及所述目标储层对应的反射检测点的上行波数据；其中，所述参考检测点和所述反射检测点是预先设置的；对所述参考检测点的下行波数据以及所述反射检测点的上行波数据进行时频分析，获得所述参考检测点的下行波数据对应的振幅谱和所述反射检测点的上行波数据对应的振幅谱；所述振幅谱为采样时间和频率的函数；根据所述参考检测点的下行波数据对应的振幅谱中参考频率对应的最大振幅值以及预设频率范围内各个频率对应的最大振幅值，分别获得所述预设频率范围内各个所述频率对应的谱均衡系数；其中，所述参考频率是所述垂直地震剖面地震数据对应的主频；根据所述预设频率范围内各个所述频率对应的谱均衡系数对所述反射检测点的上行波数据在所述预设频率范围内的振幅谱进行谱均衡处理，获得目标振幅谱；若判断获知所述目标振幅谱中在所述预设频率范围内的各个所述频率对应的振幅不稳定，则确定所述目标储层内包括流体。
本实施例提供一种计算机可读存储介质，所述计算机可读存储介质存储计算机程序，所述计算机程序使所述计算机执行上述各方法实施例所提供的方法，例如包括：获取目标储层对应的垂直地震剖面地震数据；对所述垂直地震剖面地震数据进行上下行波分离，获得所述目标储层对应的参考检测点的下行波数据以及所述目标储层对应的反射检测点的上行波数据；其中，所述参考检测点和所述反射检测点是预先设置的；对所述参考检测点的下行波数据以及所述反射检测点的上行波数据进行时频分析，获得所述参考检测点的下行波数据对应的振幅谱和所述反射检测点的上行波数据对应的振幅谱；所述振幅谱为采样时间和频率的函数；根据所述参考检测点的下行波数据对应的振幅谱中参考频率对应的最大振幅值以及预设频率范围内各个频率对应的最大振幅值，分别获得所述预设频率范围内各个所述频率对应的谱均衡系数；其中，所述参考频率是所述垂直地震剖面地震数据对应的主频；根据所述预设频率范围内各个所述频率对应的谱均衡系数对所述反射检测点的上行波数据在所述预设频率范围内的振幅谱进行谱均衡处理，获得目标振幅谱；若判断获知所述目标振幅谱中在所述预设频率范围内的各个所述频率对应的振幅不稳定，则确定所述目标储层内包括流体。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在本说明书的描述中，参考术语“一个实施例”、“一个具体实施例”、“一些实施例”、“例如”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种储层流体检测方法，其特征在于，包括：
获取目标储层对应的垂直地震剖面地震数据；
对所述垂直地震剖面地震数据进行上下行波分离，获得所述目标储层对应的参考检测点的下行波数据以及所述目标储层对应的反射检测点的上行波数据；其中，所述参考检测点和所述反射检测点是预先设置的；
对所述参考检测点的下行波数据以及所述反射检测点的上行波数据进行时频分析，获得所述参考检测点的下行波数据对应的振幅谱和所述反射检测点的上行波数据对应的振幅谱；其中，述振幅谱为采样时间和频率的函数；
根据所述参考检测点的下行波数据对应的振幅谱中参考频率对应的最大振幅值以及预设频率范围内各个频率对应的最大振幅值，分别获得所述预设频率范围内所述各个频率对应的谱均衡系数；其中，所述参考频率是所述垂直地震剖面地震数据对应的主频；
根据所述预设频率范围内所述各个频率对应的谱均衡系数对所述反射检测点的上行波数据在所述预设频率范围内的振幅谱进行谱均衡处理，获得目标振幅谱；
若判断获知所述预设频率范围内的所述各个频率对应的振幅值随着频率的增大而减小，则确定所述目标储层内包括流体。
2.根据权利要求1所述的方法，其特征在于，所述参考检测点为各个检测点中最靠近所述目标储层的上分界面的检测点，或者为各个所述检测点中最靠近所述目标储层的下分界面的检测点。
3.根据权利要求1所述的方法，其特征在于，所述根据所述参考检测点的下行波数据对应的振幅谱中参考频率对应的最大振幅值以及预设频率范围内各个频率对应的最大振幅值，分别获得所述预设频率范围内所述各个频率对应的谱均衡系数振幅谱包括：
将所述参考检测点的下行波数据对应的振幅谱中所述参考频率对应的最大振幅值除以所述预设频率范围内每个频率对应的最大振幅值的结果，作为所述预设频率范围内每个频率对应的谱均衡系数。
4.根据权利要求1至3任一项所述的方法，其特征在于，还包括：
将所述目标振幅谱以图形显示。
5.一种储层流体检测装置，其特征在于，包括：
获取单元，用于获取目标储层对应的垂直地震剖面地震数据；
分离单元，用于对所述垂直地震剖面地震数据进行上下行波分离，获得所述目标储层对应的参考检测点的下行波数据以及所述目标储层对应的反射检测点的上行波数据；其中，所述参考检测点和所述反射检测点是预先设置的；
时频分析单元，用于对所述参考检测点的下行波数据以及所述反射检测点的上行波数据进行时频分析，获得所述参考检测点的下行波数据对应的振幅谱和所述反射检测点的上行波数据对应的振幅谱；所述振幅谱为采样时间和频率的函数；
获得单元，用于根据所述参考检测点的下行波数据对应的振幅谱中参考频率对应的最大振幅值以及预设频率范围内各个频率对应的最大振幅值，分别获得所述预设频率范围内所述各个频率对应的谱均衡系数；其中，所述参考频率是所述垂直地震剖面地震数据对应的主频；
处理单元，用于根据所述预设频率范围内所述各个频率对应的谱均衡系数对所述反射检测点的上行波数据在所述预设频率范围内的振幅谱进行谱均衡处理，获得目标振幅谱；
确定单元，用于在判断获知所述预设频率范围内的所述各个频率对应的振幅值随着频率的增大而减小之后，确定所述目标储层内包括流体。
6.根据权利要求5所述的装置，其特征在于，所述参考检测点为各个检测点中最靠近所述目标储层的上分界面的检测点，或者为各个所述检测点中最靠近所述目标储层的下分界面的检测点。
7.根据权利要求5所述的装置，其特征在于，所述获得单元具体用于：
将所述参考检测点的下行波数据对应的振幅谱中所述参考频率对应的最大振幅值除以所述预设频率范围内每个频率对应的最大振幅值的结果，作为所述预设频率范围内每个频率对应的谱均衡系数。
8.根据权利要求5至7任一项所述的装置，其特征在于，还包括显示单元：
所述显示单元，用于将所述目标振幅谱以图形显示。
9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至4任一项所述方法的步骤。
10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该程序被处理器执行时实现权利要求1至4任一项所述方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910575689.3A CN112147675B (zh) | 2019-06-28 | 2019-06-28 | 一种储层流体检测方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910575689.3A CN112147675B (zh) | 2019-06-28 | 2019-06-28 | 一种储层流体检测方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112147675A CN112147675A (zh) | 2020-12-29 |
| CN112147675B true CN112147675B (zh) | 2022-07-05 |
Family
ID=73869384
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910575689.3A Active CN112147675B (zh) | 2019-06-28 | 2019-06-28 | 一种储层流体检测方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112147675B (zh) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011109839A2 (en) * | 2010-03-05 | 2011-09-09 | Vialogy Llc | Active noise injection computations for improved predictability in oil and gas reservoir discovery and characterization |
| CN102478670A (zh) * | 2010-11-24 | 2012-05-30 | 中国石油天然气集团公司 | 一种利用地震衰减属性进行储层流体性质预测的方法 |
| CN104155693A (zh) * | 2014-08-29 | 2014-11-19 | 成都理工大学 | 储层流体流度的角道集地震响应数值计算方法 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6941227B2 (en) * | 2001-04-30 | 2005-09-06 | The Regents Of The University Of California | Frequency-dependent processing and interpretation (FDPI) of seismic data for identifying, imaging and monitoring fluid-saturated underground reservoirs |
| CN102305943B (zh) * | 2011-07-22 | 2013-09-04 | 中国石油天然气股份有限公司 | 基于地震子波衰减谱的油气检测方法及装置 |
-
2019
- 2019-06-28 CN CN201910575689.3A patent/CN112147675B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2011109839A2 (en) * | 2010-03-05 | 2011-09-09 | Vialogy Llc | Active noise injection computations for improved predictability in oil and gas reservoir discovery and characterization |
| CN102478670A (zh) * | 2010-11-24 | 2012-05-30 | 中国石油天然气集团公司 | 一种利用地震衰减属性进行储层流体性质预测的方法 |
| CN104155693A (zh) * | 2014-08-29 | 2014-11-19 | 成都理工大学 | 储层流体流度的角道集地震响应数值计算方法 |
Non-Patent Citations (7)
| Title |
|---|
| Q-compensated acoustic impedance inversion of attenuated seismic data: Numerical and field-data experiments;Xintao Chai et al.;《GEOPHYSICS》;20181231;第83卷(第6期);第R553-R567页 * |
| Seismic low-frequency effects in monitoring fluid-saturated reservoirs;Valeri A. Korneev et al.;《GEOPHYSICS》;20040430;第69卷(第2期);第522-532页 * |
| 一种基于地震干涉测量的 VSP 资料波场重构及成像应用;梁上林等;《地球物理学进展》;20181231;第33卷(第3期);第1075-1080页 * |
| 地层岩性油气藏勘探;熊翥;《岩性油气藏》;20081231;第20卷(第04期);第1-8页 * |
| 基于叠前反演的流体识别与频散 AVO 反演研究;周伟等;《科学技术与工程》;20171130;第17卷(第33期);第26-35页 * |
| 基于吸收和速度频散的孔隙弹性介质流体识别研究;王永生等;《石油天然气学报》;20130331;第35卷(第03期);第72-78页 * |
| 敏感频散流体识别因子的构建与应用;丁娟;《资源信息与工程》;20180430;第33卷(第2期);第48-50、52页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112147675A (zh) | 2020-12-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US7391675B2 (en) | Microseismic event detection and location by continuous map migration | |
| US6789018B1 (en) | Mapping reservoir rocks using frequency spectral broadening and the presence of the slow-wave | |
| CA3032077C (en) | Seismic spectral balancing | |
| CA2931435A1 (en) | Method for developing a geomechanical model based on seismic data, well logs and sem analysis of horizontal and vertical drill cuttings | |
| CN110286410B (zh) | 基于绕射波能量的裂缝反演方法和装置 | |
| CN111856566A (zh) | 湖相滩坝砂体中薄储层预测方法及装置 | |
| WO2016134210A1 (en) | Microseismic sensitivity analysis and scenario modelling | |
| Germán Rubino et al. | Seismic characterization of thin beds containing patchy carbon dioxide-brine distributions: A study based on numerical simulations | |
| CN110389382A (zh) | 一种基于卷积神经网络的油气藏储层表征方法 | |
| EP3211594A1 (en) | Seismic modeling system providing seismic survey data inpainting based upon suspect region boundary comparisons and related methods | |
| CN113552624B (zh) | 孔隙度预测方法及装置 | |
| CN112147677B (zh) | 油气储层参数标签数据生成方法及装置 | |
| CN109471162A (zh) | 层间多次波处理方法、系统、电子设备及可读介质 | |
| CN112147675B (zh) | 一种储层流体检测方法及装置 | |
| CN106125133B (zh) | 一种基于气云区约束下的精细速度建模方法 | |
| US20130262070A1 (en) | System and method for subsurface reservoir characterization | |
| CN111381280A (zh) | 预测储层含烃饱和度的方法及装置 | |
| CN111239809A (zh) | 一种基于属性融合原理的优势储层确定方法及系统 | |
| CN113267809B (zh) | I类页岩储层预测方法及装置 | |
| CN108107482B (zh) | 一种少井条件下烃源岩多参数预测方法及系统 | |
| CN113050191A (zh) | 基于双参数的页岩油toc预测方法及装置 | |
| RU2789759C1 (ru) | Способ определения границ трещиноватой зоны | |
| Carcione et al. | Seismic modeling study of the Earth's deep crust | |
| CN111856584B (zh) | 一种基于超道集的薄储层烃类检测方法及装置 | |
| CN112764094B (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 |