CN112394416A - 非均质断裂控制储层预测方法及装置 - Google Patents
非均质断裂控制储层预测方法及装置 Download PDFInfo
- Publication number
- CN112394416A CN112394416A CN202011126267.7A CN202011126267A CN112394416A CN 112394416 A CN112394416 A CN 112394416A CN 202011126267 A CN202011126267 A CN 202011126267A CN 112394416 A CN112394416 A CN 112394416A
- Authority
- CN
- China
- Prior art keywords
- data
- seismic
- heterogeneous
- well
- determining
- 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 40
- 230000004044 response Effects 0.000 claims abstract description 57
- 235000019738 Limestone Nutrition 0.000 claims abstract description 43
- 239000006028 limestone Substances 0.000 claims abstract description 43
- 239000011435 rock Substances 0.000 claims abstract description 17
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 15
- 238000010276 construction Methods 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims description 36
- 238000005553 drilling Methods 0.000 claims description 35
- 238000004590 computer program Methods 0.000 claims description 15
- 238000012937 correction Methods 0.000 claims description 14
- 238000011160 research Methods 0.000 claims description 13
- 230000015572 biosynthetic process Effects 0.000 claims description 8
- 238000001514 detection method Methods 0.000 claims description 8
- 238000003860 storage Methods 0.000 claims description 8
- 239000003921 oil Substances 0.000 description 19
- 238000010586 diagram Methods 0.000 description 15
- 238000011161 development Methods 0.000 description 10
- 238000004519 manufacturing process Methods 0.000 description 8
- 238000005755 formation reaction Methods 0.000 description 6
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 description 5
- 235000009508 confectionery Nutrition 0.000 description 5
- 239000012530 fluid Substances 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000005520 cutting process Methods 0.000 description 3
- 238000013467 fragmentation Methods 0.000 description 3
- 238000006062 fragmentation reaction Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000035515 penetration Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 239000011324 bead Substances 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 235000019994 cava Nutrition 0.000 description 1
- 239000004568 cement Substances 0.000 description 1
- 238000005260 corrosion Methods 0.000 description 1
- 230000007797 corrosion Effects 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 235000021185 dessert Nutrition 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 239000010459 dolomite Substances 0.000 description 1
- 229910000514 dolomite Inorganic materials 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 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/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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/616—Data from specific type of measurement
- G01V2210/6169—Data from specific type of measurement using well-logging
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种非均质断裂控制储层预测方法及装置，其中方法包括：获得研究区的地震数据和测井数据；根据所述地震数据和测井数据，确定井震标定关系数据；对所述地震数据进行子波分解，得到致密岩性区的多子波库；根据所述多子波库，确定致密灰岩地层的背景数据；根据所述地震数据和致密灰岩地层的背景数据，确定地震响应特征数据；根据所述地震响应特征数据，井震标定关系数据和预先建立的均质阻抗模型，得到非均质阻抗模型，所述均质阻抗模型根据历史构造数据和历史测井数据预先建立；根据所述非均质阻抗模型，进行非均质断裂控制储层预测。本发明可以预测非均质断裂控制储层，有效提高储层预测精度。
Description
技术领域
本发明涉及油气勘探开发技术领域，尤其涉及非均质断裂控制储层预测方法及装置。
背景技术
在石油和天然气的勘探、开发行业中，把保存油气的圈闭称为油藏。塔里木盆地碳酸盐岩是塔里木油田油气勘探开发的重要领域，而塔里木碳酸盐岩储层非均质性强，储集层特征主要表现为三大类：裂缝-洞穴型、裂缝-孔洞型和裂缝型储层，据统计，塔中地区钻遇裂缝-孔洞型储层的井数占86.21％，塔北地区钻遇裂缝-孔洞型储层的井数占31％，进入油田开发阶段，裂缝-孔洞型储层是最主要的研究目标，优选高产稳产目标即是要找到连通性好且有一定储集规模的缝洞体系。
塔里木盆地奥陶系碳酸盐岩地层内发育着大量的走滑断裂，岩溶液沿走滑断裂破碎带和破裂带、裂缝进行多期溶蚀形成一种不规则状的三维空间展布的裂缝-孔洞集合体。经过上覆盖层封堵及致密灰岩侧向遮挡等因素形成断控岩溶缝洞体圈闭，即“断溶体”。碳酸盐岩断溶体描述技术系列注重利用多层系的高精度属性、大量钻测井资料和生产数据，综合缝洞体系量化雕刻、高精度储层反演、三维地质建模等技术方法，精细地描述了区内走滑断裂系统及其与“断溶体”立体成藏的关系。实践证实，断溶体油气藏是塔里木盆地极为特殊的一种油气藏类型，不同位置的断溶体油气藏生产效果差异显著，总结其成藏规律对下一步井位优选起着至关重要的作用。
目前这类缝洞型储层是油田产能建设的核心靶区。特别是在实际钻井过程中，我们有数口井都在地震明显断裂破碎带位置遭遇了优质储层而出现大量泥浆漏失，很多井都获得稳定的高产工业油气流，产生了灰岩地层缝洞型储层产油速度的新纪录，产油速度达到500吨每天高产井。目前油田开发阶段，这种断裂破碎带的地震响应主要是次一级反射强度和中-小规模的串珠群反射集合，也有不少与周围反射差异较大的弱反射特征。针对断裂破碎带的研究是油田管理中的一项非常重要的工作。但是，这类储层非均质性非常强，储层预测面临着大困难。
非均质断裂控制形成的缝洞型储层精细研究一直是个世界级难题，实际生产研究中，多应用常规波阻抗反演、地震属性分析等方法预测该类储层。但是，采用这些方法不能有效提升描述这类非均质储层空间发育规律，更无法有效解释与开发生产动态数据的矛盾，储层精度预测低。
因此，亟需一种可以克服上述问题的非均质断裂控制储层预测方案。
发明内容
本发明实施例提供一种非均质断裂控制储层预测方法，用以预测非均质断裂控制储层，有效提高储层预测精度，该方法包括：
获得研究区的地震数据和测井数据；
根据所述地震数据和测井数据，确定井震标定关系数据；
对所述地震数据进行子波分解，得到致密岩性区的多子波库；
根据所述多子波库，确定致密灰岩地层的背景数据；
根据所述地震数据和致密灰岩地层的背景数据，确定地震响应特征数据；
根据所述地震响应特征数据，井震标定关系数据和预先建立的均质阻抗模型，得到非均质阻抗模型，所述均质阻抗模型根据历史构造数据和历史测井数据预先建立；
根据所述非均质阻抗模型，进行非均质断裂控制储层预测。
本发明实施例提供一种非均质断裂控制储层预测装置，用以预测非均质断裂控制储层，有效提高储层预测精度，该装置包括：
数据获得模块，用于获得研究区的地震数据和测井数据；
标定关系确定模块，用于根据所述地震数据和测井数据，确定井震标定关系数据；
子波分解模块，用于对所述地震数据进行子波分解，得到致密岩性区的多子波库；
背景数据确定模块，用于根据所述多子波库，确定致密灰岩地层的背景数据；
特征数据确定模块，用于根据所述地震数据和致密灰岩地层的背景数据，确定地震响应特征数据；
模型确定模块，用于根据所述地震响应特征数据，井震标定关系数据和预先建立的均质阻抗模型，得到非均质阻抗模型，所述均质阻抗模型根据历史构造数据和历史测井数据预先建立；
预测模块，用于根据所述非均质阻抗模型，进行非均质断裂控制储层预测。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述非均质断裂控制储层预测方法。
本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述非均质断裂控制储层预测方法的计算机程序。
相比于现有的应用常规波阻抗反演、地震属性分析等方法预测非均质断裂控制储层的方案而言，本发明实施例通过获得研究区的地震数据和测井数据；根据所述地震数据和测井数据，确定井震标定关系数据；对所述地震数据进行子波分解，得到致密岩性区的多子波库；根据所述多子波库，确定致密灰岩地层的背景数据；根据所述地震数据和致密灰岩地层的背景数据，确定地震响应特征数据；根据所述地震响应特征数据，井震标定关系数据和预先建立的均质阻抗模型，得到非均质阻抗模型，所述均质阻抗模型根据历史构造数据和历史测井数据预先建立；根据所述非均质阻抗模型，进行非均质断裂控制储层预测。本发明实施例确定致密灰岩地层的背景数据之后，根据地震数据和致密灰岩地层的背景数据确定地震响应特征数据，该地震响应特征数据能够有效反映走滑断裂破碎带及储层甜点，进而结合井震标定关系数据，将地震响应特征数据嵌套到预先建立的均质阻抗模型中，根据地震响应特征数据，井震标定关系数据和均质阻抗模型得到非均质阻抗模型，由于地震响应特征数据的嵌套使得模型具备断裂破碎带、断控缝洞储层的辨识能力，从而根据非均质阻抗模型进行非均质断裂控制储层预测，有效提升断裂破碎带、断控缝洞储层描述精度，为油气藏开发方案编制提供有力依据。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：
图1为本发明实施例中非均质断裂控制储层预测方法示意图；
图2为本发明实施例中非均质断裂控制储层预测方法示意图；
图3～图6为本发明具体实施例中非均质断裂控制储层预测方法示意图；
图7为本发明实施例中非均质断裂控制储层预测装置结构图；
图8为本发明实施例中非均质断裂控制储层预测装置结构图；
图9是本发明实施例的计算机设备结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。
为了预测非均质断裂控制储层，有效提高储层预测精度，本发明实施例提供一种非均质断裂控制储层预测方法，如图1所示，该方法可以包括：
步骤101、获得研究区的地震数据和测井数据；
步骤102、根据所述地震数据和测井数据，确定井震标定关系数据；
步骤103、对所述地震数据进行子波分解，得到致密岩性区的多子波库；
步骤104、根据所述多子波库，确定致密灰岩地层的背景数据；
步骤105、根据所述地震数据和致密灰岩地层的背景数据，确定地震响应特征数据；
步骤106、根据所述地震响应特征数据，井震标定关系数据和预先建立的均质阻抗模型，得到非均质阻抗模型，所述均质阻抗模型根据历史构造数据和历史测井数据预先建立；
步骤107、根据所述非均质阻抗模型，进行非均质断裂控制储层预测。
由图1所示可以得知，本发明实施例通过获得研究区的地震数据和测井数据；根据所述地震数据和测井数据，确定井震标定关系数据；对所述地震数据进行子波分解，得到致密岩性区的多子波库；根据所述多子波库，确定致密灰岩地层的背景数据；根据所述地震数据和致密灰岩地层的背景数据，确定地震响应特征数据；根据所述地震响应特征数据，井震标定关系数据和预先建立的均质阻抗模型，得到非均质阻抗模型，所述均质阻抗模型根据历史构造数据和历史测井数据预先建立；根据所述非均质阻抗模型，进行非均质断裂控制储层预测。本发明实施例确定致密灰岩地层的背景数据之后，根据地震数据和致密灰岩地层的背景数据确定地震响应特征数据，该地震响应特征数据能够有效反映走滑断裂破碎带及储层甜点，进而结合井震标定关系数据，将地震响应特征数据嵌套到预先建立的均质阻抗模型中，根据地震响应特征数据，井震标定关系数据和均质阻抗模型得到非均质阻抗模型，由于地震响应特征数据的嵌套使得模型具备断裂破碎带、断控缝洞储层的辨识能力，从而根据非均质阻抗模型进行非均质断裂控制储层预测，有效提升断裂破碎带、断控缝洞储层描述精度，为油气藏开发方案编制提供有力依据。
具体实施时，获得研究区的地震数据和测井数据。
实施例中，所述地震数据为叠后地震数据体。
具体实施时，根据所述地震数据和测井数据，确定井震标定关系数据。
实施例中，如图2所示，图1的非均质断裂控制储层预测方法还包括：步骤108、获得研究区测井数据之后，对所述测井数据进行标准化处理及VSP曲线校正处理；根据所述地震数据和测井数据，确定井震标定关系数据包括：根据所述地震数据，和标准化处理及VSP曲线校正处理后的测井数据，确定井震标定关系数据。
本实施例中，专门选择测量VSP的斜度井、水平井作为标定井，可以提升针对断裂破碎带宽度的数据分析精度。首先对选择井的测井数据进行标准化处理与VSP曲线校正处理，通过精细井震标定确定井震标定关系数据。
实施例中，所述井震标定关系数据包括：录井气测数据(TG)，钻井液漏失数据，钻具放空数据，地震振幅阈值其中之一或任意组合。录井气测数据(TG)是从安置在振动筛前的脱气器可获得从井底返回的钻井液所携带的气体，对其进行组分和含量的检测和编录，从而判断油气层的工作，它是直接测定钻井液中可燃气体含量的一种录井方法，在钻进过程中进行的，利用气测资料能及时发现油气显示，并能预报井喷，在钻井中广泛采用。钻井液漏失数据是指在钻进、固井、测试或修井等井下作业中各种工作液(包括钻井液、水泥浆、完井液以及其他流体等)在压差作用下直接进入地层的一种井下复杂情况。井漏包括渗透性滤失、裂缝性滤失、溶洞性滤失。钻具放空是钻入裂缝、溶洞时，悬重突增，钻压突减，钻具能无阻地送入一定长度的现象。钻时指每钻进单位厚度所需要的时间(min/m)。单位纯钻进时间的钻头进尺称为机械钻速单位为m/h。松软的地层比坚硬地层钻时低，如疏松砂岩比致密砂岩钻时低:多孔的碳酸盐岩比白云岩钻时低。在无测井资料或尚未进行测井的井段，钻时曲线与岩屑录井剖面相结合，是划分层位、与邻井地层对比、修正地质预告、卡准目的层、判断油气显示层位、确定钻井取心位置的重要依据。
具体实施时，对所述地震数据进行子波分解，得到致密岩性区的多子波库。
实施例中，致密岩性区可以为无储层区。对地震数据进行子波分解，统计分析目的层中致密岩性区的多子波库。
具体实施时，根据所述多子波库，确定致密灰岩地层的背景数据。
实施例中，根据所述多子波库，确定致密灰岩地层的背景数据包括：根据所述多子波库，进行子波重构处理；根据子波重构处理的结果，确定致密灰岩地层的背景数据。
具体实施时，根据所述地震数据和致密灰岩地层的背景数据，确定地震响应特征数据。
实施例中，根据所述地震数据和致密灰岩地层的背景数据，确定地震响应特征数据：对所述地震数据和致密灰岩地层的背景数据进行求差运算或求商运算，根据运算的结果确定地震响应特征数据。
本实施例中，对所求得多子波进行重构，得到可反映致密灰岩地层的背景数据，将原始地震数据与背景数据进行求差或求商运算，最终求得反映走滑断裂破碎带及储层甜点的地震响应特征数据。
具体实施时，根据所述地震响应特征数据，井震标定关系数据和预先建立的均质阻抗模型，得到非均质阻抗模型，所述均质阻抗模型根据历史构造数据和历史测井数据预先建立。
实施例中，根据所述地震响应特征数据，井震标定关系数据和预先建立的均质阻抗模型，得到非均质阻抗模型包括：根据所述地震响应特征数据和井震标定关系数据，确定对应的均方差数据；根据所述均方差数据和预设的交汇关系比例值，确定边界值；以所述边界值作为数据嵌套边界，将所述地震响应特征数据嵌套到预先建立的均质阻抗模型中，得到非均质阻抗模型。
本实施例中，按如下公式根据地震响应特征数据计算均方差RMS数据：
其中，A为地震响应特征数据，θi为地层倾角值，Δd为地震数据的道间距，Δt＝Δd×tan(θi),(i＝1,3,5,7)，
然后根据井震标定关系数据，求取RMS数据中与储层位置对应的RMS值，并通过交汇关系比例求得与测井阻抗值相匹配的值，从而确定断裂破碎带及储层甜点的边界值。
具体实施时，根据所述非均质阻抗模型，进行非均质断裂控制储层预测。
实施例中，根据所述非均质阻抗模型，进行非均质断裂控制储层预测包括：根据所述非均质阻抗模型进行井震约束反演，确定纵波阻抗数据；根据所述纵波阻抗数据，提取沿井轨迹的波阻抗数据；将所述沿井轨迹的波阻抗数据与测井曲线进行交汇，得到相关系数，若所述相关系数小于预设阈值，则对所述非均质阻抗模型进行井震误差校正处理，并重新提取沿井轨迹的波阻抗数据，直至相关系数达到预设阈值；根据所述波阻抗数据和测井数据，进行非均质断裂控制储层预测。
本实施例中，相关系数的预设阈值可以为0.97，达到0.97以上，视为合格，否则需对非均质阻抗模型进行井震误差校正，重复提取波阻抗数据，并与测井数据交汇，直至误差合格。
实施例中，所述测井数据包括：钻井曲线，测井孔隙度曲线和远探测电阻率曲线；根据所述波阻抗数据和测井数据，进行非均质断裂控制储层预测包括：根据所述波阻抗数据，钻井曲线，测井孔隙度曲线和远探测电阻率曲线，得到孔隙度数据体；根据所述孔隙度数据体，进行非均质断裂控制储层预测。以测井孔隙度1.8％的缝洞储层的下限，对储层进行刻画。
下面给出一个具体实施例，说明本发明实施例非均质断裂控制储层预测方法的具体应用。本实施例中，应用该发明在塔里木盆地塔北塔河南地区某油田A地区进行试验研究，该区走滑断裂破碎带非常发育，断裂控制溶蚀的储层集中发育。2015年，决定钻探WELL301井，采用常规反演储层预测预估可采出原油2.4315×103m3。目前实际采出油气当量超5.1120×103m3，并且目前保持高产状态。预测与实际的巨大差异，一方面是对油藏的厚度预判是不合适的，另一方面是对油藏横向规模预判与实际有误差。图3为相同位置地震剖面、常规阻抗模型剖面、断裂破碎带及储层甜点数据剖面与非均质阻抗模型剖面对比图，图4为相同位置地震剖面、非均质断控储层预测剖面、常规储层预测剖面对比及野外地质露头照片对比，图5为WELL301井区常规储层预测与非均质断控储层预测平面结果对比图，图6为A区块常规储层预测与非均质断控储层预测平面结果对比图。如图3～图6所示，本发明纵波阻抗数据剖面，可以清楚的储层亮点，储层沿断裂破碎带集中发育的特征明显。WELL301井区的特征是一个纵向、横向较大规模的油藏，很好地解释清楚WELL301井的生产动态与早期储量预测数据不匹配的矛盾。
基于同一发明构思，本发明实施例还提供了一种非均质断裂控制储层预测装置，如下面的实施例所述。由于这些解决问题的原理与非均质断裂控制储层预测方法相似，因此装置的实施可以参见方法的实施，重复之处不再赘述。
图7为本发明实施例中非均质断裂控制储层预测装置的结构图，如图7所示，该装置包括：
数据获得模块701，用于获得研究区的地震数据和测井数据；
标定关系确定模块702，用于根据所述地震数据和测井数据，确定井震标定关系数据；
子波分解模块703，用于对所述地震数据进行子波分解，得到致密岩性区的多子波库；
背景数据确定模块704，用于根据所述多子波库，确定致密灰岩地层的背景数据；
特征数据确定模块705，用于根据所述地震数据和致密灰岩地层的背景数据，确定地震响应特征数据；
模型确定模块706，用于根据所述地震响应特征数据，井震标定关系数据和预先建立的均质阻抗模型，得到非均质阻抗模型，所述均质阻抗模型根据历史构造数据和历史测井数据预先建立；
预测模块707，用于根据所述非均质阻抗模型，进行非均质断裂控制储层预测。
一个实施例中，如图8所示，图7的非均质断裂控制储层预测装置还包括：
数据处理模块708，用于获得研究区测井数据之后，对所述测井数据进行标准化处理及VSP曲线校正处理；
标定关系确定模块702进一步用于：根据所述地震数据，和标准化处理及VSP曲线校正处理后的测井数据，确定井震标定关系数据。
一个实施例中，所述井震标定关系数据包括：录井气测数据，钻井液漏失数据，钻具放空数据，地震振幅阈值其中之一或任意组合。
一个实施例中，所述背景数据确定模块704进一步用于：
根据所述多子波库，进行子波重构处理；
根据子波重构处理的结果，确定致密灰岩地层的背景数据。
一个实施例中，所述特征数据确定模块705进一步用于：
对所述地震数据和致密灰岩地层的背景数据进行求差运算或求商运算，根据运算的结果确定地震响应特征数据。
一个实施例中，所述模型确定模706块进一步用于：
根据所述地震响应特征数据和井震标定关系数据，确定对应的均方差数据；
根据所述均方差数据和预设的交汇关系比例值，确定边界值；
以所述边界值作为数据嵌套边界，将所述地震响应特征数据嵌套到预先建立的均质阻抗模型中，得到非均质阻抗模型。
一个实施例中，所述预测模块707进一步用于：
根据所述非均质阻抗模型进行井震约束反演，确定纵波阻抗数据；
根据所述纵波阻抗数据，提取沿井轨迹的波阻抗数据；
将所述沿井轨迹的波阻抗数据与测井曲线进行交汇，得到相关系数，若所述相关系数小于预设阈值，则对所述非均质阻抗模型进行井震误差校正处理，并重新提取沿井轨迹的波阻抗数据，直至相关系数达到预设阈值；
根据所述波阻抗数据和测井数据，进行非均质断裂控制储层预测。
一个实施例中，所述测井数据包括：钻井曲线，测井孔隙度曲线和远探测电阻率曲线；
所述预测模块707进一步用于：根据所述波阻抗数据，钻井曲线，测井孔隙度曲线和远探测电阻率曲线，得到孔隙度数据体；根据所述孔隙度数据体，进行非均质断裂控制储层预测。
综上所述，本发明实施例通过获得研究区的地震数据和测井数据；根据所述地震数据和测井数据，确定井震标定关系数据；对所述地震数据进行子波分解，得到致密岩性区的多子波库；根据所述多子波库，确定致密灰岩地层的背景数据；根据所述地震数据和致密灰岩地层的背景数据，确定地震响应特征数据；根据所述地震响应特征数据，井震标定关系数据和预先建立的均质阻抗模型，得到非均质阻抗模型，所述均质阻抗模型根据历史构造数据和历史测井数据预先建立；根据所述非均质阻抗模型，进行非均质断裂控制储层预测。本发明实施例确定致密灰岩地层的背景数据之后，根据地震数据和致密灰岩地层的背景数据确定地震响应特征数据，该地震响应特征数据能够有效反映走滑断裂破碎带及储层甜点，进而结合井震标定关系数据，将地震响应特征数据嵌套到预先建立的均质阻抗模型中，根据地震响应特征数据，井震标定关系数据和均质阻抗模型得到非均质阻抗模型，由于地震响应特征数据的嵌套使得模型具备断裂破碎带、断控缝洞储层的辨识能力，从而根据非均质阻抗模型进行非均质断裂控制储层预测，有效提升断裂破碎带、断控缝洞储层描述精度，为油气藏开发方案编制提供有力依据。
基于前述发明构思，如图9所示，本发明还提出了一种计算机设备900，包括存储器910、处理器920及存储在存储器910上并可在处理器920上运行的计算机程序930，所述处理器920执行所述计算机程序930时实现前述油气钻井岩屑返出量确定方法。
基于前述发明构思，本发明提出了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现前述油气钻井岩屑返出量确定方法。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应说明的是：以上所述实施例，仅为本发明的具体实施方式，用以说明本发明的技术方案，而非对其限制，本发明的保护范围并不局限于此，尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化，或者对其中部分技术特征进行等同替换；而这些修改、变化或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。
Claims (18)
1.一种非均质断裂控制储层预测方法，其特征在于，包括：
获得研究区的地震数据和测井数据；
根据所述地震数据和测井数据，确定井震标定关系数据；
对所述地震数据进行子波分解，得到致密岩性区的多子波库；
根据所述多子波库，确定致密灰岩地层的背景数据；
根据所述地震数据和致密灰岩地层的背景数据，确定地震响应特征数据；
根据所述地震响应特征数据，井震标定关系数据和预先建立的均质阻抗模型，得到非均质阻抗模型，所述均质阻抗模型根据历史构造数据和历史测井数据预先建立；
根据所述非均质阻抗模型，进行非均质断裂控制储层预测。
2.如权利要求1所述的非均质断裂控制储层预测方法，其特征在于，还包括：获得研究区测井数据之后，对所述测井数据进行标准化处理及VSP曲线校正处理；
根据所述地震数据和测井数据，确定井震标定关系数据包括：根据所述地震数据，和标准化处理及VSP曲线校正处理后的测井数据，确定井震标定关系数据。
3.如权利要求1所述的非均质断裂控制储层预测方法，其特征在于，所述井震标定关系数据包括：录井气测数据，钻井液漏失数据，钻具放空数据，地震振幅阈值其中之一或任意组合。
4.如权利要求1所述的非均质断裂控制储层预测方法，其特征在于，根据所述多子波库，确定致密灰岩地层的背景数据包括：
根据所述多子波库，进行子波重构处理；
根据子波重构处理的结果，确定致密灰岩地层的背景数据。
5.如权利要求1所述的非均质断裂控制储层预测方法，其特征在于，根据所述地震数据和致密灰岩地层的背景数据，确定地震响应特征数据：
对所述地震数据和致密灰岩地层的背景数据进行求差运算或求商运算，根据运算的结果确定地震响应特征数据。
6.如权利要求1所述的非均质断裂控制储层预测方法，其特征在于，根据所述地震响应特征数据，井震标定关系数据和预先建立的均质阻抗模型，得到非均质阻抗模型包括：
根据所述地震响应特征数据和井震标定关系数据，确定对应的均方差数据；
根据所述均方差数据和预设的交汇关系比例值，确定边界值；
以所述边界值作为数据嵌套边界，将所述地震响应特征数据嵌套到预先建立的均质阻抗模型中，得到非均质阻抗模型。
7.如权利要求1所述的非均质断裂控制储层预测方法，其特征在于，根据所述非均质阻抗模型，进行非均质断裂控制储层预测包括：
根据所述非均质阻抗模型进行井震约束反演，确定纵波阻抗数据；
根据所述纵波阻抗数据，提取沿井轨迹的波阻抗数据；
将所述沿井轨迹的波阻抗数据与测井曲线进行交汇，得到相关系数，若所述相关系数小于预设阈值，则对所述非均质阻抗模型进行井震误差校正处理，并重新提取沿井轨迹的波阻抗数据，直至相关系数达到预设阈值；
根据所述波阻抗数据和测井数据，进行非均质断裂控制储层预测。
8.如权利要求7所述的非均质断裂控制储层预测方法，其特征在于，所述测井数据包括：钻井曲线，测井孔隙度曲线和远探测电阻率曲线；
根据所述波阻抗数据和测井数据，进行非均质断裂控制储层预测包括：根据所述波阻抗数据，钻井曲线，测井孔隙度曲线和远探测电阻率曲线，得到孔隙度数据体；根据所述孔隙度数据体，进行非均质断裂控制储层预测。
9.一种非均质断裂控制储层预测装置，其特征在于，包括：
数据获得模块，用于获得研究区的地震数据和测井数据；
标定关系确定模块，用于根据所述地震数据和测井数据，确定井震标定关系数据；
子波分解模块，用于对所述地震数据进行子波分解，得到致密岩性区的多子波库；
背景数据确定模块，用于根据所述多子波库，确定致密灰岩地层的背景数据；
特征数据确定模块，用于根据所述地震数据和致密灰岩地层的背景数据，确定地震响应特征数据；
模型确定模块，用于根据所述地震响应特征数据，井震标定关系数据和预先建立的均质阻抗模型，得到非均质阻抗模型，所述均质阻抗模型根据历史构造数据和历史测井数据预先建立；
预测模块，用于根据所述非均质阻抗模型，进行非均质断裂控制储层预测。
10.如权利要求9所述的非均质断裂控制储层预测装置，其特征在于，还包括：
数据处理模块，用于获得研究区测井数据之后，对所述测井数据进行标准化处理及VSP曲线校正处理；
标定关系确定模块进一步用于：根据所述地震数据，和标准化处理及VSP曲线校正处理后的测井数据，确定井震标定关系数据。
11.如权利要求9所述的非均质断裂控制储层预测装置，其特征在于，所述井震标定关系数据包括：录井气测数据，钻井液漏失数据，钻具放空数据，地震振幅阈值其中之一或任意组合。
12.如权利要求9所述的非均质断裂控制储层预测装置，其特征在于，所述背景数据确定模块进一步用于：
根据所述多子波库，进行子波重构处理；
根据子波重构处理的结果，确定致密灰岩地层的背景数据。
13.如权利要求9所述的非均质断裂控制储层预测装置，其特征在于，所述特征数据确定模块进一步用于：
对所述地震数据和致密灰岩地层的背景数据进行求差运算或求商运算，根据运算的结果确定地震响应特征数据。
14.如权利要求9所述的非均质断裂控制储层预测装置，其特征在于，所述模型确定模块进一步用于：
根据所述地震响应特征数据和井震标定关系数据，确定对应的均方差数据；
根据所述均方差数据和预设的交汇关系比例值，确定边界值；
以所述边界值作为数据嵌套边界，将所述地震响应特征数据嵌套到预先建立的均质阻抗模型中，得到非均质阻抗模型。
15.如权利要求9所述的非均质断裂控制储层预测装置，其特征在于，所述预测模块进一步用于：
根据所述非均质阻抗模型进行井震约束反演，确定纵波阻抗数据；
根据所述纵波阻抗数据，提取沿井轨迹的波阻抗数据；
将所述沿井轨迹的波阻抗数据与测井曲线进行交汇，得到相关系数，若所述相关系数小于预设阈值，则对所述非均质阻抗模型进行井震误差校正处理，并重新提取沿井轨迹的波阻抗数据，直至相关系数达到预设阈值；
根据所述波阻抗数据和测井数据，进行非均质断裂控制储层预测。
16.如权利要求15所述的非均质断裂控制储层预测装置，其特征在于，所述测井数据包括：钻井曲线，测井孔隙度曲线和远探测电阻率曲线；
所述预测模块进一步用于：根据所述波阻抗数据，钻井曲线，测井孔隙度曲线和远探测电阻率曲线，得到孔隙度数据体；根据所述孔隙度数据体，进行非均质断裂控制储层预测。
17.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至8任一所述方法。
18.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质存储有执行权利要求1至8任一所述方法的计算机程序。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011126267.7A CN112394416A (zh) | 2020-10-20 | 2020-10-20 | 非均质断裂控制储层预测方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011126267.7A CN112394416A (zh) | 2020-10-20 | 2020-10-20 | 非均质断裂控制储层预测方法及装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112394416A true CN112394416A (zh) | 2021-02-23 |
Family
ID=74596913
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011126267.7A Pending CN112394416A (zh) | 2020-10-20 | 2020-10-20 | 非均质断裂控制储层预测方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112394416A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115822563A (zh) * | 2022-12-30 | 2023-03-21 | 中国石油化工股份有限公司 | 一种超深层断控缝洞型油气藏高产井设计方法 |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102759745A (zh) * | 2011-04-28 | 2012-10-31 | 中国石油天然气集团公司 | 一种基于数字地质露头模型正演的碳酸盐岩储层预测方法 |
| CN104977618A (zh) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | 一种评价页岩气储层及寻找甜点区的方法 |
| WO2016032353A1 (ru) * | 2014-08-26 | 2016-03-03 | Закрытое Акционерное Общество "Технологии Обратных Задач" | Способ поиска залежей углеводородов, приуроченных к трещинно-кавернозным коллекторам |
| US20160363691A1 (en) * | 2015-06-15 | 2016-12-15 | Petrochina Company Limited | Physical simulation method and experiment device of fracture-cavity carbonate reservoir hydrocarbon charge |
| CN109061765A (zh) * | 2018-09-26 | 2018-12-21 | 西南石油大学 | 非均质薄砂岩互层油藏的圈闭评价方法 |
| CN110007344A (zh) * | 2019-04-08 | 2019-07-12 | 中国石油化工股份有限公司 | 一种断溶体储层连通性的地震识别方法及装置 |
| CN110794476A (zh) * | 2018-08-01 | 2020-02-14 | 中国石油化工股份有限公司 | 一种基于断溶体相控的反演方法 |
| US20200132869A1 (en) * | 2016-07-15 | 2020-04-30 | Hohai University | Seismic rock physics inversion method based on large area tight reservoir |
| CN111679318A (zh) * | 2019-03-11 | 2020-09-18 | 中国石油天然气股份有限公司 | 储层中断裂系统的识别方法及系统 |
-
2020
- 2020-10-20 CN CN202011126267.7A patent/CN112394416A/zh active Pending
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102759745A (zh) * | 2011-04-28 | 2012-10-31 | 中国石油天然气集团公司 | 一种基于数字地质露头模型正演的碳酸盐岩储层预测方法 |
| CN104977618A (zh) * | 2014-04-09 | 2015-10-14 | 中国石油集团东方地球物理勘探有限责任公司 | 一种评价页岩气储层及寻找甜点区的方法 |
| WO2016032353A1 (ru) * | 2014-08-26 | 2016-03-03 | Закрытое Акционерное Общество "Технологии Обратных Задач" | Способ поиска залежей углеводородов, приуроченных к трещинно-кавернозным коллекторам |
| US20160363691A1 (en) * | 2015-06-15 | 2016-12-15 | Petrochina Company Limited | Physical simulation method and experiment device of fracture-cavity carbonate reservoir hydrocarbon charge |
| US20200132869A1 (en) * | 2016-07-15 | 2020-04-30 | Hohai University | Seismic rock physics inversion method based on large area tight reservoir |
| CN110794476A (zh) * | 2018-08-01 | 2020-02-14 | 中国石油化工股份有限公司 | 一种基于断溶体相控的反演方法 |
| CN109061765A (zh) * | 2018-09-26 | 2018-12-21 | 西南石油大学 | 非均质薄砂岩互层油藏的圈闭评价方法 |
| CN111679318A (zh) * | 2019-03-11 | 2020-09-18 | 中国石油天然气股份有限公司 | 储层中断裂系统的识别方法及系统 |
| CN110007344A (zh) * | 2019-04-08 | 2019-07-12 | 中国石油化工股份有限公司 | 一种断溶体储层连通性的地震识别方法及装置 |
Non-Patent Citations (5)
| Title |
|---|
| XIANGWEN LI 等: "Heterogeneous reservoir prediction of ultra-deep strike-slip fault-damaged zone constrained with local seismic anomaly data", EARTH SCIENCE INFORMATICS, pages 1427 - 1441 * |
| 刘玲;沃玉进;张涛;孙炜;: "运用岩性因子预测四川盆地涪陵地区茅口组白云岩储层", 石油与天然气地质, no. 01, pages 144 - 156 * |
| 李相文 等: "基于断裂破碎带地质约束的超深致密灰岩非均质储层迭代反演方法", 2021年中国地球科学联合学术年会论文集（九）, pages 1353 - 1356 * |
| 陈胜 等: "匹配追踪子波分解重构技术在气层检测中的应用", 岩性油气藏, vol. 26, no. 06, pages 111 - 114 * |
| 鲜强;冯许魁;刘永雷;吕东;安海亭;蔡志东;王祖君;苑恒超;: "塔中地区碳酸盐岩缝洞型储层叠前流体识别", 石油与天然气地质, no. 01, pages 196 - 204 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115822563A (zh) * | 2022-12-30 | 2023-03-21 | 中国石油化工股份有限公司 | 一种超深层断控缝洞型油气藏高产井设计方法 |
| CN115822563B (zh) * | 2022-12-30 | 2023-06-09 | 中国石油化工股份有限公司 | 一种超深层断控缝洞型油气藏高产井设计方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU2013377864B2 (en) | Reservoir segment evaluation for well planning | |
| US8762118B2 (en) | Modeling hydraulic fractures | |
| US9836561B2 (en) | Optimizing multistage hydraulic fracturing design based on three-dimensional (3D) continuum damage mechanics | |
| AU2013384285B2 (en) | System, method and computer program product for predicting well production | |
| WO2015103582A1 (en) | Multistage oilfield design optimization under uncertainty | |
| US10495782B2 (en) | System, method and computer program product for determining placement of perforation intervals using facies, fluid boundaries, geobodies and dynamic fluid properties | |
| CA2881718A1 (en) | System and method for performing stimulation operations | |
| RU2601232C2 (ru) | Способы и системы для включения мест выделенных точек псевдоповерхности в модели скорости сейсмических волн | |
| US20170002630A1 (en) | Method of performing additional oilfield operations on existing wells | |
| US10605955B2 (en) | Multi-step subsidence inversion for modeling lithospheric layer thickness through geological time | |
| US8593462B2 (en) | Method and system of correlating a measured log to a predicted log | |
| CN104251135B (zh) | 大斜度井空间归位方法 | |
| CN112394416A (zh) | 非均质断裂控制储层预测方法及装置 | |
| AlQassab et al. | Estimating the size and orientation of hydraulic fractures using microseismic events | |
| CN110603370B (zh) | 确定岩层含量 | |
| Moradi et al. | Learning from Behavioral Frac Maps: A Montney Case Study in Integration of Modern Microseismic and Production Data Analyses | |
| CN115144929A (zh) | 油气藏储量确定方法、装置、存储介质及电子设备 | |
| US10460051B2 (en) | Computationally-efficient modeling of viscous fingering effect for enhanced oil recovery (EOR) agent injected at multiple injection concentrations | |
| CN110596760A (zh) | 确定油藏开发过程中孔吼半径的方法 | |
| WO2015153118A1 (en) | Bit performance analysis | |
| CN111260100B (zh) | 储层产能的预测方法及装置 | |
| US11680480B2 (en) | Multi-layer gas reservoir field development system and method | |
| US20230313649A1 (en) | System and Method to Develop Naturally Fractured Hydrocarbon Reservoirs Using A Fracture Density Index | |
| Tandon | Identification of productive zones in unconventional reservoirs | |
| Tamulonis | A methodology for integrating unconventional geologic and engineering data into a geocellular model |
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 |