CN111221038A - 薄储层厚度定量预测的方法和装置 - Google Patents
薄储层厚度定量预测的方法和装置 Download PDFInfo
- Publication number
- CN111221038A CN111221038A CN201811415763.7A CN201811415763A CN111221038A CN 111221038 A CN111221038 A CN 111221038A CN 201811415763 A CN201811415763 A CN 201811415763A CN 111221038 A CN111221038 A CN 111221038A
- Authority
- CN
- China
- Prior art keywords
- thin reservoir
- thickness
- thin
- reservoir
- tuned
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 230000003595 spectral effect Effects 0.000 claims abstract description 49
- 239000004215 Carbon black (E152) Substances 0.000 claims abstract description 30
- 229930195733 hydrocarbon Natural products 0.000 claims abstract description 30
- 150000002430 hydrocarbons Chemical class 0.000 claims abstract description 30
- 238000012546 transfer Methods 0.000 claims abstract description 6
- 238000006243 chemical reaction Methods 0.000 claims description 25
- 238000004364 calculation method Methods 0.000 claims description 10
- 238000001228 spectrum Methods 0.000 claims description 8
- 238000011161 development Methods 0.000 abstract description 4
- 230000006870 function Effects 0.000 description 35
- 239000007789 gas Substances 0.000 description 26
- 238000005259 measurement Methods 0.000 description 5
- 238000003860 storage Methods 0.000 description 5
- 230000015654 memory Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000009412 basement excavation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 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
- 239000003345 natural gas Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 238000011084 recovery 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/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
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
本发明公开了一种薄储层厚度定量预测的方法和装置,属于油气开发技术领域。所述方法包括:测量目标区域内的多口已勘测油气井的薄储层的实际厚度;基于每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,计算薄储层的调谐厚度;基于每口已勘测油气井的薄储层的调谐厚度和实际厚度,确定所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数;基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度;将所述待测位置的薄储层的调谐厚度,代入所述转换函数,得到所述待测位置的薄储层的实际厚度。采用本发明,可以得到较为准确的待测位置的薄储层的实际厚度。
Description
技术领域
本发明涉及油气开发技术领域,特别涉及一种薄储层厚度定量预测的方法和装置。
背景技术
石油和天然气是国家重要战略资源,是国家经济发展的命脉。而厚储层水淹现象严重,现在我国东部老油田油藏综合含水达到了90%以上,剩余油主要富集在薄储层中,然而薄储层厚度很难准确预测,会严重影响注采井网络部署和石油挖掘。因此,薄储层厚度很难准确预测成为了限制老油田采收率的一大问题。
现有技术中对薄储层厚度定量预测主要是利用薄储层的调谐厚度公式来直接估算出薄储层的实际厚度。
在实现本发明的过程中,发明人发现现有技术至少存在以下问题:
直接使用薄储层的调谐公式估算得到的薄储层的预测厚度与实际钻井所得到的薄储层的实际厚度间,误差很大。
发明内容
为了解决现有技术的问题,本发明实施例提供了一种薄储层厚度定量预测的方法和装置。所述技术方案如下:
第一方面,提供了一种薄储层厚度定量预测的方法,所述方法包括:
测量并记录目标区域内的多口已勘测油气井的薄储层的实际厚度;
获取所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率;
基于所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,计算每口已勘测油气井的薄储层的调谐厚度;
基于每口已勘测油气井的薄储层的调谐厚度和实际厚度,进行函数拟合,确定所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数;
获取所述目标区域内待测位置的薄储层处的地震波的传播速度,基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度;
将所述待测位置的薄储层的调谐厚度,代入所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数,得到所述待测位置的薄储层的实际厚度。
可选的,所述基于所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,计算每口已勘测油气井的薄储层的调谐厚度,包括:
将所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,代入薄储层调谐厚度公式:Z=v/(4f0),计算每口已勘测油气井的薄储层的调谐厚度;
其中,Z为薄储层调谐厚度,v为薄储层处的地震波的传播速度,f0为薄储层处的地震波的第一谱峰频率;
所述基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度,包括:
将所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,代入薄储层调谐厚度公式:Z=v/(4f0),计算所述目标区域内待测位置的薄储层的调谐厚度。
可选的,所述基于每口已勘测油气井的薄储层的调谐厚度和实际厚度,进行函数拟合,确定所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数,包括:
基于每口已勘测油气井的薄储层的调谐厚度和实际厚度,进行函数拟合,确定所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数:
Z实际=a×Z+b;
其中,Z实际为目标区域内薄储层的实际厚度,Z为目标区域内薄储层的调谐厚度,a和b为基于所述函数拟合确定的常数系数。
可选的,所述获取所述目标区域内待测位置的薄储层处的地震波的传播速度,基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度,包括:
基于所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度,建立所述目标区域内薄储层的地震波速度场;
基于所述目标区域内薄储层的地震波速度场,确定所述目标区域内待测位置的薄储层处的地震波的传播速度;
基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度。
可选的,所述得到所述待测位置的薄储层的实际厚度,之后还包括:
基于所述薄储层的实际厚度,确定薄储层的砂岩厚度等值图。
第二方面,提供了一种薄储层厚度定量预测的装置,所述装置包括:
测量模块,同于测量并记录目标区域内的多口已勘测油气井的薄储层的实际厚度;
获取模块,用于获取所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率;
第一计算模块,用于基于所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,计算每口已勘测油气井的薄储层的调谐厚度;
拟合模块,同于基于每口已勘测油气井的薄储层的调谐厚度和实际厚度,进行函数拟合,确定所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数;
第二计算模块,用于获取所述目标区域内待测位置的薄储层处的地震波的传播速度,基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度;
转换模块,用于将所述待测位置的薄储层的调谐厚度,代入所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数,得到所述待测位置的薄储层的实际厚度。
可选的,所述第一计算模块,用于:
将所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,代入薄储层调谐厚度公式:Z=v/(4f0),计算每口已勘测油气井的薄储层的调谐厚度;
其中,Z为薄储层调谐厚度,v为薄储层处的地震波的传播速度,f0为薄储层处的地震波的第一谱峰频率;
所述第二计算模块,用于:
将所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,代入薄储层调谐厚度公式:Z=v/(4f0),计算所述目标区域内待测位置的薄储层的调谐厚度。
可选的,所述拟合模块,用于:
基于每口已勘测油气井的薄储层的调谐厚度和实际厚度,进行函数拟合,确定所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数:
Z实际=a×Z+b;
其中,Z实际为目标区域内薄储层的实际厚度,Z为目标区域内薄储层的调谐厚度,a和b为基于所述函数拟合确定的常数系数。
可选的,所述第二计算模块,用于:
基于所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度,建立所述目标区域内薄储层的地震波速度场;
基于所述目标区域内薄储层的地震波速度场,确定所述目标区域内待测位置的薄储层处的地震波的传播速度;
基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度。
可选的,所述装置还包括:
确定模块,用于基于所述薄储层的实际厚度,确定薄储层的砂岩厚度等值图。
第三方面,提供了一种终端,所述终端包括钻井测井设备、处理器和存储器,所述存储器中存储有至少一条指令,所述至少一条指令所述处理器加载并执行以实现如上述第一方面所述的薄储层厚度定量预测的方法。
第四方面,提供了一种计算机可读存储介质,所述存储介质中存储有至少一条指令,所述至少一条指令由所述处理器加载并执行以实现如上述第一方面所述的薄储层厚度定量预测的方法。
本发明实施例提供的技术方案带来的有益效果至少包括:
本发明实施例中,首先获取目标区域内的多口已经勘测的油气井的薄储层的实际厚度,然后再基于这些已经勘探的油气井的薄储层的层速度和第一谱峰频率,计算出这些油气井的薄储层的调谐厚度。再通过函数拟合确定调谐厚度和实际厚度之间的转换函数。然后,通过获取待测位置的薄储层处的层速度和第一谱峰频率,计算出待测位置的薄储层的调谐厚度,最后将该调谐厚度代入到转换函数中,便可以得到待测位置的薄储层的实际厚度。这样,不是把调谐厚度直接作为实际厚度,而是通过拟合出调谐厚度和实际厚度之间的函数关系,再根据该函数关系求取实际厚度,能有效的减少实际厚度测量的误差,从而得到较为准确的待测位置的薄储层的实际厚度。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种薄储层厚度定量预测的方法的流程图;
图2是本发明实施例提供的一种薄储层的调谐厚度和实际厚度的拟合关系图;
图3是本发明实施例提供的一种薄储层厚度定量预测的装置的结构示意图;
图4是本发明实施例提供的一种终端的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
本发明实施例提供了一种薄储层厚度定量预测的方法,该方法可以由终端实现。其中,终端可以包括有钻井和测井设备、数据处理设备等。钻井和测井设备设置在某井区内,用于获取薄储层的层速度、地震波信息等数据。井区内要包括有已经勘测过的油气井,待测位置可以位于多个已勘测的油气井之间。
如图1所示,该方法的处理流程可以包括如下的步骤:
在步骤101中,测量并记录目标区域内的多口已勘测油气井的薄储层的实际厚度。
其中,目标区域可以是包括有待测位置和多口已勘测油气井的井区。
在实施中,技术人员对目标区域内的多口已经实际勘探的油气井的薄储层的实际厚度进行测量。再对每口井的薄储层的实际厚度进行记录,并存储。下面表1为记录有井区X中的已勘测油气井的薄储层实际厚度的表格。
表1
在步骤102中,获取多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率。
其中,薄储层处的地震波的传播速度可以称为薄储层的层速度。
在实施中,技术人员可以采用声波测井法,获得每口已勘测油气井的薄储层的声波时差,而声波时差的倒数就是该薄储层所对应的层速度。然后根据获取的地震波信息,制作出薄储层处的地震波的频谱,从频谱中可以获取地震波的第一谱峰频率。
在步骤103中,基于多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,计算每口已勘测油气井的薄储层的调谐厚度。
在实施中,将获取的多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,代入到薄储层调谐厚度公式:Z=v/(4f0),便可以计算得出每口已勘测油气井的薄储层的调谐厚度。其中,Z为薄储层调谐厚度,v为薄储层处的地震波的传播速度,f0为薄储层处的地震波的第一谱峰频率。需要说明的是,该薄储层调谐厚度公式为现有公式,其推导过程不再说明。
在步骤104中,基于每口已勘测油气井的薄储层的调谐厚度和实际厚度,进行函数拟合,确定目标区域内薄储层的调谐厚度与实际厚度之间的转换函数。
在实施中,技术人员可以利用MATLAB等具有函数拟合功能的软件,基于表2所示的每口已勘测油气井的薄储层的调谐厚度和实际厚度的数据,进行函数拟合。
表2
井号 | 调谐厚度 | 实际厚度 |
G1 | Z<sub>1</sub> | Z<sub>实际1</sub> |
G2 | Z<sub>2</sub> | Z<sub>实际2</sub> |
G3 | Z<sub>3</sub> | Z<sub>实际3</sub> |
G4 | Z<sub>4</sub> | Z<sub>实际4</sub> |
G5 | Z<sub>5</sub> | Z<sub>实际5</sub> |
...... | ...... | ...... |
可以获得如图2所示的拟合结果。从而可以确定出目标区域内薄储层的调谐厚度与实际厚度之间的转换函数:
Z实际=a×Z+b;
其中,Z实际为目标区域内薄储层的实际厚度,Z为目标区域内薄储层的调谐厚度,a和b为基于所述函数拟合确定的常数系数。下面以某油气开发区X井区的勘探为例,进行说明。对X井区内16口已经勘测的油井G101、G102、G103、G104、G105,......G116的第一谱峰频率、声波时差和实际厚度数据进行获取,可以得到,如下表3的数据。
表3
需要说明的是,表3中示例性的提供了16口已经勘测的油井中的5口油井的数据,其中的声波时差单位为μs/m,那么,通过求取声波时差的倒数来,获得测速度时,得到的层速度单位就为m/μs,在进行函数拟合时,要对该层速度的单位进行换算,换算关系为1m/μs=1×106m/s,通过换算可以得到G101的薄储层的层速度约为3012.0m/s,G102的薄储层的层速度约为2796.9m/s,G103的薄储层的层速度约为2508.8m/s,G104的薄储层的层速度为2561.5m/s,G105的薄储层的层速度约为2708.3m/s,以此类推,获得全部16口井的薄储层的层速度。然后,利用上述薄储层的调谐厚度公式Z=v/(4f0),可以得到G101的薄储层的调谐厚度约为30.12m,G102的薄储层的调谐厚度约为23.31m,G103的薄储层的调谐厚度约为20.91m,G104的薄储层的调谐厚度约为22.87m,G105的薄储层的调谐厚度约为19.34m,以此类推,获得全部16口油井的薄储层的调谐厚度。然后,将这16个油井的薄储层的调谐厚度和实际厚度输入到函数拟合软件中,可以得到薄储层的调谐厚度和实际厚度间的函数关系为:
Z实际=1.062×Z-12.54。
在步骤105中,获取目标区域内待测位置的薄储层处的地震波的传播速度,基于目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算待测位置的薄储层的调谐厚度。
在实施中,可以使用DepthTeamExpress转换工具,利用目标区域内的多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度,对井间薄储层处的地震波的传播速度进行插值,建立该目标区域内薄储层的地震波速度场,然后可以根据待测位置在该目标区域内薄储层的地震波速度场中所对应的位置,便可以获得待测位置的薄储层处的地震波的传播速度。再根据获取的地震波信息,制作出薄储层处的地震波的频谱,从频谱中可以获取地震波的第一谱峰频率。然后将得到的待测位置的薄储层处的地震波的传播速度和第一谱峰频率代入到公式:Z=v/(4f0)中,便可以得到待测位置的薄储层的调谐厚度。
在步骤106中,将待测位置的薄储层的调谐厚度,代入目标区域内薄储层的调谐厚度与实际厚度之间的转换函数,得到待测位置的薄储层的实际厚度。
在实施中,将得到的待测位置的薄储层的调谐厚度,代入到目标区域内薄储层的调谐厚度与实际厚度之间的转换函数:Z实际=a×Z+b中,便可以获得待测位置的薄储层的实际厚度。
可选的,为了反映目标区域的薄储层分布状态,可以在获取待测位置薄储层的实际厚度之后,还可以进行如下处理:基于薄储层的实际厚度,确定薄储层的砂岩厚度等值图。
在实施中,薄储层主要为砂岩岩层,在获得待测位置的薄储层的实际厚度后,可以结合已勘测的油气井的薄储层的实际厚度,标定出砂岩岩层的边界,绘制出砂岩岩层厚度等值图。
基于相同的技术构思,本发明实施例还提供了一种薄储层厚度定量预测的装置,如图3所示,该装置包括:测量模块301,获取模块302、第一计算模块303、拟合模块304、第二计算模块305、转换模块306和确定模块307。
测量模块301,同于测量并记录目标区域内的多口已勘测油气井的薄储层的实际厚度;
获取模块302,用于获取所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率;
第一计算模块303,用于基于所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,计算每口已勘测油气井的薄储层的调谐厚度;
拟合模块304,同于基于每口已勘测油气井的薄储层的调谐厚度和实际厚度,进行函数拟合,确定所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数;
第二计算模块305,用于获取所述目标区域内待测位置的薄储层处的地震波的传播速度,基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度;
转换模块306,用于将所述待测位置的薄储层的调谐厚度,代入所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数,得到所述待测位置的薄储层的实际厚度。
可选的,所述第一计算模块303,用于:
将所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,代入薄储层调谐厚度公式:Z=v/(4f0),计算每口已勘测油气井的薄储层的调谐厚度;
其中,Z为薄储层调谐厚度,v为薄储层处的地震波的传播速度,f0为薄储层处的地震波的第一谱峰频率;
所述第二计算模块305,用于:
将所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,代入薄储层调谐厚度公式:Z=v/(4f0),计算所述目标区域内待测位置的薄储层的调谐厚度。
可选的,所述拟合模块304,用于:
基于每口已勘测油气井的薄储层的调谐厚度和实际厚度,进行函数拟合,确定所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数:
Z实际=a×Z+b;
其中,Z实际为目标区域内薄储层的实际厚度,Z为目标区域内薄储层的调谐厚度,a和b为基于所述函数拟合确定的常数系数。
可选的,所述第二计算模块305,用于:
基于所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度,建立所述目标区域内薄储层的地震波速度场;
基于所述目标区域内薄储层的地震波速度场,确定所述目标区域内待测位置的薄储层处的地震波的传播速度;
基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度。
可选的,所述装置还包括:
确定模块307,用于基于所述薄储层的实际厚度,确定薄储层的砂岩厚度等值图。
关于上述实施例中的装置,其中各个模块执行操作的具体方式已经在有关该方法的实施例中进行了详细描述,此处将不做详细阐述说明。
需要说明的是:上述实施例提供的薄储层厚度定量预测的装置在定量预测薄储层厚度时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的功能模块完成,即将设备的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的薄储层厚度定量预测的装置与薄储层厚度定量预测的方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。
在示例性实施例中,还提供了一种计算机可读存储介质,存储介质中存储有至少一条指令,至少一条指令由处理器加载并执行以实现上述实施例中的识别动作类别的方法。例如,所述计算机可读存储介质可以是ROM、随机存取存储器(RAM)、CD-ROM、磁带、软盘和光数据存储设备等。
图4是本发明实施例提供的一种终端的结构示意图,该终端400可因配置或性能不同而产生比较大的差异,可以包括钻井测井设备401、一个或一个以上处理器(centralprocessing units,CPU)402和一个或一个以上的存储器403,其中,所述存储器403中存储有至少一条指令,所述至少一条指令由所述处理器401加载并执行以实现上述薄储层厚度定量预测的方法。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种薄储层厚度定量预测的方法,其特征在于,所述方法包括:
测量并记录目标区域内的多口已勘测油气井的薄储层的实际厚度;
获取所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率;
基于所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,计算每口已勘测油气井的薄储层的调谐厚度;
基于每口已勘测油气井的薄储层的调谐厚度和实际厚度,进行函数拟合,确定所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数;
获取所述目标区域内待测位置的薄储层处的地震波的传播速度,基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度;
将所述待测位置的薄储层的调谐厚度,代入所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数,得到所述待测位置的薄储层的实际厚度。
2.根据权利要求1所述的方法,其特征在于,所述基于所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,计算每口已勘测油气井的薄储层的调谐厚度,包括:
将所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,代入薄储层调谐厚度公式:Z=v/(4f0),计算每口已勘测油气井的薄储层的调谐厚度;
其中,Z为薄储层调谐厚度,v为薄储层处的地震波的传播速度,f0为薄储层处的地震波的第一谱峰频率;
所述基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度,包括:
将所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,代入薄储层调谐厚度公式:Z=v/(4f0),计算所述目标区域内待测位置的薄储层的调谐厚度。
3.根据权利要求1所述的方法,其特征在于,所述基于每口已勘测油气井的薄储层的调谐厚度和实际厚度,进行函数拟合,确定所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数,包括:
基于每口已勘测油气井的薄储层的调谐厚度和实际厚度,进行函数拟合,确定所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数:
Z实际=a×Z+b;
其中,Z实际为目标区域内薄储层的实际厚度,Z为目标区域内薄储层的调谐厚度,a和b为基于所述函数拟合确定的常数系数。
4.根据权利要求1所述方法,其特征在于,所述获取所述目标区域内待测位置的薄储层处的地震波的传播速度,基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度,包括:
基于所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度,建立所述目标区域内薄储层的地震波速度场;
基于所述目标区域内薄储层的地震波速度场,确定所述目标区域内待测位置的薄储层处的地震波的传播速度;
基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度。
5.根据权利要求1所述方法,其特征在于,所述得到所述待测位置的薄储层的实际厚度,之后还包括:
基于所述薄储层的实际厚度,确定薄储层的砂岩厚度等值图。
6.一种薄储层厚度定量预测的装置,其特征在于,所述装置包括:
测量模块,同于测量并记录目标区域内的多口已勘测油气井的薄储层的实际厚度;
获取模块,用于获取所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率;
第一计算模块,用于基于所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,计算每口已勘测油气井的薄储层的调谐厚度;
拟合模块,同于基于每口已勘测油气井的薄储层的调谐厚度和实际厚度,进行函数拟合,确定所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数;
第二计算模块,用于获取所述目标区域内待测位置的薄储层处的地震波的传播速度,基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度;
转换模块,用于将所述待测位置的薄储层的调谐厚度,代入所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数,得到所述待测位置的薄储层的实际厚度。
7.根据权利要求6所述的装置,其特征在于,所述第一计算模块,用于:
将所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度和第一谱峰频率,代入薄储层调谐厚度公式:Z=v/(4f0),计算每口已勘测油气井的薄储层的调谐厚度;
其中,Z为薄储层调谐厚度,v为薄储层处的地震波的传播速度,f0为薄储层处的地震波的第一谱峰频率;
所述第二计算模块,用于:
将所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,代入薄储层调谐厚度公式:Z=v/(4f0),计算所述目标区域内待测位置的薄储层的调谐厚度。
8.根据权利要求6所述的装置,其特征在于,所述拟合模块,用于:
基于每口已勘测油气井的薄储层的调谐厚度和实际厚度,进行函数拟合,确定所述目标区域内薄储层的调谐厚度与实际厚度之间的转换函数:
Z实际=a×Z+b;
其中,Z实际为目标区域内薄储层的实际厚度,Z为目标区域内薄储层的调谐厚度,a和b为基于所述函数拟合确定的常数系数。
9.根据权利要求6所述的装置,其特征在于,所述第二计算模块,用于:
基于所述多口已勘测油气井中每口已勘测油气井的薄储层处的地震波的传播速度,建立所述目标区域内薄储层的地震波速度场;
基于所述目标区域内薄储层的地震波速度场,确定所述目标区域内待测位置的薄储层处的地震波的传播速度;
基于所述目标区域内待测位置的薄储层处的地震波的传播速度和第一谱峰频率,计算所述待测位置的薄储层的调谐厚度。
10.根据权利要求6所述装置,其特征在于,所述装置还包括:
确定模块,用于基于所述薄储层的实际厚度,确定薄储层的砂岩厚度等值图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811415763.7A CN111221038B (zh) | 2018-11-26 | 2018-11-26 | 薄储层厚度定量预测的方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811415763.7A CN111221038B (zh) | 2018-11-26 | 2018-11-26 | 薄储层厚度定量预测的方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111221038A true CN111221038A (zh) | 2020-06-02 |
CN111221038B CN111221038B (zh) | 2022-02-01 |
Family
ID=70806738
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811415763.7A Active CN111221038B (zh) | 2018-11-26 | 2018-11-26 | 薄储层厚度定量预测的方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111221038B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113009567A (zh) * | 2021-03-31 | 2021-06-22 | 中国石油化工股份有限公司 | 泥包砂型薄储层地震特征的识别方法 |
CN113850007A (zh) * | 2021-11-30 | 2021-12-28 | 浙江中自庆安新能源技术有限公司 | 一种基于有限元分析的油藏有效油层厚度预测方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102109613A (zh) * | 2009-12-23 | 2011-06-29 | 中国石油天然气股份有限公司 | 一种复杂地质条件下目标储层有效厚度的确定方法 |
CN102478668A (zh) * | 2010-11-30 | 2012-05-30 | 中国石油天然气集团公司 | 一种应用地震多属性参数预测煤层厚度的方法 |
CN103217714A (zh) * | 2013-03-26 | 2013-07-24 | 中国石油大学(华东) | 基于频率-尺度匹配的地震沉积学解释方法 |
CN103412332A (zh) * | 2013-01-22 | 2013-11-27 | 中国地质大学(北京) | 一种定量计算薄储层厚度的方法 |
CN105005077A (zh) * | 2015-07-06 | 2015-10-28 | 成都理工大学 | 稀井条件下实钻井与虚拟井联合的薄层厚度预测方法 |
CN105974473A (zh) * | 2016-06-15 | 2016-09-28 | 中国石油天然气集团公司 | 一种确定薄储层的地震属性的方法和装置 |
CN106405647A (zh) * | 2016-12-22 | 2017-02-15 | 成都晶石石油科技有限公司 | 一种关于沉积地层厚度的调谐反演方法 |
US20180088249A1 (en) * | 2015-04-17 | 2018-03-29 | Schlumberger Technology Corporation | Thin Bed Tuning Frequency and Thickness Estimation |
-
2018
- 2018-11-26 CN CN201811415763.7A patent/CN111221038B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102109613A (zh) * | 2009-12-23 | 2011-06-29 | 中国石油天然气股份有限公司 | 一种复杂地质条件下目标储层有效厚度的确定方法 |
CN102478668A (zh) * | 2010-11-30 | 2012-05-30 | 中国石油天然气集团公司 | 一种应用地震多属性参数预测煤层厚度的方法 |
CN103412332A (zh) * | 2013-01-22 | 2013-11-27 | 中国地质大学(北京) | 一种定量计算薄储层厚度的方法 |
CN103217714A (zh) * | 2013-03-26 | 2013-07-24 | 中国石油大学(华东) | 基于频率-尺度匹配的地震沉积学解释方法 |
US20180088249A1 (en) * | 2015-04-17 | 2018-03-29 | Schlumberger Technology Corporation | Thin Bed Tuning Frequency and Thickness Estimation |
CN105005077A (zh) * | 2015-07-06 | 2015-10-28 | 成都理工大学 | 稀井条件下实钻井与虚拟井联合的薄层厚度预测方法 |
CN105974473A (zh) * | 2016-06-15 | 2016-09-28 | 中国石油天然气集团公司 | 一种确定薄储层的地震属性的方法和装置 |
CN106405647A (zh) * | 2016-12-22 | 2017-02-15 | 成都晶石石油科技有限公司 | 一种关于沉积地层厚度的调谐反演方法 |
Non-Patent Citations (6)
Title |
---|
CHARLES I. PURYEAR,等: "Layer-thickness determination and stratigraphic interpretation using spectral inversion: Theory and application", 《GEOPHYSICS》 * |
GREG PARTYKA,等: "Interpretational applications of spectral decomposition in reservoir characterization", 《THE LEADING EDGE》 * |
刘化清,等: "薄砂体预测的地震沉积学研究方法", 《岩性油气藏》 * |
周宗良,等: "用最大熵谱分解定量预测曲流河薄砂体", 《断块油气田》 * |
孙鲁平,等: "薄层地震峰值频率与厚度关系研究", 《石油地球物理勘探》 * |
张会卿,等: "高含水油藏薄砂体精细刻画技术探讨", 《天津科技》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113009567A (zh) * | 2021-03-31 | 2021-06-22 | 中国石油化工股份有限公司 | 泥包砂型薄储层地震特征的识别方法 |
CN113009567B (zh) * | 2021-03-31 | 2022-07-05 | 中国石油化工股份有限公司 | 泥包砂型薄储层地震特征的识别方法 |
CN113850007A (zh) * | 2021-11-30 | 2021-12-28 | 浙江中自庆安新能源技术有限公司 | 一种基于有限元分析的油藏有效油层厚度预测方法及系统 |
CN113850007B (zh) * | 2021-11-30 | 2022-03-11 | 浙江中自庆安新能源技术有限公司 | 一种基于有限元分析的油藏有效油层厚度预测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111221038B (zh) | 2022-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20100191514A1 (en) | Method, Program and Computer System for Scaling Hydrocarbon Reservoir Model Data | |
CN111221038B (zh) | 薄储层厚度定量预测的方法和装置 | |
WO2013158873A2 (en) | System and method for calibrating permeability for use in reservoir modeling | |
CN116430474A (zh) | 一种地热资源储量计算方法、系统及电子设备 | |
CN111241652B (zh) | 一种确定地层原油粘度的方法及装置 | |
CN111535796A (zh) | 一种测井参数的获取方法和装置 | |
CN115012903B (zh) | 判别泥页岩层理构造发育的测井评价方法 | |
WO2016093793A1 (en) | Determining non-linear petrofacies using cross-plot partitioning | |
CN111594156B (zh) | 一种天然气水合物饱和度计算方法及系统 | |
US11346968B2 (en) | Estimation of reservoir flow properties from seismic data | |
CN115327630A (zh) | 声波测井校正方法、装置、设备、存储介质及程序产品 | |
CN112347707B (zh) | 油藏流动非均质性的表征方法及装置 | |
CN111458751B (zh) | 一种综合时深标定方法及装置 | |
RU2530324C2 (ru) | Способ определения положения координат глубин маркера при построении геологической модели месторождения | |
CN109781039B (zh) | 溶洞的纵向尺度定量预测的方法 | |
CN112925018A (zh) | 一种异常压力成因判别方法及装置 | |
CN113933910B (zh) | 物性参数同时计算的方法、系统、存储介质与电子设备 | |
US20240168189A1 (en) | Systems, devices, and methods for generating average velocity maps of subsurface formations | |
CN113123780B (zh) | 单井各级烃源岩厚度的确定方法、装置及存储介质 | |
CN114063169B (zh) | 一种波阻抗反演方法、系统、设备和存储介质 | |
CN113514890B (zh) | 一种使用地震数据预测地层压力系数的方法、装置及设备 | |
CN115034151B (zh) | 一种断层垂向优势运移通道分析方法及系统 | |
CN112746839B (zh) | 一种原油体积系数的确定方法及装置 | |
Qin et al. | High‐resolution characterization of the firn layer near the West Antarctic ice sheet divide camp with active and passive seismic data | |
CN113622903A (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 |